Print evolution of state

src/one-body-sample.cc

@@ -452,6 +452,18 @@ int main(int argc, char *argv[]) {
       u4 += u4_diff;
       s4_old = *s4_new;
 
+      if (parset.get<bool>("printEvolution")) {
+        // Print displacement on frictional boundary
+        boost::format const formatter("s[%02d] = %+3e, "
+                                      "%|40t|u[%02d] = %+3e");
+        for (size_t i = 0; i < frictionalNodes.size(); ++i)
+          if (frictionalNodes[i][0]) {
+            std::cout << boost::format(formatter) % i % (*s4_new)[i] % i % u4[i]
+                      << std::endl;
+            break; // Only print the first element
+          }
+      }
+
       // Compute von Mises stress and write everything to a file
       if (parset.get<bool>("writeVTK")) {
         auto const displacement =

src/one-body-sample.parset

@@ -2,6 +2,7 @@
 timesteps = 30
 
 verbose = false
+printEvolution = true
 printFrictionalBoundary = false
 printProgress = false
 printDifference = false