Re: [Numpy-discussion] large float32 array issue

2010-11-04 Thread Vincent Schut
On 11/03/2010 03:04 PM, Bruce Southey wrote: > On 11/03/2010 06:52 AM, Pauli Virtanen wrote: >> Wed, 03 Nov 2010 12:39:08 +0100, Vincent Schut wrote: >> [clip] >>> Btw, should I file a bug on this? >> One can argue that mean() and sum() should use a numerically stabler >> algorithm, so yes, a bug c

Re: [Numpy-discussion] large float32 array issue

2010-11-03 Thread braingateway
Vincent Schut : > Hi, I'm running in this strange issue when using some pretty large > float32 arrays. In the following code I create a large array filled with > ones, and calculate mean and sum, first with a float64 version, then > with a float32 version. Note the difference between the two. NB

Re: [Numpy-discussion] large float32 array issue

2010-11-03 Thread Bruce Southey
On 11/03/2010 06:52 AM, Pauli Virtanen wrote: > Wed, 03 Nov 2010 12:39:08 +0100, Vincent Schut wrote: > [clip] >> Btw, should I file a bug on this? > One can argue that mean() and sum() should use a numerically stabler > algorithm, so yes, a bug can be filed if there is not yet one already. > This

Re: [Numpy-discussion] large float32 array issue

2010-11-03 Thread Charles R Harris
On Wed, Nov 3, 2010 at 5:52 AM, Pauli Virtanen wrote: > Wed, 03 Nov 2010 12:39:08 +0100, Vincent Schut wrote: > [clip] > > Btw, should I file a bug on this? > > One can argue that mean() and sum() should use a numerically stabler > algorithm, so yes, a bug can be filed if there is not yet one alr

Re: [Numpy-discussion] large float32 array issue

2010-11-03 Thread Warren Weckesser
On Wed, Nov 3, 2010 at 6:39 AM, Vincent Schut wrote: > > > On 11/03/2010 12:31 PM, Warren Weckesser wrote: > > > > > > On Wed, Nov 3, 2010 at 5:59 AM, Warren Weckesser > > mailto:warren.weckes...@enthought.com>> > > wrote: > > > > > > > > On Wed, Nov 3, 2010 at 3:54 AM, Vincent Schut > <

Re: [Numpy-discussion] large float32 array issue

2010-11-03 Thread Pauli Virtanen
Wed, 03 Nov 2010 12:39:08 +0100, Vincent Schut wrote: [clip] > Btw, should I file a bug on this? One can argue that mean() and sum() should use a numerically stabler algorithm, so yes, a bug can be filed if there is not yet one already. -- Pauli Virtanen ___

Re: [Numpy-discussion] large float32 array issue

2010-11-03 Thread Vincent Schut
On 11/03/2010 12:31 PM, Warren Weckesser wrote: > > > On Wed, Nov 3, 2010 at 5:59 AM, Warren Weckesser > mailto:warren.weckes...@enthought.com>> > wrote: > > > > On Wed, Nov 3, 2010 at 3:54 AM, Vincent Schut > wrote: > > Hi, I'm running in this strange

Re: [Numpy-discussion] large float32 array issue

2010-11-03 Thread Warren Weckesser
On Wed, Nov 3, 2010 at 5:59 AM, Warren Weckesser < warren.weckes...@enthought.com> wrote: > > > On Wed, Nov 3, 2010 at 3:54 AM, Vincent Schut wrote: > >> Hi, I'm running in this strange issue when using some pretty large >> float32 arrays. In the following code I create a large array filled with

Re: [Numpy-discussion] large float32 array issue

2010-11-03 Thread Warren Weckesser
On Wed, Nov 3, 2010 at 3:54 AM, Vincent Schut wrote: > Hi, I'm running in this strange issue when using some pretty large > float32 arrays. In the following code I create a large array filled with > ones, and calculate mean and sum, first with a float64 version, then > with a float32 version. Not

[Numpy-discussion] large float32 array issue

2010-11-03 Thread Vincent Schut
Hi, I'm running in this strange issue when using some pretty large float32 arrays. In the following code I create a large array filled with ones, and calculate mean and sum, first with a float64 version, then with a float32 version. Note the difference between the two. NB the float64 version is