global.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_GLOBAL_HH
#define DUNE_PDELAB_FINITEELEMENTMAP_GLOBAL_HH
 
#include <dune/pdelab/finiteelementmap/finiteelementmap.hh>
 
namespace Dune {
  namespace PDELab {
 
 
    template<class Factory>
    class GeometryFiniteElementMap
    {
      Factory& factory;
 
    public:
      typedef FiniteElementMapTraits<typename Factory::FiniteElement> Traits;
 
 
      GeometryFiniteElementMap(Factory& factory_) : factory(factory_) {}
 
 
      template<class Element>
      typename Traits::FiniteElementType find(const Element& e) const {
        return factory.make(e.geometry());
      }
    };
 
 
    template<class FEFactory, class VOFactory>
    class GeometryVertexOrderFiniteElementMap
    {
      FEFactory& feFactory;
      const VOFactory& voFactory;
 
    public:
      typedef FiniteElementMapTraits<typename FEFactory::FiniteElement> Traits;
 
 
      GeometryVertexOrderFiniteElementMap(FEFactory& feFactory_,
                                          const VOFactory & voFactory_) :
        feFactory(feFactory_), voFactory(voFactory_)
      {}
 
 
      template<class Element>
      typename Traits::FiniteElementType find(const Element& e) const {
        return feFactory.make(e.geometry(), voFactory.make(e));
      }
    };
 
  } // namespace PDELab
} // namespace Dune
 
#endif // DUNE_PDELAB_FINITEELEMENTMAP_GLOBAL_HH