ICN.h

#ifndef gridripper_amr1d_gridint_ICN_h
#define gridripper_amr1d_gridint_ICN_h

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

namespace gridripper { namespace amr1d {

class Grid;
class Grad;

namespace gridint {

using namespace std;

class ICN: public Integrator
{
private:
    Grid* tmpg1;
    Grid* tmpg2;
    Grid* resultGrid;
    Grid* richardsonResultGrid1;
    Grid* richardsonResultGrid2;
    O4Dissipation dissipationOp;

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

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

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

public:
    ~ICN();

    int getNumVelocity() const {
        return 6;
    }

    Grad& getGrad() const {
        return *gradOp;
    }

    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, Grid& g3,
                   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_ICN_h */