File rk4.hxx

class RK4Solver : public Solver

Public Functions

RK4Solver(Options *options)

virtual void resetInternalFields() override

Should wipe out internal field vector and reset from current field object data.

virtual void setMaxTimestep(BoutReal dt) override

Set a maximum internal timestep (only for explicit schemes)

inline virtual BoutReal getCurrentTimestep() override

Return the current internal timestep.

virtual int init(int nout, BoutReal tstep) override

Initialise the solver NOTE: nout and tstep should be passed to run, not init. Needed because of how the PETSc TS code works

virtual int run() override

Run the solver, calling monitors nout times, at intervals of tstep. This function is called by solve(), and is specific to each solver type

This should probably be protected, since it shouldn’t be called by users.

Private Functions

void take_step(BoutReal curtime, BoutReal dt, Array<BoutReal> &start, Array<BoutReal> &result)

Private Members

BoutReal atol

BoutReal rtol

BoutReal max_timestep

int mxstep

Array<BoutReal> f0

Array<BoutReal> f1

Array<BoutReal> f2

BoutReal out_timestep

int nsteps

BoutReal timestep

bool adaptive

int nlocal

int neq

Array<BoutReal> k1

Array<BoutReal> k2

Array<BoutReal> k3

Array<BoutReal> k4

Array<BoutReal> k5