In #127 I realized that `cs{rc}_matrix.sum` doesn’t work as documented and filed scipy/scipy#23768. Should we do that on our side or wait for scipy?