euler.hh

Go to the documentation of this file.

 
#ifndef euler_hh
#define euler_hh
 
#include "basics/typedefs.hh"
#include "function/vector.hh"
#include "timestepping/strategy.hh"
#include "timestepping/vectors.hh"
 
namespace concepts {
  template<typename F>
  class Operator;
}
 
namespace timestepping {
 
  using concepts::Real;
 
  // ****************************************************************** Euler **
 
  class Euler : public TimeStepStrategy {
  public:
    Euler(concepts::SolverFabric<Real>& fabric,
          concepts::Operator<Real>& D1,
          concepts::Operator<Real>& D0,
          timestepping::TimeVector& trhs,
          const concepts::Vector<Real>& Y0,
          Real dt);
    virtual ~Euler();
  protected:
    virtual std::ostream& info(std::ostream& os) const;
    virtual void next();
  private:
    concepts::Operator<Real> &D0_;
    TimeVector &trhs_;
    TimeVector Yn1_;    // Yn1 = Y(t_{n-1})
  };
 
  // ********************************************************** LimitingEuler **
 
  class LimitingEuler : public TimeStepStrategy {
  public:
    LimitingEuler(concepts::SolverFabric<Real>& fabric,
                  concepts::Operator<Real>& D1,
                  concepts::Operator<Real>& D0,
                  concepts::Operator<Real>& limiter,
                  timestepping::TimeVector& trhs,
                  const concepts::Vector<Real>& Y0,
                  Real dt);
    virtual ~LimitingEuler();
  protected:
    virtual std::ostream& info(std::ostream& os) const;
    virtual void next();
  private:
    concepts::Operator<Real> &D0_, &limiter_;
    TimeVector &trhs_;
    TimeVector Yn1_;    // Yn1 = Y(t_{n-1})
  };
 
}
 
#endif // euler_hh