atl_predcor_cerk4_module Module

time integration approach for local timestepping with local predictor and global corrector local predictor: Continuous Extension Runge Kutta with 4 steps (order 3) global corrector: gauss-quadrature with 2 points in time (order 3)


Uses

Used by

  • module~~atl_predcor_cerk4_module~~UsedByGraph module~atl_predcor_cerk4_module atl_predcor_cerk4_module module~atl_global_time_integration_module atl_global_time_integration_module module~atl_global_time_integration_module->module~atl_predcor_cerk4_module module~atl_container_module atl_container_module module~atl_container_module->module~atl_global_time_integration_module module~atl_program_module atl_program_module module~atl_program_module->module~atl_global_time_integration_module module~atl_program_module->module~atl_container_module module~atl_initialize_module atl_initialize_module module~atl_program_module->module~atl_initialize_module module~atl_initialize_module->module~atl_container_module program~ateles ateles program~ateles->module~atl_container_module program~ateles->module~atl_program_module program~atl_harvesting atl_harvesting program~atl_harvesting->module~atl_container_module program~atl_harvesting->module~atl_program_module program~atl_harvesting->module~atl_initialize_module

Subroutines

public subroutine atl_init_explicitLocalPredictorGlobalCorrector(me, minLevel, maxLevel, steps, statedata_list)

Routine to initialize explicit runge kutta scheme for timestepping.

Arguments

Type IntentOptional Attributes Name
type(atl_global_timestep_type), intent(inout) :: me

The datatype to initialize.
integer, intent(in) :: minLevel

The minimum of level of the mesh.
integer, intent(in) :: maxLevel

The maximum of level of the mesh.
integer, intent(in) :: steps

The number of steps in the runge kutta procedure
type(atl_statedata_type), intent(in) :: statedata_list(minLevel:maxLevel)

The state list used in your solver, for each level one entry.

private subroutine update_timestep_predcor_cerk4(me, timestepInfo)

Levelwise updating of runge kutta of order 4

Arguments

Type IntentOptional Attributes Name
class(atl_timestep_type), intent(inout) :: me

The type of your timestepping.
type(atl_local_timestep_type), intent(in) :: timestepInfo

Local timestepping information for that part of the mesh

private subroutine elemental_timestep_predcor_cerk4(me, state, cell, dof, sideFlux)

Elemental operation for timestepping of order 4.

Arguments

Type IntentOptional Attributes Name
class(atl_timestep_type), intent(inout) :: me

Description of the timestep integration method.
real(kind=rk), intent(inout) :: state(:,:,:)

The state of all cells on this level. This field will be updated ad the cell position. See kerneldata type for more explanations.
integer, intent(in) :: cell

Position of the cell to update in the state vector.
integer, intent(in) :: dof

The degree of freedom to update
real(kind=rk), intent(in) :: sideFlux(:)

The flux for one of the sides of this cell. The length of this array is the number of conservative variables of your equation.

private subroutine elemental_timestep_vec_predcor_cerk4(me, state, kerneldata)

Elemental operation for timestepping of order 4.

Arguments

Type IntentOptional Attributes Name
class(atl_timestep_type), intent(inout) :: me

Description of the timestep integration method.
real(kind=rk), intent(inout) :: state(:,:,:)

The state of all cells on this level. This field will be updated ad the cell position. See kerneldata type for more explanations.
type(atl_kerneldata_type), intent(in) :: kerneldata

Complete kerneldata to get the flux from with additional information.

private subroutine mesh_timestep_predcor_cerk4(minLevel, maxLevel, currentLevel, cubes, tree, timestep_list, nSteps, equation, general, commStateTimer, poly_proj_list)

Subroutine for timestepping with explicit runge kutta of order 4.

Read more…

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: minLevel

The minimum refinement level of the mesh.
integer, intent(in) :: maxLevel

The maximum refinement level of the mesh.
integer, intent(in) :: currentLevel

The level the timestep has to be performed for.
type(atl_cube_container_type), intent(inout) :: cubes

Container for the cubical elements.
type(treelmesh_type), intent(in) :: tree

treelm mesh
type(atl_timestep_type), intent(inout) :: timestep_list(minLevel:)

List of levelwise timestepping algorihtms
integer, intent(in) :: nSteps

The number of steps of the time stepping scheme (assumed to be 4)
type(atl_Equations_type), intent(inout) :: equation

The equation you are operating with.
type(tem_general_type), intent(inout) :: general

General treelm settings
integer, intent(inout) :: commStateTimer

Timer for measuring the communication time inside this routine.
type(ply_poly_project_type), intent(inout) :: poly_proj_list(:)

unique list for projection methods

private subroutine local_predictor_substep(mesh, kerneldata, statedata, scheme, general, poly_proj_list, timestep, equation, material, boundary)

Subroutine calculates a substep of the local predictor

Read more…

Arguments

Type IntentOptional Attributes Name
type(atl_cube_elem_type), intent(inout) :: mesh

List of mesh parts. For each level we have one.
type(atl_kerneldata_type), intent(inout) :: kerneldata

List of kerneldatas. For each level we have one
type(atl_statedata_type), intent(inout) :: statedata

List of states you want to calc the rhs for. For each level we have one.
type(atl_scheme_type), intent(inout) :: scheme

List of schemes, for each level.
type(tem_general_type), intent(inout) :: general

General treelm settings
type(ply_poly_project_type), intent(inout) :: poly_proj_list(:)

List of project, for each level.
type(atl_timestep_type), intent(inout) :: timestep

List of levelwise timestepping algorihtms
type(atl_Equations_type), intent(inout) :: equation

The equation you are operating with.
type(atl_material_type), intent(inout) :: material

Material parameter description.
type(atl_level_boundary_type), intent(in) :: boundary

private recursive subroutine global_corrector_substep(minLevel, maxLevel, currentLevel, mesh_list, tree, kerneldata_list, statedata_list, facedata_list, source, penalizationdata_list, boundary_list, bc, scheme_list, poly_proj_pos, poly_proj_list, timestep_list, equation, material_list, general)

Subroutine calculates a substep of corrector, this is the same as a usual substep of the RKDG

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: minLevel

The minimum refinement level of the mesh.
integer, intent(in) :: maxLevel

The maximum refinement level of the mesh.
integer, intent(in) :: currentLevel

The level the timestep has to be performed for.
type(atl_cube_elem_type), intent(inout) :: mesh_list(minLevel:maxLevel)

List of mesh parts. For each level we have one.
type(treelmesh_type) :: tree

treelm mesh
type(atl_kerneldata_type), intent(inout) :: kerneldata_list(minLevel:maxLevel)

List of kerneldatas. For each level we have one
type(atl_statedata_type), intent(inout) :: statedata_list(minLevel:maxLevel)

List of states you want to calc the rhs for. For each level we have one.
type(atl_facedata_type), intent(inout) :: facedata_list(minLevel:maxLevel)

List of faces states you want to calc the rhs for. For each level we have one.
type(atl_source_type), intent(inout) :: source

List of sources, for each level
type(atl_penalizationData_type), intent(inout) :: penalizationdata_list(minLevel:maxLevel)

List of penalization data, for each level.
type(atl_level_boundary_type), intent(inout) :: boundary_list(minLevel:maxLevel)

List of boundaries, for each level.
type(atl_boundary_type), intent(in) :: bc(:)

Global description of the boundaries
type(atl_scheme_type), intent(inout) :: scheme_list(minLevel:maxLevel)

List of schemes, for each level.
integer, intent(inout) :: poly_proj_pos(minLevel:maxLevel)

List of projection pointer, for each level.
type(ply_poly_project_type), intent(inout) :: poly_proj_list(:)

List of projection, for each level.
type(atl_timestep_type), intent(inout) :: timestep_list(minLevel:maxLevel)

List of levelwise timestepping algorihtms
type(atl_Equations_type), intent(inout) :: equation

The equation you are operating with.
type(atl_material_type), intent(inout) :: material_list(minlevel:maxlevel)

Material parameter description.
type(tem_general_type), intent(inout) :: general

General treelm settings