SparseArrays.jl
SparseArrays.jl copied to clipboard
JuliaSparse
SparseArrays.jl is a Julia stdlib
Just like https://github.com/JuliaSparse/SparseArrays.jl/pull/183, we should ensure CHOLMOD wrappers also zero out everything. Some examples of the cholmod crashes: https://build.julialang.org/#/builders/63/builds/6826/steps/5/logs/stdio ``` From worker 3: Please submit a bug report with steps...
`rowval` is defined on both, but curiously `getrowval` is not. I'm not sure exactly why we have two identical accessor functions.
Hi, The current implementation seems to produce a dense result vector which is then converted to sparse before returning. I ran into memory scaling issues. Perhaps something along the lines...
Calling three argument `dot` with a Symmetric view of a sparse matrix is very slow. Example: ```julia A = sprandn(1000,1000,0.1) A = A+A' #symmetric B = Symmetric(triu(A)) ``` The quadratic...
As discussed in https://stackoverflow.com/questions/70648351/getting-an-error-when-using-backslash-operator-on-sparse-matrix-and-dense-vector/70662601#70662601, `\` does not work for sparse matrices with Float32 elements. It would be nice to bring back Float32 support in sparse solvers. The issue is that...
As discussed in https://github.com/JuliaData/PooledArrays.jl/issues/76 and https://github.com/JuliaData/DataFrames.jl/issues/2957 here is the current behavior of `map` on `SparseVector`: ``` julia> x = spzeros(5) 5-element SparseVector{Float64, Int64} with 0 stored entries julia> map(_ ->...
```julia julia> a = sprand(Float64, 10, 10, 0.2); julia> a = a'*a + I; julia> Symmetric(a) \ rand(10); julia> a = sprand(Float32, 10, 10, 0.2); julia> a = a'*a +...
Just like LU in JuliaSparse/SparseArrays.jl#131, we also need these: - [ ] cholesky(a') - [ ] ldlt(a') - [ ] det(a') - [ ] inv(a') - [ ] qr(a') works,...
```julia using LinearAlgebra, SparseArrays using BenchmarkTools julia> M = rand(500, 300); julia> S = spzeros(size(M)...); julia> @btime dot($M, $S); 521.241 μs (0 allocations: 0 bytes) ``` In comparison, even materializing...
Using `opnorm` on a sparse matrix (of type `SparseMatrixCSC`) gives this error: ``` ERROR: ArgumentError: 2-norm not yet implemented for sparse matrices. Try opnorm(Array(A)) or opnorm(A, p) where p=1 or...
