Re: spherical bessel functions and utilities implemented

2009-05-17 Thread Vinzent Steinberg
2009/5/17 Alan Bromborsky abro...@verizon.net A standard technique is to search by halves until you get an interval you know you have a root in and switch to the secant, newton, or some higher order method to polish the root. Secant might fail if there are multiple roots in the interval

Re: spherical bessel functions and utilities implemented

2009-05-17 Thread Fredrik Johansson
On Sun, May 17, 2009 at 12:25 PM, Vinzent Steinberg vinzent.steinb...@googlemail.com wrote: 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? There's a tradeoff, and there are examples where 0.01 fails too...

Re: spherical bessel functions and utilities implemented

2009-05-17 Thread Vinzent Steinberg
2009/5/17 Fredrik Johansson fredrik.johans...@gmail.com On Sun, May 17, 2009 at 12:25 PM, Vinzent Steinberg vinzent.steinb...@googlemail.com wrote: 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? There's

Re: spherical bessel functions and utilities implemented

2009-05-17 Thread Ondrej Certik
On Sun, May 17, 2009 at 12:39 PM, Vinzent Steinberg vinzent.steinb...@googlemail.com wrote: 2009/5/17 Fredrik Johansson fredrik.johans...@gmail.com On Sun, May 17, 2009 at 12:25 PM, Vinzent Steinberg vinzent.steinb...@googlemail.com wrote: Yeah, this is what I changed in my patch. 0.001 was