From 86536dc56118f146364e11835c9cec85a6489821 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Wed, 26 Jan 2011 14:43:32 +0000
Subject: [PATCH] bugfix: the totalDirichletNodes field wasn't set up correctly
[[Imported from SVN: r6888]]
---
dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh b/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
index 7522f41a..e6b3edad 100644
--- a/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
+++ b/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
@@ -1225,9 +1225,10 @@ iterateWithContact(std::map<std::pair<std::string,std::string>, RigidBodyMotion<
int offset = 0;
for (int i=0; i<continuumName.size(); i++) {
const LeafBoundaryPatch<ContinuumGridType>& dirichletBoundary = complex_.continuum(continuumName[i]).dirichletBoundary_;
- for (int j=0; j<dirichletBoundary.numVertices(); j++)
+ int nGridVertices = dirichletBoundary.gridView().indexSet().size(dim);
+ for (int j=0; j<nGridVertices; j++)
totalDirichletNodes[offset + j] = dirichletBoundary.containsVertex(j);
- offset += dirichletBoundary.numVertices();
+ offset += nGridVertices;
}
--
GitLab