Print the coefficient of friction

src/one-body-sample.cc

@@ -282,6 +282,25 @@ int main(int argc, char *argv[]) {
           print_evolution<SingletonVectorType, VectorType>(
               frictionalNodes, *s4_new, u4, functions.get("neumannCondition"),
               run, h * run, octave_writer);
+
+        if (parset.get<bool>("printCoefficient"))
+          for (size_t i = 0; i < frictionalNodes.size(); ++i) {
+            if (!frictionalNodes[i][0])
+              continue;
+
+            double const V = u4_diff[i].two_norm();
+            double const state = (*s4_new)[i];
+            auto const L = parset.get<double>("boundary.friction.ruina.L");
+            auto const a = parset.get<double>("boundary.friction.ruina.a");
+            auto const b = parset.get<double>("boundary.friction.ruina.b");
+            auto const eta = parset.get<double>("boundary.friction.eta");
+            auto const mu = parset.get<double>("boundary.friction.mu");
+
+            std::cout << mu + a *std::log(V * eta) +
+                             b * (state - std::log(eta * L)) << std::endl;
+
+            break;
+          }
       }
 
       u4 += u4_diff;

src/one-body-sample.parset

@@ -2,7 +2,8 @@
 timesteps = 600
 
 verbose = false
-printEvolution = true
+printCoefficient = true
+printEvolution = false
 printFrictionalBoundary = false
 printProgress = false
 
@@ -13,7 +14,7 @@ enable = false
 iterations = 50
 
 [grid]
-refinements = 6
+refinements = 4
 
 [body]
 E = 1e8