SparseArrays.jl
SparseArrays.jl copied to clipboard
SparseArrays.jl is a Julia stdlib
```julia julia> spzeros(2) + [1.4, 3] 2-element SparseVector{Float64, Int64} with 2 stored entries: ⋮ [2] = 3.0 julia> spzeros(2,2) + [1.4 3 0 2] 2×2 Matrix{Float64}: 1.4 3.0 0.0 2.0...
Minimal working example: ```julia N = 3000 A = sprand(ComplexF64, N, N, 1/N) # A = rand(ComplexF64, N, N) x = randn(N) + 1im*randn(N) y = randn(N) + 1im*randn(N) α...
When a sparse matrix is multiplied elementwise by a dense one, the result remains sparse. I expect the same to happen when the sparse matrix is multiplied by an NxN...
Please refer to details on the discourse thread https://discourse.julialang.org/t/arnoldimethod-strange-error/108876
Reverts much of https://github.com/JuliaSparse/SparseArrays.jl/pull/471 This PR introduces a framework for installing test deps during stdlib CI in a way that respects `[compat]` entries. We can likely use this same framework...
The type-pirated `vcat` methods introduced here are leading to method ambiguities on Julia v1.10, of the form: ```julia special vcat: Error During Test at /home/jishnu/Dropbox/JuliaPackages/InfiniteArrays.jl/test/runtests.jl:721 Test threw exception Expression: [randn(2,...
(I was [told](https://github.com/JuliaLang/julia/issues/51966) to open this issue here.) This seems weird: ``` sebastienloisel@macbook-pro julia % julia _ _ _ _(_)_ | Documentation: https://docs.julialang.org (_) | (_) (_) | _ _...
It's a minor issue, but I accidentally evaluated a large sparse matrix (~500M rows) in REPL and it caused a ~5-second freeze. Not sure what's happening under the hood but...
Update: another surprising consequence of https://github.com/JuliaSparse/SparseArrays.jl/issues/190 that led to hours of debugging... SparseMatrixCSC Float64.(A) drops structural nonzeros, but float(A) doesn't: ``` julia> A = SparseArrays.sparse([1 2 0; 0 0 3;...
Very closely related to [this issue](https://github.com/JuliaGPU/CUDA.jl/issues/1855) in CUDA.jl. It was suggested that the GPU component of SuiteSparse/CHOLMOD could be implemented as a package extension to SparseArrays.jl when CUDA.jl is also...