<div class="gmail_quote">On Wed, Jun 27, 2012 at 1:25 PM, Denis Laxalde <span dir="ltr">&lt;<a href="mailto:denis@laxalde.org" target="_blank">denis@laxalde.org</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

Skipper Seabold a écrit :<br>
<div class="im">&gt;&gt;&gt;&gt;&gt;&gt;    TypeError: Cannot cast ufunc add output from dtype(&#39;complex128&#39;) to<br>
&gt;&gt;&gt;&gt;&gt;&gt;    dtype(&#39;float64&#39;) with casting rule &#39;same_kind&#39;<br>
&gt;&gt;&gt;&gt; I get &quot;nonlin.py:314: ComplexWarning: Casting complex values to real<br>
&gt;&gt;&gt;&gt; discards the imaginary part&quot; with Numpy 1.6.1.<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; There&#39;s indeed a probable bug on line 314 --- should maybe do<br>
&gt;&gt;&gt;&gt; &quot;x = x + d&quot; instead of &quot;x += d&quot;.<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; I agree. That would be consistent with what&#39;s done with line_search.<br>
&gt;&gt;<br>
&gt;&gt; Fixed in master and maintenance/0.11.x.<br>
&gt;&gt;<br>
&gt;<br>
&gt; Great, thanks for looking into this.<br>
<br>
</div>That being said, your example still fails for me because of NaNs...<br></blockquote><div><br>Hmm, I&#39;m not sure if it&#39;s a &#39;typlcal&#39; root finding problem, but I was definitely having problems with most of the solvers trying to get a solution when I don&#39;t know good starting values. It&#39;s why I ended up finding the closed form solution.<br>

<br>Skipper<br></div></div>