<br><br><div class="gmail_quote">On Wed, Mar 31, 2010 at 4:15 PM, T J <span dir="ltr"><<a href="mailto:tjhnson@gmail.com">tjhnson@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
On Wed, Mar 31, 2010 at 1:21 PM, Charles R Harris<br>
<<a href="mailto:charlesr.harris@gmail.com">charlesr.harris@gmail.com</a>> wrote:<br>
><br>
> Looks like roundoff error.<br>
><br>
<br>
So this is "expected" behavior?<br>
<br>
In [1]: np.logaddexp2(-1.5849625007211563, -53.584962500721154)<br>
Out[1]: -1.5849625007211561<br>
<br>
In [2]: np.logaddexp2(-0.5849625007211563, -53.584962500721154)<br>
Out[2]: nan<br>
<br>
In [3]: np.log2(np.exp2(-0.5849625007211563) + np.exp2(-53.584962500721154))<br>
Out[3]: -0.58496250072115608<br>
<br></blockquote><div><br>I don't see that<br><br>In [1]: np.logaddexp2(-0.5849625007211563, -53.584962500721154)<br>Out[1]: -0.58496250072115619<br><br>In [2]: np.logaddexp2(-1.5849625007211563, -53.584962500721154)<br>
Out[2]: -1.5849625007211561<br><br>What system are you running on.<br><br>Chuck<br></div></div><br>