modeladapter.hh

#ifndef __DUNE_ACFEM_MODELADAPTER_HH__
#define __DUNE_ACFEM_MODELADAPTER_HH__
 
#include "operatorparts.hh"
#include "../modelexpressionbase.hh"
 
namespace Dune {
 
  namespace ACFem {
 
    template<class OperatorParts, class GridPart>
    class OperatorPartsAdapterModel
      : public ModelExpression<OperatorPartsAdapterModel<OperatorParts, GridPart> >
    {
      typedef OperatorPartsAdapterModel ThisType;
      typedef DefaultModel<ThisType>    BaseType;
      typedef ModelInterface<ThisType>  InterfaceType;
      typedef ModelTraits<ThisType>     TraitsType;
 
     public:
      typedef OperatorParts OperatorPartsType;
      typedef typename OperatorPartsType::FunctionSpaceType FunctionSpaceType;
      typedef GridPart GridPartType;
      typedef typename BaseType::EntityType EntityType;
 
      typedef typename FunctionSpaceType::DomainType DomainType;
      typedef typename FunctionSpaceType::RangeType RangeType;
      typedef typename FunctionSpaceType::JacobianRangeType JacobianRangeType;
      typedef typename FunctionSpaceType::HessianRangeType HessianRangeType;
 
      enum {
        dimDomain = FunctionSpaceType::dimDomain,
        dimRange = FunctionSpaceType::dimRange
      };
 
      typedef typename FunctionSpaceType::DomainFieldType DomainFieldType;
      typedef typename FunctionSpaceType::RangeFieldType RangeFieldType;
 
      OperatorPartsAdapterModel(OperatorPartsType&& parts)
        : operatorParts_(parts)
      {}
 
      std::string name() const {
        return operatorParts_.name();
      }
 
      OperatorPartsType operatorParts() const {
        return operatorParts_;
      }
 
     protected:
      OperatorPartsType operatorParts_;
    };
 
    template<class OperatorParts, class GridPart>
    struct ModelTraits<OperatorPartsAdapterModel<OperatorParts, GridPart> >
      : public DefaultModelTraits<typename OperatorParts::FunctionSpaceType, GridPart>
    {
      typedef OperatorParts OperatorPartsType;
    };
 
    template<class OperatorParts, class GridPart>
    OperatorPartsAdapterModel<OperatorParts, GridPart>
    operatorPartsAdapterModel(const OperatorParts& operatorParts, const GridPart& gridPart)
    {
      return OperatorPartsAdapterModel<OperatorParts, GridPart>(operatorParts);
    }
 
  }
}
 
#endif // __DUNE_ACFEM_MODELADAPTER_HH__