rt0cube2dfem.hh

// -*- tab-width: 4; indent-tabs-mode: nil -*-
#ifndef DUNE_PDELAB_FINITEELEMENTMAP_RT0CUBE2DFEM_HH
#define DUNE_PDELAB_FINITEELEMENTMAP_RT0CUBE2DFEM_HH
 
#include<vector>
#include<dune/localfunctions/raviartthomas/raviartthomas0cube2d.hh>
#include"finiteelementmap.hh"
 
namespace Dune {
  namespace PDELab {
 
    template<typename GV, typename D, typename R>
    class RT0Cube2DLocalFiniteElementMap :
      public RTLocalFiniteElementMap<
        GV,
        Dune::RT0Cube2DLocalFiniteElement<D,R>,
        RT0Cube2DLocalFiniteElementMap<GV,D,R>,
        16>
    {
      typedef Dune::RT0Cube2DLocalFiniteElement<D,R> FE;
 
    public:
      typedef LocalFiniteElementMapTraits<FE> Traits;
 
      RT0Cube2DLocalFiniteElementMap (const GV& gv)
        : RTLocalFiniteElementMap<
          GV,
          Dune::RT0Cube2DLocalFiniteElement<D,R>,
          RT0Cube2DLocalFiniteElementMap<GV,D,R>,
          16>(gv)
      {}
 
      static constexpr bool fixedSize()
      {
        return true;
      }
 
      static constexpr bool hasDOFs(int codim)
      {
        return codim == 1;
      }
 
      static constexpr std::size_t size(GeometryType gt)
      {
        return gt == GeometryTypes::line ? 1 : 0;
      }
 
      static constexpr std::size_t maxLocalSize()
      {
        return 4;
      }
 
    };
  } // end namespace PDELab
} // end namespace Dune
 
#endif // DUNE_PDELAB_FINITEELEMENTMAP_RT0CUBE2DFEM_HH