Implement csr_binop_csr for canonical CSR matrices

Implement csr_binop_csr for canonical CSR matrices
So make a comment in the code, that the result will have indices in sorted order, and so we just need to sum over the duplicates to make it canonical.