SparseArrays.jl
SparseArrays.jl copied to clipboard
JuliaSparse
SparseArrays.jl is a Julia stdlib
The reason seems to be that in the 3-arg `ldiv!` its 2-arg-version is called, which is missing for sparse Cholesky factorizations. MWE: ``` using LinearAlgebra, SparseArrays A = sprand(100, 100,...
https://github.com/JuliaLang/julia/blob/b56a9f07948255dfbe804eef25bdbada06ec2a57/stdlib/SuiteSparse/src/spqr.jl#L207 The code should be ```julia size(A, 1) !=size(Q.factors,2) ``` as in https://github.com/JuliaLang/julia/blob/80516ca20297a67b996caa08c38786332379b6a5/stdlib/LinearAlgebra/src/qr.jl#L533
The following errors for `cholesky` factorizations: ```julia julia> using LinearAlgebra, SparseArrays julia> A = cholesky(sparse(Matrix(I,10,10))); # sparse cholesky factorization with real elements julia> z = rand(ComplexF64, 10); # complex RHS...
Moving from https://github.com/JuliaLinearAlgebra/SuiteSparse.jl/issues/4. CHOLMOD can provide an estimate of the reciprocal condition number (via cholmod_rcond) once a factorisation is available, but there is no way to access this functionality from...
A\B works if A is a sparse matrix and B is a sparse vector. A\B does not work if B is a sparse matrix. Example: A = sparse(I,100,100) B =...
If `A` is a sparse matrix, then we get ```jl julia> qr(A)' \ c ERROR: MethodError: no method matching adjoint(::SuiteSparse.SPQR.QRSparse{Float64,Int64}) ``` It would be nice to have this implemented. For...
`lq(A)` is not working with a sparse matrix `A`. It should be possible to use the sparse QR factorization available in SuiteSparse on `A'`.
``` julia> using LinearAlgebra julia> using SparseArrays julia> M = Symmetric(sparse(rand(4,4))); julia> F = LinearAlgebra.ldlt(M) SuiteSparse.CHOLMOD.Factor{Float64} type: LDLt method: simplicial maxnnz: 16 nnz: 10 success: true julia> sparse(F.D) ERROR: getindex...
I would like to know if this would be a good PR to the SparseArrays stdlib I can't workout what this operation is called. It does `A.*b` inplace for `A::SparseMatrixCSC`...
Currently, `mul!(C, A, B, α, β)` for the case `A` is a sparse matrix is defined as https://github.com/JuliaLang/julia/blob/2d4f4d26a0c5a74717d826dc7e1e62f7650a2e9d/stdlib/SparseArrays/src/linalg.jl#L34-L52 This method is not used in case `B` is `Symmetric`, sparse, etc....
