sparse-sparse matrix product BLAS

[usermg]

Hello everybody,

I have written a code in fortran capable of making the product between a large sparse matrix and another large sparse matrix (quite fast). In order to verify if my results are correct and if the code is enough fast I was wondering if some external subroutine could do the same. In this sense, BLAS offers some functions and subroutines in the context of sparse vectors and matrices. However, as far as I have been able to find, with BLAS one can make the product sparse matrix-normal matrix, sparse matrix-vector. So, do you know if with BLAS one can make the product sparse matrix-sparse-matrix?

Thank you very much in advance,

Miquel

 

[mikrom]

Hi usermg,
I found this article

http://www.dddas.org/~douglas/Preprints/pub0034.pdf

with this source

http://www.netlib.org/aicm/smmp

which should contain subroutine NUMBMM for sparse-matrix * sparse-matrix multiplication

Other package is SPARSEKIT, but like a sparse BLAS it doesn't seem to have a sparse-matrix * sparse-matrix multiplication:

http://stackoverflow.com/questions/12328482/fortran-90-95-library-for-sparse-matrices

http://people.sc.fsu.edu/~jburkardt/f_src/sparsekit/sparsekit.html