foodie_test_integrand_ladvection Module

  1. foodie_test_integrand_ladvection.f90
  2. foodie_test_integrand_ladvection
module~~foodie_test_integrand_ladvection~~UsesGraph module~foodie_test_integrand_ladvection foodie_test_integrand_ladvection wenoof wenoof wenoof->module~foodie_test_integrand_ladvection module~foodie_test_integrand_tester_object foodie_test_integrand_tester_object module~foodie_test_integrand_tester_object->module~foodie_test_integrand_ladvection module~foodie foodie module~foodie->module~foodie_test_integrand_ladvection module~foodie->module~foodie_test_integrand_tester_object penf penf penf->module~foodie_test_integrand_ladvection penf->module~foodie_test_integrand_tester_object penf->module~foodie module~foodie_integrator_adams_bashforth foodie_integrator_adams_bashforth penf->module~foodie_integrator_adams_bashforth module~foodie_integrator_ms_runge_kutta_ssp foodie_integrator_ms_runge_kutta_ssp penf->module~foodie_integrator_ms_runge_kutta_ssp module~foodie_error_codes foodie_error_codes penf->module~foodie_error_codes module~foodie_integrator_lmm_ssp foodie_integrator_lmm_ssp penf->module~foodie_integrator_lmm_ssp module~foodie_integrator_euler_explicit foodie_integrator_euler_explicit penf->module~foodie_integrator_euler_explicit module~foodie_integrator_multistep_object foodie_integrator_multistep_object penf->module~foodie_integrator_multistep_object module~foodie_integrator_runge_kutta_emd foodie_integrator_runge_kutta_emd penf->module~foodie_integrator_runge_kutta_emd module~foodie_integrator_adams_bashforth_moulton foodie_integrator_adams_bashforth_moulton penf->module~foodie_integrator_adams_bashforth_moulton module~foodie_integrator_adams_moulton foodie_integrator_adams_moulton penf->module~foodie_integrator_adams_moulton module~foodie_integrator_runge_kutta_lssp foodie_integrator_runge_kutta_lssp penf->module~foodie_integrator_runge_kutta_lssp module~foodie_integrator_lmm_ssp_vss foodie_integrator_lmm_ssp_vss penf->module~foodie_integrator_lmm_ssp_vss module~foodie_integrator_runge_kutta_ssp foodie_integrator_runge_kutta_ssp penf->module~foodie_integrator_runge_kutta_ssp module~foodie_integrand_object foodie_integrand_object penf->module~foodie_integrand_object module~foodie_integrator_runge_kutta_low_storage foodie_integrator_runge_kutta_low_storage penf->module~foodie_integrator_runge_kutta_low_storage module~foodie_integrator_multistage_object foodie_integrator_multistage_object penf->module~foodie_integrator_multistage_object module~foodie_integrator_leapfrog foodie_integrator_leapfrog penf->module~foodie_integrator_leapfrog module~foodie_integrator_multistage_multistep_object foodie_integrator_multistage_multistep_object penf->module~foodie_integrator_multistage_multistep_object module~foodie_integrator_object foodie_integrator_object penf->module~foodie_integrator_object module~foodie_integrator_backward_differentiation_formula foodie_integrator_backward_differentiation_formula penf->module~foodie_integrator_backward_differentiation_formula iso_fortran_env iso_fortran_env iso_fortran_env->module~foodie_test_integrand_ladvection iso_fortran_env->module~foodie iso_fortran_env->module~foodie_integrator_multistep_object iso_fortran_env->module~foodie_integrator_multistage_object iso_fortran_env->module~foodie_integrator_multistage_multistep_object iso_fortran_env->module~foodie_integrator_object flap flap flap->module~foodie_test_integrand_ladvection flap->module~foodie_test_integrand_tester_object module~foodie_integrator_adams_bashforth->module~foodie module~foodie_integrator_adams_bashforth->module~foodie_integrator_adams_bashforth_moulton module~foodie_integrator_ms_runge_kutta_ssp->module~foodie module~foodie_error_codes->module~foodie module~foodie_error_codes->module~foodie_integrator_adams_bashforth module~foodie_error_codes->module~foodie_integrator_ms_runge_kutta_ssp module~foodie_error_codes->module~foodie_integrator_lmm_ssp module~foodie_error_codes->module~foodie_integrator_euler_explicit module~foodie_error_codes->module~foodie_integrator_runge_kutta_emd module~foodie_error_codes->module~foodie_integrator_adams_bashforth_moulton module~foodie_error_codes->module~foodie_integrator_adams_moulton module~foodie_error_codes->module~foodie_integrator_runge_kutta_lssp module~foodie_error_codes->module~foodie_integrator_lmm_ssp_vss module~foodie_error_codes->module~foodie_integrator_runge_kutta_ssp module~foodie_error_codes->module~foodie_integrator_runge_kutta_low_storage module~foodie_error_codes->module~foodie_integrator_leapfrog module~foodie_error_codes->module~foodie_integrator_backward_differentiation_formula module~foodie_integrator_lmm_ssp->module~foodie module~foodie_integrator_euler_explicit->module~foodie module~foodie_integrator_multistep_object->module~foodie module~foodie_integrator_multistep_object->module~foodie_integrator_adams_bashforth module~foodie_integrator_multistep_object->module~foodie_integrator_lmm_ssp module~foodie_integrator_multistep_object->module~foodie_integrator_adams_bashforth_moulton module~foodie_integrator_multistep_object->module~foodie_integrator_adams_moulton module~foodie_integrator_multistep_object->module~foodie_integrator_lmm_ssp_vss module~foodie_integrator_multistep_object->module~foodie_integrator_leapfrog module~foodie_integrator_multistep_object->module~foodie_integrator_backward_differentiation_formula module~foodie_integrator_runge_kutta_emd->module~foodie module~foodie_integrator_adams_bashforth_moulton->module~foodie module~foodie_integrator_adams_moulton->module~foodie module~foodie_integrator_adams_moulton->module~foodie_integrator_adams_bashforth_moulton module~foodie_integrator_runge_kutta_lssp->module~foodie module~foodie_integrator_lmm_ssp_vss->module~foodie module~foodie_integrator_runge_kutta_ssp->module~foodie module~foodie_integrand_object->module~foodie module~foodie_integrand_object->module~foodie_integrator_adams_bashforth module~foodie_integrand_object->module~foodie_integrator_ms_runge_kutta_ssp module~foodie_integrand_object->module~foodie_integrator_lmm_ssp module~foodie_integrand_object->module~foodie_integrator_euler_explicit module~foodie_integrand_object->module~foodie_integrator_multistep_object module~foodie_integrand_object->module~foodie_integrator_runge_kutta_emd module~foodie_integrand_object->module~foodie_integrator_adams_bashforth_moulton module~foodie_integrand_object->module~foodie_integrator_adams_moulton module~foodie_integrand_object->module~foodie_integrator_runge_kutta_lssp module~foodie_integrand_object->module~foodie_integrator_lmm_ssp_vss module~foodie_integrand_object->module~foodie_integrator_runge_kutta_ssp module~foodie_integrand_object->module~foodie_integrator_runge_kutta_low_storage module~foodie_integrand_object->module~foodie_integrator_multistage_object module~foodie_integrand_object->module~foodie_integrator_leapfrog module~foodie_integrand_object->module~foodie_integrator_multistage_multistep_object module~foodie_integrand_object->module~foodie_integrator_backward_differentiation_formula module~foodie_integrator_runge_kutta_low_storage->module~foodie module~foodie_integrator_multistage_object->module~foodie module~foodie_integrator_multistage_object->module~foodie_integrator_euler_explicit module~foodie_integrator_multistage_object->module~foodie_integrator_runge_kutta_emd module~foodie_integrator_multistage_object->module~foodie_integrator_runge_kutta_lssp module~foodie_integrator_multistage_object->module~foodie_integrator_runge_kutta_ssp module~foodie_integrator_multistage_object->module~foodie_integrator_runge_kutta_low_storage module~foodie_integrator_leapfrog->module~foodie module~foodie_integrator_multistage_multistep_object->module~foodie module~foodie_integrator_multistage_multistep_object->module~foodie_integrator_ms_runge_kutta_ssp module~foodie_integrator_object->module~foodie module~foodie_integrator_object->module~foodie_integrator_adams_bashforth module~foodie_integrator_object->module~foodie_integrator_ms_runge_kutta_ssp module~foodie_integrator_object->module~foodie_integrator_lmm_ssp module~foodie_integrator_object->module~foodie_integrator_euler_explicit module~foodie_integrator_object->module~foodie_integrator_multistep_object module~foodie_integrator_object->module~foodie_integrator_runge_kutta_emd module~foodie_integrator_object->module~foodie_integrator_adams_bashforth_moulton module~foodie_integrator_object->module~foodie_integrator_adams_moulton module~foodie_integrator_object->module~foodie_integrator_runge_kutta_lssp module~foodie_integrator_object->module~foodie_integrator_lmm_ssp_vss module~foodie_integrator_object->module~foodie_integrator_runge_kutta_ssp module~foodie_integrator_object->module~foodie_integrator_runge_kutta_low_storage module~foodie_integrator_object->module~foodie_integrator_multistage_object module~foodie_integrator_object->module~foodie_integrator_leapfrog module~foodie_integrator_object->module~foodie_integrator_multistage_multistep_object module~foodie_integrator_object->module~foodie_integrator_backward_differentiation_formula module~foodie_integrator_backward_differentiation_formula->module~foodie
Help

Define integrand_ladvection, the 1D linear advection PDE test field that is a concrete extension of the abstract integrand type.

Used By

module~~foodie_test_integrand_ladvection~~UsedByGraph module~foodie_test_integrand_ladvection foodie_test_integrand_ladvection module~foodie_test_object foodie_test_object module~foodie_test_integrand_ladvection->module~foodie_test_object program~foodie_tester foodie_tester module~foodie_test_object->program~foodie_tester
Help

Derived Types

Functions

Subroutines


Derived Types

type, public, extends(integrand_tester_object) :: integrand_ladvection

1D linear advection field.

Components

TypeVisibility AttributesNameInitial
real(kind=R_P), public :: CFL =0._R_P

CFL value.
real(kind=R_P), public :: Dx =0._R_P

Space step.
integer(kind=I_P), public :: Ng =0

Ghost cells number.
integer(kind=I_P), public :: Ni =0

Space dimension.
real(kind=R_P), public :: a =0._R_P

Advection coefficient.
character(len=99), public :: initial_state

Initial state.
class(interpolator_object), public, allocatable:: interpolator

WENO interpolator.
real(kind=R_P), public :: length =0._R_P

Domain length.
real(kind=R_P), public, allocatable:: u(:)

Integrand (state) variable.
character(len=99), public :: w_scheme =''

WENO Scheme used.
real(kind=R_P), public :: weno_eps =0._R_P

WENO epsilon to avoid division by zero, default value.
integer(kind=I_P), public :: weno_order =0

WENO reconstruction order.

Type-Bound Procedures

generic, public :: add_fast => integrand_add_integrand_fast

Overloading add_fast method.

procedure, public, pass(lhs) :: assign_integrand

= operator.

procedure, public, pass(lhs) :: assign_real

= real operator.

generic, public :: assignment(=) => assign_integrand, assign_real

Overloading = assignament.

procedure, public, pass(self) :: compute_dx

Compute the space step by means of CFL condition.

procedure, public, pass(self) :: description

Return an informative description of the test.

procedure, public, pass(self) :: destroy

Destroy field.

procedure, public, pass(self) :: dt => compute_dt

Compute the current time step by means of CFL condition.

procedure, public, pass(self) :: error

Return error.

procedure, public, pass(self) :: exact_solution

Return exact solution.

procedure, public, pass(self) :: export_tecplot

Export integrand to Tecplot file.

procedure, private, pass(self) :: impose_boundary_conditions

Impose boundary conditions.

procedure, public, pass(self) :: initialize

Initialize integrand.

procedure, public, pass(lhs) :: integrand_add_integrand

+ operator.

procedure, public, pass(opr) :: integrand_add_integrand_fast

+ fast operator.

procedure, public, pass(lhs) :: integrand_add_real

+ real operator.

procedure, public, pass(self) :: integrand_dimension

Return integrand dimension.

procedure, public, pass(lhs) :: integrand_multiply_integrand

* operator.

procedure, public, pass(opr) :: integrand_multiply_integrand_fast

* fast operator.

procedure, public, pass(lhs) :: integrand_multiply_real

* real operator.

procedure, public, pass(lhs) :: integrand_multiply_real_scalar

* real_scalar operator.

procedure, public, pass(opr) :: integrand_multiply_real_scalar_fast

* real_scalar fast operator.

procedure, public, pass(lhs) :: integrand_sub_integrand

- operator.

procedure, public, pass(lhs) :: integrand_sub_real

- real operator.

procedure, public, pass(opr) :: integrand_subtract_integrand_fast

- fast operator.

procedure, public, pass(lhs) :: local_error

||integrand_ladvection - integrand_ladvection|| operator.

generic, public :: multiply_fast => integrand_multiply_integrand_fast, integrand_multiply_real_scalar_fast

Overloading multiply_fast method.

generic, public :: operator(*) => integrand_multiply_integrand, integrand_multiply_real, real_multiply_integrand, integrand_multiply_real_scalar, real_scalar_multiply_integrand

Overloading * operator.

generic, public :: operator(+) => integrand_add_integrand, integrand_add_real, real_add_integrand

Overloading + operator.

generic, public :: operator(-) => integrand_sub_integrand, integrand_sub_real, real_sub_integrand

Overloading - operator.

generic, public :: operator(.lterror.) => local_error

Estimate local truncation error.

procedure, public, pass(self) :: output

Extract integrand state field.

procedure, public, pass(self) :: parse_cli

Initialize from command line interface.

procedure, public, pass(rhs) :: real_add_integrand

real + operator.

procedure, public, pass(rhs) :: real_multiply_integrand

real * operator.

procedure, public, pass(rhs) :: real_scalar_multiply_integrand

real_scalar * operator.

procedure, public, pass(rhs) :: real_sub_integrand

real - operator.

procedure, private, pass(self) :: reconstruct_interfaces

Reconstruct interface states.

procedure, public, nopass :: set_cli

Set command line interface.

procedure, private, pass(self) :: set_sin_wave_initial_state

Set initial state as a sin wave.

procedure, private, pass(self) :: set_square_wave_initial_state

Set initial state as a square wave.

generic, public :: subtract_fast => integrand_subtract_integrand_fast

Overloading subtract_fast method.

procedure, public, pass(self) :: t => dU_dt

Time derivative, residuals.

procedure, public, pass(self) :: t_fast

Time derivative, residuals, fast mode.

Functions

private pure function compute_dt(self, final_time, t) result(Dt)

Compute the current time step by means of CFL condition.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(in) :: self

Advection field.
real(kind=R_P), intent(in) :: final_time

Maximum integration time.
real(kind=R_P), intent(in), optional :: t

Time.

Return Value real(kind=R_P)

Time step.

private pure function compute_dx(self, Dt) result(Dx)

Compute the space step step by means of CFL condition.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(in) :: self

Advection field.
real(kind=R_P), intent(in) :: Dt

Time step.

Return Value real(kind=R_P)

Space step.

private function dU_dt(self, t) result(dState_dt)

Time derivative of advection field, the residuals function.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(in) :: self

Advection field.
real(kind=R_P), intent(in), optional :: t

Time.

Return Value real(kind=R_P), allocatable, (:)

Advection field time derivative.

private pure function description(self, prefix) result(desc)

Return informative integrator description.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(in) :: self

Integrand.
character(len=*), intent(in), optional :: prefix

Prefixing string.

Return Value character(len=:), allocatable

Description.

private pure function error(self, t, t0, U0)

Return error.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(in) :: self

Integrand.
real(kind=R_P), intent(in) :: t

Time.
real(kind=R_P), intent(in), optional :: t0

Initial time.
class(integrand_object), intent(in), optional :: U0

Initial conditions.

Return Value real(kind=R_P), allocatable, (:)

Error.

private pure function exact_solution(self, t, t0, U0) result(exact)

Return exact solution.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(in) :: self

Integrand.
real(kind=R_P), intent(in) :: t

Time.
real(kind=R_P), intent(in), optional :: t0

Initial time.
class(integrand_object), intent(in), optional :: U0

Initial conditions.

Return Value real(kind=R_P), allocatable, (:)

Exact solution.

private pure function integrand_add_integrand(lhs, rhs) result(opr)

+ operator.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(in) :: lhs

Left hand side.
class(integrand_object), intent(in) :: rhs

Right hand side.

Return Value real(kind=R_P), allocatable, (:)

Operator result.

private pure function integrand_add_real(lhs, rhs) result(opr)

+ real operator.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(in) :: lhs

Left hand side.
real(kind=R_P), intent(in) :: rhs(1:)

Right hand side.

Return Value real(kind=R_P), allocatable, (:)

Operator result.

private pure function integrand_dimension(self)

return integrand dimension.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(in) :: self

integrand.

Return Value integer(kind=I_P)

integrand dimension.

private pure function integrand_multiply_integrand(lhs, rhs) result(opr)

* operator.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(in) :: lhs

Left hand side.
class(integrand_object), intent(in) :: rhs

Right hand side.

Return Value real(kind=R_P), allocatable, (:)

Operator result.

private pure function integrand_multiply_real(lhs, rhs) result(opr)

* real_scalar operator.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(in) :: lhs

Left hand side.
real(kind=R_P), intent(in) :: rhs(1:)

Right hand side.

Return Value real(kind=R_P), allocatable, (:)

Operator result.

private pure function integrand_multiply_real_scalar(lhs, rhs) result(opr)

* real_scalar operator.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(in) :: lhs

Left hand side.
real(kind=R_P), intent(in) :: rhs

Right hand side.

Return Value real(kind=R_P), allocatable, (:)

Operator result.

private pure function integrand_sub_integrand(lhs, rhs) result(opr)

- operator.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(in) :: lhs

Left hand side.
class(integrand_object), intent(in) :: rhs

Right hand side.

Return Value real(kind=R_P), allocatable, (:)

Operator result.

private pure function integrand_sub_real(lhs, rhs) result(opr)

- real operator.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(in) :: lhs

Left hand side.
real(kind=R_P), intent(in) :: rhs(1:)

Right hand side.

Return Value real(kind=R_P), allocatable, (:)

Operator result.

private function local_error(lhs, rhs) result(error)

Estimate local truncation error between 2 advection approximations.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(in) :: lhs

Left hand side.
class(integrand_object), intent(in) :: rhs

Right hand side.

Return Value real(kind=R_P)

Error estimation.

private pure function output(self) result(state)

Output the advection field state.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(in) :: self

Advection field.

Return Value real(kind=R_P), allocatable, (:)

Advection state

private pure function real_add_integrand(lhs, rhs) result(opr)

real + operator.

Arguments

Type IntentOptional AttributesName
real(kind=R_P), intent(in) :: lhs(1:)

Left hand side.
class(integrand_ladvection), intent(in) :: rhs

Left hand side.

Return Value real(kind=R_P), allocatable, (:)

Operator result.

private pure function real_multiply_integrand(lhs, rhs) result(opr)

real_scalar * operator.

Arguments

Type IntentOptional AttributesName
real(kind=R_P), intent(in) :: lhs(1:)

Left hand side.
class(integrand_ladvection), intent(in) :: rhs

Right hand side.

Return Value real(kind=R_P), allocatable, (:)

Operator result.

private pure function real_scalar_multiply_integrand(lhs, rhs) result(opr)

real_scalar * operator.

Arguments

Type IntentOptional AttributesName
real(kind=R_P), intent(in) :: lhs

Left hand side.
class(integrand_ladvection), intent(in) :: rhs

Right hand side.

Return Value real(kind=R_P), allocatable, (:)

Operator result.

private pure function real_sub_integrand(lhs, rhs) result(opr)

real - operator.

Arguments

Type IntentOptional AttributesName
real(kind=R_P), intent(in) :: lhs(1:)

Left hand side.
class(integrand_ladvection), intent(in) :: rhs

Left hand side.

Return Value real(kind=R_P), allocatable, (:)

Operator result.


Subroutines

private subroutine assign_integrand(lhs, rhs)

= operator.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(inout) :: lhs

Left hand side.
class(integrand_object), intent(in) :: rhs

Right hand side.

private pure subroutine assign_real(lhs, rhs)

Assign one real to an advection field.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(inout) :: lhs

Left hand side.
real(kind=R_P), intent(in) :: rhs(1:)

Right hand side.

private subroutine destroy(self)

Destroy field.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(inout) :: self

Advection field.

private subroutine export_tecplot(self, file_name, t, scheme, close_file, with_exact_solution, U0)

Export integrand to Tecplot file.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(in) :: self

Advection field.
character(len=*), intent(in), optional :: file_name

File name.
real(kind=R_P), intent(in), optional :: t

Time.
character(len=*), intent(in), optional :: scheme

Scheme used to integrate integrand.
logical, intent(in), optional :: close_file

Flag for closing file.
logical, intent(in), optional :: with_exact_solution

Flag for export also exact solution.
class(integrand_object), intent(in), optional :: U0

Initial conditions.

private pure subroutine impose_boundary_conditions(self, u)

Impose boundary conditions.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(in) :: self

Advection field.
real(kind=R_P), intent(inout) :: u(1-self%Ng:)

Conservative variables.

private subroutine initialize(self, Dt)

Initialize integrand.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(inout) :: self

Integrand.
real(kind=R_P), intent(in) :: Dt

Time step.

private subroutine parse_cli(self, cli)

Initialize from command line interface.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(inout) :: self

Advection field.
type(command_line_interface), intent(inout) :: cli

Command line interface handler.

private subroutine reconstruct_interfaces(self, conservative, r_conservative)

Reconstruct interfaces states.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(in) :: self

Advection field.
real(kind=R_P), intent(in) :: conservative(1-self%Ng:)

Conservative variables.
real(kind=R_P), intent(inout) :: r_conservative(1:,0:)

Reconstructed conservative vars.

private subroutine set_cli(cli)

Set command line interface.

Arguments

Type IntentOptional AttributesName
type(command_line_interface), intent(inout) :: cli

Command line interface handler.

private subroutine set_sin_wave_initial_state(self)

Set initial state as a sin wave.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(inout) :: self

Advection field.

private subroutine set_square_wave_initial_state(self)

Set initial state as a square wave.

Arguments

Type IntentOptional AttributesName
class(integrand_ladvection), intent(inout) :: self

Advection field.

All Modules