On the applicability of constrained symplectic integrators in general   relativity
  The purpose of this note is to point out that a naive application of symplectic integration schemes for Hamiltonian systems with constraints such as SHAKE or RATTLE which preserve holonomic constraints encounters difficulties when applied to the numerical treatment of the equations of general relativity. 
J. Frauendiener