RK2.h

Go to the documentation of this file.

#ifndef RK2_H
#define RK2_H

#include "timeInt.h"
#include "hybrid.h"


class RK2 : public TimeIntegrator
{
  public:

    // Need some work arrays
    double *p1cons, *p1prims, *p1aux, *args1, *args2;


    RK2(Data * data, Model * model, Bcs * bcs, FluxMethod * fluxMethod, ModelExtension * modelExtension = NULL);

    virtual ~RK2();


    void predictorStep(double * cons, double * prims, double * aux, double dt);

    // Corrector
    void correctorStep(double * cons, double * prims, double * aux, double dt);


    void step(double * cons, double * prims, double * aux, double dt=0);

};


#endif