Matrix-matrix and matrix-vector product operations are the core of linear algebra.
For dense matrices, the underlying implementation uses OpenBLAS and it hugely benefits from using multithreading.
For sparse matrices we have only sequential Julia functions, implementing matmat or matvec.
Are there any plans to develop parallel (multithreaded/multiworker) matmat and matvec operations for sparse matrices?
I believe it is an important part and should be included into stdlib/sparse.
Is there some kind of recommended replacement to stick to in the mean-time?