SparseArrays.jl
SparseArrays.jl copied to clipboard
JuliaSparse
SparseArrays.jl is a Julia stdlib
On Julia 1.6 and master, broadcasting the scaled sum of 5-6 complex sparse arrays sporadically produces eltypes of type `Any`. Examples: ```julia julia> using SparseArrays julia> A, B, C, D,...
`SparseVector` is designed assuming that all values that are not stored explicitly are equal to `zero`. I have a use case in which the non-stored value should be `one(::MyType)`, rather...
Concatenation of a sparse and a dense array usually makes a sparse array. This fails completely when the resulting eltype cannot be used for a sparse array: ```julia julia> hcat([1,2,3],...
```julia promote_type(typeof(sprand(10,10,0.1)),typeof(sprand(ComplexF64,10,10,0.1))) SparseMatrixCSC{Tv,Int64} where Tv ``` vs in the dense case: ```julia promote_type(typeof(rand(10,10)),typeof(rand(ComplexF64,10,10))) Array{Complex{Float64},2} ```
Possibly related to JuliaLang/julia#21370: On 0.5.1: ```julia julia> A = sprand(10,10,0.2); julia> @benchmark 3.0*A BenchmarkTools.Trial: memory estimate: 672 bytes allocs estimate: 4 -------------- minimum time: 139.775 ns (0.00% GC) median...
Calling `getindex` on an adjoint of a sparse array falls back (tested on 1.2 and on master built from source today, both on Linux) to a general `AbstractArray` method, making...
https://github.com/JuliaLang/julia/pull/42647 improves sparse matrix slicing when the matrix isn't too sparse and the number of requested indices isn't too large relative to the average number of non-zeros in a column...
If `A` is a sparse matrix, `B` is a sparse or dense vector of non-zeros with appropriate size, then `C = A ./ B` and `C = A ./ B'`...
`nonzeros` is not defined for several structured matrix types (`Symmetric`, `Tridiagonal`, etc. as shown below). It also has behavior that is accurately described by the docs, but counter to what...
In JuliaLang/julia#19518, `broadcast` for sparse matrices underwent a major revision, overall a move in the right direction, IMHO. But it also resulted in a regression (in my eyes) if type...
