Implement csr_binop_csr for canonical CSR matrices

According to the comment [here](https://github.com/scipy/scipy/blob/master/scipy/sparse/sparsetools/csr.h#L740) we only need to sum duplicates to make sure that resulting matrix is in canonical format.