don't fill! on a SparseMatrix

JuliaDiff / SparseDiffTools.jl

Fast jacobian computation through sparsity exploitation and matrix coloring

MIT License

240 stars 42 forks source link

don't fill! on a SparseMatrix #179

Closed ChrisRackauckas closed 2 years ago

ChrisRackauckas commented 2 years ago

julia> x = sprand(2,2,0.5)
2×2 SparseMatrixCSC{Float64, Int64} with 1 stored entry:
 0.507603   ⋅ 
  ⋅         ⋅ 

julia> fill!(x, true)
2×2 SparseMatrixCSC{Float64, Int64} with 4 stored entries:
 1.0  1.0
 1.0  1.0

Seems like an issue for Base.

ChrisRackauckas commented 2 years ago

https://github.com/JuliaLang/julia/issues/17670

ChrisRackauckas commented 2 years ago

Most now supports CUDA https://github.com/JuliaGPU/CUDA.jl/pull/1367