SuiteSparseGraphBLAS.jl

SuiteSparseGraphBLAS.jl copied to clipboard

I expect the following script should use 16 threads:

using SparseArrays
using SuiteSparseGraphBLAS
using Arpack

gbset(:nthreads, 16)
A = GBMatrix(sprand(100000, 100000, 0.001))
H = A' + A

@info "matrix constructed"

eigs(H, nev=6, which=:SR, maxiter=1000)

However, the CPU usage keeps < 100% during eigs, for which I believe most of the time should be doing matrix vector multiplication:

This is falling back on a generic matmul right now. I will try to figure out how to deal with that.