SparseArrays.jl
SparseArrays.jl copied to clipboard
JuliaSparse
SparseArrays.jl is a Julia stdlib
Exactly the same issue as https://github.com/JuliaSparse/SparseArrays.jl/pull/93 since the other one was closed :/ ```julia julia> using SparseArrays julia> direction = spzeros(3); direction[1] = 2 2 julia> argmin(direction) === nothing true...
```julia julia> A = sparse([rand() < 0.01 ? 1 : 0 for _ in 1:50, _ in 1:50]) 50×50 SparseMatrixCSC{Int64,Int64} with 23 stored entries: [...] julia> sum(A, dims=1) 1×50 Array{Int64,2}:...
```Julia julia> using LinearAlgebra julia> using SparseArrays julia> 0 == [0;;] false julia> fill(0,1,1) == fill([0;;],1,1) false julia> sparse(fill(0,1,1)) == sparse(fill([0;;],1,1)) true ``` The values `0` and `[0;;]` (a 1x1...
The expression "sparse \ sparse" returns a dense matrix, while the equivalent "sparse / sparse" returns a sparse matrix: ```Julia julia> using LinearAlgebra julia> using SparseArrays julia> A = sprandn(2,2,...
It was [pointed out on discourse](https://discourse.julialang.org/t/show-and-displaysize/86126/5?u=stevengj) that the Braille display of sparse matrices unconditionally limits the output to `displaysize` here: https://github.com/JuliaSparse/SparseArrays.jl/blob/3b1ce54fe4575913140cc4f3b001d456da6091e9/src/sparsematrix.jl#L304-L305 Arguably, this should be conditional on whether `get(io, :limit,...
This is a fix for https://github.com/JuliaLang/julia/issues/46355, where `SymTridiagonal` matrices backed by sparse arrays couldn't be summed/subtracted/etc. The root cause is that adding `@view`s of SparseArrays produces a dense array. https://github.com/JuliaLang/julia/blob/ec98087cbf19bac26f6be05df9282746b0cffe78/stdlib/LinearAlgebra/src/tridiag.jl#L207...
Consider ``` julia> using SymPy julia> A = spzeros(Sym, 2,2) 2×2 SparseMatrixCSC{Sym, Int64} with 0 stored entries: ⋅ ⋅ ⋅ ⋅ julia> B = sparse(Sym[1 0;0 1]) 2×2 SparseMatrixCSC{Sym, Int64}...
Hi, it would be nice to have [dmperm](https://github.com/DrTimothyAldenDavis/SuiteSparse/blob/master/CXSparse/Source/cs_dmperm.c) implemented as well. I coudn't find a julia implementation and am using a python workaround for now.
We have made a significant number of improvements over the past few months thanks to @SobhanMP, @Wimmerer , and @dkarrasch. We should make sure to submit a NEWS update for...
We should audit the UMFPACK build and configuration to ensure those issues are not causing https://github.com/JuliaSparse/SparseArrays.jl/issues/147 https://github.com/DrTimothyAldenDavis/SuiteSparse/blob/master/SuiteSparse_config/SuiteSparse_config.mk#L272 > -D'LONGBLAS=long' or -DLONGBLAS='long long' defines the integers used by LAPACK and the...
