Make printing of the difference optional

src/one-body-sample.cc

@@ -419,20 +419,22 @@ int main(int argc, char *argv[]) {
     }
     std::cout << std::endl;
 
-    VectorType diff2 = u2;
-    diff2 -= u1;
-    std::cout << "sup |u1 - u2| = " << diff2.infinity_norm() << ", "
-              << "|u1 - u2| = " << diff2.two_norm() << std::endl;
-
-    VectorType diff3 = u3;
-    diff3 -= u1;
-    std::cout << "sup |u1 - u3| = " << diff3.infinity_norm() << ", "
-              << "|u1 - u3| = " << diff3.two_norm() << std::endl;
-
-    VectorType diff4 = u4;
-    diff4 -= u1;
-    std::cout << "sup |u1 - u4| = " << diff4.infinity_norm() << ", "
-              << "|u1 - u4| = " << diff4.two_norm() << std::endl;
+    if (parset.get<bool>("printDifference")) {
+      VectorType diff2 = u2;
+      diff2 -= u1;
+      std::cout << "sup |u1 - u2| = " << diff2.infinity_norm() << ", "
+                << "|u1 - u2| = " << diff2.two_norm() << std::endl;
+
+      VectorType diff3 = u3;
+      diff3 -= u1;
+      std::cout << "sup |u1 - u3| = " << diff3.infinity_norm() << ", "
+                << "|u1 - u3| = " << diff3.two_norm() << std::endl;
+
+      VectorType diff4 = u4;
+      diff4 -= u1;
+      std::cout << "sup |u1 - u4| = " << diff4.infinity_norm() << ", "
+                << "|u1 - u4| = " << diff4.two_norm() << std::endl;
+    }
 
     if (parset.get<bool>("printFrictionalBoundary")) {
       // Print displacement on frictional boundary

src/one-body-sample.parset

@@ -4,6 +4,7 @@ timesteps = 100
 verbose = false
 printFrictionalBoundary = false
 printProgress = true
+printDifference = false
 
 useNonlinearGS = false # true
 useGS = false