<br><br><div class="gmail_quote">On Tue, Mar 24, 2009 at 1:12 PM, Pauli Virtanen <span dir="ltr">&lt;<a href="mailto:pav@iki.fi">pav@iki.fi</a>&gt;</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;">
<div class="im">Tue, 24 Mar 2009 13:15:21 -0400, Darren Dale wrote:<br>
&gt; I just performed an svn update, deleted my old build/ and<br>
&gt; site-packages/numpy*, reinstalled, and I see a new test failure on a 64<br>
&gt; bit linux machine:<br>
&gt;<br>
&gt; ======================================================================<br>
&gt; FAIL: test_umath.TestComplexFunctions.test_loss_of_precision_longcomplex<br>
&gt; ----------------------------------------------------------------------<br>
</div>[clip]<br>
<div class="im">&gt;   &quot;/usr/lib64/python2.6/site-packages/numpy/core/tests/test_umath.py&quot;,<br>
&gt; line 497, in check<br>
&gt;     &#39;arctanh&#39;)<br>
&gt; AssertionError: (135, 3.4039637354191726288e-09,<br>
&gt; 3.9031278209478159624e-18, &#39;arctanh&#39;)<br>
<br>
</div>I can reproduce this (on another 64-bit machine). This time around, it&#39;s<br>
the real function that is faulty:<br>
<br>
&gt;&gt;&gt; x = np.longdouble(3e-9)<br>
&gt;&gt;&gt; np.arctanh(x+0j).real - x<br>
9.0876776281460559983e-27<br>
&gt;&gt;&gt; np.arctanh(x).real - x<br>
0.0<br>
&gt;&gt;&gt; np.finfo(np.longdouble).eps * x<br>
3.2526065174565132804e-28<br>
<br>
So, the system atanhl is ~ 30 relative eps away from the correct answer:<br>
</blockquote><div><br>I see this also. The compiler is gcc version 4.3.0 20080428 (Red Hat 4.3.0-8) (GCC). Maybe we should ping the compiler folks? I could also open a Fedora bug for this.<br><br>Chuck  <br><br></div></div>
<br>