From cc41c30910741763d558b789ffe950a0460b70b7 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Wed, 30 Jan 2013 18:24:39 +0100
Subject: [PATCH] Control damping through a parameter
---
src/one-body-sample.cc | 2 +-
src/one-body-sample.parset | 1 +
2 files changed, 2 insertions(+), 1 deletion(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index fee77c9e..77520b29 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 c17a57ad..5b6a7433 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
--
GitLab