Alfredo Braunstein
Alfredo Braunstein
`sprand` accepts a function argument `rfn` that is used to generate the non zero values of a sparse random matrix or vector. I find the calling convention with rfn to...
See https://github.com/JuliaGraphs/Graphs.jl/pull/150#issuecomment-1374486801: To compare I've used ```julia for N in (500, 5000), is_directed in (true, false), p in (10/N, 0.5, 1-10/N) @btime erdos_renyi($N, $p; is_directed=$is_directed) end ``` That gives ```...
I was thinking, in the spirit of this package, maybe `rand(Normal(),SparseMatrixCSC,p,m,n) ` could be better expressed as `rand(Bernoulli(p,Normal()),SparseMatrixCSC,m,n)` where `Bernoulli(p, Normal())` would be the "Gauss-Bernoulli" or "Spike-and-Slab" mixture distribution P(x)...