Nuke useless option

src/one-body-sample.cc

@@ -268,11 +268,9 @@ int main(int argc, char *argv[]) {
             u4_diff[i][0] /= refinement_factor;
           }
 
-        auto const fixed_point_iterations =
-            parset.get<bool>("boundary.friction.state.evolve")
-                ? parset.get<size_t>("solver.tnnmg.fixed_point_iterations")
-                : 1;
-        for (size_t state_fpi = 0; state_fpi < fixed_point_iterations;
+        for (size_t state_fpi = 0;
+             state_fpi <
+                 parset.get<size_t>("solver.tnnmg.fixed_point_iterations");
              ++state_fpi) {
           auto myGlobalNonlinearity =
               assemble_nonlinearity<MatrixType, VectorType>(
@@ -291,9 +289,6 @@ int main(int argc, char *argv[]) {
               false); // absolute error
           overallSolver.solve();
 
-          if (!parset.get<bool>("boundary.friction.state.evolve"))
-            continue;
-
           for (size_t i = 0; i < frictionalNodes.size(); ++i) {
             if (frictionalNodes[i][0]) {
               double const unorm = u4_diff[i].two_norm();

src/one-body-sample.parset

@@ -71,7 +71,6 @@ eta = 1
 model = Exponential
 
 [boundary.friction.state]
-evolve = true
 # log(alpha(0))
 initial = -2.23805
 model = Dieterich # Ruina