diff --git a/src/myblockproblem.hh b/src/myblockproblem.hh
index 283f303f636b7c35a4be90b7b1172fb392f9be2d..4881f59ba352709b1b89fdd7cf6432438dec42a4 100644
--- a/src/myblockproblem.hh
+++ b/src/myblockproblem.hh
@@ -108,9 +108,9 @@ class MyBlockProblem<ConvexProblemTypeTEMPLATE>::IterateObject {
         if (j == m)
           localA = &(*it); // localA = &A[m][m]
 
-        (*it).umv(u[j], localb); // localb += A[m][j] * u[j]
+        it->umv(u[j], localb); // localb += A[m][j] * u[j]
       }
-      localb -= problem.f[m];
+      localb -= problem.f[m]; // localb -= b[m]
       assert(localA != NULL);
 
       // FIXME: Hardcoding a fixed function here for now