2 #ifndef DUNE_PDELAB_CONSTRAINTS_P0GHOST_HH
3 #define DUNE_PDELAB_CONSTRAINTS_P0GHOST_HH
5 #include "../common/geometrywrapper.hh"
6 #include<dune/grid/common/gridenums.hh>
8 #include <dune/typetree/typetree.hh>
19 :
public TypeTree::LeafNode
34 template<
typename P,
typename EG,
typename LFS,
typename T>
35 void volume (
const P& param,
const EG& eg,
const LFS& lfs, T& trafo)
const
38 if (eg.entity().partitionType()==Dune::InteriorEntity)
42 else if (eg.entity().partitionType()==Dune::GhostEntity){
43 typename T::RowType empty;
44 typedef typename LFS::Traits::SizeType size_type;
45 for (size_type i=0; i<lfs.size(); i++){
46 trafo[lfs.dofIndex(i)] = empty;
57 #endif // DUNE_PDELAB_CONSTRAINTS_P0GHOST_HH