2009/5/17 Fredrik Johansson <fredrik.johans...@gmail.com> > > On Sun, May 17, 2009 at 11:00 AM, Ondrej Certik <ond...@certik.cz> wrote: > > Thanks for all the suggestions and thanks Vinzent for looking into > > that. Let's fix the findroot() in mpmath, e..g create a newton method, > > that only needs one point, close to the root and it converges by > > itself. > > For a well-behaving function (spherical Bessel functions qualify), > the secant method with a small initial step size like 0.001 > should converge and is faster than Newton's method, > provided that the starting point is good. >
Yeah, this is what I changed in my patch. 0.001 was too small, but 0.01 did the job. Shouldn't this be the default instead of 0.25? Vinzent --~--~---------~--~----~------------~-------~--~----~ You received this message because you are subscribed to the Google Groups "sympy-patches" group. To post to this group, send email to sympy-patches@googlegroups.com To unsubscribe from this group, send email to sympy-patches+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sympy-patches?hl=en -~----------~----~----~----~------~----~------~--~---