diff --git a/Project.toml b/Project.toml
index c39dee1c..d4be4163 100644
--- a/Project.toml
+++ b/Project.toml
@@ -59,11 +59,11 @@ MuladdMacro = "0.2"
 OpenCL = "0.10"
 Parameters = "0.12, 0.13"
 RecursiveArrayTools = "4"
-SciMLBase = "3"
+SciMLBase = "3.1"
 Setfield = "1"
 SimpleDiffEq = "1.11"
 SimpleNonlinearSolve = "2"
-StaticArrays = "1.9"
+StaticArrays = "1.9.14"
 TOML = "1"
 ZygoteRules = "0.2.7"
 julia = "1.10"
diff --git a/test/Project.toml b/test/Project.toml
index 558c0313..630ebcaa 100644
--- a/test/Project.toml
+++ b/test/Project.toml
@@ -41,7 +41,7 @@ ForwardDiff = "1"
 GPUArraysCore = "0.2"
 JLArrays = "0.3"
 KernelAbstractions = "0.9"
-ModelingToolkit = "11.17.0"
+ModelingToolkit = "11.23.0"
 OpenCL = "0.10"
 Optimization = "5"
 OptimizationOptimisers = "0.3"
@@ -50,8 +50,8 @@ OrdinaryDiffEqRosenbrock = "2"
 OrdinaryDiffEqSDIRK = "2"
 OrdinaryDiffEqStabilizedRK = "2"
 SafeTestsets = "0.1"
-SciMLBase = "3"
-StaticArrays = "1.9"
+SciMLBase = "3.1"
+StaticArrays = "1.9.14"
 Statistics = "1"
 StochasticDiffEq = "7"
 TOML = "1"
diff --git a/test/gpu_kernel_de/stiff_ode/gpu_ode_mass_matrix.jl b/test/gpu_kernel_de/stiff_ode/gpu_ode_mass_matrix.jl
index e35fd652..9db8fffc 100644
--- a/test/gpu_kernel_de/stiff_ode/gpu_ode_mass_matrix.jl
+++ b/test/gpu_kernel_de/stiff_ode/gpu_ode_mass_matrix.jl
@@ -35,8 +35,11 @@ monteprob = EnsembleProblem(prob, safetycopy = false)
 
 alg = GPURosenbrock23()
 
+# dt must match the Float32 eltype of u0/p. A Float64 dt makes calc_W build a
+# Float64 StaticWOperator that OrdinaryDiffEqRosenbrock's JacReuseState (typed
+# from the Float32 state) cannot store, erroring in the cached_W assignment.
 bench_sol = solve(
-    prob, Rosenbrock23(), dt = 0.1, abstol = 1.0f-5, reltol = 1.0f-5,
+    prob, Rosenbrock23(), dt = 0.1f0, abstol = 1.0f-5, reltol = 1.0f-5,
     initializealg = BrownFullBasicInit()
 )