diff --git a/src/NonlinearSolve.jl b/src/NonlinearSolve.jl
index 15e1ee3a9..903a1f33d 100644
--- a/src/NonlinearSolve.jl
+++ b/src/NonlinearSolve.jl
@@ -17,11 +17,9 @@ using LazyArrays: LazyArrays, ApplyArray, cache
 using LinearAlgebra: LinearAlgebra, ColumnNorm, Diagonal, I, LowerTriangular, Symmetric,
                      UpperTriangular, axpy!, cond, diag, diagind, dot, issuccess, istril,
                      istriu, lu, mul!, norm, pinv, tril!, triu!
-using LineSearch: LineSearch, AbstractLineSearchAlgorithm, AbstractLineSearchCache,
-                  NoLineSearch, RobustNonMonotoneLineSearch, BackTracking, LineSearchesJL,
-                  LiFukushimaLineSearch
-using LinearSolve: LinearSolve, LUFactorization, QRFactorization,
-                   needs_concrete_A, AbstractFactorization,
+using LineSearch: LineSearch, AbstractLineSearchCache, LineSearchesJL, NoLineSearch,
+                  RobustNonMonotoneLineSearch, BackTracking,  LiFukushimaLineSearch
+using LinearSolve: LinearSolve, QRFactorization, needs_concrete_A, AbstractFactorization,
                    DefaultAlgorithmChoice, DefaultLinearSolver
 using MaybeInplace: @bb
 using Printf: @printf
@@ -110,10 +108,8 @@ include("default.jl")
         NewtonRaphson(),
         TrustRegion(),
         LevenbergMarquardt(),
-        # PseudoTransient(),
         Broyden(),
         Klement(),
-        # DFSane(),
         nothing
     )
 
@@ -139,9 +135,6 @@ include("default.jl")
         LevenbergMarquardt(),
         GaussNewton(),
         TrustRegion(),
-        # LevenbergMarquardt(; linsolve = LUFactorization()),
-        # GaussNewton(; linsolve = LUFactorization()),
-        # TrustRegion(; linsolve = LUFactorization()),
         nothing
     )
 
diff --git a/test/misc/bruss_tests.jl b/test/misc/bruss_tests.jl
index 80b1bd540..32a3fff68 100644
--- a/test/misc/bruss_tests.jl
+++ b/test/misc/bruss_tests.jl
@@ -52,11 +52,6 @@
     sol = solve(prob_brusselator_2d_sparse, NewtonRaphson(); abstol = 1e-8)
     @test norm(sol.resid, Inf) < 1e-8
 
-    # Deprecated
-    sol = solve(prob_brusselator_2d,
-        NewtonRaphson(autodiff = AutoSparse(AutoFiniteDiff())); abstol = 1e-8)
-    @test norm(sol.resid, Inf) < 1e-8
-
     f! = (du, u) -> brusselator_2d_loop(du, u, p)
     du0 = similar(u0)
     jac_prototype = ADTypes.jacobian_sparsity(f!, du0, u0, TracerSparsityDetector())