[Help-gsl] convergence check in odeiv_rk4imp ?

[Forest Yang]

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
property,
it has better conservation of the first integral.


Best regards,
Forest.