begin_integration_rk_variable_step_class Subroutine

private subroutine begin_integration_rk_variable_step_class(me)

Begin a rk_variable_step_class integration.

Type Bound

rk_variable_step_class

Arguments

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

Calls

proc~~begin_integration_rk_variable_step_class~~CallsGraph proc~begin_integration_rk_variable_step_class rk_variable_step_class%begin_integration_rk_variable_step_class begin begin proc~begin_integration_rk_variable_step_class->begin proc~destroy_fsal_cache rk_variable_step_class%destroy_fsal_cache proc~begin_integration_rk_variable_step_class->proc~destroy_fsal_cache proc~stepsize_class_reset stepsize_class%stepsize_class_reset proc~begin_integration_rk_variable_step_class->proc~stepsize_class_reset

Called by

proc~~begin_integration_rk_variable_step_class~~CalledByGraph proc~begin_integration_rk_variable_step_class rk_variable_step_class%begin_integration_rk_variable_step_class proc~integrate_to_event_variable_step rk_variable_step_class%integrate_to_event_variable_step proc~integrate_to_event_variable_step->proc~begin_integration_rk_variable_step_class proc~integrate_variable_step rk_variable_step_class%integrate_variable_step proc~integrate_variable_step->proc~begin_integration_rk_variable_step_class program~rklib_example rklib_example program~rklib_example->proc~integrate_variable_step

Source Code

    subroutine begin_integration_rk_variable_step_class(me)
    class(rk_variable_step_class),intent(inout) :: me

    call me%begin() ! base

    ! variable step params:
    me%num_rejected_steps = 0
    me%last_accepted_step_size = zero

    ! report at a fixed step size using the dense output (only if reporting):
    me%dense_reporting = me%dense_output_step > zero .and. &
                         associated(me%report) .and. me%report_rate > 0
    me%dense_output_k = 0
    call me%destroy_fsal_cache()
    call me%stepsize_method%reset()

    end subroutine begin_integration_rk_variable_step_class