[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
Fri, 3 Sep 2010 13:22:10 +0400
On 9/3/10, Tuomo Keskitalo <address@hidden> wrote:
> > Right. Fixed version (with jacobian, attached below) seems to be working.
> > Actually, it works only for eps_abs ~ round off error. Already for
> > eps_abs = 1e-10 it's easy to see small (linear?) drift in hamiltonian
> > from the constant.
> Sorry, I'm not familiar with your problem. Can you please elaborate?