CSR matrix-matrix products algorithm

Question

I found a lot of examples and algorithms on CSR matrix creation, visualize, product with vector. But I didn't found any algorithm on CSR-matrix * CSR-matrix product.

I found a nice link in which is explained the matrix dot vector product.

Could somebody help me?

Answer 1

You can look into the references and source code of some libraries implemented this operation

• SuiteSparse has SSMULT
• PETSc is a library collection and has MatMatMult
• SPARSEKIT has AMUB

Or even easier you could just use one of these libraries.