The dynamical motion of mechanical systems possesses underlying geometric structures, and preserving these structures in numerical integration improves the qualitative accuracy and reduces the long-time error of the simulation. For a single mechanical system, structure preservation can be achieved by adopting the variational integrator construction. This construction has been generalized to mor...