simplemodel.hh

#ifndef SIMPLE_MODEL_HH
#define SIMPLE_MODEL_HH
 
#include <cassert>
#include <cmath>
 
#include <dune/acfem/models/modelinterface.hh>
 
namespace Dune {
 
namespace ACFem {
 
 
template<class FunctionSpace, class GridPart>
class FluxModel
  : public DefaultModel<FluxModel<FunctionSpace, GridPart> >
{
  typedef DefaultModel<FluxModel<FunctionSpace, GridPart> > BaseType;
  typedef ModelInterface<FluxModel<FunctionSpace, GridPart> > InterfaceType;
  typedef FluxModel ThisType;
 public:
  typedef typename BaseType::FunctionSpaceType FunctionSpaceType;
  typedef typename BaseType::GridPartType GridPartType;
  typedef typename BaseType::EntityType EntityType;
 
  typedef typename BaseType::DomainType DomainType;
  typedef typename BaseType::RangeType RangeType;
  typedef typename BaseType::JacobianRangeType JacobianRangeType;
  typedef typename BaseType::HessianRangeType HessianRangeType;
 
  typedef typename BaseType::DomainFieldType DomainFieldType;
  typedef typename BaseType::RangeFieldType RangeFieldType;
 
  FluxModel(const RangeFieldType a = RangeFieldType(0.)):a_(a)
  {}
 
  template<class Point>
  void flux(const EntityType &entity,
            const Point &x,
            const RangeType &value,
            const JacobianRangeType &gradient,
            JacobianRangeType &flux ) const
  {
    //Default is  0
    flux = JacobianRangeType(a_) ;
  }
 
 private:
  RangeFieldType a_;
};
 
template<class FunctionSpace, class GridPart>
struct ModelTraits<FluxModel<FunctionSpace, GridPart> >
  : public DefaultModelTraits<FunctionSpace, GridPart>
{
 private:
  typedef DefaultModelTraits<FunctionSpace, GridPart> BaseType;
 public:
};
 
} // namespace ACFem
 
}//Namespace Dune
 
#endif // #ifndef SIMPLE_MODEL_HH