SparseArrays.jl
SparseArrays.jl copied to clipboard
SparseArrays.jl is a Julia stdlib
I often find myself needing to operate on a subset of a sparse array, typically a few adjacent columns. Creating a new `SparseMatrix` is relatively expensive, so being able to...
``` julia> vcat(1, sparse([1 2])) ERROR: DimensionMismatch: tried to assign 1-element array to 1×2 destination Stacktrace: [1] throw_setindex_mismatch(X::SparseVector{Int64, Int64}, I::Tuple{Int64, Int64}) @ Base ./indices.jl:193 [2] setindex_shape_check @ ./indices.jl:248 [inlined] [3]...
`nonzeros` appears to be giving the wrong answer for this case. See example: > julia> a = sparse([1 0 2; 0 3 0]) 2×3 SparseMatrixCSC{Int64, Int64} with 3 stored entries:...
This PR uses the features of julia https://github.com/JuliaLang/julia/pull/48861 This PR may not be merged before the above PR is enabled.
An example code ```julia A = Symmetric(sprand(20000, 20000, 0.001)) B = rand(100, 20000) C = B * A ``` After executing this code, the type of `C` is a sparse...
I am not exactly sure how to formulate this. I'm not exactly experiencing a bug, but behaviour that reduces the user experience. An exception is thrown When I try to...
It seems like this should be supported: ``` julia> using SparseArrays, LinearAlgebra julia> A = sprand(10,10,0.1); A = A'A + I; julia> F = cholesky(A) SuiteSparse.CHOLMOD.Factor{Float64} type: LLt method: simplicial...
Following a discussion on slack about sparse matrix format interfaces (https://julialang.slack.com/archives/C6A044SQH/p1715862033197389) I would like to bring up the discussion again here on GitHub, also for future reference on the topic....
CSR interface for #538 . Don't mind the typo in the branch name and commit. ## TODOs * [ ] Increase test coverage * [ ] Docs * [ ]...
Creating a brand-new `conda` environment and doing `conda install python`, `conda install R`, and `conda install julia`. This installed Python 3.12.4, R 4.3.3, and Julia 1.10.4: ``` julia> versioninfo() Julia...