localfunctionsetadapter.hh

#ifndef DUNE_FEM_FUNCTION_LOCALFUNCTION_LOCALFUNCTIONSETADAPTER_HH
#define DUNE_FEM_FUNCTION_LOCALFUNCTION_LOCALFUNCTIONSETADAPTER_HH
 
#include <cassert>
#include <cstddef>
 
#include <dune/fem/common/coordinate.hh>
 
namespace Dune
{
 
  namespace Fem
  {
 
    // LocalFunctionSetAdapter
    // -----------------------
 
    template< class Entity, class FunctionSet >
    struct LocalFunctionSetAdapter
    {
      typedef Entity EntityType;
      typedef FunctionSet FunctionSetType;
 
      typedef typename FunctionSet::FunctionSpaceType FunctionSpaceType;
 
      typedef typename FunctionSpaceType::DomainType DomainType;
      typedef typename FunctionSpaceType::RangeType RangeType;
      typedef typename FunctionSpaceType::JacobianRangeType JacobianRangeType;
      typedef typename FunctionSpaceType::HessianRangeType HessianRangeType;
 
      explicit LocalFunctionSetAdapter ( const FunctionSet &functionSet = FunctionSet() )
        : functionSet_( functionSet )
      {}
 
      explicit LocalFunctionSetAdapter ( const Entity &entity, const FunctionSet &functionSet = FunctionSet() )
        : entity_( &entity ), functionSet_( functionSet )
      {}
 
      int order () const { return functionSet_.order(); }
 
      const EntityType &entity () const
      {
        assert( entity_ );
        return *entity_;
      }
 
      std::size_t size () const { return functionSet().size(); }
 
      template< class Point, class Functor >
      void evaluateEach ( const Point &x, Functor functor ) const
      {
        DomainType y = geometry().global( coordinate( x ) );
        functionSet().evaluateEach( y, functor );
      }
 
      template< class Point, class Functor >
      void jacobianEach ( const Point &x, Functor functor ) const
      {
        DomainType y = geometry().global( coordinate( x ) );
        functionSet().jacobianEach( y, functor );
      }
 
      template< class Point, class Functor >
      void hessianEach ( const Point &x, Functor functor ) const
      {
        DomainType y = geometry().global( coordinate( x ) );
        functionSet().hessianEach( y, functor );
      }
 
 
      // Non-interface methods
      // ---------------------
 
      const FunctionSetType functionSet () const { return functionSet_; }
 
    private:
      typename EntityType::Geometry geometry () const { return entity().geometry(); }
 
      const EntityType *entity_ = nullptr;
      FunctionSetType functionSet_;
    };
 
  } // namespace Fem
 
} // namespace Dune
 
#endif // #ifndef DUNE_FEM_FUNCTION_LOCALFUNCTION_LOCALFUNCTIONSETADAPTER_HH