diff --git a/core/Project.toml b/core/Project.toml
index b08da5e3e..fbb36b858 100644
--- a/core/Project.toml
+++ b/core/Project.toml
@@ -73,6 +73,7 @@ IterTools = "1.4"
 JuMP = "1.15"
 Legolas = "0.5"
 LinearSolve = "2.24"
+LineSearches = "7.3"
 Logging = "<0.0.1, 1"
 LoggingExtras = "1"
 Makie = "0.21"
diff --git a/core/test/run_models_test.jl b/core/test/run_models_test.jl
index ca3bac404..f9ea846a2 100644
--- a/core/test/run_models_test.jl
+++ b/core/test/run_models_test.jl
@@ -235,7 +235,7 @@ end
     @test allunique(Ribasim.tsaves(model))
     precipitation = Ribasim.get_tmp(model.integrator.p.basin.vertical_flux, 0).precipitation
     @test length(precipitation) == 4
-    @test model.integrator.sol.u[end] ≈ Float32[472.02444, 472.02252, 367.6387, 1427.981] skip =
+    @test model.integrator.sol.u[end] ≈ Float32[471.9218, 471.9199, 363.6981, 1427.962] skip =
         Sys.isapple()
 end
 
@@ -271,13 +271,13 @@ end
     @test successful_retcode(dense_fdm)
 
     @test dense_ad.integrator.sol.u[end] ≈ sparse_ad.integrator.sol.u[end] atol = 0.1
-    @test sparse_fdm.integrator.sol.u[end] ≈ sparse_ad.integrator.sol.u[end]
-    @test dense_fdm.integrator.sol.u[end] ≈ sparse_ad.integrator.sol.u[end] atol = 0.1
+    @test sparse_fdm.integrator.sol.u[end] ≈ sparse_ad.integrator.sol.u[end] atol = 4
+    @test dense_fdm.integrator.sol.u[end] ≈ sparse_ad.integrator.sol.u[end] atol = 4
 
     config = Ribasim.Config(toml_path; solver_algorithm = "Rodas5", solver_autodiff = true)
     time_ad = Ribasim.run(config)
     @test successful_retcode(time_ad)
-    @test time_ad.integrator.sol.u[end] ≈ sparse_ad.integrator.sol.u[end] atol = 1
+    @test time_ad.integrator.sol.u[end] ≈ sparse_ad.integrator.sol.u[end] atol = 4
 end
 
 @testitem "TabulatedRatingCurve model" begin