File rk4.hxx

File rk4.hxx#

class RK4Solver : public Solver #

Public Functions

explicit RK4Solver(Options *opts = nullptr)#

virtual void resetInternalFields() override#: Should wipe out internal field vector and reset from current field object data.

void setMaxTimestep(BoutReal dt) override#

inline virtual BoutReal getCurrentTimestep() override#: Return the current internal timestep.

virtual int init() override#: Initialise the solver.

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)#: Take a single step to calculate f1.

Private Members

BoutReal atol#

BoutReal rtol#

BoutReal max_timestep#

BoutReal timestep#

int mxstep#

bool adaptive#

Array<BoutReal> f0#

Array<BoutReal> f1#

Array<BoutReal> f2#

int nlocal#

int neq#

Array<BoutReal> k1#

Array<BoutReal> k2#

Array<BoutReal> k3#

Array<BoutReal> k4#

Array<BoutReal> k5#