From 6384b92de0a888e1687fa449cba831f5fc26b8fa Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Sun, 25 May 2014 04:41:20 +0000
Subject: [PATCH] [bugfix] The solver must always use the global Dirichlet
vector
Even if there is only one processor. Because even then the local
and global vectors differ in their ordering.
[[Imported from SVN: r9757]]
---
dune/gfe/riemanniantrsolver.cc | 5 +----
1 file changed, 1 insertion(+), 4 deletions(-)
diff --git a/dune/gfe/riemanniantrsolver.cc b/dune/gfe/riemanniantrsolver.cc
index d3c02028..592e9545 100644
--- a/dune/gfe/riemanniantrsolver.cc
+++ b/dune/gfe/riemanniantrsolver.cc
@@ -105,10 +105,7 @@ setup(const GridType& grid,
MonotoneMGStep<MatrixType, CorrectionType>* mmgStep = new MonotoneMGStep<MatrixType, CorrectionType>;
mmgStep->setMGType(mu, nu1, nu2);
- if (mpiHelper.size()==1)
- mmgStep->ignoreNodes_ = &dirichletNodes;
- else
- mmgStep->ignoreNodes_ = globalDirichletNodes;
+ mmgStep->ignoreNodes_ = globalDirichletNodes;
mmgStep->basesolver_ = baseSolver;
mmgStep->setSmoother(presmoother, postsmoother);
mmgStep->obstacleRestrictor_= new MandelObstacleRestrictor<CorrectionType>();
--
GitLab