On Wed, Jun 1, 2011 at 10:44, Craig Yoshioka <crai...@me.com> wrote: > would anyone object to fixing the numpy mean and stdv functions, so that they > always used a 64-bit value to track sums, or so that they used a running > calculation. That way > > np.mean(np.zeros([4000,4000],'f4')+500) > > would not equal 511.493408?
Yes, I object. You can set the accumulator dtype explicitly if you need it: np.mean(arr, dtype=np.float64) -- Robert Kern "I have come to believe that the whole world is an enigma, a harmless enigma that is made terrible by our own mad attempt to interpret it as though it had an underlying truth." -- Umberto Eco _______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion