Merge pull request #1048 from prbzrg/new-AlgorithmInterpretation

update for new `AlgorithmInterpretation`

Author: ChrisRackauckas

Project.toml

@@ -34,9 +34,9 @@ NonlinearSolve = "3"
 OrdinaryDiffEq = "6.53"
 RecursiveArrayTools = "3"
 Reexport = "1.0"
-SciMLBase = "2"
+SciMLBase = "2.51"
 SteadyStateDiffEq = "2"
-StochasticDiffEq = "6.61"
+StochasticDiffEq = "6.69"
 Sundials = "4.19"
 julia = "1.9"
 

src/DifferentialEquations.jl

@@ -21,6 +21,7 @@ using LinearAlgebra
 import DiffEqBase: solve
 @reexport using LinearSolve
 using NonlinearSolve
+import SciMLBase
 
 include("default_solve.jl")
 include("default_arg_parsing.jl")

src/sde_default_alg.jl

@@ -13,16 +13,17 @@ function default_algorithm(
     end
 
     is_stiff = :stiff ∈ alg_hints
-    is_stratonovich = :Stratonovich ∈ alg_hints
+    is_stratonovich = :stratonovich ∈ alg_hints
     if is_stiff || prob.f.mass_matrix !== I
         alg = ImplicitRKMil(autodiff = false)
     end
 
     if is_stratonovich
         if is_stiff || prob.f.mass_matrix !== I
-            alg = ImplicitRKMil(autodiff = false, interpretation = :Stratonovich)
+            alg = ImplicitRKMil(autodiff = false,
+                interpretation = SciMLBase.AlgorithmInterpretation.Stratonovich)
         else
-            alg = RKMil(interpretation = :Stratonovich)
+            alg = RKMil(interpretation = SciMLBase.AlgorithmInterpretation.Stratonovich)
         end
     end
 

test/default_sde_alg_test.jl

@@ -1,4 +1,5 @@
 using DifferentialEquations, Test
+import SciMLBase
 
 f_additive(u, p, t) = @. p[2] / sqrt(1 + t) - u / (2 * (1 + t))
 σ_additive(u, p, t) = @. p[1] * p[2] / sqrt(1 + t)
@@ -17,8 +18,9 @@ sol = solve(prob, dt = 1 / 2^(3))
 sol = solve(prob, dt = 1 / 2^(3), alg_hints = [:additive])
 @test sol.alg isa SOSRA
 
-sol = solve(prob, dt = 1 / 2^(3), alg_hints = [:Stratonovich])
-@test StochasticDiffEq.alg_interpretation(sol.alg) == :Stratonovich
+sol = solve(prob, dt = 1 / 2^(3), alg_hints = [:stratonovich])
+@test SciMLBase.alg_interpretation(sol.alg) ==
+      SciMLBase.AlgorithmInterpretation.Stratonovich
 @test sol.alg isa RKMil
 
 f = (du, u, p, t) -> du .= 1.01u
@@ -43,5 +45,5 @@ sol = solve(prob, dt = 1 / 2^(3), alg_hints = [:stiff])
 sol = solve(prob, dt = 1 / 2^(3), alg_hints = [:additive])
 @test sol.alg isa SOSRA
 
-sol = solve(prob, dt = 1 / 2^(3), alg_hints = [:Stratonovich])
+sol = solve(prob, dt = 1 / 2^(3), alg_hints = [:stratonovich])
 @test sol.alg isa LambaEulerHeun