Bump compat for OrdinaryDiffEq v7 / SciMLBase v3 ecosystem

Cascade major-version compat relaxations for the OrdinaryDiffEq v7
ecosystem release (SciML/OrdinaryDiffEq.jl#3562 + #3565):

Project.toml / docs/Project.toml / test/extensions/Project.toml:
- OrdinaryDiffEqBDF:        "1" → "1, 2"
- OrdinaryDiffEqCore:   "3.22" → "3.22, 4"
- OrdinaryDiffEqDefault:    "1" → "1, 2"
- OrdinaryDiffEqRosenbrock: "1" → "1, 2"
- OrdinaryDiffEqSDIRK:      "1" → "1, 2"  (docs only)
- OrdinaryDiffEqTsit5:      "1" → "1, 2"
- OrdinaryDiffEqVerner:     "1" → "1, 2"
- SciMLBase:             "2.84" → "2.84, 3"
- StochasticDiffEq:     "6.101" → "6.101, 7"

Source migrations for removed/changed APIs:
- test/simulation_and_solving/hybrid_models.jl: remove `verbose = false`
  from `solve` calls (Bool no longer accepted; use ODEVerbosity in v7)
- test/simulation_and_solving/simulate_SDEs.jl: iterate
  `cat_sol.u` / `hand_sol.u` instead of the EnsembleSolution directly
  (RAT v4: iterating EnsembleSolution now yields scalar elements, not
  trajectory solutions)

Co-Authored-By: Chris Rackauckas <accounts@chrisrackauckas.com>

Author: ChrisRackauckas

Project.toml

@@ -65,14 +65,19 @@ MacroTools = "0.5.5"
 Makie = "0.22.1, 0.23, 0.24"
 ModelingToolkitBase = "1.17"
 NetworkLayout = "0.4.7"
-OrdinaryDiffEqCore = "3.22"
+OrdinaryDiffEqBDF = "1, 2"
+OrdinaryDiffEqCore = "3.22, 4"
+OrdinaryDiffEqDefault = "1, 2"
+OrdinaryDiffEqRosenbrock = "1, 2"
+OrdinaryDiffEqTsit5 = "1, 2"
+OrdinaryDiffEqVerner = "1, 2"
 Parameters = "0.12"
 Reexport = "1.0"
 RuntimeGeneratedFunctions = "0.5.12"
-SciMLBase = "2.84"
+SciMLBase = "2.84, 3"
 SciMLPublic = "1"
 Setfield = "1"
-StochasticDiffEq = "6.101"
+StochasticDiffEq = "6.101, 7"
 StructuralIdentifiability = "0.5.11"
 SymbolicIndexingInterface = "0.3"
 SymbolicUtils = "4.9.1"

docs/Project.toml

@@ -82,24 +82,24 @@ OptimizationBase = "4, 5.0"
 OptimizationNLopt = "0.3"
 OptimizationOptimJL = "0.4"
 OptimizationOptimisers = "0.3"
-OrdinaryDiffEqBDF = "1"
-OrdinaryDiffEqDefault = "1"
-OrdinaryDiffEqRosenbrock = "1"
-OrdinaryDiffEqSDIRK = "1"
-OrdinaryDiffEqTsit5 = "1"
-OrdinaryDiffEqVerner = "1"
+OrdinaryDiffEqBDF = "1, 2"
+OrdinaryDiffEqDefault = "1, 2"
+OrdinaryDiffEqRosenbrock = "1, 2"
+OrdinaryDiffEqSDIRK = "1, 2"
+OrdinaryDiffEqTsit5 = "1, 2"
+OrdinaryDiffEqVerner = "1, 2"
 PEtab = "5"
 Plots = "1.40"
 QuasiMonteCarlo = "0.3"
 ReactionNetworkImporters = "1.1.0"
 SBMLImporter = "4.0.0"
-SciMLBase = "2.46"
+SciMLBase = "2.46, 3"
 SciMLSensitivity = "7.60"
 SpecialFunctions = "2.4"
 StaticArrays = "1.9"
 StatsPlots = "0.15.8"
 SteadyStateDiffEq = "2.2"
-StochasticDiffEq = "6.65"
+StochasticDiffEq = "6.65, 7"
 StructuralIdentifiability = "0.5.11"
 SymbolicIndexingInterface = "0.3"
 Symbolics = "7.13"

test/extensions/Project.toml

@@ -17,3 +17,8 @@ StableRNGs = "860ef19b-820b-49d6-a774-d7a799459cd3"
 SteadyStateDiffEq = "9672c7b4-1e72-59bd-8a11-6ac3964bc41f"
 StructuralIdentifiability = "220ca800-aa68-49bb-acd8-6037fa93a544"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
+
+[compat]
+OrdinaryDiffEqDefault = "1, 2"
+OrdinaryDiffEqTsit5 = "1, 2"
+OrdinaryDiffEqVerner = "1, 2"

test/simulation_and_solving/hybrid_models.jl

@@ -188,7 +188,7 @@ let
 
     # Solves and checks values of solution (do this before integrator mutation test
     # since integrator mutation affects shared parameter state).
-    sol = solve(prob, Tsit5(); maxiters = 10, verbose = false)
+    sol = solve(prob, Tsit5(); maxiters = 10)
     @test sol[:X][1] == 1.0
     @test sol[:X2][1] == 0.0
     @test sol[:Y][1] == 0.0
@@ -215,7 +215,7 @@ let
     @test prob2.ps[:d] == 0.5   # not changed
 
     # Test that we can solve the remade problem and verify solution values.
-    sol2 = solve(prob2, Tsit5(); maxiters = 10, verbose = false)
+    sol2 = solve(prob2, Tsit5(); maxiters = 10)
     @test sol2[:X][1] == 3.0
     @test sol2[:X2][1] == 2.0
     @test sol2[:Y][1] == 1.0

test/simulation_and_solving/simulate_SDEs.jl

@@ -618,8 +618,8 @@ let
     cat_sol = solve(cat_ensemble, EM(); dt = 0.001, trajectories = N, saveat = 20.0)
     hand_sol = solve(hand_ensemble, EM(); dt = 0.001, trajectories = N, saveat = 20.0)
 
-    cat_endpoints = [sol[1, end] for sol in cat_sol]
-    hand_endpoints = [sol[1, end] for sol in hand_sol]
+    cat_endpoints = [sol[1, end] for sol in cat_sol.u]
+    hand_endpoints = [sol[1, end] for sol in hand_sol.u]
 
     cat_mean = mean(cat_endpoints)
     hand_mean = mean(hand_endpoints)