RK2.h

#ifndef gridripper_amr1d_gridint_RK2_h
#define gridripper_amr1d_gridint_RK2_h

#include <gridripper/amr1d/Integrator.h>
#include <gridripper/amr1d/Dissipation.h>

namespace gridripper { namespace amr1d {

class Grid;
class Grad;
class Sigma;

namespace gridint {

using namespace std;

class RK2: public Integrator
{
private:
    Grid* tmpg;
    Grid* resultGrid;
    Grid* richardsonResultGrid1;
    Grid* richardsonResultGrid2;
    O4Dissipation dissipationOp;

public:
    RK2(PDE& pde, const tvalarray<string>* opts, Sigma* sigma);

protected:
    RK2(const RK2& gi, Grid& g);

private:
    RK2(Grid* g, PDE& pde, Grad* d, Grad* d2, int deftype, Sigma* sigma);

public:
    ~RK2();

    int getNumVelocity() const {
        return 4;
    }

    void integrate(int ir);

    void richardson(const AMError& amerr, int ir);

    Integrator* createIntegrator(Grid& g, const Sigma* s) const;

    Integrator* cloneIntegrator(Grid& g) const;

    void getResult(int i, tvalarray<GReal_t>& w) const;

private:
    void integrate(const Grid& g0, Grid& g1, Grid& g2,
                   Grad& grad, Dissipation* dissip, GReal_t dt, int ir);

    void step(GReal_t t, GReal_t dt, const Grid& g0, const Grid& gin,
              Grid& gout, Grad& grad, Dissipation* dissip, int ir);
};

} } } // namespace gridripper::amr1d::gridint

#endif /* gridripper_amr1d_gridint_RK2_h */