  We construct the exponential map associated to a nonholonomic system that
allows us to define an exact discrete nonholonomic constraint submanifold. We
reproduce the continuous nonholonomic flow as a discrete flow on this discrete
constraint submanifold deriving an exact discrete version of the nonholonomic
equations. Finally, we derive a general family of nonholonomic integrators.
