Skip to content
Snippets Groups Projects
Commit c730cead authored by Omar Awile's avatar Omar Awile
Browse files

Activated schemes in ODE create procedure

parent e8e284bc
No related branches found
No related tags found
No related merge requests found
......@@ -61,6 +61,21 @@ subroutine ode_create(this,scheme,fields,rhsfunc,rhs_fields_discr,info,options,k
allocate(ppm_t_sts::this%integrator,STAT=info)
call this%integrator%create(fields,rhsfunc,rhs_fields_discr,info,options)
or_fail("Creating STS failed")
case(ppm_param_ode_scheme_tvdrk2)
! allocate changes array
allocate(ppm_t_tvdrk2::this%integrator,STAT=info)
call this%integrator%create(fields,rhsfunc,rhs_fields_discr,info,options)
or_fail("Creating tvd RK2 failed")
case(ppm_param_ode_scheme_midrk2)
! allocate changes array
allocate(ppm_t_midrk2::this%integrator,STAT=info)
call this%integrator%create(fields,rhsfunc,rhs_fields_discr,info,options)
or_fail("Creating mid RK2 failed")
case(ppm_param_ode_scheme_rk4)
! allocate changes array
allocate(ppm_t_rk4::this%integrator,STAT=info)
call this%integrator%create(fields,rhsfunc,rhs_fields_discr,info,options)
or_fail("Creating RK4 failed")
case default
ppm_fail("Integrator not implemented")
end select
......@@ -86,6 +101,21 @@ subroutine ode_create(this,scheme,fields,rhsfunc,rhs_fields_discr,info,options,k
allocate(ppm_t_sts::this%kickoff,STAT=info)
call this%kickoff%create(fields,rhsfunc,rhs_fields_discr,info,options)
or_fail("Creating STS failed")
case(ppm_param_ode_scheme_tvdrk2)
! allocate changes array
allocate(ppm_t_tvdrk2::this%kickoff,STAT=info)
call this%kickoff%create(fields,rhsfunc,rhs_fields_discr,info,options)
or_fail("Creating tvd RK2 failed")
case(ppm_param_ode_scheme_midrk2)
! allocate changes array
allocate(ppm_t_midrk2::this%kickoff,STAT=info)
call this%kickoff%create(fields,rhsfunc,rhs_fields_discr,info,options)
or_fail("Creating mid RK2 failed")
case(ppm_param_ode_scheme_rk4)
! allocate changes array
allocate(ppm_t_rk4::this%kickoff,STAT=info)
call this%kickoff%create(fields,rhsfunc,rhs_fields_discr,info,options)
or_fail("Creating RK4 failed")
case default
ppm_fail("Integrator not implemented")
end select
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment