diff --git a/src/one-body-sample.org b/src/one-body-sample.org
index ab8164bf39528973208cae9143143c23af2639f9..6055ddd29492560d496e866573c0a1b90b6af9d1 100644
--- a/src/one-body-sample.org
+++ b/src/one-body-sample.org
@@ -485,49 +485,48 @@
       auto const state_fpi_max
         = parset.get<size_t>("solver.tnnmg.fixed_point_iterations");
       for (size_t run = 1; run <= timesteps; ++run) {
+        stateUpdater->nextTimeStep();
+        timeSteppingScheme->nextTimeStep();
+  
         double const time = tau * run;
-        {
-          stateUpdater->nextTimeStep();
-          timeSteppingScheme->nextTimeStep();
-  
-          VectorType ell(finestSize);
-          assemble_neumann<GridType, GridView, SmallVector, P1Basis>
-            (leafView, p1Basis, neumannNodes, ell, neumannFunction, time);
-          ell += gravityFunctional;
-  
-          MatrixType problem_A;
-          VectorType problem_rhs(finestSize);
-          VectorType problem_iterate(finestSize);
-  
-          stateUpdater->setup(tau);
-          timeSteppingScheme->setup(ell, tau, time,
-                                    problem_rhs, problem_iterate, problem_A);
-  
-          VectorType u_saved;
-          for (size_t state_fpi = 1; state_fpi <= state_fpi_max; ++state_fpi) {
-            <<setupAndSolveProblem>>;
-  
-            timeSteppingScheme->postProcess(problem_iterate);
-            timeSteppingScheme->extractDisplacement(u);
-            timeSteppingScheme->extractVelocity(ud);
-  
-            stateUpdater->solve(ud);
-            stateUpdater->extractState(alpha);
-  
-            <<printInnerProgress>>;
-            if (state_fpi > 1 && energyNorm.diff(u_saved, u)
-                < parset.get<double>("solver.tnnmg.fixed_point_tolerance"))
-              break;
-            else
-              u_saved = u;
-  
-            if (state_fpi == state_fpi_max)
-              std::cerr << "[ref = " << refinements
-                        << "]: FPI did not converge after "
-                        << state_fpi_max << " iterations" << std::endl;
-          }
-          <<printOuterProgress>>;
+  
+        VectorType ell(finestSize);
+        assemble_neumann<GridType, GridView, SmallVector, P1Basis>
+          (leafView, p1Basis, neumannNodes, ell, neumannFunction, time);
+        ell += gravityFunctional;
+  
+        MatrixType problem_A;
+        VectorType problem_rhs(finestSize);
+        VectorType problem_iterate(finestSize);
+  
+        stateUpdater->setup(tau);
+        timeSteppingScheme->setup(ell, tau, time,
+                                  problem_rhs, problem_iterate, problem_A);
+  
+        VectorType u_saved;
+        for (size_t state_fpi = 1; state_fpi <= state_fpi_max; ++state_fpi) {
+          <<setupAndSolveProblem>>;
+  
+          timeSteppingScheme->postProcess(problem_iterate);
+          timeSteppingScheme->extractDisplacement(u);
+          timeSteppingScheme->extractVelocity(ud);
+  
+          stateUpdater->solve(ud);
+          stateUpdater->extractState(alpha);
+  
+          <<printInnerProgress>>;
+          if (state_fpi > 1 && energyNorm.diff(u_saved, u)
+              < parset.get<double>("solver.tnnmg.fixed_point_tolerance"))
+            break;
+          else
+            u_saved = u;
+  
+          if (state_fpi == state_fpi_max)
+            std::cerr << "[ref = " << refinements
+                      << "]: FPI did not converge after "
+                      << state_fpi_max << " iterations" << std::endl;
         }
+        <<printOuterProgress>>;
   
         <<writeData>>;
         <<assembleAndWriteStress>>;