Numerical integration in Celestial Mechanics: a case for contact geometry