rewritten extraction of dirichletNodes in DirichletBC

src/amdis/DirichletBC.inc.hpp

 #pragma once
 
+#include <dune/functions/functionspacebases/boundarydofs.hh>
 #include <dune/functions/functionspacebases/interpolate.hh>
 #include <dune/functions/functionspacebases/subspacebasis.hh>
 #include <amdis/LinearAlgebra.hpp>
@@ -13,8 +14,13 @@ namespace AMDiS
   void DirichletBC<WorldVector, Range>::initImpl(
       RowBasis const& rowBasis, Dune::TypeTree::LeafNodeTag)
   {
-    using Dune::Functions::interpolate;
-    interpolate(rowBasis, dirichletNodes_, predicate_);
+    using Dune::Functions::forEachBoundaryDOF;
+    using LV = typename RowBasis::LocalView;
+    using I = typename RowBasis::GridView::Intersection;
+    dirichletNodes_.resize(rowBasis.dimension());
+    forEachBoundaryDOF(rowBasis, [&](int localIndex, LV const& localView, I const& intersection) {
+      dirichletNodes_[localView.index(localIndex)] = predicate_(intersection.geometry().center());
+    });
   }
 
   template <class WorldVector, class Range>