diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 40425882ac6a92566f96a7c8e34a8c17fa8c4c17..1242952618edff01d63d8fe7d8225d839a088885 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -342,29 +342,6 @@ int main(int argc, char *argv[]) {
         MyConvexProblemType myConvexProblem(stiffnessMatrix,
                                             *myGlobalNonlinearity, b4, u4_diff);
 
-        // FIXME: This is sort of nasty
-        // Choose an initial iterate that is admissable
-        // and non-zero in as many places as possible
-        if (run == 1)
-          for (size_t i = 0; i < ignoreNodes.size(); ++i) {
-            size_t const count = ignoreNodes[i].count();
-
-            if (count == dim)
-              continue;
-
-            if (count == 0) {
-              u4_diff[i] = 1e-4;
-              continue;
-            }
-
-            assert(count == 1);
-            assert(ignoreNodes[i][1]); // Only the y coordinate is allowed to be
-                                       // fixed
-            for (int j = 0; j < dim; ++j)
-              if (j != 1)
-                u4_diff[i][j] = 1e-4;
-          }
-
         MyBlockProblemType myBlockProblem(parset, myConvexProblem);
         multigridStep.setProblem(u4_diff, myBlockProblem);