From 47090a3578b93ded8944484a63fa95b3668b5cf5 Mon Sep 17 00:00:00 2001
From: Simon Praetorius <simon.praetorius@tu-dresden.de>
Date: Thu, 6 Jun 2019 17:49:19 +0200
Subject: [PATCH] do not overwrite dirichlet node classification when iterating
over intersections
---
src/amdis/DirichletBC.inc.hpp | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/src/amdis/DirichletBC.inc.hpp b/src/amdis/DirichletBC.inc.hpp
index 6f285d24..0a1a576c 100644
--- a/src/amdis/DirichletBC.inc.hpp
+++ b/src/amdis/DirichletBC.inc.hpp
@@ -20,8 +20,9 @@ init(RB const& rowBasis, CB const& colBasis)
using LV = typename CB::LocalView;
using IS = typename CB::GridView::Intersection;
dirichletNodes_.resize(colBasis.dimension());
+ std::fill(dirichletNodes_.begin(), dirichletNodes_.end(), false);
forEachBoundaryDOF(colBasis, [&](int localIndex, LV const& localView, IS const& intersection) {
- dirichletNodes_[localView.index(localIndex)] = onBoundary(intersection);
+ dirichletNodes_[localView.index(localIndex)] = dirichletNodes_[localView.index(localIndex)] || onBoundary(intersection);
});
}
--
GitLab