<br><br><div class="gmail_quote">On Thu, Mar 19, 2009 at 1:01 PM, Charles R Harris <span dir="ltr">&lt;<a href="mailto:charlesr.harris@gmail.com">charlesr.harris@gmail.com</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;">
<br><br><div class="gmail_quote">2009/3/19 Robert Pyle <span dir="ltr">&lt;<a href="mailto:rpyle@post.harvard.edu" target="_blank">rpyle@post.harvard.edu</a>&gt;</span><div class="im"><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">

<div><br>
On Mar 19, 2009, at 1:38 PM, Pauli Virtanen wrote:<br>
<br>
Thanks for tracking this! Can you check what your platform gives for:<br>
<br>
&gt; import numpy as np<br>
&gt; info = np.finfo(np.longcomplex)<br>
&gt; print &quot;eps:&quot;, info.eps, info.eps.dtype<br>
&gt; print &quot;tiny:&quot;, info.tiny, info.tiny.dtype<br>
&gt; print &quot;log10:&quot;, np.log10(info.tiny), np.log10(info.tiny/info.eps)<br>
<br>
</div>eps: 1.3817869701e-76 float128<br>
tiny: -1.08420217274e-19 float128<br>
log10: nan nan<br>
</blockquote></div><div><br>The log of a negative number is nan, so part of the problem is the value of tiny. The size of the values also look suspect to me. On my machine<br><br>In [8]: finfo(longcomplex).eps<br>Out[8]: 1.084202172485504434e-19<br>

<br>In [9]: finfo(float128).tiny<br>Out[9]: array(3.3621031431120935063e-4932, dtype=float128)<br><br>So at a minimum eps and tiny are reversed.<br></div></div></blockquote><div><br>I started to look at the code for this but my eyes rolled up in my head and I passed out. It could use some improvements...<br>
<br>Chuck  <br></div></div><br>