pub fn sparse_mul(a: &SparseCsr, b: &SparseCsr) -> Result<SparseCsr, String>
Element-wise multiplication (Hadamard product) of two CSR matrices. Only positions where BOTH matrices have non-zeros produce non-zeros.