#include <RKPlus.h>

Inheritance diagram for RK2B:

Collaboration diagram for RK2B:

[legend]

Public Member Functions
	RK2B (Data data, Model model, Bcs bcs, FluxMethod fluxMethod, ModelExtension *modelExtension=NULL)

virtual	~RK2B ()

void	stage1 (double cons, double prims, double *aux, double dt)

void	stage2 (double cons, double prims, double *aux, double dt)

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

Public Member Functions inherited from RKPlus
	RKPlus (Data data, Model model, Bcs bcs, FluxMethod fluxMethod, ModelExtension *modelExtension=NULL)

virtual	~RKPlus ()

virtual void	rhs (double cons, double prims, double aux, double rhsVec)

Public Member Functions inherited from TimeIntegrator
	TimeIntegrator (Data data, Model model, Bcs bcs, FluxMethod fluxMethod, ModelExtension *modelExtension=NULL)

virtual	~TimeIntegrator ()

void	finalise (double cons, double prims, double *aux)

Public Member Functions inherited from TimeIntegratorBase
	TimeIntegratorBase (Data data, Model model, Bcs bcs, FluxMethod fluxMethod, ModelExtension *modelExtension=NULL)

virtual	~TimeIntegratorBase ()

Public Attributes

double *	u1cons

double *	u1prims

double *	u1aux

double *	rhs1

double *	rhs2

Public Attributes inherited from RKPlus
double *	fluxCont

Public Attributes inherited from TimeIntegratorBase
Data *	data

Model *	model

Bcs *	bcs

FluxMethod *	fluxMethod

ModelExtension *	modelExtension

Detailed Description

Second order RK.

: SSPRK(2,2) from Gottlieb 2009. This can run the AdvectionSingleFluid test with cfl=0.6

The step function performs the following stages. Stage 1:

\begin{align} U^{(1)} = U^n + \Delta t \mathcal{L}(U^n). \end{align}

Stage 2:

\begin{align} U^{n+1} = \frac{1}{2} U^n + \frac{1}{2} U^{(1)} + \frac{1}{2} \Delta t \mathcal{L}(U^{(1)}) \end{align}

Definition at line 97 of file RKPlus.h.

Constructor & Destructor Documentation

◆ RK2B()

RK2B::RK2B	(	Data *	data,
		Model *	model,
		Bcs *	bcs,
		FluxMethod *	fluxMethod,
		ModelExtension *	modelExtension = `NULL`
	)

Constructor.

Constructor requires simulation data and the flux and source functions from the model class. Stores the necessary pointer.

Parameters

[in]	*cons	pointer to conserved vector work array. Size is \(N_{cons} \times N_x \times N_y \times N_z\)
[in]	*prims	pointer to primitive vector work array. Size is \(N_{prims} \times N_x \times N_y \times N_z\)
[in]	*aux	pointer to auxiliary vector work array. Size is \(N_{aux} \times N_x \times N_y \times N_z\)
	dt	the step size desired to move by. Defaults to the value in the Data class

See also: TimeIntegrator::TimeIntegrator; RKPlus::RKPlus

◆ ~RK2B()

virtual RK2B::~RK2B ( )

virtual

Member Function Documentation

◆ stage1()

void RK2B::stage1	(	double *	cons,
		double *	prims,
		double *	aux,
		double	dt
	)

Stage 1 result.

Compute and store stage 1 result.

Parameters

[in]	*cons	pointer to conserved vector work array. Size is \(N_{cons} \times N_x \times N_y \times N_z\)
[in]	*prims	pointer to primitive vector work array. Size is \(N_{prims} \times N_x \times N_y \times N_z\)
[in]	*aux	pointer to auxiliary vector work array. Size is \(N_{aux} \times N_x \times N_y \times N_z\)
	dt	the step size desired to move by. Defaults to the value in the Data class

See also: RK2B::RK2B

◆ stage2()

void RK2B::stage2	(	double *	cons,
		double *	prims,
		double *	aux,
		double	dt
	)

Stage 2 result.

Compute and store stage 2 result.

Parameters

[in]	*cons	pointer to conserved vector work array. Size is \(N_{cons} \times N_x \times N_y \times N_z\)
[in]	*prims	pointer to primitive vector work array. Size is \(N_{prims} \times N_x \times N_y \times N_z\)
[in]	*aux	pointer to auxiliary vector work array. Size is \(N_{aux} \times N_x \times N_y \times N_z\)
	dt	the step size desired to move by. Defaults to the value in the Data class

See also: RK2B::RK2B

◆ step()

void RK2B::step	(	double *	cons,
		double *	prims,
		double *	aux,
		double	dt = `0`
	)

virtual

Performs a single time step.

The timestep will use the current values of the conserved, primitive and auxiliary variables at t=t0 and compute the values of all of them at time \(t=t_0 + dt\). I.e. the conserved vector is evolved forward, and the corresponding prims and aux vars are found.

Parameters

[in]	*cons	pointer to conserved vector work array. Size is \(N_{cons} \times N_x \times N_y \times N_z\)
[in]	*prims	pointer to primitive vector work array. Size is \(N_{prims} \times N_x \times N_y \times N_z\)
[in]	*aux	pointer to auxiliary vector work array. Size is \(N_{aux} \times N_x \times N_y \times N_z\)
	dt	the step size desired to move by. Defaults to the value in the Data class

See also: TimeIntegrator::step

Implements TimeIntegratorBase.