Closed stefano-meschiari closed 10 years ago
Ended up keeping leapfrog for simplicity's sake — using Kepler's equation does not save a lot of computational time (if at all), and now that collisions are detected, there is no need for variable-step integration.
For systems without planet-planet interactions, use Keplerian orbital elements to compute cartesian coordinates of the planets — cuts down on calculations, enables fixed time steps, avoids precession for highly eccentric objects and makes it easier to detect collisions with the central star.