diff --git a/benchmark/Project.toml b/benchmark/Project.toml
index e7d2caf0..f8803a23 100644
--- a/benchmark/Project.toml
+++ b/benchmark/Project.toml
@@ -1,4 +1,5 @@
 [deps]
 BenchmarkTools = "6e4b80f9-dd63-53aa-95a3-0cdb28fa8baf"
+NonlinearSolveFirstOrder = "5959db7a-ea39-4486-b5fe-2dd0bf03d60d"
 ODEInterface = "54ca160b-1b9f-5127-a996-1867f4bc2a2c"
 OrdinaryDiffEq = "1dea7af3-3e70-54e6-95c3-0bf5283fa5ed"
diff --git a/benchmark/benchmarks.jl b/benchmark/benchmarks.jl
index 91f07cb9..2fdec47e 100644
--- a/benchmark/benchmarks.jl
+++ b/benchmark/benchmarks.jl
@@ -1,5 +1,5 @@
 using BenchmarkTools
-using BoundaryValueDiffEq, OrdinaryDiffEq
+using BoundaryValueDiffEq, OrdinaryDiffEq, NonlinearSolveFirstOrder
 
 include("simple_pendulum.jl")
 
diff --git a/lib/BoundaryValueDiffEqFIRK/test/runtests.jl b/lib/BoundaryValueDiffEqFIRK/test/runtests.jl
index 1a91fd7b..ad766138 100644
--- a/lib/BoundaryValueDiffEqFIRK/test/runtests.jl
+++ b/lib/BoundaryValueDiffEqFIRK/test/runtests.jl
@@ -15,4 +15,5 @@ const RETESTITEMS_NWORKER_THREADS = parse(Int,
 
 ReTestItems.runtests(
     BoundaryValueDiffEqFIRK; tags = (GROUP == "all" ? nothing : [Symbol(GROUP)]),
-    nworkers = RETESTITEMS_NWORKERS, nworker_threads = RETESTITEMS_NWORKER_THREADS)
+    nworkers = RETESTITEMS_NWORKERS,
+    nworker_threads = RETESTITEMS_NWORKER_THREADS, testitem_timeout = 120 * 60)
diff --git a/test/misc/non_vector_input_tests.jl b/test/misc/non_vector_input_tests.jl
index b4a967f3..5c6cbf1f 100644
--- 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
+    using LinearAlgebra, NonlinearSolveFirstOrder
 
     for order in (2, 3, 4, 5, 6)
         s = Symbol("MIRK$(order)")
diff --git a/test/wrappers/odeinterface_tests.jl b/test/wrappers/odeinterface_tests.jl
index 1f4153dd..e3af08c1 100644
--- a/test/wrappers/odeinterface_tests.jl
+++ b/test/wrappers/odeinterface_tests.jl
@@ -45,7 +45,7 @@ end
 
 # Just test that it runs. BVPSOL only works with linearly separable BCs.
 @testitem "BVPSOL" setup=[ODEInterfaceWrapperTestSetup] begin
-    using ODEInterface, RecursiveArrayTools
+    using ODEInterface, RecursiveArrayTools, NonlinearSolveFirstOrder
 
     tpprob = TwoPointBVProblem(ex7_f!, (ex7_2pbc1!, ex7_2pbc2!), u0, tspan,
         p; bcresid_prototype = (zeros(1), zeros(1)))