On 24/12/14 13:23, Julian Taylor wrote: > hm this is a brute force search, probably fast enough but slower than > scipy's code (if it also were cython)
That was what I tought as well when I wrote it. But it turned out that regular numbers are so close and abundant that was damn fast, even in Python :) > I also ported it to C a while back, so that could be used for numpy if > speed is an issue. Code attached. Very nice :) Sturla _______________________________________________ NumPy-Discussion mailing list [email protected] http://mail.scipy.org/mailman/listinfo/numpy-discussion
