Numerical long-term integration of geodesic equations of motion
Numerical long-term integration of geodesic equations of motion
Numerical long-term integration of geodesic equations of motion
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
Structure preserving integrator for Geodesic <strong>equations</strong>Structure preserving step size algorithmsVariable step size and efficient symplectic<strong>integration</strong> are not compatible,cf. [St<strong>of</strong>fer, 1988]◮ Reversibility + symmetry yield conservation properties similar tosymplecticity, cf. [Hairer et al., 2006]⇒ Try to preserve symmetry and reversibility!◮ First attempts for classical celestial mechanics, cf.[Hairer and Söderlind, 2005].◮ These fail for Geodesic <strong>equations</strong> (see experiments below)10/26 | Jonathan Seyrich (<strong>Numerical</strong> <strong>long</strong>-<strong>term</strong> <strong>integration</strong>) Tübingen | January 21, 2013