diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index fee77c9e912d1ba336bad0040d3050b6429ba978..77520b29ed614ad7350f3d2469a4cf518b274aca 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -392,6 +392,7 @@ int main(int argc, char *argv[]) {
       VectorType u_saved;
       double const fixedPointTolerance =
           parset.get<double>("solver.tnnmg.fixed_point_tolerance");
+      double const damping = parset.get<double>("solver.damping");
       for (size_t state_fpi = 1; state_fpi <= state_fpi_max; ++state_fpi) {
         stateUpdater->solve(ud);
         if (state_fpi == 1)
@@ -399,7 +400,6 @@ int main(int argc, char *argv[]) {
         else {
           SingletonVectorType computed_state;
           stateUpdater->extractState(computed_state);
-          double const damping = 0.5;
           alpha *= damping;
           alpha.axpy(1.0 - damping, computed_state);
         }
diff --git a/src/one-body-sample.parset b/src/one-body-sample.parset
index c17a57ad00054777fdf1b5d29366884b38fc216a..5b6a743344a638bee26c49aa97a606a18ee55ed5 100644
--- a/src/one-body-sample.parset
+++ b/src/one-body-sample.parset
@@ -23,6 +23,7 @@ width = 5
 
 [solver]
 tolerance = 1e-10
+damping = 0.0
 
 [solver.tnnmg]
 maxiterations = 1000000