On Wed, Mar 31, 2010 at 11:25 PM, <josef.p...@gmail.com> wrote: > On Thu, Apr 1, 2010 at 1:22 AM, <josef.p...@gmail.com> wrote: > > On Thu, Apr 1, 2010 at 1:17 AM, Charles R Harris > > <charlesr.har...@gmail.com> wrote: > >> > >> > >> On Wed, Mar 31, 2010 at 6:08 PM, <josef.p...@gmail.com> wrote: > >>> > >>> On Wed, Mar 31, 2010 at 7:37 PM, Warren Weckesser > >>> <warren.weckes...@enthought.com> wrote: > >>> > T J wrote: > >>> >> On Wed, Mar 31, 2010 at 1:21 PM, Charles R Harris > >>> >> <charlesr.har...@gmail.com> wrote: > >>> >> > >>> >>> Looks like roundoff error. > >>> >>> > >>> >>> > >>> >> > >>> >> So this is "expected" behavior? > >>> >> > >>> >> In [1]: np.logaddexp2(-1.5849625007211563, -53.584962500721154) > >>> >> Out[1]: -1.5849625007211561 > >>> >> > >>> >> In [2]: np.logaddexp2(-0.5849625007211563, -53.584962500721154) > >>> >> Out[2]: nan > >>> >> > >>> > > >>> > Is any able to reproduce this? I don't get 'nan' in either 1.4.0 or > >>> > 2.0.0.dev8313 (32 bit Mac OSX). In an earlier email T J reported > using > >>> > 1.5.0.dev8106. > >>> > >>> > >>> > >>> >>> np.logaddexp2(-0.5849625007211563, -53.584962500721154) > >>> nan > >>> >>> np.logaddexp2(-1.5849625007211563, -53.584962500721154) > >>> -1.5849625007211561 > >>> > >>> >>> np.version.version > >>> '1.4.0' > >>> > >>> WindowsXP 32 > >>> > >> > >> What compiler? Mingw? > > > > yes, mingw 3.4.5. , official binaries release 1.4.0 by David > > sse2 Pentium M > > Can you try the exp2/log2 functions with the problem data and see if something goes wrong?
Chuck
_______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion
