Time reversible and symplectic integrators for molecular dynamics simulations of rigid molecules.