This seems incorrect? Shouldn't it be using the equivalent of `crsum` from http://bugs.python.org/file10357/msum4.py ?