From 59547d328b3a89753151b07ecaa7807355ae2ba8 Mon Sep 17 00:00:00 2001
From: podlesny <podlesny@zedat.fu-berlin.de>
Date: Sun, 14 Mar 2021 17:24:31 +0100
Subject: [PATCH] only do 1 iteration step, add debug print
---
dune/tectonic/spatial-solving/fixedpointiterator.cc | 9 ++++++++-
1 file changed, 8 insertions(+), 1 deletion(-)
diff --git a/dune/tectonic/spatial-solving/fixedpointiterator.cc b/dune/tectonic/spatial-solving/fixedpointiterator.cc
index 586272ff..f21aa580 100644
--- a/dune/tectonic/spatial-solving/fixedpointiterator.cc
+++ b/dune/tectonic/spatial-solving/fixedpointiterator.cc
@@ -124,6 +124,8 @@ FixedPointIterator<Factory, NBodyAssembler, Updaters, ErrorNorms>::run(
Vector total_v;
nBodyAssembler_.nodalToTransformed(velocityIterates, total_v);
+ print(velocityIterates, "velocityIterates start:");
+
for (fixedPointIteration = 0; fixedPointIteration < fixedPointMaxIterations_;
++fixedPointIteration) {
@@ -148,8 +150,11 @@ FixedPointIterator<Factory, NBodyAssembler, Updaters, ErrorNorms>::run(
updaters.state_->extractAlpha(newAlpha);
nBodyAssembler_.postprocess(total_v, velocityIterates);
+ //Rprint(velocityIterates, "velocityIterates loop:");
updaters.rate_->postProcess(velocityIterates);
- bool breakCriterion = displacementCriterion(updaters, last_u); //stateCriterion(alpha, newAlpha);
+ bool breakCriterion = true; //displacementCriterion(updaters, last_u); //stateCriterion(alpha, newAlpha);
+
+ printRegularityTruncation(globalFriction_, total_v);
if (lambda_ < 1e-12 or breakCriterion) {
fixedPointIteration++;
@@ -171,6 +176,8 @@ FixedPointIterator<Factory, NBodyAssembler, Updaters, ErrorNorms>::run(
nBodyAssembler_.postprocess(total_v, velocityIterates);
updaters.rate_->postProcess(velocityIterates);
+ print(velocityIterates, "velocityIterates end:");
+
// Cannot use return { fixedPointIteration, multigridIterations };
// with gcc 4.9.2, see also http://stackoverflow.com/a/37777814/179927
FixedPointIterationCounter ret;
--
GitLab