diff --git a/integration/VODE/_parameters b/integration/VODE/_parameters index 7bfebf611..1b2165b85 100644 --- a/integration/VODE/_parameters +++ b/integration/VODE/_parameters @@ -3,3 +3,8 @@ # for the step rejection logic on mass fractions, we only consider # species that are > X_reject_buffer * atol_spec X_reject_buffer real 1.0 + +# set the max integration order for VODE +# For the backward differentiation formula (BDF) integration +# the maximum order should be no greater than 5. +VODE_MAXORD int 5 diff --git a/integration/VODE/vode_type.H b/integration/VODE/vode_type.H index 3c246bfb5..15218e495 100644 --- a/integration/VODE/vode_type.H +++ b/integration/VODE/vode_type.H @@ -22,10 +22,7 @@ constexpr amrex::Real HMIN = 0.0_rt; constexpr amrex::Real vode_increase_change_factor = 4.0_rt; constexpr amrex::Real vode_decrease_change_factor = 0.25_rt; -// For the backward differentiation formula (BDF) integration -// the maximum order should be no greater than 5. -constexpr int VODE_MAXORD = 5; -constexpr int VODE_LMAX = VODE_MAXORD + 1; +int VODE_LMAX = integrator_rp::VODE_MAXORD + 1; // How many timesteps should pass before refreshing the Jacobian constexpr int max_steps_between_jacobian_evals = 50;