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 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 straightforward to implement. I can try to submit a
 patch for it.

 Cheers,
 Daniele
 ___
 NumPy-Discussion mailing list
 NumPy-Discussion@scipy.org
 http://mail.scipy.org/mailman/listinfo/numpy-discussion



-- 
Nathaniel J. Smith
Postdoctoral researcher - Informatics - University of Edinburgh
http://vorpus.org
___
NumPy-Discussion mailing list
NumPy-Discussion@scipy.org
http://mail.scipy.org/mailman/listinfo/numpy-discussion

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:
http://www.python.org/dev/peps/pep-0450/

On 02.12.2013 19:07, Nathaniel Smith wrote:
 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 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 straightforward to implement. I can try to submit a
 patch for it.

 Cheers,
 Daniele
 ___
 NumPy-Discussion mailing list
 NumPy-Discussion@scipy.org
 http://mail.scipy.org/mailman/listinfo/numpy-discussion
 
 
 

___
NumPy-Discussion mailing list
NumPy-Discussion@scipy.org
http://mail.scipy.org/mailman/listinfo/numpy-discussion