Czy scipy obsługuje wielowątkowość dla rzadkiego mnożenia macierzy podczas używania MKL BLAS?
Zgodnie z dokumentacją MKL BLAS „Wszystkie operacje macierzowo-macierzowe (poziom 3) są gwintowane zarówno dla gęstych, jak i rzadkich BLAS”.http://software.intel.com/en-us/articles/parallelism-in-the-intel-math-kernel-library
Zbudowałem Scipy z MKL BLAS. Używając poniższego kodu testowego, widzę oczekiwane przyspieszenie wielowątkowe dla gęstego, ale nie rzadkiego mnożenia macierzy. Czy są jakieś zmiany w Scipy, aby włączyć wielowątkowe rzadkie operacje?
# test dense matrix multiplication
from numpy import *
import time
x = random.random((10000,10000))
t1 = time.time()
foo = dot(x.T, x)
print time.time() - t1
# test sparse matrix multiplication
from scipy import sparse
x = sparse.rand(10000,10000)
t1 = time.time()
foo = dot(x.T, x)
print time.time() - t1