2 #ifndef DUNE_PDELAB_MFDCONSTRAINTS_HH
3 #define DUNE_PDELAB_MFDCONSTRAINTS_HH
5 #include<dune/geometry/referenceelements.hh>
6 #include<dune/geometry/type.hh>
31 template<
typename B,
typename IG,
typename LFS,
typename T>
32 void boundary(
const B& b,
const IG&
ig,
const LFS& lfs, T& trafo)
const
34 static const unsigned int dimIntersection = IG::dimension - 1;
35 typedef typename IG::ctype ctype;
37 GeometryType gt = ig.intersection().type();
39 Dune::FieldVector<ctype,dimIntersection> center
40 = ReferenceElements<ctype,dimIntersection>::general(gt).position(0,0);
41 if(b.isDirichlet(ig, center) )
43 typename T::RowType empty;
44 trafo[ig.intersectionIndex()] = empty;
Mimetic Constraints.
Definition: mfd.hh:16
const IG & ig
Definition: common/constraints.hh:146
void boundary(const B &b, const IG &ig, const LFS &lfs, T &trafo) const
boundary constraints
Definition: mfd.hh:32