diff --git a/Project.toml b/Project.toml
index be4775003..315a75108 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,6 +1,6 @@
 name = "NonlinearSolve"
 uuid = "8913a72c-1f9b-4ce2-8d82-65094dcecaec"
-version = "4.20.1"
+version = "4.20.2"
 authors = ["SciML"]
 
 [deps]
@@ -91,7 +91,7 @@ LeastSquaresOptim = "0.8.5"
 LineSearch = "0.1.9"
 LineSearches = "7.3"
 LinearAlgebra = "1.10"
-LinearSolve = "3.48"
+LinearSolve = "3.48, 4"
 MINPACK = "1.2"
 MPI = "0.20.22"
 NLSolvers = "0.5"
diff --git a/docs/Project.toml b/docs/Project.toml
index 21d8fe609..22a2fb313 100644
--- a/docs/Project.toml
+++ b/docs/Project.toml
@@ -77,7 +77,7 @@ DocumenterInterLinks = "1.0.0"
 IncompleteLU = "0.2"
 InteractiveUtils = "<0.0.1, 1"
 LineSearch = "0.1"
-LinearSolve = "2, 3"
+LinearSolve = "2, 3, 4"
 MPI = "0.20"
 NonlinearSolve = "4"
 NonlinearSolveBase = "1.15, 2"
diff --git a/lib/NonlinearSolveBase/Project.toml b/lib/NonlinearSolveBase/Project.toml
index b0a94fa85..1be11c77d 100644
--- a/lib/NonlinearSolveBase/Project.toml
+++ b/lib/NonlinearSolveBase/Project.toml
@@ -82,7 +82,7 @@ FunctionWrappersWrappers = "1"
 InteractiveUtils = "<0.0.1, 1"
 LineSearch = "0.1.4"
 LinearAlgebra = "1.10"
-LinearSolve = "3.48"
+LinearSolve = "3.48, 4"
 LogExpFunctions = "0.3.29, 1"
 Markdown = "1.10"
 MaybeInplace = "0.1.4"
diff --git a/lib/NonlinearSolveFirstOrder/Project.toml b/lib/NonlinearSolveFirstOrder/Project.toml
index 0a217019f..47a213370 100644
--- a/lib/NonlinearSolveFirstOrder/Project.toml
+++ b/lib/NonlinearSolveFirstOrder/Project.toml
@@ -1,7 +1,7 @@
 name = "NonlinearSolveFirstOrder"
 uuid = "5959db7a-ea39-4486-b5fe-2dd0bf03d60d"
 authors = ["Avik Pal <avikpal@mit.edu> and contributors"]
-version = "2.1.2"
+version = "2.1.3"
 
 [deps]
 ADTypes = "47edcb42-4c32-4615-8424-f2b9edc5f35b"
@@ -42,7 +42,7 @@ InteractiveUtils = "<0.0.1, 1"
 LineSearch = "0.1.4"
 LineSearches = "7.3.0"
 LinearAlgebra = "1.10"
-LinearSolve = "2.36.1, 3"
+LinearSolve = "2.36.1, 3, 4"
 MaybeInplace = "0.1.4"
 NonlinearProblemLibrary = "0.1.2"
 NonlinearSolveBase = "2.20"
diff --git a/lib/NonlinearSolveQuasiNewton/Project.toml b/lib/NonlinearSolveQuasiNewton/Project.toml
index c129a6abc..09286ca63 100644
--- a/lib/NonlinearSolveQuasiNewton/Project.toml
+++ b/lib/NonlinearSolveQuasiNewton/Project.toml
@@ -1,7 +1,7 @@
 name = "NonlinearSolveQuasiNewton"
 uuid = "9a2c21bd-3a47-402d-9113-8faf9a0ee114"
 authors = ["Avik Pal <avikpal@mit.edu> and contributors"]
-version = "1.13.2"
+version = "1.13.3"
 
 [deps]
 ArrayInterface = "4fba245c-0d91-5ea0-9b3e-6abc04ee57a9"
@@ -40,7 +40,7 @@ InteractiveUtils = "<0.0.1, 1"
 LineSearch = "0.1.4"
 LineSearches = "7.3.0"
 LinearAlgebra = "1.10"
-LinearSolve = "2.36.1, 3"
+LinearSolve = "2.36.1, 3, 4"
 MaybeInplace = "0.1.4"
 NonlinearProblemLibrary = "0.1.2"
 NonlinearSolveBase = "2.20"
diff --git a/test/gpu/Project.toml b/test/gpu/Project.toml
index 109ec0363..b5c608fe6 100644
--- a/test/gpu/Project.toml
+++ b/test/gpu/Project.toml
@@ -29,7 +29,7 @@ ADTypes = "1.9"
 BracketingNonlinearSolve = "1.6"
 CUDA = "6"
 LinearAlgebra = "1.10"
-LinearSolve = "3.48"
+LinearSolve = "3.48, 4"
 NonlinearSolve = "4.19"
 NonlinearSolveBase = "2.20"
 NonlinearSolveFirstOrder = "2"