Control damping through a parameter

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);
         }

src/one-body-sample.parset

@@ -23,6 +23,7 @@ width = 5
 
 [solver]
 tolerance = 1e-10
+damping = 0.0
 
 [solver.tnnmg]
 maxiterations = 1000000