pub fn sparse_matmul(a: &SparseCsr, b: &SparseCsr) -> Result<SparseCsr, String>Expand description
Sparse matrix-matrix multiplication (SpGEMM): C = A * B. Uses row-wise accumulation with BTreeMap for deterministic column ordering. All floating-point reductions use binned summation.