[Numpy-discussion] math.fsum() like ufunc for numpy

2013-12-02 Thread Daniele Nicolodi
Hello, there would be interest in adding a floating point accurate summation function like Python's math.fsum() in the form of an ufunc to NumPy? I had a look at the algorithm (http://code.activestate.com/recipes/393090-binary-floating-point-summation-accurate-to-full-p/) and it looks quite

Re: [Numpy-discussion] math.fsum() like ufunc for numpy

2013-12-02 Thread Nathaniel Smith
I think that would be great. Technically what you'd want is a gufunc. -n On Mon, Dec 2, 2013 at 9:44 AM, Daniele Nicolodi dani...@grinta.net wrote: Hello, there would be interest in adding a floating point accurate summation function like Python's math.fsum() in the form of an ufunc to

Re: [Numpy-discussion] math.fsum() like ufunc for numpy

2013-12-02 Thread Julian Taylor
related this PR attempts to improve the accuracy of summation: https://github.com/numpy/numpy/pull/3685 but math.fsum gives the exact result so it would a valuable ufunc even when that PR is merged. python3.4 will have yet another accurate summation in the statistics module: