Merge pull request #481 from ChrisRackauckas-Claude/bump-ordinarydiffeq-v7-ecosystem

Bump compat for OrdinaryDiffEq v7 / SciMLBase v3 ecosystem

Author: ChrisRackauckas

Project.toml

@@ -42,18 +42,18 @@ BoundaryValueDiffEqMIRK = "1"
 BoundaryValueDiffEqMIRKN = "1"
 BoundaryValueDiffEqShooting = "1"
 DiffEqBase = "6.183, 7"
-DiffEqDevTools = "2.48"
+DiffEqDevTools = "2.48, 3"
 FastClosures = "0.3.2"
 ForwardDiff = "0.10.38, 1"
 Hwloc = "3.3"
 InteractiveUtils = "<0.0.1, 1"
 LinearAlgebra = "1.10"
 LinearSolve = "3.12"
 NonlinearSolveBase = "2"
-NonlinearSolveFirstOrder = "1.3"
+NonlinearSolveFirstOrder = "1.3, 2"
 ODEInterface = "0.5"
-OrdinaryDiffEqLowOrderRK = "1.2"
-OrdinaryDiffEqTsit5 = "1.1"
+OrdinaryDiffEqLowOrderRK = "1.2, 2"
+OrdinaryDiffEqTsit5 = "1.1, 2"
 OptimizationBase = "5"
 Pkg = "1.10.0"
 Random = "1.10"

docs/Project.toml

@@ -46,7 +46,7 @@ LineSearch = "0.1.4"
 LinearAlgebra = "1.10"
 LinearSolve = "2.36.2, 3"
 OptimizationIpopt = "0.2.0, 1.0"
-OrdinaryDiffEq = "6.90.1"
+OrdinaryDiffEq = "7"
 Plots = "1"
 SciMLBase = "2.60.0, 3.3"
 SimpleBoundaryValueDiffEq = "1.1.0"

lib/BoundaryValueDiffEqAscher/Project.toml

@@ -24,7 +24,7 @@ AlmostBlockDiagonals = "0.1.10"
 Aqua = "0.8"
 BoundaryValueDiffEqCore = "2.0"
 ConcreteStructs = "0.2.3"
-DiffEqDevTools = "2.44"
+DiffEqDevTools = "2.44, 3"
 DifferentiationInterface = "0.7"
 FastClosures = "0.3.2"
 ForwardDiff = "0.10.38, 1"

lib/BoundaryValueDiffEqCore/src/utils.jl

@@ -149,8 +149,8 @@ end
         resid::Tuple, ::TwoPointBVProblem, (bca!, bcb!)::BC,
         sol::SciMLBase.ODESolution, p, t
     ) where {BC}
-    ua = first(sol)
-    ub = last(sol)
+    ua = sol.u[1]
+    ub = sol.u[end]
     bca!(resid[1], ua, p)
     bcb!(resid[2], ub, p)
     return resid
@@ -659,7 +659,7 @@ function __build_solution(prob::AbstractBVProblem, odesol, nlsol::SciMLBase.Nonl
 end
 function __build_solution(prob::AbstractBVProblem, odesol, optsol::SciMLBase.OptimizationSolution)
     retcode = ifelse(SciMLBase.successful_retcode(optsol), odesol.retcode, optsol.retcode)
-    return SciMLBase.solution_new_original_retcode(odesol, optsol, retcode, zeros(length(first(odesol)))) # Need a patch in SciMLBase
+    return SciMLBase.solution_new_original_retcode(odesol, optsol, retcode, zeros(length(first(odesol.u)))) # Need a patch in SciMLBase
 end
 
 # Fix3

lib/BoundaryValueDiffEqFIRK/Project.toml

@@ -34,7 +34,7 @@ ArrayInterface = "7.18"
 BandedMatrices = "1.7.5"
 BoundaryValueDiffEqCore = "2.2"
 ConcreteStructs = "0.2.3"
-DiffEqDevTools = "2.44"
+DiffEqDevTools = "2.44, 3"
 DifferentiationInterface = "0.7"
 Enzyme = "0.13.33"
 FastAlmostBandedMatrices = "0.1.4"
@@ -44,8 +44,8 @@ Hwloc = "3"
 InteractiveUtils = "<0.0.1, 1"
 LinearAlgebra = "1.10"
 LinearSolve = "2.36.2, 3"
-Mooncake = "0.4.146"
-OrdinaryDiffEqRosenbrock = "1.15.1"
+Mooncake = "0.4.146, 0.5"
+OrdinaryDiffEqRosenbrock = "1.15.1, 2"
 PreallocationTools = "0.4, 1"
 PrecompileTools = "1.2"
 Preferences = "1.4"

lib/BoundaryValueDiffEqMIRK/Project.toml

@@ -33,7 +33,7 @@ ArrayInterface = "7.18"
 BandedMatrices = "1.7.5"
 BoundaryValueDiffEqCore = "2.2"
 ConcreteStructs = "0.2.3"
-DiffEqDevTools = "2.44"
+DiffEqDevTools = "2.44, 3"
 DifferentiationInterface = "0.7"
 Enzyme = "0.13.33"
 FastAlmostBandedMatrices = "0.1.4"
@@ -43,9 +43,9 @@ Hwloc = "3"
 InteractiveUtils = "<0.0.1, 1"
 LinearAlgebra = "1.10"
 LinearSolve = "2.36.2, 3"
-Mooncake = "0.4"
+Mooncake = "0.4, 0.5"
 OptimizationIpopt = "1"
-OrdinaryDiffEqRosenbrock = "1"
+OrdinaryDiffEqRosenbrock = "1, 2"
 PreallocationTools = "0.4, 1"
 PrecompileTools = "1.2"
 Preferences = "1.4"

lib/BoundaryValueDiffEqMIRKN/Project.toml

@@ -33,7 +33,7 @@ ArrayInterface = "7.18"
 BandedMatrices = "1.7.5"
 BoundaryValueDiffEqCore = "2.0"
 ConcreteStructs = "0.2.3"
-DiffEqDevTools = "2.44"
+DiffEqDevTools = "2.44, 3"
 DifferentiationInterface = "0.7"
 FastAlmostBandedMatrices = "0.1.4"
 FastClosures = "0.3.2"

lib/BoundaryValueDiffEqShooting/Project.toml

@@ -32,17 +32,17 @@ ArrayInterface = "7.18"
 BandedMatrices = "1.11.0"
 BoundaryValueDiffEqCore = "2.0"
 ConcreteStructs = "0.2.3"
-DiffEqDevTools = "2.44"
+DiffEqDevTools = "2.44, 3"
 DifferentiationInterface = "0.7"
 FastClosures = "0.3.2"
 ForwardDiff = "0.10.38, 1"
 Hwloc = "3"
 InteractiveUtils = "<0.0.1, 1"
 LinearAlgebra = "1.10"
-OrdinaryDiffEqLowOrderRK = "1"
-OrdinaryDiffEqRosenbrock = "1"
-OrdinaryDiffEqTsit5 = "1"
-OrdinaryDiffEqVerner = "1"
+OrdinaryDiffEqLowOrderRK = "1, 2"
+OrdinaryDiffEqRosenbrock = "1, 2"
+OrdinaryDiffEqTsit5 = "1, 2"
+OrdinaryDiffEqVerner = "1, 2"
 PrecompileTools = "1.2"
 Preferences = "1.5"
 Pkg = "1.10.0"

lib/BoundaryValueDiffEqShooting/src/multiple_shooting.jl

@@ -43,7 +43,7 @@ function SciMLBase.__solve(
         M = length(bcresid_prototype)
     end
 
-    internal_ode_kwargs = (; verbose, kwargs..., odesolve_kwargs..., save_end = true)
+    internal_ode_kwargs = (; kwargs..., odesolve_kwargs..., save_end = true)
 
     solve_internal_odes! = @closure (
         resid_nodes,

lib/BoundaryValueDiffEqShooting/src/single_shooting.jl

@@ -30,7 +30,7 @@ function SciMLBase.__solve(
     resid_prototype = __vec(bcresid_prototype)
 
     # Construct the residual function
-    actual_ode_kwargs = (; kwargs..., verbose, odesolve_kwargs...)
+    actual_ode_kwargs = (; kwargs..., odesolve_kwargs...)
     # For TwoPointBVPs we don't need to save every step
     if prob.problem_type isa TwoPointBVProblem
         ode_kwargs = (; save_everystep = false, actual_ode_kwargs...)