Commit 47090a35 authored by Praetorius, Simon's avatar Praetorius, Simon
Browse files

do not overwrite dirichlet node classification when iterating over intersections

parent c2cb20a7
......@@ -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);
});
}
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment