Rework frictional boundary printing

src/one-body-sample.cc

@@ -357,6 +357,12 @@ int main(int argc, char *argv[]) {
         }
       }
 
+      if (parset.get<bool>("printFrictionalBoundary")) {
+        for (size_t i = 0; i < frictionalNodes.size(); ++i)
+          if (frictionalNodes[i][0])
+            std::cout << u_diff[i][0] << " ";
+        std::cout << std::endl;
+      }
       u += u_diff;
       alpha_old = alpha;
 
@@ -378,14 +384,6 @@ int main(int argc, char *argv[]) {
       std::cout << std::endl << "Making " << timesteps << " time steps took "
                 << timer.elapsed() << "s" << std::endl;
 
-    if (parset.get<bool>("printFrictionalBoundary")) {
-      // Print displacement on frictional boundary
-      boost::format const formatter("u[%02d] = %+3e");
-      for (size_t i = 0; i < frictionalNodes.size(); ++i)
-        if (frictionalNodes[i][0])
-          std::cout << (boost::format(formatter) % i % u[i]) << std::endl;
-    }
-
     octave_writer.close();
     coefficient_writer.close();
     velocity_stepping_writer.close();