[SciPy-User] Broyden with complex numbers [Was Re: ANN: SciPy 0.11.0 beta 1 released]

[Skipper Seabold]

On Wed, Jun 27, 2012 at 1:25 PM, Denis Laxalde <denis at laxalde.org> wrote:

> Skipper Seabold a écrit :
> >>>>>>    TypeError: Cannot cast ufunc add output from dtype('complex128')
> to
> >>>>>>    dtype('float64') with casting rule 'same_kind'
> >>>> I get "nonlin.py:314: ComplexWarning: Casting complex values to real
> >>>> discards the imaginary part" with Numpy 1.6.1.
> >>>>
> >>>> There's indeed a probable bug on line 314 --- should maybe do
> >>>> "x = x + d" instead of "x += d".
> >>>
> >>> I agree. That would be consistent with what's done with line_search.
> >>
> >> Fixed in master and maintenance/0.11.x.
> >>
> >
> > Great, thanks for looking into this.
>
> That being said, your example still fails for me because of NaNs...
>

Hmm, I'm not sure if it's a 'typlcal' root finding problem, but I was
definitely having problems with most of the solvers trying to get a
solution when I don't know good starting values. It's why I ended up
finding the closed form solution.

Skipper
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20120627/b376afec/attachment.html>