diff --git a/Project.toml b/Project.toml
index 5c51b73d0..ce67ddcac 100644
--- a/Project.toml
+++ b/Project.toml
@@ -25,6 +25,7 @@ PrecompileTools = "aea7be01-6a6a-4083-8856-8a6e6704d82a"
 Preferences = "21216c6a-2e73-6563-6e65-726566657250"
 Reexport = "189a3867-3050-52da-a836-e630ba90ab69"
 SciMLBase = "0bca4576-84f4-4d90-8ffe-ffa030f20462"
+SciMLJacobianOperators = "19f34311-ddf3-4b8b-af20-060888a46c0e"
 SimpleNonlinearSolve = "727e6d20-b764-4bd8-a329-72de5adea6c7"
 SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
 SparseMatrixColorings = "0a514795-09f3-496d-8182-132a7b665d35"
diff --git a/lib/NonlinearSolveFirstOrder/Project.toml b/lib/NonlinearSolveFirstOrder/Project.toml
index 1df55eaa8..942a6928d 100644
--- a/lib/NonlinearSolveFirstOrder/Project.toml
+++ b/lib/NonlinearSolveFirstOrder/Project.toml
@@ -56,6 +56,7 @@ Setfield = "1.1.1"
 SparseArrays = "1.10"
 SparseConnectivityTracer = "0.6.8"
 SparseMatrixColorings = "0.4.5"
+Sparspak = "0.3.9"
 StableRNGs = "1"
 StaticArrays = "1.9.8"
 StaticArraysCore = "1.4.3"
@@ -81,10 +82,11 @@ ReTestItems = "817f1d60-ba6b-4fd5-9520-3cf149f6a823"
 SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
 SparseConnectivityTracer = "9f842d2f-2579-4b1d-911e-f412cf18a3f5"
 SparseMatrixColorings = "0a514795-09f3-496d-8182-132a7b665d35"
+Sparspak = "e56a9233-b9d6-4f03-8d0f-1825330902ac"
 StableRNGs = "860ef19b-820b-49d6-a774-d7a799459cd3"
 StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 Zygote = "e88e6eb3-aa80-5325-afca-941959d7151f"
 
 [targets]
-test = ["Aqua", "BandedMatrices", "BenchmarkTools", "ForwardDiff", "Enzyme", "ExplicitImports", "Hwloc", "InteractiveUtils", "LineSearch", "LineSearches", "NonlinearProblemLibrary", "Pkg", "Random", "ReTestItems", "SparseArrays", "SparseConnectivityTracer", "SparseMatrixColorings", "StableRNGs", "StaticArrays", "Test", "Zygote"]
+test = ["Aqua", "BandedMatrices", "BenchmarkTools", "ForwardDiff", "Enzyme", "ExplicitImports", "Hwloc", "InteractiveUtils", "LineSearch", "LineSearches", "NonlinearProblemLibrary", "Pkg", "Random", "ReTestItems", "SparseArrays", "SparseConnectivityTracer", "SparseMatrixColorings", "Sparspak", "StableRNGs", "StaticArrays", "Test", "Zygote"]
diff --git a/lib/NonlinearSolveFirstOrder/test/rootfind_tests.jl b/lib/NonlinearSolveFirstOrder/test/rootfind_tests.jl
index 3aa61b256..fcbb70a7e 100644
--- a/lib/NonlinearSolveFirstOrder/test/rootfind_tests.jl
+++ b/lib/NonlinearSolveFirstOrder/test/rootfind_tests.jl
@@ -393,7 +393,7 @@ end
 end
 
 @testitem "Simple Sparse AutoDiff" setup=[CoreRootfindTesting] tags=[:core] begin
-    using ADTypes, SparseConnectivityTracer, SparseMatrixColorings
+    using ADTypes, SparseConnectivityTracer, SparseMatrixColorings, Sparspak
 
     @testset for ad in (AutoForwardDiff(), AutoFiniteDiff(), AutoZygote(), AutoEnzyme())
         @testset for u0 in ([1.0, 1.0], 1.0)
diff --git a/lib/NonlinearSolveFirstOrder/test/sparsity_tests.jl b/lib/NonlinearSolveFirstOrder/test/sparsity_tests.jl
index 3e4df284f..22a38b690 100644
--- a/lib/NonlinearSolveFirstOrder/test/sparsity_tests.jl
+++ b/lib/NonlinearSolveFirstOrder/test/sparsity_tests.jl
@@ -1,6 +1,6 @@
 @testitem "Brusselator 2D" tags=[:core] begin
     using LinearAlgebra, SparseArrays, SparseConnectivityTracer, ADTypes,
-          SparseMatrixColorings
+          SparseMatrixColorings, Sparspak
 
     const N = 32
     const xyd_brusselator = range(0, stop = 1, length = N)
@@ -72,7 +72,7 @@ end
 
 @testitem "Structured Jacobians" tags=[:core] begin
     using SparseConnectivityTracer, BandedMatrices, LinearAlgebra, SparseArrays,
-          SparseMatrixColorings
+          SparseMatrixColorings, Sparspak
 
     N = 16
     p = rand(N)