p0.hh

// -*- tab-width: 4; indent-tabs-mode: nil -*-
#ifndef DUNE_PDELAB_CONSTRAINTS_P0_HH
#define DUNE_PDELAB_CONSTRAINTS_P0_HH
 
#include "../common/geometrywrapper.hh"
#include <dune/typetree/typetree.hh>
 
namespace Dune {
  namespace PDELab {
 
 
    class P0ParallelConstraints
      : public TypeTree::LeafNode
    {
    public:
      enum{doBoundary=false};
      enum{doProcessor=true};
      enum{doSkeleton=false};
      enum{doVolume=false};
 
 
      template<typename I, typename LFS, typename T>
      void processor (const Dune::PDELab::IntersectionGeometry<I>& ig,
                      const LFS& lfs, T& trafo) const
      {
        typename T::RowType empty;
        typedef typename LFS::Traits::SizeType size_type;
        for (size_type i=0; i<lfs.size(); i++){
          trafo[lfs.dofIndex(i)] = empty;
        }
      }
    };
 
  }
}
 
#endif // DUNE_PDELAB_CONSTRAINTS_P0_HH