edges0.5fem.hh

// -*- tab-width: 8; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=8 sw=2 sts=2:
 
#ifndef DUNE_PDELAB_FINITEELEMENTMAP_EDGES0_5FEM_HH
#define DUNE_PDELAB_FINITEELEMENTMAP_EDGES0_5FEM_HH
 
#include <dune/localfunctions/whitney/edges0.5.hh>
 
#include <dune/pdelab/finiteelementmap/global.hh>
 
namespace Dune {
  namespace PDELab {
 
 
    template<class Geometry, class VertexOrderFactory, class RF>
    class EdgeS0_5FiniteElementMap :
      public GeometryVertexOrderFiniteElementMap<
        EdgeS0_5FiniteElementFactory<Geometry, RF>, VertexOrderFactory
        >
    {
      typedef EdgeS0_5FiniteElementFactory<Geometry, RF> FEFactory;
      typedef GeometryVertexOrderFiniteElementMap<
        FEFactory, VertexOrderFactory
        > Base;
 
      static FEFactory &feFactory() {
        static FEFactory feFactory_;
        return feFactory_;
      }
 
    public:
 
      static constexpr int dimension = Geometry::mydimension;
 
      EdgeS0_5FiniteElementMap(const VertexOrderFactory &voFactory) :
        Base(feFactory(), voFactory)
      { }
 
      static constexpr bool fixedSize()
      {
        return true;
      }
 
      static constexpr bool hasDOFs(int codim)
      {
        return Geometry::mydimension - codim == 1;
      }
 
      static constexpr std::size_t size(GeometryType gt)
      {
        return gt == GeometryTypes::line ? 1 : 0;
      }
 
      static constexpr std::size_t maxLocalSize()
      {
        // This is a bit of a hack. It only provides an upper bound which is not
        // sharp for non-cubes. However, a better compile time guess is hard to
        // obtain from latest dune-localfunctions...
        return 2 * dimension;
      }
 
    };
  }
}
 
#endif // DUNE_PDELAB_FINITEELEMENTMAP_EDGES0_5FEM_HH