From fd7cb708bb07bb1c612a99c043ca14ed677dcf36 Mon Sep 17 00:00:00 2001 From: Christopher Rackauckas Date: Wed, 11 Sep 2024 10:19:18 +0200 Subject: [PATCH] Update lib/OrdinaryDiffEqStabilizedIRK/src/irkc_perform_step.jl --- lib/OrdinaryDiffEqStabilizedIRK/src/irkc_perform_step.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) 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))