From 3917d6ac256d84a1035e0a044c698a9cfa94ba8e Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Wed, 12 Sep 2012 22:43:43 +0200
Subject: [PATCH] Revert "Use ud as the accuracy criterion"
Use u again since with absolute errors, it's independent of refinement
in time (whereas ud is not)
---
src/one-body-sample.cc | 10 +++++-----
1 file changed, 5 insertions(+), 5 deletions(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 0e0ce17a..b44d6687 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -338,8 +338,7 @@ int main(int argc, char *argv[]) {
timeSteppingScheme->setup(ell, tau, time, problem_rhs, problem_iterate,
problem_A);
- VectorType ud_saved(finestSize);
- ud_saved = 0.0; // FIXME
+ VectorType u_saved;
auto const state_fpi_max =
parset.get<size_t>("solver.tnnmg.fixed_point_iterations");
for (size_t state_fpi = 1; state_fpi <= state_fpi_max; ++state_fpi) {
@@ -388,11 +387,12 @@ int main(int argc, char *argv[]) {
std::cerr << '.';
std::cerr.flush();
}
- if (energyNorm.diff(ud_saved, ud) <
- parset.get<double>("solver.tnnmg.fixed_point_tolerance"))
+ if (state_fpi > 1 &&
+ energyNorm.diff(u_saved, u) <
+ parset.get<double>("solver.tnnmg.fixed_point_tolerance"))
break;
else
- ud_saved = ud;
+ u_saved = u;
if (state_fpi == state_fpi_max)
std::cerr << "[ref = " << refinements
--
GitLab