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>The new integratorNew step size control algorithm◮ Use constant underlying step size ɛ < 1.◮ Variable step size in step y n → y n+1 :◮ Resulting scheme:h n (ɛ, Y i ) =12[∥∂f∂yɛ](Y 1 ) +y n+1 = y n + h n (ɛ, Y i )Y i = y n + h n (ɛ, Y i )[ ]∂f∂y(Y s ) ∥s∑b i f (Y i ),i=1s∑a ij f (Y j ),j=111/26 | Jonathan Seyrich (<strong>Numerical</strong> <strong>long</strong>-<strong>term</strong> <strong>integration</strong>) Tübingen | January 21, 2013