Skip to content

Commit

Permalink
Update lib/OrdinaryDiffEqStabilizedIRK/src/irkc_perform_step.jl
Browse files Browse the repository at this point in the history
  • Loading branch information
ChrisRackauckas authored Sep 11, 2024
1 parent aeafc10 commit fd7cb70
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion lib/OrdinaryDiffEqStabilizedIRK/src/irkc_perform_step.jl
Original file line number Diff line number Diff line change
Expand Up @@ -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))
Expand Down

0 comments on commit fd7cb70

Please sign in to comment.