On Wed, Mar 31, 2010 at 11:25 PM, <josef.pktd@gmail.com> wrote:
On Thu, Apr 1, 2010 at 1:22 AM,  <josef.pktd@gmail.com> wrote:
> On Thu, Apr 1, 2010 at 1:17 AM, Charles R Harris
> <charlesr.harris@gmail.com> wrote:
>>
>>
>> On Wed, Mar 31, 2010 at 6:08 PM, <josef.pktd@gmail.com> wrote:
>>>
>>> On Wed, Mar 31, 2010 at 7:37 PM, Warren Weckesser
>>> <warren.weckesser@enthought.com> wrote:
>>> > T J wrote:
>>> >> On Wed, Mar 31, 2010 at 1:21 PM, Charles R Harris
>>> >> <charlesr.harris@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