Cleanup

src/myblockproblem.hh

@@ -109,10 +109,6 @@ class MyBlockProblem<ConvexProblemTypeTEMPLATE>::IterateObject {
           localA = &(*it); // localA = &A[m][m]
 
         (*it).umv(u[j], localb); // localb += A[m][j] * u[j]
-        LocalMatrixType foo = problem.A[j][m];
-        foo -= problem.A[m][j];
-        // if (foo.infinity_norm() != 0)
-        //   std::cout << "Indices: " << col << ", " << m << std::endl;
       }
       localb -= problem.f[m];
       assert(localA != NULL);
@@ -124,11 +120,10 @@ class MyBlockProblem<ConvexProblemTypeTEMPLATE>::IterateObject {
 
       LocalVectorType ui_copy = ui;
 
-      // std::cout << "Norm before: " << ui.two_norm() << std::endl;
       LocalVectorType correction;
       for (size_t i = 1; i <= 5; ++i) {
         Dune::minimise(localJ, ui_copy, correction);
-        ui_copy += correction; // This should be correct
+        ui_copy += correction;
       }
 
       ui = ui_copy;