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>
31 template<
typename P,
typename EG,
typename LFS,
typename T>
32 void volume (
const P& param,
const EG& eg,
const LFS& lfs, T& trafo)
const
35 if (eg.entity().partitionType()==Dune::InteriorEntity)
39 else if (eg.entity().partitionType()==Dune::GhostEntity){
40 typename T::RowType empty;
41 typedef typename LFS::Traits::SizeType size_type;
42 for (size_type i=0; i<lfs.size(); i++){
43 trafo[lfs.dofIndex(i)] = empty;
54 #endif // DUNE_PDELAB_CONSTRAINTS_P0GHOST_HH