pub fn matmul<F>( a: &CsrMatrix<F>, b: &CsrMatrix<F>, ) -> SparseResult<CsrMatrix<F>>where F: Float + NumAssign + Sum + SparseElement + 'static + Debug,
Matrix multiplication