How to regularize a symplectic-energy-momentum integrator