diff --git a/src/cmfd_solver.cpp b/src/cmfd_solver.cpp
index 472107d0cb..11d6d97b5d 100644
--- a/src/cmfd_solver.cpp
+++ b/src/cmfd_solver.cpp
@@ -133,7 +133,6 @@ int cmfd_linsolver_1g(const double* A_data, const double* b, double* x,
 
     // Check convergence
     err = std::sqrt(err / cmfd::dim);
-    std::cout << err << "\n";
     if (err < tol)
       return igs;
 
@@ -220,7 +219,7 @@ int cmfd_linsolver_2g(const double* A_data, const double* b, double* x,
 
         // Compute residual and update error
         double res = (tmpx[irow] - x[irow]) / tmpx[irow];
-        err = res * res;
+        err += res * res;
       }
     }