timestepparameter.hh

Go to the documentation of this file.

 
#ifndef __DUNE_ACFEM_TIMESTEP_PARAMETER_HH__
#define __DUNE_ACFEM_TIMESTEP_PARAMETER_HH__
 
#include <dune/fem/solver/timeprovider.hh>
#include "../expressions/parameterinterface.hh"
#include "timeview.hh"
 
namespace Dune {
 
  namespace ACFem {
 
    template<class TimeProvider>
    class TimeStepParameter
      : public DefaultParameter<TimeStepParameter<TimeProvider> >
    {
      typedef TimeStepParameter            ThisType;
     public:
      typedef typename ParameterSpaceTraits<double>::ParameterSpaceType ParameterSpaceType;
      typedef typename ParameterSpaceType::ScalarParameterSpaceType ScalarParameterSpaceType;
      enum { dimension = ParameterSpaceType::dimension };
      typedef typename ParameterSpaceType::FieldType FieldType;
      typedef typename ParameterSpaceType::ValueType ValueType;
      typedef TimeProvider TimeProviderType;
 
      TimeStepParameter(const TimeProviderType& tp) : timeProvider_(tp) {}
 
      ValueType value() const { return timeProvider_.deltaT(); }
 
     protected:
      const TimeProviderType& timeProvider_;
    };
 
    template<class TimeProvider>
    struct ParameterTraits<TimeStepParameter<TimeProvider> >
    {
      typedef typename ParameterSpaceTraits<double>::ParameterSpaceType ParameterSpaceType;
    };
 
    template<class TimeProvider>
    class InverseTimeStepParameter
      : public DefaultParameter<InverseTimeStepParameter<TimeProvider> >
    {
      typedef InverseTimeStepParameter     ThisType;
     public:
      typedef typename ParameterSpaceTraits<double>::ParameterSpaceType ParameterSpaceType;
      typedef typename ParameterSpaceType::ScalarParameterSpaceType ScalarParameterSpaceType;
      enum { dimension = ParameterSpaceType::dimension };
      typedef typename ParameterSpaceType::FieldType FieldType;
      typedef typename ParameterSpaceType::ValueType ValueType;
      typedef TimeProvider TimeProviderType;
 
      InverseTimeStepParameter(const TimeProviderType& tp) : timeProvider_(tp) {}
 
      ValueType value() const {
        return timeProvider_.inverseDeltaT();
      }
 
     protected:
      const TimeProviderType& timeProvider_;
    };
 
    template<class TimeProvider>
    struct ParameterTraits<InverseTimeStepParameter<TimeProvider> >
    {
      typedef typename ParameterSpaceTraits<double>::ParameterSpaceType ParameterSpaceType;
    };
 
    template<class TimeProvider>
    class TimeViewParameter
      : public DefaultParameter<TimeViewParameter<TimeProvider> >
    {
      typedef TimeViewParameter ThisType;
     public:
      typedef typename ParameterSpaceTraits<double>::ParameterSpaceType ParameterSpaceType;
      typedef typename ParameterSpaceType::ScalarParameterSpaceType ScalarParameterSpaceType;
      enum { dimension = ParameterSpaceType::dimension };
      typedef typename ParameterSpaceType::FieldType FieldType;
      typedef typename ParameterSpaceType::ValueType ValueType;
      typedef TimeProvider TimeProviderType;
      typedef TimeView<TimeProviderType> TimeViewType;
 
      TimeViewParameter(const TimeProviderType& tp, double theta = 0.0)
        : timeView_(tp, theta)
      {}
 
      ValueType value() const {
        return timeView().time();
      }
 
      TimeViewType& timeView()
      {
        return timeView_;
      }
 
      const TimeViewType& timeView() const
      {
        return timeView_;
      }
     protected:
      TimeViewType timeView_;
    };
 
    template<class TimeProvider>
    struct ParameterTraits<TimeViewParameter<TimeProvider> >
    {
      typedef typename ParameterSpaceTraits<double>::ParameterSpaceType ParameterSpaceType;
    };
 
 
 
  }
 
}
 
#endif // __DUNE_ACFEM_TIMESTEP_PARAMETER_HH__