13.07.2015 Views

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

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!