Same question here: is it really worth having cython code for this thing? Doesn't this work very fast in  a numpy-based python implementation?