From 71f6cc1bb552d34ca25f5fe8410430f35931159d Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Mon, 13 Jun 2011 16:48:42 +0000
Subject: [PATCH] Bugfix for the case a dof is Dirichlet in a few directions,
and Neumann in others
[[Imported from SVN: r7428]]
---
dune/gfe/riemanniantrsolver.cc | 7 ++++---
1 file changed, 4 insertions(+), 3 deletions(-)
diff --git a/dune/gfe/riemanniantrsolver.cc b/dune/gfe/riemanniantrsolver.cc
index 0f477377..f8028222 100644
--- a/dune/gfe/riemanniantrsolver.cc
+++ b/dune/gfe/riemanniantrsolver.cc
@@ -226,10 +226,11 @@ void RiemannianTrustRegionSolver<GridType,TargetSpace>::solve()
for (; cIt!=cEndIt; ++cIt) {
for (int k=0; k<blocksize; k++) {
- if (ignoreNodes_->operator[](j)[k])
+ if (ignoreNodes_->operator[](j)[k]) {
(*cIt)[k] = 0;
- if (j==cIt.index())
- (*cIt)[k][k] = 1;
+ if (j==cIt.index())
+ (*cIt)[k][k] = 1;
+ }
}
}
--
GitLab