tvd_rk2.hh

Go to the documentation of this file.

 
#ifndef tvd_rk2_hh
#define tvd rk2_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;
 
  // class TimeVector;
 
  // ****************************************************** TVD Runge Kutta **
  class TvdRK2 : public TimeStepStrategy {
  public:
    TvdRK2(concepts::SolverFabric<Real>& fabric,
           concepts::Operator<Real>& D1,
           concepts::Operator<Real>& D0,
           TimeVector& trhs,
           const concepts::Vector<Real>& Y0,
           Real dt);
    virtual ~TvdRK2();
  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})
  };
 
  class LimitingTvdRK2 : public TimeStepStrategy {
  public:
    LimitingTvdRK2(concepts::SolverFabric<Real>& fabric,
                   concepts::Operator<Real>& D1,
                   concepts::Operator<Real>& D0,
                   concepts::Operator<Real>& limiter,
                   TimeVector& trhs,
                   const concepts::Vector<Real>& Y0,
                   Real dt);
    virtual ~LimitingTvdRK2();
  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 // tvd_rk2