From 0638ec57c2026e9761c58b8b9d3d1e29fe0f949f Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Fri, 22 Aug 2014 14:34:40 +0000
Subject: [PATCH] [bugfix] Use setRhs to properly reset the rhs of a monotone
multigrid step
[[Imported from SVN: r9857]]
---
dune/gfe/mixedriemanniantrsolver.cc | 9 +++------
1 file changed, 3 insertions(+), 6 deletions(-)
diff --git a/dune/gfe/mixedriemanniantrsolver.cc b/dune/gfe/mixedriemanniantrsolver.cc
index aaa7648d..0fd01f85 100644
--- a/dune/gfe/mixedriemanniantrsolver.cc
+++ b/dune/gfe/mixedriemanniantrsolver.cc
@@ -450,16 +450,13 @@ void MixedRiemannianTrustRegionSolver<GridType,Basis0,TargetSpace0,Basis1,Target
{
std::cout << "Iteration " << ii << std::endl;
residual0 = rhs_global0;
- residual1 = rhs_global1;
-
stiffnessMatrix01.mmv(corr_global1, residual0);
- mmgStep0->setProblem(stiffnessMatrix00, corr_global0, residual0);
-
+ mmgStep0->setRhs(residual0);
mmgStep0->iterate();
+ residual1 = rhs_global1;
stiffnessMatrix10.mmv(corr_global0, residual1);
- mmgStep1->setProblem(stiffnessMatrix11, corr_global1, residual1);
-
+ mmgStep1->setRhs(residual1);
mmgStep1->iterate();
}
--
GitLab