diff --git a/src/default.jl b/src/default.jl index 13697bd14..48e8a9b1f 100644 --- a/src/default.jl +++ b/src/default.jl @@ -159,8 +159,7 @@ function defaultalg(A, b, assump::OperatorAssumptions{Bool}) (__conditioning(assump) === OperatorCondition.IllConditioned || __conditioning(assump) === OperatorCondition.WellConditioned) if length(b) <= 10 - # DefaultAlgorithmChoice.GenericLUFactorization - DefaultAlgorithmChoice.RFLUFactorization + DefaultAlgorithmChoice.GenericLUFactorization elseif appleaccelerate_isavailable() DefaultAlgorithmChoice.AppleAccelerateLUFactorization elseif (length(b) <= 100 || (isopenblas() && length(b) <= 500) || diff --git a/src/factorization.jl b/src/factorization.jl index 60fe39c66..51431647b 100644 --- a/src/factorization.jl +++ b/src/factorization.jl @@ -1020,12 +1020,12 @@ function init_cacheval(alg::NormalCholeskyFactorization, A, b, u, Pl, Pr, return ArrayInterface.cholesky_instance(Symmetric(MType{eltype(A), 2}(undef,0,0)), alg.pivot) end -const PREALLOCATED_NORMALCHOLESKY = ArrayInterface.cholesky_instance( +const PREALLOCATED_NORMALCHOLESKY_SYMMETRIC = ArrayInterface.cholesky_instance( Symmetric(rand(1, 1)), NoPivot()) -function init_cacheval(alg::NormalCholeskyFactorization, A::Matrix, b, u, Pl, Pr, +function init_cacheval(alg::NormalCholeskyFactorization, A::Matrix{Float64}, b, u, Pl, Pr, maxiters::Int, abstol, reltol, verbose::Bool, assumptions::OperatorAssumptions) - return PREALLOCATED_NORMALCHOLESKY + return PREALLOCATED_NORMALCHOLESKY_SYMMETRIC end function init_cacheval(alg::NormalCholeskyFactorization,