diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 1be2924711df161322b12d8d100f0939fe350f83..f7dedf7750f30841c45b1b79b00fbdba2e11c4b7 100644
--- a/src/one-body-sample.cc
+++ b/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();
diff --git a/src/one-body-sample.parset b/src/one-body-sample.parset
index 21267d844609c718fbc7b5a68c6b587057f07557..2702a3e490a35e91eb1386fb65555b1ab90098e7 100644
--- a/src/one-body-sample.parset
+++ b/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