diff --git a/Project.toml b/Project.toml
index a0d35dc89..a9869c8b9 100644
--- a/Project.toml
+++ b/Project.toml
@@ -19,7 +19,7 @@ SymbolicIndexingInterface = "2efcf032-c050-4f8e-a9bb-153293bab1f5"
 [compat]
 FastBroadcast = "0.3.5"
 LinearAlgebra = "1"
-LinearSolve = "2.32"
+LinearSolve = "2.32, 3"
 MuladdMacro = "0.2.4"
 OrdinaryDiffEqCore = "1.16"
 Reexport = "1.2.2"
diff --git a/docs/Project.toml b/docs/Project.toml
index 597c40722..1ceb3c4ca 100644
--- a/docs/Project.toml
+++ b/docs/Project.toml
@@ -19,7 +19,7 @@ DiffEqCallbacks = "4"
 Documenter = "1"
 InteractiveUtils = "1"
 LinearAlgebra = "1"
-LinearSolve = "2.32"
+LinearSolve = "2.32, 3"
 OrdinaryDiffEqRosenbrock = "1.4"
 Pkg = "1"
 Plots = "1"
diff --git a/examples/Project.toml b/examples/Project.toml
index 12fd9d33a..bbe7c8183 100644
--- a/examples/Project.toml
+++ b/examples/Project.toml
@@ -14,7 +14,7 @@ SymbolicIndexingInterface = "2efcf032-c050-4f8e-a9bb-153293bab1f5"
 [compat]
 BenchmarkTools = "1"
 DiffEqDevTools = "2"
-LinearSolve = "2"
+LinearSolve = "2, 3"
 OrdinaryDiffEq = "6"
 Plots = "1"
 PositiveIntegrators = "0.1"
diff --git a/test/Project.toml b/test/Project.toml
index a24c9dc5e..2723606f5 100644
--- a/test/Project.toml
+++ b/test/Project.toml
@@ -10,6 +10,7 @@ OrdinaryDiffEqSDIRK = "2d112036-d095-4a1e-ab9a-08536f3ecdbf"
 OrdinaryDiffEqTsit5 = "b1df2697-797e-41e3-8120-5422d3b24e4a"
 OrdinaryDiffEqVerner = "79d7bb75-1356-48c1-b8c0-6832512096c2"
 Plots = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
+RecursiveFactorization = "f2c3362d-daeb-58d1-803e-2bc74f2840b4"
 SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
 StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
@@ -21,13 +22,14 @@ ADTypes = "1.12"
 Aqua = "0.7, 0.8"
 ExplicitImports = "1.0.1"
 LinearAlgebra = "1"
-LinearSolve = "2.32"
+LinearSolve = "2.32, 3"
 OrdinaryDiffEqLowOrderRK = "1.2"
 OrdinaryDiffEqRosenbrock = "1.4"
 OrdinaryDiffEqSDIRK = "1.2"
 OrdinaryDiffEqTsit5 = "1.1"
 OrdinaryDiffEqVerner = "1.1"
 Plots = "1.25.11"
+RecursiveFactorization = "0.2.21"
 SparseArrays = "1"
 StaticArrays = "1.9.7"
 Statistics = "1"
diff --git a/test/runtests.jl b/test/runtests.jl
index b47d833bc..f1e9f0afe 100644
--- a/test/runtests.jl
+++ b/test/runtests.jl
@@ -15,6 +15,8 @@ using OrdinaryDiffEqTsit5: Tsit5
 using OrdinaryDiffEqVerner: Vern9
 using PositiveIntegrators
 
+# load RecursiveFactorization to get RFLUFactorization
+using RecursiveFactorization: RecursiveFactorization
 using LinearSolve: RFLUFactorization, LUFactorization, KrylovJL_GMRES
 
 using Aqua: Aqua