diff --git a/src/main/c/src/math/ode/conjugate_gradient.c b/src/main/c/src/math/ode/conjugate_gradient.c
index 8d46705e..279e3265 100644
--- a/src/main/c/src/math/ode/conjugate_gradient.c
+++ b/src/main/c/src/math/ode/conjugate_gradient.c
@@ -230,27 +230,19 @@ float solver_conjugate_gradient_iterate_serial(float * phi, float * phi0, float
         printf("Divide by 0! %f \n", denominator);
         fflush(stdout);
     }
-     alpha = convergence / denominator;
+    alpha = convergence / denominator;
+    r_new_dot = 0;
     for(k=1; k<DIM-1; k++){
         for(j=1; j<DIM-1; j++){
             for(i = 1; i < DIM-1; i++){
                 phi[ode_index(i,j,k,DIM)] = phi[ode_index(i,j,k,DIM)] + alpha * p[ode_index(i,j,k,DIM)];
                 r[ode_index(i,j,k,DIM)] = r[ode_index(i,j,k,DIM)] - alpha * A[ode_index(i,j,k,DIM)];
-            }
-        }
-    }
-
-    r_new_dot = 0;
-    for(k=1; k<DIM-1; k++){
-        for(j=1; j<DIM-1; j++){
-            for(i = 1; i < DIM-1; i++){
                 r_new_dot = r_new_dot + r[ode_index(i,j,k,DIM)] * r[ode_index(i,j,k,DIM)];
             }
         }
     }
 
     beta = r_new_dot / convergence;
-
     for(k=1; k<DIM-1; k++){
         for(j=1; j<DIM-1; j++){
             for(i = 1; i < DIM-1; i++){