Does parallelism in the scipy.sparse solvers only come from multi-threaded BLAS? I.e., the solvers themselves aren't inherently parallel other than through a multi-threaded BLAS? Follow-up question: Numpy/Scipy are using the intel mkl, will scipy.sparse matrices use intel mkl's sparse blas? -Chris