Skip to content
Snippets Groups Projects
Commit 38118ee7 authored by Praetorius, Simon's avatar Praetorius, Simon
Browse files

Merge branch 'issue/dirichlet_nodes_overwrite' into 'master'

do not overwrite dirichlet node classification when iterating over intersections

See merge request !62
parents c4f8fac9 47090a35
No related branches found
No related tags found
1 merge request!62do not overwrite dirichlet node classification when iterating over intersections
...@@ -20,8 +20,9 @@ init(RB const& rowBasis, CB const& colBasis) ...@@ -20,8 +20,9 @@ init(RB const& rowBasis, CB const& colBasis)
using LV = typename CB::LocalView; using LV = typename CB::LocalView;
using IS = typename CB::GridView::Intersection; using IS = typename CB::GridView::Intersection;
dirichletNodes_.resize(colBasis.dimension()); dirichletNodes_.resize(colBasis.dimension());
std::fill(dirichletNodes_.begin(), dirichletNodes_.end(), false);
forEachBoundaryDOF(colBasis, [&](int localIndex, LV const& localView, IS const& intersection) { 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);
}); });
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment