Merge branch 'master' into qqy/more_AD

ErikQQY · web-flow · commit a0dd726f7ce0 · 2025-05-01T16:33:07.000+08:00
diff --git a/Project.toml b/Project.toml
@@ -1,6 +1,6 @@
 name = "BoundaryValueDiffEq"
 uuid = "764a87c0-6b3e-53db-9096-fe964310641d"
-version = "5.17.0"
+version = "5.16.1"
 
 [deps]
 ADTypes = "47edcb42-4c32-4615-8424-f2b9edc5f35b"
@@ -43,15 +43,16 @@ BoundaryValueDiffEqShooting = "1"
 DiffEqBase = "6.167"
 DiffEqDevTools = "2.48"
 FastClosures = "0.3.2"
-ForwardDiff = "0.10.38"
+ForwardDiff = "0.10.38, 1"
 Hwloc = "3.3"
 InteractiveUtils = "<0.0.1, 1"
 JET = "0.9.18"
 LinearAlgebra = "1.10"
 LinearSolve = "3.9"
 NonlinearSolveFirstOrder = "1.3"
 ODEInterface = "0.5"
-OrdinaryDiffEq = "6.93"
+OrdinaryDiffEqLowOrderRK = "1.2"
+OrdinaryDiffEqTsit5 = "1.1"
 Pkg = "1.10.0"
 Random = "1.10"
 RecursiveArrayTools = "3.31.2"
@@ -72,7 +73,8 @@ JET = "c3a54625-cd67-489e-a8e7-0a5a0ff4e31b"
 LinearSolve = "7ed4a6bd-45f5-4d41-b270-4a48e9bafcae"
 NonlinearSolveFirstOrder = "5959db7a-ea39-4486-b5fe-2dd0bf03d60d"
 ODEInterface = "54ca160b-1b9f-5127-a996-1867f4bc2a2c"
-OrdinaryDiffEq = "1dea7af3-3e70-54e6-95c3-0bf5283fa5ed"
+OrdinaryDiffEqLowOrderRK = "1344f307-1e59-4825-a18e-ace9aa3fa4c6"
+OrdinaryDiffEqTsit5 = "b1df2697-797e-41e3-8120-5422d3b24e4a"
 Pkg = "44cfe95a-1eb2-52ea-b672-e2afdf69b78f"
 Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 ReTestItems = "817f1d60-ba6b-4fd5-9520-3cf149f6a823"
@@ -82,4 +84,4 @@ StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 
 [targets]
-test = ["Aqua", "DiffEqDevTools", "Hwloc", "InteractiveUtils", "JET", "LinearSolve", "NonlinearSolveFirstOrder", "ODEInterface", "OrdinaryDiffEq", "Pkg", "Random", "ReTestItems", "RecursiveArrayTools", "Sparspak", "StaticArrays", "Test"]
+test = ["Aqua", "DiffEqDevTools", "Hwloc", "InteractiveUtils", "JET", "LinearSolve", "NonlinearSolveFirstOrder", "ODEInterface", "OrdinaryDiffEqLowOrderRK", "OrdinaryDiffEqTsit5", "Pkg", "Random", "ReTestItems", "RecursiveArrayTools", "Sparspak", "StaticArrays", "Test"]
diff --git a/lib/BoundaryValueDiffEqAscher/Project.toml b/lib/BoundaryValueDiffEqAscher/Project.toml
@@ -1,7 +1,7 @@
 name = "BoundaryValueDiffEqAscher"
 uuid = "7227322d-7511-4e07-9247-ad6ff830280e"
 authors = ["Qingyu Qu <erikqqy123@gmail.com>"]
-version = "1.5.0"
+version = "1.5.1"
 
 [deps]
 ADTypes = "47edcb42-4c32-4615-8424-f2b9edc5f35b"
@@ -29,7 +29,7 @@ DiffEqBase = "6.167"
 DiffEqDevTools = "2.44"
 DifferentiationInterface = "0.6.42"
 FastClosures = "0.3.2"
-ForwardDiff = "0.10.38"
+ForwardDiff = "0.10.38, 1"
 Hwloc = "3"
 InteractiveUtils = "<0.0.1, 1"
 JET = "0.9.18"
diff --git a/lib/BoundaryValueDiffEqCore/Project.toml b/lib/BoundaryValueDiffEqCore/Project.toml
@@ -1,7 +1,7 @@
 name = "BoundaryValueDiffEqCore"
 uuid = "56b672f2-a5fe-4263-ab2d-da677488eb3a"
 authors = ["Qingyu Qu <erikqqy123@gmail.com>"]
-version = "1.8.0"
+version = "1.8.1"
 
 [deps]
 ADTypes = "47edcb42-4c32-4615-8424-f2b9edc5f35b"
@@ -30,7 +30,7 @@ Aqua = "0.8"
 ArrayInterface = "7.18"
 ConcreteStructs = "0.2.3"
 DiffEqBase = "6.167"
-ForwardDiff = "0.10.38"
+ForwardDiff = "0.10.38, 1"
 InteractiveUtils = "<0.0.1, 1"
 JET = "0.9.18"
 LineSearch = "0.1.4"
diff --git a/lib/BoundaryValueDiffEqFIRK/Project.toml b/lib/BoundaryValueDiffEqFIRK/Project.toml
@@ -1,6 +1,6 @@
 name = "BoundaryValueDiffEqFIRK"
 uuid = "85d9eb09-370e-4000-bb32-543851f73618"
-version = "1.6.0"
+version = "1.6.1"
 
 [deps]
 ADTypes = "47edcb42-4c32-4615-8424-f2b9edc5f35b"
@@ -39,14 +39,14 @@ DifferentiationInterface = "0.6.42"
 Enzyme = "0.13.33"
 FastAlmostBandedMatrices = "0.1.4"
 FastClosures = "0.3.2"
-ForwardDiff = "0.10.38"
+ForwardDiff = "0.10.38, 1"
 Hwloc = "3"
 InteractiveUtils = "<0.0.1, 1"
 JET = "0.9.18"
 LinearAlgebra = "1.10"
 LinearSolve = "2.36.2, 3"
 Mooncake = "0.4.108"
-OrdinaryDiffEq = "6.93"
+OrdinaryDiffEqRosenbrock = "1"
 PreallocationTools = "0.4.24"
 PrecompileTools = "1.2"
 Preferences = "1.4"
@@ -70,12 +70,12 @@ InteractiveUtils = "b77e0a4c-d291-57a0-90e8-8db25a27a240"
 JET = "c3a54625-cd67-489e-a8e7-0a5a0ff4e31b"
 LinearSolve = "7ed4a6bd-45f5-4d41-b270-4a48e9bafcae"
 Mooncake = "da2b9cff-9c12-43a0-ae48-6db2b0edb7d6"
-OrdinaryDiffEq = "1dea7af3-3e70-54e6-95c3-0bf5283fa5ed"
+OrdinaryDiffEqRosenbrock = "43230ef6-c299-4910-a778-202eb28ce4ce"
 Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 ReTestItems = "817f1d60-ba6b-4fd5-9520-3cf149f6a823"
 RecursiveArrayTools = "731186ca-8d62-57ce-b412-fbd966d074cd"
 StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 
 [targets]
-test = ["Aqua", "DiffEqDevTools", "Enzyme", "Hwloc", "InteractiveUtils", "JET", "LinearSolve", "Mooncake", "OrdinaryDiffEq", "Random", "ReTestItems", "RecursiveArrayTools", "StaticArrays", "Test"]
+test = ["Aqua", "DiffEqDevTools", "Enzyme", "Hwloc", "InteractiveUtils", "JET", "LinearSolve", "OrdinaryDiffEqRosenbrock", "Mooncake", "Random", "ReTestItems", "RecursiveArrayTools", "StaticArrays", "Test"]
diff --git a/lib/BoundaryValueDiffEqFIRK/test/expanded/vectorofvector_initials_tests.jl b/lib/BoundaryValueDiffEqFIRK/test/expanded/vectorofvector_initials_tests.jl
@@ -1,5 +1,5 @@
 @testitem "VectorOfVector Initial Condition" begin
-    using BoundaryValueDiffEqFIRK, OrdinaryDiffEq
+    using BoundaryValueDiffEqFIRK, OrdinaryDiffEqRosenbrock
     #System Constants
     ss = 1 #excitatory parameter
     sj = 0 #inhibitory parameter
diff --git a/lib/BoundaryValueDiffEqFIRK/test/nested/vectorofvector_initials_tests.jl b/lib/BoundaryValueDiffEqFIRK/test/nested/vectorofvector_initials_tests.jl
@@ -1,5 +1,5 @@
 @testitem "VectorOfVector Initial Condition" begin
-    using BoundaryValueDiffEqFIRK, OrdinaryDiffEq
+    using BoundaryValueDiffEqFIRK, OrdinaryDiffEqRosenbrock
     #System Constants
     ss = 1 #excitatory parameter
     sj = 0 #inhibitory parameter
diff --git a/lib/BoundaryValueDiffEqMIRK/Project.toml b/lib/BoundaryValueDiffEqMIRK/Project.toml
@@ -1,6 +1,6 @@
 name = "BoundaryValueDiffEqMIRK"
 uuid = "1a22d4ce-7765-49ea-b6f2-13c8438986a6"
-version = "1.6.0"
+version = "1.6.1"
 
 [deps]
 ADTypes = "47edcb42-4c32-4615-8424-f2b9edc5f35b"
@@ -39,14 +39,14 @@ DifferentiationInterface = "0.6.42"
 Enzyme = "0.13.33"
 FastAlmostBandedMatrices = "0.1.4"
 FastClosures = "0.3.2"
-ForwardDiff = "0.10.38"
+ForwardDiff = "0.10.38, 1"
 Hwloc = "3"
 InteractiveUtils = "<0.0.1, 1"
 JET = "0.9.18"
 LinearAlgebra = "1.10"
 LinearSolve = "2.36.2, 3"
 Mooncake = "0.4.108"
-OrdinaryDiffEq = "6.93"
+OrdinaryDiffEqRosenbrock = "1"
 PreallocationTools = "0.4.24"
 PrecompileTools = "1.2"
 Preferences = "1.4"
@@ -70,12 +70,12 @@ InteractiveUtils = "b77e0a4c-d291-57a0-90e8-8db25a27a240"
 JET = "c3a54625-cd67-489e-a8e7-0a5a0ff4e31b"
 LinearSolve = "7ed4a6bd-45f5-4d41-b270-4a48e9bafcae"
 Mooncake = "da2b9cff-9c12-43a0-ae48-6db2b0edb7d6"
-OrdinaryDiffEq = "1dea7af3-3e70-54e6-95c3-0bf5283fa5ed"
+OrdinaryDiffEqRosenbrock = "43230ef6-c299-4910-a778-202eb28ce4ce"
 Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 ReTestItems = "817f1d60-ba6b-4fd5-9520-3cf149f6a823"
 RecursiveArrayTools = "731186ca-8d62-57ce-b412-fbd966d074cd"
 StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 
 [targets]
-test = ["Aqua", "DiffEqDevTools", "Enzyme", "Hwloc", "InteractiveUtils", "JET", "LinearSolve", "Mooncake", "OrdinaryDiffEq", "Random", "ReTestItems", "RecursiveArrayTools", "StaticArrays", "Test"]
+test = ["Aqua", "DiffEqDevTools", "Enzyme", "Hwloc", "InteractiveUtils", "JET", "LinearSolve", "Mooncake", "OrdinaryDiffEqRosenbrock", "Random", "ReTestItems", "RecursiveArrayTools", "StaticArrays", "Test"]
diff --git a/lib/BoundaryValueDiffEqMIRK/test/vectorofvector_initials_tests.jl b/lib/BoundaryValueDiffEqMIRK/test/vectorofvector_initials_tests.jl
@@ -1,5 +1,5 @@
 @testitem "VectorOfVector Initial Condition" begin
-    using BoundaryValueDiffEqMIRK, OrdinaryDiffEq
+    using BoundaryValueDiffEqMIRK, OrdinaryDiffEqRosenbrock
     #System Constants
     ss = 1 #excitatory parameter
     sj = 0 #inhibitory parameter
diff --git a/lib/BoundaryValueDiffEqMIRKN/Project.toml b/lib/BoundaryValueDiffEqMIRKN/Project.toml
@@ -1,7 +1,7 @@
 name = "BoundaryValueDiffEqMIRKN"
 uuid = "9255f1d6-53bf-473e-b6bd-23f1ff009da4"
 authors = ["Qingyu Qu <erikqqy123@gmail.com>"]
-version = "1.5.0"
+version = "1.5.1"
 
 [deps]
 ADTypes = "47edcb42-4c32-4615-8424-f2b9edc5f35b"
@@ -39,12 +39,11 @@ DiffEqDevTools = "2.44"
 DifferentiationInterface = "0.6.42"
 FastAlmostBandedMatrices = "0.1.4"
 FastClosures = "0.3.2"
-ForwardDiff = "0.10.38"
+ForwardDiff = "0.10.38, 1"
 Hwloc = "3"
 InteractiveUtils = "<0.0.1, 1"
 JET = "0.9.18"
 LinearAlgebra = "1.10"
-OrdinaryDiffEq = "6.93"
 PreallocationTools = "0.4.24"
 PrecompileTools = "1.2"
 Preferences = "1.4"
@@ -65,12 +64,11 @@ DiffEqDevTools = "f3b72e0c-5b89-59e1-b016-84e28bfd966d"
 Hwloc = "0e44f5e4-bd66-52a0-8798-143a42290a1d"
 InteractiveUtils = "b77e0a4c-d291-57a0-90e8-8db25a27a240"
 JET = "c3a54625-cd67-489e-a8e7-0a5a0ff4e31b"
-OrdinaryDiffEq = "1dea7af3-3e70-54e6-95c3-0bf5283fa5ed"
 Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 ReTestItems = "817f1d60-ba6b-4fd5-9520-3cf149f6a823"
 RecursiveArrayTools = "731186ca-8d62-57ce-b412-fbd966d074cd"
 StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 
 [targets]
-test = ["Aqua", "DiffEqDevTools", "Hwloc", "InteractiveUtils", "JET", "OrdinaryDiffEq", "Random", "ReTestItems", "RecursiveArrayTools", "StaticArrays", "Test"]
+test = ["Aqua", "DiffEqDevTools", "Hwloc", "InteractiveUtils", "JET", "Random", "ReTestItems", "RecursiveArrayTools", "StaticArrays", "Test"]
diff --git a/lib/BoundaryValueDiffEqShooting/Project.toml b/lib/BoundaryValueDiffEqShooting/Project.toml
@@ -1,6 +1,6 @@
 name = "BoundaryValueDiffEqShooting"
 uuid = "ed55bfe0-3725-4db6-871e-a1dc9f42a757"
-version = "1.6.0"
+version = "1.6.1"
 
 [deps]
 ADTypes = "47edcb42-4c32-4615-8424-f2b9edc5f35b"
@@ -38,12 +38,15 @@ DiffEqDevTools = "2.44"
 DifferentiationInterface = "0.6.42"
 FastAlmostBandedMatrices = "0.1.4"
 FastClosures = "0.3.2"
-ForwardDiff = "0.10.38"
+ForwardDiff = "0.10.38, 1"
 Hwloc = "3"
 InteractiveUtils = "<0.0.1, 1"
 JET = "0.9.18"
 LinearAlgebra = "1.10"
-OrdinaryDiffEq = "6.93"
+OrdinaryDiffEqLowOrderRK = "1"
+OrdinaryDiffEqRosenbrock = "1"
+OrdinaryDiffEqTsit5 = "1"
+OrdinaryDiffEqVerner = "1"
 Pkg = "1.10.0"
 PreallocationTools = "0.4.24"
 PrecompileTools = "1.2"
@@ -65,7 +68,10 @@ DiffEqDevTools = "f3b72e0c-5b89-59e1-b016-84e28bfd966d"
 Hwloc = "0e44f5e4-bd66-52a0-8798-143a42290a1d"
 InteractiveUtils = "b77e0a4c-d291-57a0-90e8-8db25a27a240"
 JET = "c3a54625-cd67-489e-a8e7-0a5a0ff4e31b"
-OrdinaryDiffEq = "1dea7af3-3e70-54e6-95c3-0bf5283fa5ed"
+OrdinaryDiffEqLowOrderRK = "1344f307-1e59-4825-a18e-ace9aa3fa4c6"
+OrdinaryDiffEqTsit5 = "b1df2697-797e-41e3-8120-5422d3b24e4a"
+OrdinaryDiffEqVerner = "79d7bb75-1356-48c1-b8c0-6832512096c2"
+OrdinaryDiffEqRosenbrock = "43230ef6-c299-4910-a778-202eb28ce4ce"
 Pkg = "44cfe95a-1eb2-52ea-b672-e2afdf69b78f"
 Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 ReTestItems = "817f1d60-ba6b-4fd5-9520-3cf149f6a823"
@@ -74,4 +80,4 @@ StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 
 [targets]
-test = ["Aqua", "DiffEqDevTools", "Hwloc", "InteractiveUtils", "JET", "OrdinaryDiffEq", "Pkg", "Random", "ReTestItems", "RecursiveArrayTools", "StaticArrays", "Test"]
+test = ["Aqua", "DiffEqDevTools", "Hwloc", "InteractiveUtils", "JET", "OrdinaryDiffEqLowOrderRK", "OrdinaryDiffEqRosenbrock", "OrdinaryDiffEqTsit5", "OrdinaryDiffEqVerner", "Pkg", "Random", "ReTestItems", "RecursiveArrayTools", "StaticArrays", "Test"]
diff --git a/lib/BoundaryValueDiffEqShooting/test/basic_problems_tests.jl b/lib/BoundaryValueDiffEqShooting/test/basic_problems_tests.jl
@@ -1,5 +1,5 @@
 @testitem "Basic Shooting" begin
-    using BoundaryValueDiffEqShooting, LinearAlgebra, OrdinaryDiffEq, JET
+    using BoundaryValueDiffEqShooting, LinearAlgebra, OrdinaryDiffEqTsit5, JET
 
     SOLVERS = [Shooting(Tsit5(), NewtonRaphson(; autodiff = AutoFiniteDiff())),
         Shooting(Tsit5(), NewtonRaphson(; autodiff = AutoForwardDiff(; chunksize = 2))),
@@ -126,7 +126,7 @@
 end
 
 @testitem "Shooting with Complex Values" begin
-    using BoundaryValueDiffEqShooting, OrdinaryDiffEq, LinearAlgebra
+    using BoundaryValueDiffEqShooting, OrdinaryDiffEqVerner, LinearAlgebra
 
     SOLVERS = [
         Shooting(Vern7(), NewtonRaphson(; autodiff = AutoFiniteDiff())), Shooting(Vern7()),
@@ -160,7 +160,8 @@ end
 end
 
 @testitem "Flow in a Channel" begin
-    using BoundaryValueDiffEqShooting, OrdinaryDiffEq, LinearAlgebra
+    using BoundaryValueDiffEqShooting, OrdinaryDiffEqTsit5, OrdinaryDiffEqRosenbrock,
+          LinearAlgebra
 
     function flow_in_a_channel!(du, u, p, t)
         R, P = p
@@ -209,7 +210,8 @@ end
 end
 #FIXME: MultipleShooting fails for large out-of-place BVP systems
 @testitem "Ray Tracing" begin
-    using BoundaryValueDiffEqShooting, OrdinaryDiffEq, LinearAlgebra
+    using BoundaryValueDiffEqShooting, OrdinaryDiffEqVerner, OrdinaryDiffEqRosenbrock,
+          LinearAlgebra
 
     @inline v(x, y, z, p) = 1 / (4 + cos(p[1] * x) + sin(p[2] * y) - cos(p[3] * z))
     @inline ux(x, y, z, p) = -p[1] * sin(p[1] * x)
diff --git a/lib/BoundaryValueDiffEqShooting/test/nlls_tests.jl b/lib/BoundaryValueDiffEqShooting/test/nlls_tests.jl
@@ -1,5 +1,5 @@
 @testitem "Overconstrained BVP" begin
-    using BoundaryValueDiffEqShooting, OrdinaryDiffEq, LinearAlgebra, JET
+    using BoundaryValueDiffEqShooting, OrdinaryDiffEqTsit5, LinearAlgebra, JET
 
     SOLVERS = [
         Shooting(Tsit5(), NewtonRaphson(),
diff --git a/lib/BoundaryValueDiffEqShooting/test/orbital_tests.jl b/lib/BoundaryValueDiffEqShooting/test/orbital_tests.jl
@@ -1,5 +1,5 @@
 @testitem "Lambert's Problem" begin
-    using BoundaryValueDiffEqShooting, OrdinaryDiffEq, LinearAlgebra
+    using BoundaryValueDiffEqShooting, OrdinaryDiffEqLowOrderRK, LinearAlgebra
 
     y0 = [-4.7763169762853989E+06, -3.8386398704441520E+05, -5.3500183933132319E+06,
         -5528.612564911408, 1216.8442360202787, 4845.114446429901]
diff --git a/test/misc/manifolds_tests.jl b/test/misc/manifolds_tests.jl
@@ -1,5 +1,5 @@
 @testitem "Manifolds.jl Integration" begin
-    using LinearAlgebra, OrdinaryDiffEq
+    using LinearAlgebra, OrdinaryDiffEqTsit5
 
     struct EmbeddedTorus
         R::Float64
diff --git a/test/misc/non_vector_input_tests.jl b/test/misc/non_vector_input_tests.jl
@@ -1,5 +1,5 @@
 @testitem "Non-Vector Inputs" begin
-    using LinearAlgebra, NonlinearSolveFirstOrder, OrdinaryDiffEq
+    using LinearAlgebra, NonlinearSolveFirstOrder, OrdinaryDiffEqTsit5
 
     for order in (2, 3, 4, 5, 6)
         s = Symbol("MIRK$(order)")
diff --git a/test/misc/type_stability_tests.jl b/test/misc/type_stability_tests.jl
@@ -1,5 +1,5 @@
 @testitem "Type Stability" begin
-    using LinearAlgebra, BoundaryValueDiffEq, OrdinaryDiffEq
+    using LinearAlgebra, BoundaryValueDiffEq, OrdinaryDiffEqTsit5
 
     f(u, p, t) = [p[1] * u[1] - p[2] * u[1] * u[2], p[3] * u[1] * u[2] - p[4] * u[2]]
     function f!(du, u, p, t)
diff --git a/test/wrappers/odeinterface_tests.jl b/test/wrappers/odeinterface_tests.jl
@@ -45,13 +45,13 @@ end
 
 # Just test that it runs. BVPSOL only works with linearly separable BCs.
 @testitem "BVPSOL" setup=[ODEInterfaceWrapperTestSetup] begin
-    using ODEInterface, OrdinaryDiffEq, RecursiveArrayTools, NonlinearSolveFirstOrder
+    using ODEInterface, OrdinaryDiffEqTsit5, RecursiveArrayTools, NonlinearSolveFirstOrder
 
     tpprob = TwoPointBVProblem(ex7_f!, (ex7_2pbc1!, ex7_2pbc2!), u0, tspan,
         p; bcresid_prototype = (zeros(1), zeros(1)))
 
     # Just generate a solution for bvpsol
-    sol_ms = solve(tpprob, MultipleShooting(10, DP5(), NewtonRaphson());
+    sol_ms = solve(tpprob, MultipleShooting(10, Tsit5(), NewtonRaphson());
         dt = π / 20, abstol = 1e-5, maxiters = 1000, adaptive = false)
 
     initial_u0 = [sol_ms(t) .+ rand() for t in tspan[1]:(π / 20):tspan[2]]
