powerfem.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_POWERFEM_HH
#define DUNE_PDELAB_FINITEELEMENTMAP_POWERFEM_HH
 
#include <cstddef>
 
#include <dune/localfunctions/meta/power.hh>
 
#include <dune/pdelab/finiteelementmap/finiteelementmap.hh>
 
namespace Dune {
  namespace PDELab {
 
 
    template<class BackendFEM, std::size_t dimR>
    class PowerFiniteElementMap
    {
      typedef typename BackendFEM::Traits::FiniteElementType BackendFE;
      typedef PowerFiniteElementFactory<BackendFE, dimR> Factory;
 
      const BackendFEM &backend;
      static const Factory factory;
 
    public:
      typedef FiniteElementMapTraits<typename Factory::FiniteElement> Traits;
 
      // This is not backwards compatible, but I don't know of anyone who is using
      // this thing, so it's probably not much of a problem.
      static constexpr int dimension = BackendFEM::dimension;
 
 
      PowerFiniteElementMap(const BackendFEM& backend_) : backend(backend_) { }
 
 
      template<class Element>
      typename Traits::FiniteElementType find(const Element& e) const
      { return factory.make(backend.find(e)); }
 
      bool fixedSize() const
      {
        return backend.fixedSize();
      }
 
      bool hasDOFs(int codim) const
      {
        return backend.hasDOFs(codim);
      }
 
      std::size_t size(GeometryType gt) const
      {
        return dimR * backend.size(gt);
      }
 
      std::size_t maxLocalSize() const
      {
        return dimR * backend.maxLocalSize();
      }
 
    };
 
    template<class BackendFEM, std::size_t dimR>
    const typename PowerFiniteElementMap<BackendFEM, dimR>::Factory
    PowerFiniteElementMap<BackendFEM, dimR>::factory = Factory();
 
  } // namespace PDELab
} // namespace Dune
 
#endif // DUNE_PDELAB_FINITEELEMENTMAP_POWERFEM_HH