Function spgemm 

pub fn spgemm(a: &CsrMatrix, b: &CsrMatrix) -> CsrMatrix

Sparse-sparse matrix multiply: C = A * B.

A is m × k, B is k × n, result C is m × n. Uses a dense accumulator row for each output row.