Skip to content

Commit

Permalink
more factorization init good stuff
Browse files Browse the repository at this point in the history
  • Loading branch information
ChrisRackauckas committed Dec 21, 2021
1 parent 7b2d57f commit b58b8c0
Showing 1 changed file with 1 addition and 4 deletions.
5 changes: 1 addition & 4 deletions src/factorization.jl
Original file line number Diff line number Diff line change
Expand Up @@ -225,10 +225,6 @@ function init_cacheval(alg::Union{GenericFactorization,GenericFactorization{type
end
end

function init_cacheval(alg, A::Union{Hermitian,Symmetric}, b, u, Pl, Pr, maxiters, abstol, reltol, verbose)
BunchKaufman(A.data, Array(1:size(A,1)), A.uplo, true, false, 0)
end

# Fallback, tries to make nonsingular and just factorizes
# Try to never use it.
function init_cacheval(alg::Union{QRFactorization,SVDFactorization,GenericFactorization}, A, b, u, Pl, Pr, maxiters, abstol, reltol, verbose)
Expand All @@ -241,3 +237,4 @@ end

RFLUFactorization() = GenericFactorization(;fact_alg=RecursiveFactorization.lu!)
init_cacheval(alg::GenericFactorization{typeof(RecursiveFactorization.lu!)}, A, b, u, Pl, Pr, maxiters, abstol, reltol, verbose) = ArrayInterface.lu_instance(A)
init_cacheval(alg::GenericFactorization{typeof(RecursiveFactorization.lu!)}, A::StridedMatrix{<:LinearAlgebra.BlasFloat}, b, u, Pl, Pr, maxiters, abstol, reltol, verbose) = ArrayInterface.lu_instance(A)

0 comments on commit b58b8c0

Please sign in to comment.