diff --git a/lib/OrdinaryDiffEqStabilizedIRK/src/irkc_perform_step.jl b/lib/OrdinaryDiffEqStabilizedIRK/src/irkc_perform_step.jl index f74ccfc062..5c3820701a 100644 --- a/lib/OrdinaryDiffEqStabilizedIRK/src/irkc_perform_step.jl +++ b/lib/OrdinaryDiffEqStabilizedIRK/src/irkc_perform_step.jl @@ -156,7 +156,7 @@ function perform_step!(integrator, cache::IRKCCache, repeat_step = false) # The the number of degree for Chebyshev polynomial #maxm = max(2,int(floor(sqrt(integrator.opts.internalnorm(integrator.opts.reltol,t)/(10 *eps(integrator.opts.internalnorm(uprev,t))))))) maxm = 50 - mdeg = 1 + floor(Int, sqrt(1.54 * abs(dt) * integrator.eigen_est + 1)) + mdeg = 1 + Int(floor(sqrt(1.54 * abs(dt) * integrator.eigen_est + 1))) mdeg = clamp(mdeg, minm, maxm) ω₀ = 1 + 2 / (13 * (mdeg^2))