  A new procedure to construct symplectic methods for constrained mechanical
systems is developed in this paper. The definition of a map coming from the
notion of retraction maps allows to adapt the continuous problem to the
discretization rule rather than viceversa. As a result, the constraint
submanifold is exactly preserved by the discrete flow and the extension of the
methods to the case of non-linear configuration spaces is straightforward.
