Runge-Kutta Finite Differences which insert the density floor after each stage. More...

Inheritance diagram for 20_additional-solver-steps.RungeKuttaWithDensityFloor:

[legend]

Collaboration diagram for 20_additional-solver-steps.RungeKuttaWithDensityFloor:

[legend]

Public Member Functions
	construct_output (self, exahype2_project, peano_project, time_step)
	Add the data to the Peano 4 project, generate the kernels and add the tasks of one Runge-Kutta step to time_step.

Public Member Functions inherited from exahype2.api.rkfd4.CentralDifferencesWithRungeKutta.CentralDifferencesWithRungeKutta
	__init__ (self, Solver solver, int patch_size, int rk_order, float time_step_size=None, float cfl_constant=0.5, "CentralDifferencesWithRungeKutta.RightHandSideKernel" right_hand_side_kernel=RightHandSideKernel.CentralDifferences, advection_velocity=None, advected_quantities=None, float min_volume_size=None, float max_volume_size=None, "toolbox.blockstructured.api.InterpolateHangingFaces.Scheme" interpolation=toolbox.blockstructured.api.InterpolateHangingFaces.Scheme.Linear, "toolbox.blockstructured.api.RestrictHangingFaces.Scheme" restriction=toolbox.blockstructured.api.RestrictHangingFaces.Scheme.Averaging)
	The constructor memorises the attributes and creates the data structures.

int	overlap (self)
	Width of the halo, i.e.

int	rk_order (self)

bool	has_adaptive_time_step_size (self)

str	right_hand_side_kernel_name (self)
	Name of the generated right-hand side kernel.

str	max_eigenvalue_kernel_name (self)

Solver.MeshAdaptation	create_default_mesh_adaptation (self)
	The finite differences adapt the mesh and balance the load through exahype2.api.rkfd4.AdaptMeshAndBalanceLoad by default.

"toolbox.blockstructured.api.RestrictHangingFaces"	create_restriction_of_hanging_faces (self)
	Action set which restricts the projected fine data from the hanging faces onto the coarse faces.

Protected Attributes
	_solver

Protected Attributes inherited from exahype2.api.rkfd4.CentralDifferencesWithRungeKutta.CentralDifferencesWithRungeKutta
	_solver
	Project the outcome onto the faces.

Protected Attributes inherited from exahype2.api.Solver.Solver.NumericalScheme
	_solver

Additional Inherited Members
Data Fields inherited from exahype2.api.rkfd4.CentralDifferencesWithRungeKutta.CentralDifferencesWithRungeKutta
	interpolation

	restriction

	right_hand_side_kernel

	advection_velocity

	advected_quantities

	HALO_SIZE

	patch_size

	min_volume_size

	max_volume_size

	butcher_tableau

	time_step_size

	cfl_constant

	kernel_namespace

	max_eigenvalue_kernel_call

	right_hand_side_kernel_name

	max_eigenvalue_kernel_name

	rk_order

Static Public Attributes inherited from exahype2.api.rkfd4.CentralDifferencesWithRungeKutta.CentralDifferencesWithRungeKutta
int	HALO_SIZE = 2
	Width of the halo.

Protected Member Functions inherited from exahype2.api.rkfd4.CentralDifferencesWithRungeKutta.CentralDifferencesWithRungeKutta
	_right_hand_side_action_sets (self, stage, cell_guard)
	Action sets of the task which evaluates the right-hand side of one stage.

	_patches (self)
	All patches of this numerical scheme, i.e.

Detailed Description

Runge-Kutta Finite Differences which insert the density floor after each stage.

Definition at line 148 of file 20_additional-solver-steps.py.

Member Function Documentation

◆ construct_output()

20_additional-solver-steps.RungeKuttaWithDensityFloor.construct_output	(	self,
		exahype2_project,
		peano_project,
		time_step
	)

Add the data to the Peano 4 project, generate the kernels and add the tasks of one Runge-Kutta step to time_step.

Per stage k, the time step holds the following tasks (see the class documentation for the data flow):

{SolverName}_Apply Dirichlet BCs stage k: sets the ghost grid cells along the domain boundary for Dirichlet boundary conditions. The task does not exist if all axes are periodic, or for Sommerfeld boundary conditions, which advance the ghost grid cells within the roll over of the projection onto the faces.
{SolverName}_Reconstruct stage k: reconstructs the stage's input with its halo.
{SolverName}_Evaluate right hand side stage k.
{SolverName}_Prepare stage k+1, or {SolverName}_Complete time step after the last stage: computes the linear combination of the estimates. Prepare stage k+1 is tagged with exahype2.api.Solver.IntermediateSolutionUpdateTag, Complete time step with exahype2.api.Solver.SolutionUpdateTaskTag, i.e. postprocessing such as the enforcement of constraints can be inserted directly after them (see Manipulate an existing solver's steps (e.g. inserting additional sweeps)).
{SolverName}_Project stage k+1 onto faces, or {SolverName}_Project solution onto faces after the last stage: projects the outcome of the linear combination onto the faces. After the last stage, it also computes the maximum eigenvalue of the new solution.

After the last stage, {SolverName}_Compute time step size determines the time step size of the next time step and advances the solver's time stamp. {SolverName}_Dynamic mesh refinement finally evaluates the adaptivity criterion if the solver has one (exahype2.api.blockstructured.create_dynamic_mesh_refinement_action_set.create_dynamic_mesh_refinement_action_set()), switches the load balancing on again and informs it that the time step has finished.

Reimplemented from exahype2.api.rkfd4.CentralDifferencesWithRungeKutta.CentralDifferencesWithRungeKutta.

Definition at line 151 of file 20_additional-solver-steps.py.

References 20_additional-solver-steps.RungeKuttaWithDensityFloor.construct_output(), and 20_additional-solver-steps.insert_density_floor().

Referenced by 20_additional-solver-steps.RungeKuttaWithDensityFloor.construct_output().

Here is the call graph for this function: