[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] [bug #30540] please, add convergence checks in rk4imp/rk2i
Sergey B Kirpichev
Re: [Bug-gsl] [bug #30540] please, add convergence checks in rk4imp/rk2imp
Sat, 4 Sep 2010 16:28:32 +0400
On 9/4/10, Tuomo Keskitalo <address@hidden> wrote:
> Thanks for the link. That is a special code meant for second order ODEs.
> Apparently your case requires precise solution, and therefore you should
> modify the Newton iteration code if you want to use GSL for your problem.
> However, I think your case is special. At least I do not know of a generally
> accepted all-purpose tolerance for the Newton iteration.
I just want to use, e.g. the implicit RK method of order 2 (midpoint
rule) and rely on the proved properties of the stepper function
(symmetric, symplectic). Either the method shouldn't be named
"rk2imp". Hope, this opinion is clear enough.