Long-term Integrations of the Solar System
University of Toronto
I'll discuss the differences between symplectic and non-symplectic integrators in the context of long term orbit simulations of the Solar System. In particular, I'll present our recent implementation of the non-symplectic Gauss-Radau scheme, IAS15, as well as the symplectic integrator WHFast. These algorithms are the most accurate and fastest algorithms running on standard computer architectures. Both of them achieve sub-linear energy error growth (Brouwer's law) over a billion dynamical times (more than 10^10 timesteps). I'll talk about the tricks we used to write code that is both fast but also ensures that all important floating-point operations are unbiased. Our integrators are freely available and come with an easy to use python wrapper. I will give an live demo on how to use them at the end of my talk.
Date: Tuesday, 13 October 2015 Time: 15:30 Where: McGill University Ernest Rutherford Physics Building, R.E. Bell Conference Room (room 103) Contact: Robert Rutledge