diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 9ffdeb08b1fa68ab4c55e962190398380f9222aa..d745f2566e63052be7d909ddd695be6ca20bdadf 100644
--- a/src/one-body-sample.cc
+++ b/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
diff --git a/src/one-body-sample.parset b/src/one-body-sample.parset
index 4a931af64e1686cf9c9f03e05e78d97b2fc11f61..5b0a5ac00a157b9b2edc6bfa70ec425b85d020cf 100644
--- a/src/one-body-sample.parset
+++ b/src/one-body-sample.parset
@@ -4,6 +4,7 @@ timesteps = 100
 verbose = false
 printFrictionalBoundary = false
 printProgress = true
+printDifference = false
 
 useNonlinearGS = false # true
 useGS = false