From 0c830092d05a7944beb4e05ec427db4f66118b7c Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Wed, 1 Oct 2014 15:53:23 +0000
Subject: [PATCH] [bugfix] Use the proper ignoreNodes fields for gradient norm
computations
[[Imported from SVN: r9909]]
---
dune/gfe/riemanniantrsolver.cc | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/dune/gfe/riemanniantrsolver.cc b/dune/gfe/riemanniantrsolver.cc
index 7ef0bc72..2ab1da2b 100644
--- a/dune/gfe/riemanniantrsolver.cc
+++ b/dune/gfe/riemanniantrsolver.cc
@@ -331,7 +331,7 @@ void RiemannianTrustRegionSolver<GridType,TargetSpace>::solve()
CorrectionType gradient = rhs_global;
for (size_t j=0; j<gradient.size(); j++)
for (int k=0; k<gradient[j].size(); k++)
- if ((*ignoreNodes_)[j][k])
+ if ((*mgStep->ignoreNodes_)[j][k]) // global Dirichlet nodes set
gradient[j][k] = 0;
if (this->verbosity_ == Solver::FULL and rank==0)
--
GitLab