I would like to bring the following improvements to
julia> n = 10000 julia> Random.seed!(0); A = sprandn(n, n, 0.01); b = randn(n); nnz(A) 1000598 julia> SA = Symmetric(A); julia> @btime B = SparseMatrixCSC(SA); 3.049 s (35 allocations: 40.03 MiB) julia> using SparseWrappers julia> @btime B = SparseMatrixCSC(SA); 10.060 ms (7 allocations: 15.32 MiB)
I would probably feel better if somebody welcomed this work. The implementation exists in my package SparseWrappers but needs to be converted into a PR for julia.