File rk4.hxx

class RK4Solver : public Solver

Public Functions

RK4Solver(Options *options)
void resetInternalFields()

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

void setMaxTimestep(BoutReal dt)

Set a maximum internal timestep (only for explicit schemes)

BoutReal getCurrentTimestep()

Return the current internal timestep.

int init(int nout, BoutReal tstep)

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

int run()

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