diff --git a/lib/OptimizationAuglag/Project.toml b/lib/OptimizationAuglag/Project.toml
index bac3fd160..d03e4ddc6 100644
--- a/lib/OptimizationAuglag/Project.toml
+++ b/lib/OptimizationAuglag/Project.toml
@@ -30,8 +30,9 @@ LinearAlgebra = "1.10"
 MLUtils = "0.4.8"
 OptimizationBase = "5"
 OptimizationOptimisers = "0.3.16"
-Reexport = "1.2"
+Reexport = "1.2.2"
 SciMLBase = "2.122.1, 3"
+Random = "1.10"
 Test = "1.10.0"
 julia = "1.10"
 
diff --git a/lib/OptimizationAuglag/test/qa/qa.jl b/lib/OptimizationAuglag/test/qa/qa.jl
index f02133836..9932f62a5 100644
--- a/lib/OptimizationAuglag/test/qa/qa.jl
+++ b/lib/OptimizationAuglag/test/qa/qa.jl
@@ -2,7 +2,11 @@ using OptimizationAuglag, Aqua, JET
 using Test
 
 @testset "Aqua" begin
-    Aqua.test_all(OptimizationAuglag)
+    # OptimizationOptimisers is used in tests as the inner solver, not in src.
+    Aqua.test_all(
+        OptimizationAuglag;
+        stale_deps = (; ignore = [:OptimizationOptimisers])
+    )
 end
 
 @testset "JET static analysis" begin
diff --git a/lib/OptimizationBBO/Project.toml b/lib/OptimizationBBO/Project.toml
index 487d1b570..c79d67019 100644
--- a/lib/OptimizationBBO/Project.toml
+++ b/lib/OptimizationBBO/Project.toml
@@ -18,13 +18,15 @@ Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 OptimizationBase = {path = "../OptimizationBase"}
 
 [compat]
-Pkg = "1"
-SafeTestsets = "0.1"
-julia = "1.10"
 BlackBoxOptim = "0.6.3"
 OptimizationBase = "5"
-SciMLBase = "2.122.1, 3"
+Pkg = "1"
+Random = "1.10"
 Reexport = "1.2"
+SafeTestsets = "0.1"
+SciMLBase = "2.122.1, 3"
+Test = "1.10"
+julia = "1.10"
 
 [targets]
 test = ["Random", "Test", "SafeTestsets", "Pkg"]
diff --git a/lib/OptimizationBase/Project.toml b/lib/OptimizationBase/Project.toml
index a4270668e..61755fce5 100644
--- a/lib/OptimizationBase/Project.toml
+++ b/lib/OptimizationBase/Project.toml
@@ -10,7 +10,6 @@ DifferentiationInterface = "a0c0ee7d-e4b9-4e03-894e-1c5f64a51d63"
 DocStringExtensions = "ffbed154-4ef7-542d-bbb7-c09d3a79fcae"
 FastClosures = "9aa1b823-49e4-5ca5-8b0f-3971ec8bab6a"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
-PDMats = "90014a1f-27ba-587c-ab20-58faa44d9150"
 PrecompileTools = "aea7be01-6a6a-4083-8856-8a6e6704d82a"
 Reexport = "189a3867-3050-52da-a836-e630ba90ab69"
 SciMLBase = "0bca4576-84f4-4d90-8ffe-ffa030f20462"
@@ -56,7 +55,7 @@ BenchmarkTools = "1"
 ChainRulesCore = "1"
 ComponentArrays = ">= 0.13.9"
 DifferentiationInterface = "0.7.13"
-DocStringExtensions = "0.9"
+DocStringExtensions = "0.9.5"
 Enzyme = "0.13.2"
 FastClosures = "0.3"
 FiniteDiff = "2.12"
@@ -66,13 +65,13 @@ LinearAlgebra = "1.9, 1.10"
 Lux = "1.12"
 MLDataDevices = "1"
 MLUtils = "0.4"
-Manifolds = "0.10"
+Mooncake = "0.4.138, 0.5"
+Manifolds = "0.10, 0.11"
 ModelingToolkit = "11"
 Optim = ">= 1.4.1"
 Optimisers = ">= 0.2.5"
 OptimizationLBFGSB = "1.1"
 OptimizationManopt = "1.1"
-PDMats = "0.11"
 Pkg = "1"
 PrecompileTools = "1"
 Random = "1.10"
@@ -120,4 +119,4 @@ Tracker = "9f7883ad-71c0-57eb-9f7f-b5c9e6d3789c"
 Zygote = "e88e6eb3-aa80-5325-afca-941959d7151f"
 
 [targets]
-test = ["ADTypes", "Aqua", "BenchmarkTools", "ChainRulesCore", "ComponentArrays", "DifferentiationInterface", "Enzyme", "FiniteDiff", "ForwardDiff", "IterTools", "LinearAlgebra", "Lux", "MLUtils", "Manifolds", "ModelingToolkit", "Optim", "Optimisers", "OptimizationLBFGSB", "OptimizationManopt", "Pkg", "Random", "ReverseDiff", "SafeTestsets", "SciMLBase", "SparseArrays", "SparseConnectivityTracer", "SymbolicAnalysis", "Symbolics", "Test", "Tracker", "Zygote"]
+test = ["Aqua", "BenchmarkTools", "ChainRulesCore", "ComponentArrays", "Enzyme", "FiniteDiff", "ForwardDiff", "IterTools", "Lux", "MLUtils", "Manifolds", "ModelingToolkit", "Optim", "Optimisers", "OptimizationLBFGSB", "OptimizationManopt", "Pkg", "Random", "ReverseDiff", "SafeTestsets", "SymbolicAnalysis", "Symbolics", "Test", "Tracker", "Zygote"]
diff --git a/lib/OptimizationBase/src/OptimizationDISparseExt.jl b/lib/OptimizationBase/src/OptimizationDISparseExt.jl
index 3af95e9ed..392d82882 100644
--- a/lib/OptimizationBase/src/OptimizationDISparseExt.jl
+++ b/lib/OptimizationBase/src/OptimizationDISparseExt.jl
@@ -142,25 +142,25 @@ function instantiate_function(
         cons = let f = f, p = p
             (res, θ, p_call = p) -> f.cons(res, θ, p_call)
         end
+    end
 
-        function cons_oop(x)
-            _res = zeros(eltype(x), num_cons)
-            f.cons(_res, x, p)
-            return _res
-        end
+    function cons_oop(x)
+        _res = zeros(eltype(x), num_cons)
+        f.cons(_res, x, p)
+        return _res
+    end
 
-        function cons_oop(x, i)
-            _res = zeros(eltype(x), num_cons)
-            f.cons(_res, x, p)
-            return _res[i]
-        end
+    function cons_oop(x, i)
+        _res = zeros(eltype(x), num_cons)
+        f.cons(_res, x, p)
+        return _res[i]
+    end
 
-        function lagrangian(θ, σ, λ, p)
-            if eltype(θ) <: SparseConnectivityTracer.AbstractTracer || !iszero(θ)
-                return σ * f.f(θ, p) + dot(λ, cons_oop(θ))
-            else
-                return dot(λ, cons_oop(θ))
-            end
+    function lagrangian(θ, σ, λ, p)
+        if eltype(θ) <: SparseConnectivityTracer.AbstractTracer || !iszero(θ)
+            return σ * f.f(θ, p) + dot(λ, cons_oop(θ))
+        else
+            return dot(λ, cons_oop(θ))
         end
     end
 
@@ -452,10 +452,10 @@ function instantiate_function(
         cons = let f = f, p = p
             (x, p_call = p) -> f.cons(x, p_call)
         end
+    end
 
-        function lagrangian(θ, σ, λ, p)
-            return σ * f.f(θ, p) + dot(λ, f.cons(θ, p))
-        end
+    function lagrangian(θ, σ, λ, p)
+        return σ * f.f(θ, p) + dot(λ, f.cons(θ, p))
     end
 
     cons_jac_prototype = f.cons_jac_prototype
diff --git a/lib/OptimizationBase/src/function.jl b/lib/OptimizationBase/src/function.jl
index b86ed29ab..a2696f024 100644
--- a/lib/OptimizationBase/src/function.jl
+++ b/lib/OptimizationBase/src/function.jl
@@ -273,6 +273,7 @@ function instantiate_function(
     else
         min(open_nrmlbrkt_ind, open_squigllybrkt_ind)
     end
-    adpkg = adtypestr[strtind:(open_brkt_ind - 1)]
+    lastidx = isnothing(open_brkt_ind) ? lastindex(adtypestr) : (open_brkt_ind - 1)
+    adpkg = adtypestr[strtind:lastidx]
     throw(ArgumentError("The passed automatic differentiation backend choice is not available. Please load the corresponding AD package $adpkg."))
 end
diff --git a/lib/OptimizationBase/test/qa/qa.jl b/lib/OptimizationBase/test/qa/qa.jl
index d6abb0ed5..468e4dd1b 100644
--- a/lib/OptimizationBase/test/qa/qa.jl
+++ b/lib/OptimizationBase/test/qa/qa.jl
@@ -2,7 +2,18 @@ using OptimizationBase, Aqua, JET
 using Test
 
 @testset "Aqua" begin
-    Aqua.test_all(OptimizationBase)
+    # OptimizationBase defines the Optimization solve interface itself, so it
+    # extends SciMLBase's __init/__solve and CommonSolve's init/solve/solve! (the
+    # latter reexported through SciMLBase) on SciMLBase.OptimizationProblem /
+    # AbstractOptimizationCache. Those are our *own* interface functions, so mark
+    # them as own for the piracy check rather than flagging the SciML types.
+    SB = OptimizationBase.SciMLBase
+    Aqua.test_all(
+        OptimizationBase;
+        piracies = (
+            treat_as_own = [SB.__init, SB.__solve, SB.init, SB.solve, SB.solve!],
+        )
+    )
 end
 
 @testset "JET static analysis" begin
diff --git a/lib/OptimizationEvolutionary/Project.toml b/lib/OptimizationEvolutionary/Project.toml
index a7a4fd824..bcccd9ff7 100644
--- a/lib/OptimizationEvolutionary/Project.toml
+++ b/lib/OptimizationEvolutionary/Project.toml
@@ -18,13 +18,15 @@ Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 OptimizationBase = {path = "../OptimizationBase"}
 
 [compat]
+Evolutionary = "0.11"
+OptimizationBase = "5"
 Pkg = "1"
+Random = "1.10"
+Reexport = "1.2"
 SafeTestsets = "0.1"
-julia = "1.10"
-OptimizationBase = "5"
-Evolutionary = "0.11"
 SciMLBase = "2.122.1, 3"
-Reexport = "1.2"
+Test = "1.10"
+julia = "1.10"
 
 [targets]
 test = ["Random", "Test", "SafeTestsets", "Pkg"]
diff --git a/lib/OptimizationEvolutionary/test/qa/qa.jl b/lib/OptimizationEvolutionary/test/qa/qa.jl
index 6aaf53bf8..16a740881 100644
--- a/lib/OptimizationEvolutionary/test/qa/qa.jl
+++ b/lib/OptimizationEvolutionary/test/qa/qa.jl
@@ -1,8 +1,29 @@
 using OptimizationEvolutionary, Aqua, JET
 using Test
+using Evolutionary
 
 @testset "Aqua" begin
-    Aqua.test_all(OptimizationEvolutionary)
+    # SciML trait/interface methods are our own, not piracy — mark them as such.
+    # The Evolutionary.trace! override IS genuine piracy (changes Evolutionary's
+    # tracing globally); mark the piracy test broken until it's replaced.
+    SB = OptimizationEvolutionary.SciMLBase
+    Aqua.test_all(
+        OptimizationEvolutionary;
+        piracies = (
+            broken = true,
+            treat_as_own = [
+                SB.__solve,
+                SB.allowsbounds,
+                SB.allowscallback,
+                SB.allowsconstraints,
+                SB.has_init,
+                SB.requiresconshess,
+                SB.requiresconsjac,
+                SB.requiresgradient,
+                SB.requireshessian,
+            ],
+        )
+    )
 end
 
 @testset "JET static analysis" begin
diff --git a/lib/OptimizationGCMAES/Project.toml b/lib/OptimizationGCMAES/Project.toml
index 310f90064..5bbc03bbd 100644
--- a/lib/OptimizationGCMAES/Project.toml
+++ b/lib/OptimizationGCMAES/Project.toml
@@ -18,13 +18,15 @@ Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 OptimizationBase = {path = "../OptimizationBase"}
 
 [compat]
+ForwardDiff = "0.10, 1"
+GCMAES = "0.1.34"
+OptimizationBase = "5.1"
 Pkg = "1"
+Reexport = "1.2"
 SafeTestsets = "0.1"
-julia = "1.10"
-OptimizationBase = "5.1"
 SciMLBase = "2.122.1, 3"
-Reexport = "1.2"
-GCMAES = "0.1.34"
+Test = "1.10"
+julia = "1.10"
 
 [targets]
 test = ["ForwardDiff", "Test", "SafeTestsets", "Pkg"]
diff --git a/lib/OptimizationGCMAES/src/OptimizationGCMAES.jl b/lib/OptimizationGCMAES/src/OptimizationGCMAES.jl
index 3447b345f..dc33bf30e 100644
--- a/lib/OptimizationGCMAES/src/OptimizationGCMAES.jl
+++ b/lib/OptimizationGCMAES/src/OptimizationGCMAES.jl
@@ -72,11 +72,9 @@ function SciMLBase.__solve(cache::OptimizationCache{O}) where {O <: GCMAESOpt}
         return x[1]
     end
 
-    if !isnothing(cache.f.grad)
-        g = function (θ)
-            cache.f.grad(G, θ)
-            return G
-        end
+    g = function (θ)
+        cache.f.grad(G, θ)
+        return G
     end
 
     maxiters = OptimizationBase._check_and_convert_maxiters(cache.solver_args.maxiters)
diff --git a/lib/OptimizationIpopt/Project.toml b/lib/OptimizationIpopt/Project.toml
index b474e720c..6ed49aa80 100644
--- a/lib/OptimizationIpopt/Project.toml
+++ b/lib/OptimizationIpopt/Project.toml
@@ -12,16 +12,21 @@ SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
 SymbolicIndexingInterface = "2efcf032-c050-4f8e-a9bb-153293bab1f5"
 
 [compat]
-Pkg = "1"
-SafeTestsets = "0.1"
+Aqua = "0.8"
 Ipopt = "1.10.3"
 LinearAlgebra = "1.10.0"
 ModelingToolkit = "11"
 OptimizationBase = "5.1"
+Pkg = "1"
+Random = "1.10"
 Reexport = "1.2"
+ReverseDiff = "1"
+SafeTestsets = "0.1"
 SciMLBase = "2.122.1, 3"
 SparseArrays = "1.10.0"
-SymbolicIndexingInterface = "0.3.40"
+SymbolicIndexingInterface = "0.3.46"
+Symbolics = "7"
+Test = "1.10"
 Zygote = "0.7"
 julia = "1.10"
 
diff --git a/lib/OptimizationLBFGSB/Project.toml b/lib/OptimizationLBFGSB/Project.toml
index 0297327b3..62e75b8a4 100644
--- a/lib/OptimizationLBFGSB/Project.toml
+++ b/lib/OptimizationLBFGSB/Project.toml
@@ -30,6 +30,7 @@ MLUtils = "0.4.8"
 OptimizationBase = "5"
 Reexport = "1.2"
 SciMLBase = "2.122.1, 3"
+Test = "1.10"
 Zygote = "0.7.10"
 julia = "1.10"
 
diff --git a/lib/OptimizationLBFGSB/src/OptimizationLBFGSB.jl b/lib/OptimizationLBFGSB/src/OptimizationLBFGSB.jl
index 100a8029e..abf62ffc0 100644
--- a/lib/OptimizationLBFGSB/src/OptimizationLBFGSB.jl
+++ b/lib/OptimizationLBFGSB/src/OptimizationLBFGSB.jl
@@ -160,6 +160,7 @@ function SciMLBase.__solve(cache::OptimizationCache{O}) where {O <: LBFGSB}
         end
 
         opt_ret = ReturnCode.MaxIters
+        res = (cache.u0, nothing)
         n = length(cache.u0)
 
         if cache.lb === nothing
diff --git a/lib/OptimizationMOI/Project.toml b/lib/OptimizationMOI/Project.toml
index 950913362..bd095cc48 100644
--- a/lib/OptimizationMOI/Project.toml
+++ b/lib/OptimizationMOI/Project.toml
@@ -31,7 +31,7 @@ Ipopt = "1.10.4"
 Ipopt_jll = "300.1400"
 Juniper = "0.9"
 LinearAlgebra = "1"
-MathOptInterface = "1.40.2"
+MathOptInterface = "1.45"
 ModelingToolkit = "11"
 ModelingToolkitBase = "1"
 NLopt = "1"
@@ -41,7 +41,7 @@ ReverseDiff = "1"
 SciMLBase = "2.130, 3"
 SciMLStructures = "1"
 SparseArrays = "1.6"
-SymbolicIndexingInterface = "0.3"
+SymbolicIndexingInterface = "0.3.46"
 SymbolicUtils = "4.10"
 Symbolics = "7"
 Test = "1.6"
diff --git a/lib/OptimizationMOI/test/qa/qa.jl b/lib/OptimizationMOI/test/qa/qa.jl
index 67241f196..3d7cfdac8 100644
--- a/lib/OptimizationMOI/test/qa/qa.jl
+++ b/lib/OptimizationMOI/test/qa/qa.jl
@@ -1,15 +1,33 @@
 using OptimizationMOI, Aqua, JET
 using Test
+using MathOptInterface
 
 @testset "Aqua" begin
-    # OptimizationMOI deliberately extends SciMLBase's solver-trait interface onto the
-    # MathOptInterface optimizer types, so those methods are intentional, not piracy.
+    # OptimizationMOI implements the SciML optimization interface for
+    # MathOptInterface, so the trait/interface methods it adds extend SciML's *own*
+    # functions rather than committing type piracy — mark those functions as own.
+    SB = OptimizationMOI.SciMLBase
+    OB = OptimizationMOI.OptimizationBase
     Aqua.test_all(
         OptimizationMOI;
         piracies = (
             treat_as_own = [
-                OptimizationMOI.MOI.AbstractOptimizer,
-                OptimizationMOI.MOI.OptimizerWithAttributes,
+                SB.__init,
+                SB.__solve,
+                SB.allowsbounds,
+                SB.allowscallback,
+                SB.allowsconstraints,
+                SB.get_observed,
+                SB.get_p,
+                SB.get_paramsyms,
+                SB.get_syms,
+                SB.has_init,
+                SB.requiresconshess,
+                SB.requiresconsjac,
+                SB.requiresgradient,
+                SB.requireshessian,
+                SB.supports_opt_cache_interface,
+                OB.supports_sense,
             ],
         )
     )
diff --git a/lib/OptimizationMadNLP/Project.toml b/lib/OptimizationMadNLP/Project.toml
index 727c48dc1..24b218329 100644
--- a/lib/OptimizationMadNLP/Project.toml
+++ b/lib/OptimizationMadNLP/Project.toml
@@ -14,9 +14,9 @@ SciMLBase = "0bca4576-84f4-4d90-8ffe-ffa030f20462"
 SymbolicIndexingInterface = "2efcf032-c050-4f8e-a9bb-153293bab1f5"
 
 [compat]
-Pkg = "1"
-SafeTestsets = "0.1"
-DifferentiationInterface = "0.7"
+ADTypes = "1.14"
+Aqua = "0.8"
+DifferentiationInterface = "0.7.13"
 ForwardDiff = "1.2.1"
 LinearAlgebra = "1.10.0"
 MadNLP = "0.9.1, 0.10"
@@ -24,9 +24,16 @@ ModelingToolkit = "11"
 NLPModels = "0.21.5"
 OptimizationBase = "5.1"
 OptimizationNLPModels = "1.3"
+Pkg = "1"
+Random = "1.10"
 Reexport = "1.2"
-SciMLBase = "2.122.1, 3"
-SymbolicIndexingInterface = "0.3.40"
+ReverseDiff = "1"
+SafeTestsets = "0.1"
+SciMLBase = "2.130, 3"
+SparseArrays = "1.10"
+SymbolicIndexingInterface = "0.3.46"
+Symbolics = "7"
+Test = "1.10"
 Zygote = "0.7"
 julia = "1.10"
 
diff --git a/lib/OptimizationMadNLP/test/qa/qa.jl b/lib/OptimizationMadNLP/test/qa/qa.jl
index e0566ce0f..846c95e3f 100644
--- a/lib/OptimizationMadNLP/test/qa/qa.jl
+++ b/lib/OptimizationMadNLP/test/qa/qa.jl
@@ -2,7 +2,14 @@ using OptimizationMadNLP, Aqua, JET
 using Test
 
 @testset "Aqua" begin
-    Aqua.test_all(OptimizationMadNLP)
+    # `@reexport using OptimizationBase` brings in SciMLBase's `solve!` while
+    # `using MadNLP` also pulls a `solve!` into scope; the clash leaves
+    # OptimizationMadNLP's `solve!` export pointing at neither binding.
+    # Mark broken until the reexport is restructured.
+    Aqua.test_all(
+        OptimizationMadNLP;
+        undefined_exports = (; broken = true)
+    )
 end
 
 @testset "JET static analysis" begin
diff --git a/lib/OptimizationManopt/Project.toml b/lib/OptimizationManopt/Project.toml
index d4079d768..3cb44a202 100644
--- a/lib/OptimizationManopt/Project.toml
+++ b/lib/OptimizationManopt/Project.toml
@@ -6,10 +6,10 @@ version = "1.3.3"
 Dates = "ade2ca70-3891-5945-98fb-dc099432e06a"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 ManifoldDiff = "af67fdf4-a580-4b9f-bbec-742ef357defd"
-Manifolds = "1cead3c2-87b3-11e9-0ccd-23c62b72b94e"
 ManifoldsBase = "3362f125-f0bb-47a3-aa74-596ffd7ef2fb"
 Manopt = "0fc0a36d-df90-57f3-8f93-d78a9fc72bb5"
 OptimizationBase = "bca83a33-5cc9-4baa-983d-23429ab6bcbb"
+Manifolds = "1cead3c2-87b3-11e9-0ccd-23c62b72b94e"
 Reexport = "189a3867-3050-52da-a836-e630ba90ab69"
 SciMLBase = "0bca4576-84f4-4d90-8ffe-ffa030f20462"
 
@@ -34,8 +34,17 @@ OptimizationBase = {path = "../OptimizationBase"}
 
 [compat]
 Aqua = "0.8"
+Pkg = "1"
+QuadraticModels = "0.9"
+Random = "1.10"
+ReverseDiff = "1.14"
+RipQP = "0.7"
+SafeTestsets = "0.1"
+Test = "1.10"
+Zygote = "0.6.67, 0.7"
+julia = "1.10"
 Dates = "1.10"
-DifferentiationInterface = "0.7"
+DifferentiationInterface = "0.7.13"
 Enzyme = "0.13"
 FiniteDiff = "2"
 ForwardDiff = "0.10, 1"
@@ -46,17 +55,8 @@ Manifolds = "0.10, 0.11"
 ManifoldsBase = "1, 2"
 Manopt = "0.5.25"
 OptimizationBase = "5.1"
-Pkg = "1"
-QuadraticModels = "0.9"
-Random = "1.10"
 Reexport = "1.2"
-ReverseDiff = "1"
-RipQP = "0.6"
-SafeTestsets = "0.1"
 SciMLBase = "2.122.1, 3"
-Test = "1.10"
-Zygote = "0.6, 0.7"
-julia = "1.10"
 
 [targets]
 test = ["DifferentiationInterface", "Enzyme", "ForwardDiff", "FiniteDiff", "QuadraticModels", "Random", "ReverseDiff", "RipQP", "Test", "Zygote", "SafeTestsets", "Pkg"]
diff --git a/lib/OptimizationMetaheuristics/Project.toml b/lib/OptimizationMetaheuristics/Project.toml
index a440a6945..a1b9a1b24 100644
--- a/lib/OptimizationMetaheuristics/Project.toml
+++ b/lib/OptimizationMetaheuristics/Project.toml
@@ -25,6 +25,8 @@ OptimizationBase = "5"
 Metaheuristics = "3.3"
 SciMLBase = "2.122.1, 3"
 Reexport = "1.2"
+Random = "1.10"
+Test = "1.10"
 
 [targets]
 test = ["Random", "Test", "SafeTestsets", "Pkg"]
diff --git a/lib/OptimizationMetaheuristics/src/OptimizationMetaheuristics.jl b/lib/OptimizationMetaheuristics/src/OptimizationMetaheuristics.jl
index 596d29a5a..a92e0ffbb 100644
--- a/lib/OptimizationMetaheuristics/src/OptimizationMetaheuristics.jl
+++ b/lib/OptimizationMetaheuristics/src/OptimizationMetaheuristics.jl
@@ -107,6 +107,7 @@ function SciMLBase.__solve(cache::OptimizationCache{O}) where {
         end
     end
 
+    opt_bounds = nothing
     if !isnothing(cache.lb) & !isnothing(cache.ub)
         opt_bounds = [cache.lb cache.ub]'
     end
diff --git a/lib/OptimizationMetaheuristics/test/qa/qa.jl b/lib/OptimizationMetaheuristics/test/qa/qa.jl
index f3a5bbd8b..78663cab8 100644
--- a/lib/OptimizationMetaheuristics/test/qa/qa.jl
+++ b/lib/OptimizationMetaheuristics/test/qa/qa.jl
@@ -1,8 +1,28 @@
 using OptimizationMetaheuristics, Aqua, JET
 using Test
+using Metaheuristics
 
 @testset "Aqua" begin
-    Aqua.test_all(OptimizationMetaheuristics)
+    # OptimizationMetaheuristics implements the SciML optimization interface for
+    # Metaheuristics, so the trait/interface methods it adds extend SciML's *own*
+    # functions rather than committing type piracy — mark those functions as own.
+    # `@reexport using Metaheuristics` exports `solve!`, which clashes with SciMLBase's
+    # `solve!` brought in transitively via OptimizationBase; mark broken until restructured.
+    SB = OptimizationMetaheuristics.SciMLBase
+    Aqua.test_all(
+        OptimizationMetaheuristics;
+        piracies = (
+            treat_as_own = [
+                SB.__init,
+                SB.__solve,
+                SB.allowsbounds,
+                SB.allowscallback,
+                SB.has_init,
+                SB.requiresbounds,
+            ],
+        ),
+        undefined_exports = (; broken = true)
+    )
 end
 
 @testset "JET static analysis" begin
diff --git a/lib/OptimizationMultistartOptimization/Project.toml b/lib/OptimizationMultistartOptimization/Project.toml
index 4ee7d4a33..e9d9ddbbc 100644
--- a/lib/OptimizationMultistartOptimization/Project.toml
+++ b/lib/OptimizationMultistartOptimization/Project.toml
@@ -28,6 +28,10 @@ OptimizationBase = "5"
 MultistartOptimization = "0.3"
 SciMLBase = "2.122.1, 3"
 Reexport = "1.2"
+ForwardDiff = "0.10, 1"
+OptimizationNLopt = "0.3"
+ReverseDiff = "1"
+Test = "1.10"
 
 [targets]
 test = ["ForwardDiff", "OptimizationNLopt", "ReverseDiff", "Pkg", "Test", "SafeTestsets"]
diff --git a/lib/OptimizationMultistartOptimization/test/qa/qa.jl b/lib/OptimizationMultistartOptimization/test/qa/qa.jl
index cad8a7cd8..2eff8b8b8 100644
--- a/lib/OptimizationMultistartOptimization/test/qa/qa.jl
+++ b/lib/OptimizationMultistartOptimization/test/qa/qa.jl
@@ -1,8 +1,28 @@
 using OptimizationMultistartOptimization, Aqua, JET
 using Test
+using MultistartOptimization
 
 @testset "Aqua" begin
-    Aqua.test_all(OptimizationMultistartOptimization)
+    # OptimizationMultistartOptimization implements the SciML optimization interface
+    # for MultistartOptimization, so the trait/interface methods it adds extend
+    # SciML's *own* functions rather than committing type piracy — mark those
+    # functions as own.
+    SB = OptimizationMultistartOptimization.SciMLBase
+    Aqua.test_all(
+        OptimizationMultistartOptimization;
+        piracies = (
+            treat_as_own = [
+                SB.__init,
+                SB.__solve,
+                SB.allowsbounds,
+                SB.allowscallback,
+                SB.has_init,
+                SB.requiresbounds,
+            ],
+        ),
+        # OptimizationNLopt is used in tests as the inner solver, not in src.
+        stale_deps = (; ignore = [:OptimizationNLopt])
+    )
 end
 
 @testset "JET static analysis" begin
diff --git a/lib/OptimizationNLPModels/Project.toml b/lib/OptimizationNLPModels/Project.toml
index 5657313d6..78433dfd4 100644
--- a/lib/OptimizationNLPModels/Project.toml
+++ b/lib/OptimizationNLPModels/Project.toml
@@ -18,14 +18,22 @@ OptimizationMOI = {path = "../OptimizationMOI"}
 OptimizationOptimJL = {path = "../OptimizationOptimJL"}
 
 [compat]
-Pkg = "1"
-SafeTestsets = "0.1"
 ADTypes = "1.18"
+Ipopt = "1"
 NLPModels = "0.21"
+NLPModelsTest = "0.10"
 OptimizationBase = "5"
+OptimizationLBFGSB = "1.2"
+OptimizationMOI = "1"
+OptimizationOptimJL = "0.4.10"
+Pkg = "1"
 Reexport = "1.2"
-SciMLBase = "2.122.1, 3"
+ReverseDiff = "1"
+SafeTestsets = "0.1"
+SciMLBase = "2.130, 3"
 SparseArrays = "1.10.0"
+Test = "1.10"
+Zygote = "0.6, 0.7"
 julia = "1.10"
 
 [extras]
diff --git a/lib/OptimizationNLPModels/test/qa/qa.jl b/lib/OptimizationNLPModels/test/qa/qa.jl
index b53d7d363..5eff5628d 100644
--- a/lib/OptimizationNLPModels/test/qa/qa.jl
+++ b/lib/OptimizationNLPModels/test/qa/qa.jl
@@ -1,8 +1,21 @@
 using OptimizationNLPModels, Aqua, JET
 using Test
+using NLPModels
 
 @testset "Aqua" begin
-    Aqua.test_all(OptimizationNLPModels)
+    # This package defines SciMLBase.OptimizationFunction / SciMLBase.OptimizationProblem
+    # constructors for NLPModels.AbstractNLPModel. Those constructors extend SciML's
+    # *own* types, so mark those (not the NLPModels type) as own for the piracy check.
+    SB = OptimizationNLPModels.SciMLBase
+    Aqua.test_all(
+        OptimizationNLPModels;
+        piracies = (
+            treat_as_own = [SB.OptimizationFunction, SB.OptimizationProblem],
+        ),
+        # Optimization sublibraries used to construct OptimizationProblems in tests
+        # but not `using`d in src; Aqua's static analysis can't see test-only usage.
+        stale_deps = (; ignore = [:OptimizationLBFGSB, :OptimizationMOI, :OptimizationOptimJL])
+    )
 end
 
 @testset "JET static analysis" begin
diff --git a/lib/OptimizationNLopt/Project.toml b/lib/OptimizationNLopt/Project.toml
index a513cf9c3..b482d5e71 100644
--- a/lib/OptimizationNLopt/Project.toml
+++ b/lib/OptimizationNLopt/Project.toml
@@ -28,6 +28,10 @@ OptimizationBase = "5.1"
 NLopt = "1.1"
 SciMLBase = "2.122.1, 3"
 Reexport = "1.2"
+Random = "1.10"
+ReverseDiff = "1"
+Test = "1.10"
+Zygote = "0.6, 0.7"
 
 [targets]
 test = ["ReverseDiff", "Test", "Zygote", "SafeTestsets", "Pkg"]
diff --git a/lib/OptimizationNLopt/src/OptimizationNLopt.jl b/lib/OptimizationNLopt/src/OptimizationNLopt.jl
index d0da6da12..c3b2767c3 100644
--- a/lib/OptimizationNLopt/src/OptimizationNLopt.jl
+++ b/lib/OptimizationNLopt/src/OptimizationNLopt.jl
@@ -196,14 +196,9 @@ function SciMLBase.__solve(cache::OptimizationCache{O}) where {
         eqinds = map((y) -> y[1] == y[2], zip(cache.lcons, cache.ucons))
         ineqinds = map((y) -> y[1] != y[2], zip(cache.lcons, cache.ucons))
         cons_cache = zeros(eltype(cache.u0), sum(eqinds) + sum(ineqinds))
-        thetacache = rand(size(cache.u0))
-        Jthetacache = rand(size(cache.u0))
         Jcache = zeros(eltype(cache.u0), sum(ineqinds) + sum(eqinds), length(cache.u0))
         evalcons = function (θ, ineqoreq)
-            if thetacache != θ
-                cache.f.cons(cons_cache, θ)
-                thetacache = copy(θ)
-            end
+            cache.f.cons(cons_cache, θ)
             if ineqoreq == :eq
                 return @view(cons_cache[eqinds])
             else
@@ -212,11 +207,7 @@ function SciMLBase.__solve(cache::OptimizationCache{O}) where {
         end
 
         evalconj = function (θ, ineqoreq)
-            if Jthetacache != θ
-                cache.f.cons_j(Jcache, θ)
-                Jthetacache = copy(θ)
-            end
-
+            cache.f.cons_j(Jcache, θ)
             if ineqoreq == :eq
                 return @view(Jcache[eqinds, :])'
             else
diff --git a/lib/OptimizationNLopt/test/qa/qa.jl b/lib/OptimizationNLopt/test/qa/qa.jl
index 7c01565ea..ffb0c53df 100644
--- a/lib/OptimizationNLopt/test/qa/qa.jl
+++ b/lib/OptimizationNLopt/test/qa/qa.jl
@@ -1,8 +1,34 @@
 using OptimizationNLopt, Aqua, JET
 using Test
+using NLopt
 
 @testset "Aqua" begin
-    Aqua.test_all(OptimizationNLopt)
+    # OptimizationNLopt implements the SciML optimization interface for NLopt, so
+    # the trait/interface methods it adds extend SciML's *own* functions rather
+    # than committing type piracy — mark those functions as own for Aqua's piracy
+    # check. NLopt.Algorithm is also kept for the `(::NLopt.Algorithm)()`
+    # normalization method, which extends NLopt's type directly and has no SciML
+    # function to attribute it to.
+    SB = OptimizationNLopt.SciMLBase
+    OB = OptimizationNLopt.OptimizationBase
+    Aqua.test_all(
+        OptimizationNLopt;
+        piracies = (
+            treat_as_own = [
+                SB.__init,
+                SB.__solve,
+                SB.allowsbounds,
+                SB.allowscallback,
+                SB.allowsconstraints,
+                SB.has_init,
+                SB.requiresconsjac,
+                SB.requiresgradient,
+                SB.requireshessian,
+                OB.supports_sense,
+                NLopt.Algorithm,
+            ],
+        )
+    )
 end
 
 @testset "JET static analysis" begin
diff --git a/lib/OptimizationNOMAD/Project.toml b/lib/OptimizationNOMAD/Project.toml
index a6b586d69..8e705e4c4 100644
--- a/lib/OptimizationNOMAD/Project.toml
+++ b/lib/OptimizationNOMAD/Project.toml
@@ -22,8 +22,9 @@ SafeTestsets = "0.1"
 julia = "1.10"
 OptimizationBase = "5"
 NOMAD = "2.4.1"
-SciMLBase = "2.58, 3"
+SciMLBase = "2.122.1, 3"
 Reexport = "1.2"
+Test = "1.10"
 
 [targets]
 test = ["Test", "SafeTestsets", "Pkg"]
diff --git a/lib/OptimizationNOMAD/src/OptimizationNOMAD.jl b/lib/OptimizationNOMAD/src/OptimizationNOMAD.jl
index aaf107a48..04a401f34 100644
--- a/lib/OptimizationNOMAD/src/OptimizationNOMAD.jl
+++ b/lib/OptimizationNOMAD/src/OptimizationNOMAD.jl
@@ -74,16 +74,27 @@ function SciMLBase.__solve(
         return first(x)
     end
 
+    bb = nothing
+    bounds = (;)
+    if !isnothing(prob.lb)
+        bounds = (; bounds..., lower_bound = prob.lb)
+    end
+
+    if !isnothing(prob.ub)
+        bounds = (; bounds..., upper_bound = prob.ub)
+    end
+
     if prob.f.cons === nothing
-        function bb(x)
+        bb = function (x)
             l = _loss(x)
             success = !isnan(l) && !isinf(l)
             count_eval = true
             return (success, count_eval, [l])
         end
+        opt_setup = NOMAD.NomadProblem(length(prob.u0), 1, ["OBJ"], bb; bounds...)
     else
         eqinds = findall(i -> prob.lcons[i] == prob.ucons[i], 1:length(prob.ucons))
-        function bbcons(x)
+        bbcons = function (x)
             l = _loss(x)
             c = zeros(eltype(x), length(prob.ucons))
             prob.f.cons(c, x, prob.p)
@@ -95,20 +106,6 @@ function SciMLBase.__solve(
             count_eval = true
             return (success, count_eval, vcat(l, c))
         end
-    end
-
-    bounds = (;)
-    if !isnothing(prob.lb)
-        bounds = (; bounds..., lower_bound = prob.lb)
-    end
-
-    if !isnothing(prob.ub)
-        bounds = (; bounds..., upper_bound = prob.ub)
-    end
-
-    if prob.f.cons === nothing
-        opt_setup = NOMAD.NomadProblem(length(prob.u0), 1, ["OBJ"], bb; bounds...)
-    else
         opt_setup = NOMAD.NomadProblem(
             length(prob.u0), 1 + length(prob.ucons),
             vcat("OBJ", fill(strcnsmethod(cons_method), length(prob.ucons))),
diff --git a/lib/OptimizationNOMAD/test/qa/qa.jl b/lib/OptimizationNOMAD/test/qa/qa.jl
index f65df9629..5a01837b7 100644
--- a/lib/OptimizationNOMAD/test/qa/qa.jl
+++ b/lib/OptimizationNOMAD/test/qa/qa.jl
@@ -2,7 +2,12 @@ using OptimizationNOMAD, Aqua, JET
 using Test
 
 @testset "Aqua" begin
-    Aqua.test_all(OptimizationNOMAD)
+    # `@reexport using OptimizationBase` exports `solve!`, which clashes with NOMAD's
+    # `solve!` brought in via `using NOMAD`; mark broken until the reexport is restructured.
+    Aqua.test_all(
+        OptimizationNOMAD;
+        undefined_exports = (; broken = true)
+    )
 end
 
 @testset "JET static analysis" begin
diff --git a/lib/OptimizationODE/Project.toml b/lib/OptimizationODE/Project.toml
index 6371eaefb..58cad9179 100644
--- a/lib/OptimizationODE/Project.toml
+++ b/lib/OptimizationODE/Project.toml
@@ -27,23 +27,26 @@ Sundials = "c3572dad-4567-51f8-b174-8c6c989267f4"
 SafeTestsets = "1bc83da4-3b8d-516f-aca4-4fe02f6d838f"
 
 [compat]
-Pkg = "1"
-NonlinearSolve = "4.11"
+ADTypes = "1.18"
+DiffEqBase = "7"
 ForwardDiff = "0.10, 1"
+LinearAlgebra = "1.10"
+NonlinearSolve = "4.11"
 OptimizationBase = "5"
+OrdinaryDiffEqLowOrderRK = "2"
+OrdinaryDiffEqNonlinearSolve = "2"
+OrdinaryDiffEqRosenbrock = "2"
+OrdinaryDiffEqStabilizedRK = "2"
 OrdinaryDiffEqTsit5 = "2"
 OrdinaryDiffEqVerner = "2"
+Pkg = "1"
+Reexport = "1"
+SafeTestsets = "0.1"
 SciMLBase = "3"
 SteadyStateDiffEq = "2.5"
-OrdinaryDiffEqLowOrderRK = "2"
-OrdinaryDiffEqStabilizedRK = "2"
+Sundials = "5"
+Test = "1.10"
 julia = "1.10"
-ADTypes = "1.18"
-OrdinaryDiffEqNonlinearSolve = "2"
-OrdinaryDiffEqRosenbrock = "2"
-DiffEqBase = "7"
-SafeTestsets = "0.1"
-Reexport = "1"
 
 [targets]
 test = ["Sundials", "Test", "SafeTestsets", "Pkg"]
diff --git a/lib/OptimizationOptimJL/Project.toml b/lib/OptimizationOptimJL/Project.toml
index c578c84de..65b9ae195 100644
--- a/lib/OptimizationOptimJL/Project.toml
+++ b/lib/OptimizationOptimJL/Project.toml
@@ -28,6 +28,12 @@ SparseArrays = "1.6"
 Optim = "2"
 Reexport = "1.2"
 SciMLBase = "2.122.1, 3"
+ForwardDiff = "0.10, 1"
+ModelingToolkit = "11"
+Random = "1.10"
+ReverseDiff = "1"
+Test = "1.10"
+Zygote = "0.6, 0.7"
 
 [sources]
 OptimizationBase = {path = "../OptimizationBase"}
diff --git a/lib/OptimizationOptimJL/test/qa/qa.jl b/lib/OptimizationOptimJL/test/qa/qa.jl
index d4b5566f2..294f8b08f 100644
--- a/lib/OptimizationOptimJL/test/qa/qa.jl
+++ b/lib/OptimizationOptimJL/test/qa/qa.jl
@@ -1,8 +1,33 @@
 using OptimizationOptimJL, Aqua, JET
 using Test
+using Optim
 
 @testset "Aqua" begin
-    Aqua.test_all(OptimizationOptimJL)
+    # OptimizationOptimJL implements the SciML optimization interface for Optim,
+    # so the trait/interface methods it adds extend SciML's *own* functions rather
+    # than committing type piracy — mark those functions as own.
+    SB = OptimizationOptimJL.SciMLBase
+    OB = OptimizationOptimJL.OptimizationBase
+    Aqua.test_all(
+        OptimizationOptimJL;
+        piracies = (
+            treat_as_own = [
+                SB.__init,
+                SB.__solve,
+                SB.allowsbounds,
+                SB.allowscallback,
+                SB.allowsconstraints,
+                SB.allowsfg,
+                SB.has_init,
+                SB.requiresbounds,
+                SB.requiresconshess,
+                SB.requiresconsjac,
+                SB.requiresgradient,
+                SB.requireshessian,
+                OB.supports_sense,
+            ],
+        )
+    )
 end
 
 @testset "JET static analysis" begin
diff --git a/lib/OptimizationOptimisers/Project.toml b/lib/OptimizationOptimisers/Project.toml
index fad63136c..b799ed0c7 100644
--- a/lib/OptimizationOptimisers/Project.toml
+++ b/lib/OptimizationOptimisers/Project.toml
@@ -22,6 +22,15 @@ OptimizationBase = "5"
 Reexport = "1.2.2"
 SciMLBase = "2.122.1, 3"
 julia = "1.10"
+ComponentArrays = "0.15"
+ForwardDiff = "0.10, 1"
+Lux = "1"
+MLDataDevices = "1"
+MLUtils = "0.4"
+Printf = "1.10"
+Random = "1.10"
+Test = "1.10"
+Zygote = "0.6, 0.7"
 
 [extras]
 Pkg = "44cfe95a-1eb2-52ea-b672-e2afdf69b78f"
diff --git a/lib/OptimizationOptimisers/src/OptimizationOptimisers.jl b/lib/OptimizationOptimisers/src/OptimizationOptimisers.jl
index 0c945e3d6..dfbbf7f6d 100644
--- a/lib/OptimizationOptimisers/src/OptimizationOptimisers.jl
+++ b/lib/OptimizationOptimisers/src/OptimizationOptimisers.jl
@@ -61,6 +61,7 @@ function SciMLBase.__solve(cache::OptimizationCache{O}) where {O <: AbstractRule
     G = copy(θ)
 
     local x, min_err, min_θ
+    x = cache.f(cache.u0, first(data))
     min_err = typemax(eltype(real(cache.u0))) #dummy variables
     min_opt = 1
     min_θ = cache.u0
diff --git a/lib/OptimizationOptimisers/test/qa/qa.jl b/lib/OptimizationOptimisers/test/qa/qa.jl
index c2f8d62e4..77d29b63c 100644
--- a/lib/OptimizationOptimisers/test/qa/qa.jl
+++ b/lib/OptimizationOptimisers/test/qa/qa.jl
@@ -1,8 +1,25 @@
 using OptimizationOptimisers, Aqua, JET
 using Test
+using Optimisers
 
 @testset "Aqua" begin
-    Aqua.test_all(OptimizationOptimisers)
+    # OptimizationOptimisers implements the SciML optimization interface for
+    # Optimisers, so the trait/interface methods it adds extend SciML's *own*
+    # functions rather than committing type piracy — mark those functions as own.
+    SB = OptimizationOptimisers.SciMLBase
+    Aqua.test_all(
+        OptimizationOptimisers;
+        piracies = (
+            treat_as_own = [
+                SB.__init,
+                SB.__solve,
+                SB.allowscallback,
+                SB.allowsfg,
+                SB.has_init,
+                SB.requiresgradient,
+            ],
+        )
+    )
 end
 
 @testset "JET static analysis" begin
diff --git a/lib/OptimizationPRIMA/Project.toml b/lib/OptimizationPRIMA/Project.toml
index a921fdc7d..9b65875e7 100644
--- a/lib/OptimizationPRIMA/Project.toml
+++ b/lib/OptimizationPRIMA/Project.toml
@@ -17,13 +17,17 @@ Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 ModelingToolkit = "961ee093-0014-501f-94e3-6117800e7a78"
 
 [compat]
-Pkg = "1"
-SafeTestsets = "0.1"
-julia = "1.10"
+ForwardDiff = "0.10, 1"
+ModelingToolkit = "11"
 OptimizationBase = "5.1"
+Pkg = "1"
 PRIMA = "0.2.0"
-SciMLBase = "2.122.1, 3"
 Reexport = "1"
+ReverseDiff = "1"
+SafeTestsets = "0.1"
+SciMLBase = "2.122.1, 3"
+Test = "1.10"
+julia = "1.10"
 
 [sources]
 OptimizationBase = {path = "../OptimizationBase"}
diff --git a/lib/OptimizationPRIMA/src/OptimizationPRIMA.jl b/lib/OptimizationPRIMA/src/OptimizationPRIMA.jl
index e4a47789c..a043a66a5 100644
--- a/lib/OptimizationPRIMA/src/OptimizationPRIMA.jl
+++ b/lib/OptimizationPRIMA/src/OptimizationPRIMA.jl
@@ -80,6 +80,8 @@ function get_solve_func(opt::PRIMASolvers)
         return PRIMA.lincoa
     elseif opt isa COBYLA
         return PRIMA.cobyla
+    else
+        error("Unknown PRIMA solver type: $(typeof(opt))")
     end
 end
 
diff --git a/lib/OptimizationPolyalgorithms/Project.toml b/lib/OptimizationPolyalgorithms/Project.toml
index d0116bffd..4f849ac1b 100644
--- a/lib/OptimizationPolyalgorithms/Project.toml
+++ b/lib/OptimizationPolyalgorithms/Project.toml
@@ -21,14 +21,16 @@ OptimizationOptimisers = {path = "../OptimizationOptimisers"}
 OptimizationOptimJL = {path = "../OptimizationOptimJL"}
 
 [compat]
-Pkg = "1"
-SafeTestsets = "0.1"
-julia = "1.10"
+ForwardDiff = "0.10, 1"
 OptimizationBase = "5"
 OptimizationOptimisers = "0.3.16"
-SciMLBase = "2.58, 3"
-Reexport = "1.2"
 OptimizationOptimJL = "0.4.10"
+Pkg = "1"
+Reexport = "1.2.2"
+SafeTestsets = "0.1"
+SciMLBase = "2.122.1, 3"
+Test = "1.10"
+julia = "1.10"
 
 [targets]
 test = ["ForwardDiff", "Test", "SafeTestsets", "Pkg"]
diff --git a/lib/OptimizationPyCMA/Project.toml b/lib/OptimizationPyCMA/Project.toml
index 611ef9ec1..8c45c28be 100644
--- a/lib/OptimizationPyCMA/Project.toml
+++ b/lib/OptimizationPyCMA/Project.toml
@@ -5,7 +5,6 @@ version = "1.3.2"
 [deps]
 OptimizationBase = "bca83a33-5cc9-4baa-983d-23429ab6bcbb"
 CondaPkg = "992eb4ea-22a4-4c89-a5bb-47a3300528ab"
-Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 SciMLBase = "0bca4576-84f4-4d90-8ffe-ffa030f20462"
 Reexport = "189a3867-3050-52da-a836-e630ba90ab69"
 PythonCall = "6099a3de-0909-46bc-b1f4-468b9a2dfc0d"
@@ -27,6 +26,7 @@ OptimizationBase = {path = "../OptimizationBase"}
 [extras]
 Pkg = "44cfe95a-1eb2-52ea-b672-e2afdf69b78f"
 SafeTestsets = "1bc83da4-3b8d-516f-aca4-4fe02f6d838f"
+Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 
 [targets]
-test = ["SafeTestsets", "Pkg"]
+test = ["SafeTestsets", "Pkg", "Test"]
diff --git a/lib/OptimizationQuadDIRECT/Project.toml b/lib/OptimizationQuadDIRECT/Project.toml
index 73006b5dc..34de5029d 100644
--- a/lib/OptimizationQuadDIRECT/Project.toml
+++ b/lib/OptimizationQuadDIRECT/Project.toml
@@ -19,12 +19,16 @@ Pkg = "1"
 julia = "1.10"
 OptimizationBase = "5"
 QuadDIRECT = "0.1.2"
-SciMLBase = "2.58, 3"
+SciMLBase = "2.122.1, 3"
 SafeTestsets = "0.1"
 Reexport = "1.2"
+Test = "1.10"
 
 [targets]
 test = ["Pkg", "Test", "SafeTestsets"]
 
 [sources.OptimizationBase]
 path = "../OptimizationBase"
+
+[sources.QuadDIRECT]
+url = "https://github.com/timholy/QuadDIRECT.jl.git"
diff --git a/lib/OptimizationSciPy/Project.toml b/lib/OptimizationSciPy/Project.toml
index ec0dc77a0..87566f549 100644
--- a/lib/OptimizationSciPy/Project.toml
+++ b/lib/OptimizationSciPy/Project.toml
@@ -19,13 +19,19 @@ Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 ModelingToolkit = "961ee093-0014-501f-94e3-6117800e7a78"
 
 [compat]
+ForwardDiff = "0.10, 1"
+ModelingToolkit = "11"
+OptimizationBase = "5.1"
 Pkg = "1"
+PythonCall = "0.9.24"
+Random = "1.10"
+Reexport = "1.2"
+ReverseDiff = "1"
 SafeTestsets = "0.1"
-julia = "1.10"
-OptimizationBase = "5.1"
 SciMLBase = "2.122.1, 3"
-Reexport = "1.2"
-PythonCall = "0.9.24"
+Test = "1.10"
+Zygote = "0.6, 0.7"
+julia = "1.10"
 
 [sources]
 OptimizationBase = {path = "../OptimizationBase"}
diff --git a/lib/OptimizationSciPy/src/OptimizationSciPy.jl b/lib/OptimizationSciPy/src/OptimizationSciPy.jl
index 7d83cb742..00e9ca657 100644
--- a/lib/OptimizationSciPy/src/OptimizationSciPy.jl
+++ b/lib/OptimizationSciPy/src/OptimizationSciPy.jl
@@ -29,7 +29,7 @@ function safe_get_message(result)
     if pyisinstance(msg, pybuiltins.list) || pyisinstance(msg, pybuiltins.tuple)
         return join(pyconvert(Vector{String}, msg), ", ")
     end
-    return string(pytypeof(msg))
+    return string(pytype(msg))
 end
 
 # Squash any kind of numeric object down to a Julia Float64.
diff --git a/lib/OptimizationSophia/src/OptimizationSophia.jl b/lib/OptimizationSophia/src/OptimizationSophia.jl
index 69ce63560..4fc7fb1ff 100644
--- a/lib/OptimizationSophia/src/OptimizationSophia.jl
+++ b/lib/OptimizationSophia/src/OptimizationSophia.jl
@@ -114,6 +114,7 @@ function SciMLBase.__solve(cache::OptimizationCache{O}) where {O <: Sophia}
     gₜ = zero(θ)
     mₜ = zero(θ)
     hₜ = zero(θ)
+    x = cache.f(θ, first(data))
     for epoch in 1:maxiters
         for (i, d) in enumerate(data)
             if cache.f.fg !== nothing && dataiterate
diff --git a/lib/OptimizationSpeedMapping/Project.toml b/lib/OptimizationSpeedMapping/Project.toml
index 9111c6519..eb1446b83 100644
--- a/lib/OptimizationSpeedMapping/Project.toml
+++ b/lib/OptimizationSpeedMapping/Project.toml
@@ -15,13 +15,15 @@ ForwardDiff = "f6369f11-7733-5829-9624-2563aa707210"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 
 [compat]
+ForwardDiff = "0.10, 1"
+OptimizationBase = "5"
 Pkg = "1"
+Reexport = "1.2"
 SafeTestsets = "0.1"
-julia = "1.10"
-OptimizationBase = "5"
-SpeedMapping = "0.3"
 SciMLBase = "2.122.1, 3"
-Reexport = "1.2"
+SpeedMapping = "0.3"
+Test = "1.10"
+julia = "1.10"
 
 [sources]
 OptimizationBase = {path = "../OptimizationBase"}
diff --git a/lib/SimpleOptimization/Project.toml b/lib/SimpleOptimization/Project.toml
index f2149e76e..ec6b254f7 100644
--- a/lib/SimpleOptimization/Project.toml
+++ b/lib/SimpleOptimization/Project.toml
@@ -3,21 +3,6 @@ version = "1.2.2"
 uuid = "510db2f7-4254-4e34-bbda-04240c91ca8f"
 authors = ["Utkarsh <rajpututkarsh530@gmail.com> and contributors"]
 
-[targets]
-test = ["Test", "ForwardDiff", "SafeTestsets", "Pkg"]
-
-[compat]
-SimpleNonlinearSolve = "2"
-Pkg = "1"
-julia = "1.10"
-ForwardDiff = "0.10, 1"
-ADTypes = "1"
-OptimizationBase = "5"
-SciMLBase = "2.122.1, 3"
-SafeTestsets = "0.1"
-Enzyme = "0.13"
-Reexport = "1.2"
-
 [weakdeps]
 ForwardDiff = "f6369f11-7733-5829-9624-2563aa707210"
 Enzyme = "7da242da-08ed-463a-9acd-ee780be4f1d9"
@@ -42,3 +27,20 @@ path = "../OptimizationBase"
 [extensions]
 SimpleOptimizationEnzymeExt = "Enzyme"
 SimpleOptimizationForwardDiffExt = "ForwardDiff"
+
+[compat]
+ADTypes = "1.14"
+Enzyme = "0.13"
+ForwardDiff = "0.10, 1"
+LinearAlgebra = "1.10"
+OptimizationBase = "5"
+Pkg = "1"
+Reexport = "1.2"
+SafeTestsets = "0.1"
+SciMLBase = "2.122.1, 3"
+SimpleNonlinearSolve = "1, 2"
+Test = "1.10"
+julia = "1.10"
+
+[targets]
+test = ["Test", "ForwardDiff", "SafeTestsets", "Pkg"]
diff --git a/lib/SimpleOptimization/src/SimpleOptimization.jl b/lib/SimpleOptimization/src/SimpleOptimization.jl
index b1046a2e5..b94230dbf 100644
--- a/lib/SimpleOptimization/src/SimpleOptimization.jl
+++ b/lib/SimpleOptimization/src/SimpleOptimization.jl
@@ -155,9 +155,15 @@ function instantiate_gradient(f, adtype::ADTypes.AbstractADType)
     strtind = isnothing(_strtind) ? 5 : _strtind + 5
     open_nrmlbrkt_ind = findfirst('(', adtypestr)
     open_squigllybrkt_ind = findfirst('{', adtypestr)
-    open_brkt_ind = isnothing(open_squigllybrkt_ind) ? open_nrmlbrkt_ind :
+    open_brkt_ind = if isnothing(open_squigllybrkt_ind)
+        open_nrmlbrkt_ind
+    elseif isnothing(open_nrmlbrkt_ind)
+        open_squigllybrkt_ind
+    else
         min(open_nrmlbrkt_ind, open_squigllybrkt_ind)
-    adpkg = adtypestr[strtind:(open_brkt_ind - 1)]
+    end
+    lastidx = isnothing(open_brkt_ind) ? lastindex(adtypestr) : (open_brkt_ind - 1)
+    adpkg = adtypestr[strtind:lastidx]
     throw(ArgumentError("The passed automatic differentiation backend choice is not available. Please load the corresponding AD package $adpkg."))
 end