mus_interpolate_verify_module Module

Interpolation scheme tools

For an overview over implemented interpolation methods, see Interpolation methods



Subroutines

public subroutine mus_testInterpolation(scheme, tree, general, fac, iLevel, minLevel, maxLevel, pdf)

Call tests to determine the actual error from the interpolation routines on the ghost elements. Compare against the analytical solution, which is given in terms of the initial conditions. Call this routine after the initial values are set and the ghost elements have been filled once, but no computation was started -> after fillHelperElements in the mus_aux_module

Arguments

Type IntentOptional Attributes Name
type(mus_scheme_type), intent(inout) :: scheme
type(treelmesh_type), intent(in) :: tree

treelm tree

type(tem_general_type), intent(in) :: general
type(mus_convertFac_type), intent(in) :: fac
integer, intent(in) :: iLevel

Level counter variable

integer, intent(in) :: minLevel
integer, intent(in) :: maxLevel
type(pdf_data_type), intent(inout) :: pdf(tree%global%minlevel:tree%global%maxlevel)

private subroutine mus_intp_error(pdf, level, tree, general, fac, scheme, layout, intp, ind, varSys, derVarPos, eType)

Determine the numerical error of the interpolated quantities to the given initial conditions

Arguments

Type IntentOptional Attributes Name
type(pdf_data_type), intent(in) :: pdf

global fluid parameters

integer, intent(in) :: level

level

type(treelmesh_type), intent(in) :: tree

treelm tree

type(tem_general_type), intent(in) :: general
type(mus_convertFac_type), intent(in) :: fac
type(mus_scheme_type), intent(in) :: scheme

the current scheme

type(mus_scheme_layout_type), intent(in) :: layout

the layout used

type(mus_interpolation_type), intent(in) :: intp

interpolation method info

type(grw_intarray_type), intent(in) :: ind

indirection list

type(tem_varSys_type), intent(in) :: varSys

global variable system

type(mus_derVarPos_type), intent(in) :: derVarPos(:)

required variable system which maps to global varsys

integer, intent(in) :: eType

element type, eT_ghostFromFiner, eT_ghostFromCoarser