slightly better with test

Author: ChrisRackauckas

src/default.jl

@@ -79,7 +79,9 @@ end
     function defaultalg(A::AbstractSparseMatrixCSC{<:Union{Float64, ComplexF64}, Ti}, b,
         assump::OperatorAssumptions) where {Ti}
         if assump.issq
-            if length(b) <= 10_000 && length(nonzeros(A)) / length(A) < 1e-4
+            @show length(b) <= 10_000
+            @show length(nonzeros(A)) / length(A)
+            if length(b) <= 10_000 && length(nonzeros(A)) / length(A) < 2e-4
                 DefaultLinearSolver(DefaultAlgorithmChoice.KLUFactorization)
             else
                 DefaultLinearSolver(DefaultAlgorithmChoice.UMFPACKFactorization)

test/default_algs.jl

@@ -21,7 +21,7 @@ solve(prob)
     LinearSolve.OperatorAssumptions(false)).alg ===
       LinearSolve.DefaultAlgorithmChoice.QRFactorization
 
-@test LinearSolve.defaultalg(sprand(1000, 1000, 1e-6) + I, zeros(1000)).alg ===
+@test LinearSolve.defaultalg(sprand(10^4, 10^4, 1e-5) + I, zeros(1000)).alg ===
       LinearSolve.DefaultAlgorithmChoice.KLUFactorization
 prob = LinearProblem(sprand(1000, 1000, 0.5), zeros(1000))
 solve(prob)