[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Help-gsl] convergence check in odeiv_rk4imp ?

From: Forest Yang
Subject: [Help-gsl] convergence check in odeiv_rk4imp ?
Date: Tue, 29 Dec 2009 13:35:15 -0500

Hi All,

    Is there any plan for implementing convergence check in implicit
RK4 method ?
It seems the 2 stages Gauss points are only iterate 3 steps, without
checking whether it is converged or not.
Another possible way could be using the Jacobian to solve for "Y_i",
the fixed point.

   The reason I care about this is, Implicit Gauss RK4 is symplectic,
which for Hamiltonian Dynamics simulation is a very important
it has better conservation of the first integral.

Best regards,

reply via email to

[Prev in Thread] Current Thread [Next in Thread]