diff --git a/src/ipglasma.cpp b/src/ipglasma.cpp
index a683f5dec1b0de7281004ca4da6d004c6c5bcd67..4da94cc5cca4752798e594c1f79c175aa8c2fb39 100644
--- a/src/ipglasma.cpp
+++ b/src/ipglasma.cpp
@@ -630,22 +630,16 @@ double IPGlasma::SumAmplitudes()
     WilsonLine quark;
     int thrownAway = 0;
     int iWasHere = 0;
-    #pragma omp parallel for reduction(+:sum,thrownAway,iWasHere) collapse(4)
+    #pragma omp parallel for reduction(+:sum)
     for(int i =0; i < L; i += 1){
         for(int j =0; j < L; j += 1){
-            
+            int q1[2] = {i,j};
+            quark = GetWilsonLine(WilsonLineCoordinate(i, j))
+            double r1 = DistanceToOrigin(q1);
             //#pragma omp parallel for reduction(+:sum) collapse(2)
             for(int k =0; k < L; k +=1){
-//                cout << "# Thread " << omp_get_thread_num() << " is processing (" << i << ", " << j << ", " << k <<")" << endl;
                 for(int l =0; l < L; l +=1){
-                    cout << "# Hello from inside the l loop!" << endl;
-            //TODO: pointer instead of copy?
-                    //TODO: faster to get these two layers above?
-                    iWasHere += 1;
-                    int q1[2] = {i,j};
-                    //quark = GetWilsonLine(WilsonLineCoordinate(i, j));
                     int q2[2] = {k,l};
-		    double r1 = DistanceToOrigin(q1);
                     double r2 = DistanceToOrigin(q2);
                     if (r1 > 1. && r2 > 1.){
                        // cout << "#!!Throwing away large dipole separation!" << endl;