adam_flume_cpu_object

ADAM, FLUME CPU backend object.

Implements the forest contract on the host (MPI + OpenMP). The dispatch procedure pointers are type components, bound in initialize_flume by exhaustive select case with a fatal case default (issue #35, D-10). Space operator: characteristic (or conservative) WENO flux splitting, face fluxes then flux difference (issue #35, section 3.4); the per-model loops live in adam_flume_cpu_<model>_kernels, selected here by select case on the physical model, never inside a loop (issue #41, section 4).

Source: src/app/flume/cpu/adam_flume_cpu_object.F90

Dependencies

Contents

• flume_cpu_object
• accumulate_seam_fluxes
• apply_floors
• allocate_cpu
• check_glm_ch
• compute_conservation
• compute_divb_history
• compute_q_aux
• initialize_flume
• save_residuals
• save_simulation_data
• set_boundary_conditions
• set_initial_conditions
• update_ghost
• advance_one_step_forest
• apply_reflux_to_stage_forest
• begin_stage_forest
• close_step_forest
• compute_local_dt_forest
• end_stage_forest
• fill_seam_from_peer_forest
• finalize_forest
• initialize_forest
• is_done_forest
• open_step_forest
• post_step_forest
• compute_face_mirror_indexes
• compute_flux_difference
• compute_rk_ssp_residual
• compute_residuals_weno
• integrate_rk_ls
• integrate_rk_ssp
• stages_per_step_forest

Derived Types

flume_cpu_object

FLUME CPU backend object.

Inheritance

Extends: flume_common_object

Components

Name	Type	Attributes	Description
io	type(io_object)		IO handler.
amr	type(amr_object)		AMR marker handler.
slices	type(slices_object)		Slices handler.
blanesmoan	type(blanesmoan_object)		Blanes-Moan integrator.
cfm	type(cfm_object)		Commutator-Free Magnus integrator.
leapfrog	type(leapfrog_object)		Leapfrog integrator.
flail	type(flail_object)		Linear algebra methods handler.
weno	type(weno_object)		WENO reconstructor.
ib	type(ib_object)		Immersed boundary.
rk	type(rk_object)		Runge-Kutta integrator.
adam	type(adam_object)		ADAM, grid + tree + field + maps orchestrator.
fdv_scheme	character(len=:)	allocatable	FDV scheme, fd/fv.
fdv_order	integer(kind=I4P)		Order of finite difference/volume schemes, general order.
fdv_half_stencil	integer(kind=I4P)		Half stencil length of finite difference/volume schemes.
fdv_half_stencils	integer(kind=I4P)		Half stencil length of fdv schemes for each derivative up to 6.
ngc	integer(kind=I4P)	pointer	Number of ghost cells.
ni	integer(kind=I4P)	pointer	Number of cells in i direction.
nj	integer(kind=I4P)	pointer	Number of cells in j direction.
nk	integer(kind=I4P)	pointer	Number of cells in k direction.
nb	integer(kind=I4P)	pointer	Total blocks number for MPI.
blocks_number	integer(kind=I4P)	pointer	Actual blocks number.
nv	integer(kind=I4P)	pointer	Number of variables in q vector.
realm_index	integer(kind=I4P)		Inter-realm seam coupling — currently-published stage buffer key.
stage_active	integer(kind=I4P)		Procedure pointer TBPs for FDV operators (set at initialization by backend).
compute_block_total_variation	procedure(compute_block_total_variation_interface)	pass(self), pointer	Compute TV.
compute_curl	procedure(compute_curl_interface)	pass(self), pointer	Compute curl.
compute_curl_extended	procedure(compute_curl_interface)	pass(self), pointer	Compute curl on extended support.
compute_derivative1	procedure(compute_derivative1_interface)	pass(self), pointer	Compute deriv1.
compute_derivative2	procedure(compute_derivative2_interface)	pass(self), pointer	Compute deriv2.
compute_derivative4	procedure(compute_derivative4_interface)	pass(self), pointer	Compute deriv4.
compute_divergence	procedure(compute_divergence_interface)	pass(self), pointer	Compute dive.
compute_gradient	procedure(compute_gradient_interface)	pass(self), pointer	Compute grad.
compute_gradient_extended	procedure(compute_gradient_interface)	pass(self), pointer	Compute gradient on extended support.
compute_laplacian	procedure(compute_laplacian_interface)	pass(self), pointer	Compute laplac.
amr_locked_	logical		Runtime AMR locked after initialization.
save_auxiliary_fields	logical		Save the auxiliary variables with the fields.
glm_speed_reported	real(kind=R8P)		Largest wave speed above c_h reported so far.
divb_reported	real(kind=R8P)		Largest max|div B| above divb_tol reported so far.
divb_seam	integer(kind=I4P)	allocatable	Seam faces flags of the div(B) history [nb, 6].
q	real(kind=R8P)	allocatable	Conservative variables [nv, 1-ngc:ni+ngc, ..., nb].
dq	real(kind=R8P)	allocatable	Residuals [nv, 1-ngc:ni+ngc, ..., nb].
q_aux	real(kind=R8P)	allocatable	Auxiliary variables [nv_aux, 1-ngc:ni+ngc, ..., nb].
q_name	type(string)	allocatable	Conservative variables names.
dq_name	type(string)	allocatable	Residuals names.
q_aux_name	type(string)	allocatable	Auxiliary variables names.
bc	type(flume_bc_object)		Boundary conditions.
diagnostics	type(flume_diagnostics_object)		Diagnostics.
ic	type(flume_ic_object)		Initial conditions.
numerics	type(flume_numerics_object)		Numerics.
physics	type(flume_physics_object)		Physics.
time	type(flume_time_object)		Time handler.
flx_f	real(kind=R8P)	allocatable	X-face fluxes [nv, 0:ni, 1:nj, 1:nk, nb], face i = i+1/2.
fly_f	real(kind=R8P)	allocatable	Y-face fluxes [nv, 1:ni, 0:nj, 1:nk, nb], face j = j+1/2.
flz_f	real(kind=R8P)	allocatable	Z-face fluxes [nv, 1:ni, 1:nj, 0:nk, nb], face k = k+1/2.
compute_residuals	procedure(compute_residuals_interface)	pass(self), pointer	Space operator.
integrate	procedure(integrate_interface)	pass(self), pointer	Time operator.

Type-Bound Procedures

Name	Attributes	Description
load_fdv_from_file	pass(self)	Load FDV config from file.
finalize_mpi_forest	pass(self)	Process-global MPI finalize; forest calls it ONCE after all.
after_topology_build_forest	pass(self)	Backend hook invoked once after the forest builds the seam maps.
close_block_xh5f	nopass	Close XH5F file block.
close_file_xh5f	nopass	Close XH5F file.
open_block_xh5f	pass(self)	Open block file XH5F.
open_file_xh5f	pass(self)	Open file XH5F.
save_q_xh5f	pass(self)	Save in XH5F (XDMF/HDF5) format.
compute_curl_fd_extended	pass(self)	Compute curl of vector field, finite difference, even for half gcs.
compute_gradient_fd_extended	pass(self)	Compute gradient of scalar field, finite difference, even for half
amr_update	pass(self)	Do AMR update (initialization-time only).
mark_by_geometry	pass(self)	Mark blocks to be refined by a primitive geometric box.
mark_by_gradient	pass(self)	Mark blocks by the gradient of a conservative or auxiliary variable.
mark_by_solid	pass(self)	Mark blocks crossed by the surface of an immersed solid.
accumulate_seam_skin	pass(self)	Route one weighted seam face skin to the forest's flux register.
allocate_common	pass(self)	Allocate common data.
compute_fields_number	pass(self)	Compute the block-sized fields allocated per block.
compute_phi	pass(self)	Compute the immersed solids distance function (host).
destroy_common	pass(self)	Free common data.
glm_lambda	pass(self)	Return the GLM bound of the local dt, c_h max sum_d 1/dx_d.
initialize	pass(self)	Initialize the common data.
null_freeze	pass(self)	Return the variable each null direction freezes.
report_divb	pass(self)	Reduce and save the div(B) norms, apply the divb_tol monitor.
report_glm_speed	pass(self)	Check c_h against the fastest wave (warning or stop).
load_restart_files	pass(self)	Load restart files.
save_restart_files	pass(self)	Save restart files.
save_slices	pass(self)	Save the slices on their cadence.
save_xh5f	pass(self)	Save fields in XH5F format.
set_divb_seam	pass(self)	Set the seam faces flags of the div(B) history.
set_glm_damping	pass(self)	Set the GLM damping once the grid exists.
coupling_descriptor_forest	pass(self)	Return the realm coupling descriptor.
accumulate_seam_fluxes	pass(self)	Accumulate the weighted seam face fluxes of one stage.
apply_floors	pass(self)	Apply the MHD positivity floors of a stage.
allocate_cpu	pass(self)	Allocate CPU data.
check_glm_ch	pass(self)	Check the GLM c_h against the fastest wave.
compute_conservation	pass(self)	Compute and save the conservation integrals.
compute_divb_history	pass(self)	Compute and save the div(B) norms (MHD).
compute_q_aux	pass(self)	Compute the auxiliary variables.
initialize_flume	pass(self)	Initialize the CPU backend.
save_residuals	pass(self)	Save residuals history.
save_simulation_data	pass(self)	Save fields, restart and diagnostics on their cadence.
set_boundary_conditions	pass(self)	Set boundary conditions on the crown maps.
set_initial_conditions	pass(self)	Set initial conditions.
update_ghost	pass(self)	Update ghost cells: local, MPI, boundary conditions.
advance_one_step_forest	pass(self)	Advance one full step (fast path).
apply_reflux_to_stage_forest	pass(self)	Apply the reflux correction.
begin_stage_forest	pass(self)	Begin an integrator stage (staged path).
close_step_forest	pass(self)	Close a step (staged path).
compute_local_dt_forest	pass(self)	Compute the local stability-limited time step.
end_stage_forest	pass(self)	End an integrator stage (staged path).
fill_seam_from_peer_forest	pass(self)	Fill inter-realm seam ghosts from a peer.
finalize_forest	pass(self)	Finalize the realm.
initialize_forest	pass(self)	Initialize the realm.
is_done_forest	pass(self)	Return true if the realm is done.
open_step_forest	pass(self)	Open a step (staged path).
post_step_forest	pass(self)	Post-step IO and diagnostics.
stages_per_step_forest	pass(self)	Return the integrator stages per step.

Subroutines

accumulate_seam_fluxes

Accumulate the seam face fluxes of stage s, weighted by its SSP coefficient, into the forest's flux register.

Every stage contributes (issue #35, P5): the register holds sum_s beta_s F_s, the flux of the committed step.

subroutine accumulate_seam_fluxes(self, s, flux_register)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.
s	integer(kind=I4P)	in		Runge-Kutta stage.
flux_register	class(flux_register_object)	inout		Forest's flux register.

Call graph

apply_floors

Apply the MHD positivity floors to the interior of a stage state, before its ghost exchange (issue #41, 3.8).

Euler has no floors (return before any work). MHD: the floored cells of the stage are logged by rank 0 when any; a non-positive density or pressure with the floors disabled (both zero) is fatal, reported with the global minimum density and pressure.

subroutine apply_floors(self, q)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.
q	real(kind=R8P)	inout		Conservative variables.

Call graph

allocate_cpu

Allocate CPU data.

subroutine allocate_cpu(self)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.

Call graph

check_glm_ch

Check the GLM cleaning speed against the fastest wave of the committed state (issue #41, section 3.5); a no-op without GLM.

subroutine check_glm_ch(self)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.

Call graph

compute_conservation

Compute the volume integrals of the conservative variables and save them on the diagnostics cadence.

The cell volume includes the null directions: the tree splits them too, so a refined block's cells are smaller along them, and a volume without them overweights the fine cells (issue #37).

subroutine compute_conservation(self)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.

Call graph

compute_divb_history

Compute and save the div(B) norms of the committed state on the diagnostics cadence (issue #41, D-10); a no-op for Euler. The ghost cells are refreshed first and, with sibling realms, the inter-realm seam ghosts are refilled from the peers (update_ghost does not fill them, issue #31): the stencils of the seam-local cells read them.

subroutine compute_divb_history(self, realm)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.
realm	class(realm_object)	inout	target, optional	Sibling realms.

Call graph

compute_q_aux

Compute the auxiliary variables on every cell, ghost cells included.

subroutine compute_q_aux(self, q)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.
q	real(kind=R8P)	in		Conservative variables.

Call graph

initialize_flume

Initialize the CPU backend: MPI, memory budget, common data, CPU data, dispatch.

subroutine initialize_flume(self, filename, realms_number)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.
filename	character(len=*)	in		Input file name.
realms_number	integer(kind=I4P)	in	optional	Realm count; divides the per-process budget.

Call graph

save_residuals

Save residuals history (L2 norm of dq, MPI-reduced, rank 0 writes).

subroutine save_residuals(self)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.

Call graph

save_simulation_data

Save fields, restart, slices and conservation history, each on its own cadence (one predicate per output).

subroutine save_simulation_data(self)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.

Call graph

set_boundary_conditions

Set boundary conditions on the crown maps (face ghosts by kind, edge and corner ghosts by extrapolation).

Periodic faces have no crown rows: their ghosts are filled by the ghost exchange. Edge and corner ghosts (fec > 6) are never read by the directional stencils; they are extrapolated so that every ghost holds a finite, deterministic state for the auxiliary variables computation.

subroutine set_boundary_conditions(self, q)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.
q	real(kind=R8P)	inout		Conservative variables.

Call graph

set_initial_conditions

Set initial conditions.

subroutine set_initial_conditions(self)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.

Call graph

update_ghost

Update ghost cells: intra-realm local copies, MPI exchange, boundary conditions.

Inter-realm seam ghosts are filled by the forest (fill_seam_from_peer_forest), not here.

subroutine update_ghost(self, q)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.
q	real(kind=R8P)	inout		Conservative variables.

Call graph

advance_one_step_forest

Advance one full step of size dt (fast path: single realm, no AMR seam faces).

subroutine advance_one_step_forest(self, dt)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.
dt	real(kind=R8P)	in		Time step from the forest's global reduction.

Call graph

apply_reflux_to_stage_forest

Apply the Berger-Colella reflux correction to the committed q (the forest calls it once per step, after close_step_forest, with the final stage).

For every register face whose coarse side this realm and this rank own, the coarse skin cells receive sgn dt / dx_coarse (F_coarse - F_fine_sum): the register holds the step fluxes sum_s beta_s F_s of both sides (accumulate_seam_fluxes), so the coarse flux the step used is replaced by the restricted fine one, and the coarse-fine interface conserves to round-off.

The scale uses the step the update used, time%dt: on the last step of a time-driven run the realm caps it to land on time_max, and the forest's dt argument is the uncapped value (with it the correction was off by the ratio of the two, 1.3e-9 of the mass on sod-amr; issue #37).

subroutine apply_reflux_to_stage_forest(self, stage, dt, flux_register)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.
stage	integer(kind=I4P)	in		Integrator stage.
dt	real(kind=R8P)	in		Forest time step (unused: see the scale below).
flux_register	class(flux_register_object)	in		Forest's flux register.

Call graph

begin_stage_forest

Begin integrator stage k (staged path): publish the stage and compute its state.

subroutine begin_stage_forest(self, k, K_total, dt, realm)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.
k	integer(kind=I4P)	in		Stage index (1..K_total).
K_total	integer(kind=I4P)	in		Forest-wide stage count for this step.
dt	real(kind=R8P)	in		Time step from the forest.
realm	class(realm_object)	inout	target, optional	Sibling realms (contract parity).

Call graph

close_step_forest

Close a step (staged path): assemble q, save residuals, advance time, clear the active stage.

subroutine close_step_forest(self, dt)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.
dt	real(kind=R8P)	in		Time step from the forest (the local capped value is time%dt).

Call graph

compute_local_dt_forest

Compute the local stability-limited time step, dt = CFL / max(sum_d (|u_d| + a) / dx_d) (no MPI reduction); with GLM, also dt <= CFL / (c_h max sum_d 1 / dx_d) (glm_lambda, issue #41, section 3.5).

The auxiliary variables are recomputed from the committed q (not read from q_aux, which holds the last stage state); null directions do not contribute.

subroutine compute_local_dt_forest(self, dt_local)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	in		The equation.
dt_local	real(kind=R8P)	out		Local stability-limited time step.

Call graph

end_stage_forest

End integrator stage k (staged path): residuals on the stage state, then stage assignment.

subroutine end_stage_forest(self, k, K_total, dt, realm, flux_register)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.
k	integer(kind=I4P)	in		Stage index (1..K_total).
K_total	integer(kind=I4P)	in		Forest-wide stage count.
dt	real(kind=R8P)	in		Time step from the forest.
realm	class(realm_object)	inout	target, optional	Sibling realms (parity only).
flux_register	class(flux_register_object)	inout	optional	Forest's flux register.

Call graph

fill_seam_from_peer_forest

Fill this realm's seam ghosts for peer slot p_idx from the peer's interior, on the active buffers (q when stage_active == 0, else the active stage of q_rk) of both realms.

subroutine fill_seam_from_peer_forest(self, peer, p_idx)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.
peer	class(realm_object)	in	target	Peer realm.
p_idx	integer(kind=I4P)	in		Peer slot.

Call graph

finalize_forest

Finalize the realm: close the output files and free the data (MPI is finalized once by the forest).

subroutine finalize_forest(self)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.

Call graph

initialize_forest

Initialize the realm (issue #35, section 6.1): backend init, IC (or restart), initial AMR, ghost update, initial output, output files open, AMR lock.

subroutine initialize_forest(self, filename, realms_number, memory_avail, nv, verbose)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.
filename	character(len=*)	in		Input file name.
realms_number	integer(kind=I4P)	in	optional	Realm count; divides the per-process budget.
memory_avail	real(kind=R8P)	in	optional	Unused: the budget comes from the MPI handler.
nv	integer(kind=I4P)	in	optional	Unused: nv is decided by the physics.
verbose	logical	in	optional	Unused: initialization is always verbose.

Call graph

is_done_forest

Return true if the realm has reached its end.

subroutine is_done_forest(self, done)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	in		The equation.
done	logical	out		True if the realm is done.

Call graph

open_step_forest

Open a step (staged path): time bookkeeping and Runge-Kutta stages initialization.

subroutine open_step_forest(self, dt)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.
dt	real(kind=R8P)	in		Time step from the forest.

Call graph

post_step_forest

Post-step work: the GLM c_h check, the div(B) history (MHD), fields, restart and conservation history on their cadence.

subroutine post_step_forest(self, dt, t, it, do_save_state, do_save_residuals, do_save_restart, do_amr, realm)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.
dt	real(kind=R8P)	in		Time step just advanced.
t	real(kind=R8P)	in		Time after the advance.
it	integer(kind=I4P)	in		Iteration after the advance.
do_save_state	logical	in	optional	Unused: cadence is internal.
do_save_residuals	logical	in	optional	Unused: cadence is internal.
do_save_restart	logical	in	optional	Unused: cadence is internal.
do_amr	logical	in	optional	Unused: AMR is init-time only.
realm	class(realm_object)	inout	target, optional	Sibling realms.

Call graph

compute_face_mirror_indexes

Return the donor indexes mirrored across a boundary face.

subroutine compute_face_mirror_indexes(face, ni, nj, nk, i_gc, j_gc, k_gc, idelta, jdelta, kdelta, i_d, j_d, k_d)

Arguments

Name	Type	Intent	Attributes	Description
face	integer(kind=I4P)	in		Face index, 1 to 6.
ni	integer(kind=I4P)	in		Grid dimensions.
nj	integer(kind=I4P)	in		Grid dimensions.
nk	integer(kind=I4P)	in		Grid dimensions.
i_gc	integer(kind=I4P)	in		Ghost cell indexes.
j_gc	integer(kind=I4P)	in		Ghost cell indexes.
k_gc	integer(kind=I4P)	in		Ghost cell indexes.
idelta	integer(kind=I4P)	in		Inward steps.
jdelta	integer(kind=I4P)	in		Inward steps.
kdelta	integer(kind=I4P)	in		Inward steps.
i_d	integer(kind=I4P)	out		Mirrored donor indexes.
j_d	integer(kind=I4P)	out		Mirrored donor indexes.
k_d	integer(kind=I4P)	out		Mirrored donor indexes.

Call graph

compute_flux_difference

Compute the residuals from the face fluxes, dq = -sum_d (F_{d,i+1/2} - F_{d,i-1/2}) / dx_d.

A null direction weighs zero, and the residual of the variable it freezes, freeze(d) (0: none), is zero: the normal momentum for Euler (CHASE semantics, issue #35, section 3.4), none for MHD (issue #41, M2-P3). With immersed solids (phi present, its last slot the all-solids summary), the spacing of a fluid cell is cut by the solid surface (ib_cut_spacing, CHASE semantics, D-9).

subroutine compute_flux_difference(nv, ni, nj, nk, ngc, blocks_number, is_null, freeze, dxyz, flx, fly, flz, dq, phi)

Arguments

Name	Type	Intent	Attributes	Description
nv	integer(kind=I4P)	in		Conservative variables number.
ni	integer(kind=I4P)	in		Grid dimensions.
nj	integer(kind=I4P)	in		Grid dimensions.
nk	integer(kind=I4P)	in		Grid dimensions.
ngc	integer(kind=I4P)	in		Grid dimensions.
blocks_number	integer(kind=I4P)	in		Actual blocks number.
is_null	logical	in		Null directions.
freeze	integer(kind=I4P)	in		Variable frozen by each null direction.
dxyz	real(kind=R8P)	in		Blocks space steps [3, nb].
flx	real(kind=R8P)	in		X-face fluxes.
fly	real(kind=R8P)	in		Y-face fluxes.
flz	real(kind=R8P)	in		Z-face fluxes.
dq	real(kind=R8P)	inout		Residuals.
phi	real(kind=R8P)	in	optional	Immersed solids distance function.

Call graph

compute_rk_ssp_residual

Compute the residual of a strong stability preserving step, dq = sum_s beta_s dq_s, from the stored stages.

Used with immersed solids only: the library update_q computes the step residual on its unmasked path but not on the masked one. Same summation order as the FNL kernel compute_rk_ssp_residual_dev.

subroutine compute_rk_ssp_residual(self)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.

Call graph

compute_residuals_weno

Compute the residuals with the WENO space operator: ghost update, auxiliary variables, face fluxes of the active directions, flux difference, and with GLM the psi damping source.

The face fluxes of a null direction are never computed: they keep their zero initialization. On the staged path (AMR seam faces), the seam face fluxes of every stage are accumulated into the forest's flux register. With immersed solids, the eikonal extrapolation fills the solid cells first (n_eikonal Jacobi sweeps, each followed by a ghost exchange, then the wall inversion), and the flux difference uses the spacing cut by the surface.

subroutine compute_residuals_weno(self, q, dq, s, flux_register)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.
q	real(kind=R8P)	inout		Conservative variables.
dq	real(kind=R8P)	inout		Residuals.
s	integer(kind=I4P)	in	optional	Runge-Kutta stage.
flux_register	class(flux_register_object)	inout	optional	Forest's flux register for reflux.

Call graph

integrate_rk_ls

Integrate one time step with a low-storage Runge-Kutta scheme.

subroutine integrate_rk_ls(self)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.

Call graph

integrate_rk_ssp

Integrate one time step with a strong stability preserving Runge-Kutta scheme.

subroutine integrate_rk_ssp(self)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	inout		The equation.

Call graph

Functions

stages_per_step_forest

Return the integrator stages per step: only SSP schemes are stage-splittable (staged path).

Returns: integer(kind=I4P)

function stages_per_step_forest(self) result(K)

Arguments

Name	Type	Intent	Attributes	Description
self	class(flume_cpu_object)	in		The equation.

Call graph