  Using a Newtonian model of the Solar System with all 8 planets, we perform
extensive tests on various symplectic integrators of high orders, searching for
the best splitting scheme for long term studies in the Solar System. These
comparisons are made in Jacobi and Heliocentric coordinates and the
implementation of the algorithms is fully detailed for practical use. We
conclude that high order integrators should be privileged, with a preference
for the new $(10,6,4)$ method of (Blanes et al., 2012)
