I'm fine for the A*B part. I'm hung up on the B*C part - as I understand it, the mkl_?csrmm functions take the left matrix as the sparse one, while the right is dense. I'm OK enough with linear algebra to remember that matrix multiplication is not commutative. How can I then get the product of B*C using the mkl_?csrmm function?
C := alpha*A'*B + beta*C