diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 42b96fc913e4ed39b63b6bca3eeaf9d34605bc32..637549e990f24dd941f21339ae331443718cdd27 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -478,30 +478,29 @@ int main(int argc, char *argv[]) {
             }
           }
         }
+
+        if (parset.get<bool>("printEvolution")) {
+          // Find the first node that belongs to the frictional boundary
+          for (size_t i = 0; i < frictionalNodes.size(); ++i) {
+            if (!frictionalNodes[i][0])
+              continue;
+
+            boost::format const formatter("s[%03d] = %+3e, "
+                                          "%|40t|u[%03d] = %+3e");
+            std::cout << boost::format(formatter) % run % (*s4_new)[i] % run %
+                             u4[i] << std::endl;
+            double out;
+            functions.get("sampleFunction").evaluate(h * run, out);
+            octave_writer << (*s4_new)[i] << " " << u4[i][0] * 1e6 << " " << out
+                          << std::endl;
+            break;
+          }
+        }
       }
 
       u4 += u4_diff;
       s4_old = *s4_new;
 
-      if (!parset.get<bool>("printEvolution"))
-        continue;
-
-      // Find the first node that belongs to the frictional boundary
-      for (size_t i = 0; i < frictionalNodes.size(); ++i) {
-        if (!frictionalNodes[i][0])
-          continue;
-
-        boost::format const formatter("s[%03d] = %+3e, "
-                                      "%|40t|u[%03d] = %+3e");
-        std::cout << boost::format(formatter) % run % (*s4_new)[i] % run % u4[i]
-                  << std::endl;
-        double out;
-        functions.get("sampleFunction").evaluate(h * run, out);
-        octave_writer << (*s4_new)[i] << " " << u4[i][0] * 1e6 << " " << out
-                      << std::endl;
-        break;
-      }
-
       // Compute von Mises stress and write everything to a file
       if (parset.get<bool>("writeVTK")) {
         auto const displacement =