O4Relax.h

#ifndef gridripper_odesolver_relax_method_O4Relax_h
#define gridripper_odesolver_relax_method_O4Relax_h

#include <gridripper/odesolver/relax/ODERelax.h>

namespace gridripper { namespace odesolver {

class ODE;

namespace relax {

class O4Relax: public ODERelax::Method
{
private:
    ODE* theODE;

    int maxi;

public:
    O4Relax(ODE* ode, int nx) throw(ODERelax::RelaxationException&);

    int getNumSubsteps() const;

    void set(tvalarray<GReal_t>& fieldValues, tvalarray<GReal_t>& matrix,
             tvalarray<GReal_t>& rhsvec,
             const tvalarray<GReal_t>& deltaField);
};

} } } // namespace gridripper::odesolver::relax

#endif /* gridripper_odesolver_relax_method_O4Relax_h */