Hi all
the current version of mkl_dcsrmm implements the operations
C=AB+C and C=A'B+C, where in both cases A is the csr sparse matrix.
Is there any way/alternative such that B is the csr sparse matrix??
In my specific application C and A are large dense matrices of dimension eg 3,000,000:70. B is of dimension 70:70 with up to 80% zero coefficients. Ignoring the sparse structure of B would lead to dgemm, but this would cause a large overhead for multiplications with zeros. Also mkl_dcsrmm would be feasible with A and C being transposed before the operation, however, due to their size this is not possible with regard to speed (mkl_dcsrmm is called up to 10,000 times).
Any suggestions appreciated.
Cheers