constraints: import interface function from PDELAb

src/constraints.hh

+#ifndef PACXX_PROJECTSEMINAR_2019_SRC_CONSTRAINTS_HH
+#define PACXX_PROJECTSEMINAR_2019_SRC_CONSTRAINTS_HH
+
+namespace PPS {
+  template<typename P, typename GFS, typename CG>
+  void constraints(const P& p, const GFS& gfs, CG& cg,
+                   const bool verbose = false)
+  {
+    // clear global constraints
+    cg.clear();
+    Dune::PDELab::ConstraintsAssemblerHelper<P, GFS, CG, Dune::PDELab::IsGridFunction<P>::value>::assemble(p,gfs,cg,verbose);
+  }
+}
+
+#endif // PACXX_PROJECTSEMINAR_2019_SRC_CONSTRAINTS_HH

src/driver.hh

@@ -47,6 +47,7 @@
 #endif
 //== }
 #include"problem.hh"
+#include "constraints.hh"
 #include "gridoperator.hh"
 #include "interpolate.hh"
 #include "linearsolver.hh"
@@ -92,7 +93,7 @@ void driver (const GV& gv, Dune::ParameterTree& ptree)
   typedef typename GFS::template
     ConstraintsContainer<RF>::Type CC;
   CC cc;
-  Dune::PDELab::constraints(b,gfs,cc); // assemble constraints
+  PPS::constraints(b,gfs,cc); // assemble constraints
   std::cout << "constrained dofs=" << cc.size() << " of "
             << gfs.globalSize() << std::endl;
   std::vector<std::size_t> cc_plain;