Variable time step integration with symplectic methods

From MaRDI portal





It is a fundamental problem in the symplectic integration of Hamiltonian systems that the favorable properties of such methods get lost, if one uses variable stepsizes. The author proposes to perform a time reparametrization in the Hamiltonian system instead of changing the stepsize during the numerical integration. A naive transformation of the time variable will not preserve the Hamiltonian character of the system. But the author shows how one can define for given initial data a new Hamiltonian with a transformed time variable such that the solutions of the two Hamiltonian systems for these special initial data coincide. A backward error analysis of the approach is given. The author discusses several concrete choices of time reparametrization and demonstrates their performances in numerical experiments with the Kepler problem.




Cited in
(38)








This page was built for publication: Variable time step integration with symplectic methods

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1372703)