Symplectic methods for the approximation of the exponential map and the Newton iteration on Riemannian submanifolds

We use a Hamiltonian approach and symplectic methods to compute the geodesics on a Riemannian submanifold.