Is your feature request related to a problem? Please describe.

I am not sure if I am doing something wrong but I am using scipy.sparse.csr_matrix object and contract it with a np.ndarray object using the method dot: scipy.sparse.csr_matrix.dot(np.ndarray).
I noticed that the operation is not parallelized. Is there any possibility to parallelize the sparse matrix operations (like dot/kron) in scipy?

PS: I test the execution with numpy.dot at it runs in multicore. Both my scipy and numpy are installed with conda and linked with blas/lapack

Describe the solution you'd like.

The solution I would like is the parallelization of the dot/kron products with sparse matrices.

Describe alternatives you've considered.

numpy.dot without relying in sparse matrix

Additional context (e.g. screenshots, GIFs)

No response