Skip to content

Commit

Permalink
Merge pull request #900 from SciML/ChrisRackauckas-patch-1
Browse files Browse the repository at this point in the history
Change default Rodas5 to Rodas5P
  • Loading branch information
ChrisRackauckas authored Sep 19, 2022
2 parents 9323fcc + 37dae2a commit d030563
Show file tree
Hide file tree
Showing 2 changed files with 3 additions and 3 deletions.
2 changes: 1 addition & 1 deletion src/ode_default_alg.jl
Original file line number Diff line number Diff line change
Expand Up @@ -70,7 +70,7 @@ function default_algorithm(prob::DiffEqBase.AbstractODEProblem{uType, tType, inp
if length(prob.u0) > 100
alg = AutoVern9(KenCarp47(autodiff = false), lazy = !callbacks)
else
alg = AutoVern9(Rodas5(autodiff = false), lazy = !callbacks)
alg = AutoVern9(Rodas5P(autodiff = false), lazy = !callbacks)
end
elseif tol_level == :low_tol
if length(prob.u0) > 500
Expand Down
4 changes: 2 additions & 2 deletions test/default_ode_alg_test.jl
Original file line number Diff line number Diff line change
Expand Up @@ -24,7 +24,7 @@ sol = solve(prob_ode_2Dlinear; reltol = 1e-7)
sol = solve(prob_ode_2Dlinear; reltol = 1e-10)

@test typeof(sol.alg.algs[1]) <: Vern9
@test typeof(sol.alg.algs[2]) <: Rodas5
@test typeof(sol.alg.algs[2]) <: Rodas5P

sol = solve(prob_ode_2Dlinear; alg_hints = [:stiff])

Expand All @@ -44,7 +44,7 @@ prob_ode_bigfloat2Dlinear = ODEProblem(f, map(BigFloat, rand(4, 2)) .* ones(4, 2

sol = solve(prob_ode_bigfloat2Dlinear; dt = 1 // 2^(4))
@test typeof(sol.alg.algs[1]) <: Vern9
@test typeof(sol.alg.algs[2]) <: Rodas5
@test typeof(sol.alg.algs[2]) <: Rodas5P

default_algorithm(prob_ode_bigfloat2Dlinear; alg_hints = [:stiff])

Expand Down

0 comments on commit d030563

Please sign in to comment.