diff --git a/src/myblockproblem.hh b/src/myblockproblem.hh index 47d518406b1d30a6df8ae922cc4376c7af55adf1..6e3bb6083bfae286e85fb693944b6702b9b56c55 100644 --- a/src/myblockproblem.hh +++ b/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;