Toward making the constraint hypersurface an attractor in free evolution

When constructing numerical solutions to systems of evolution equations subject to a constraint, one must decide what role the constraint equations will play in the evolution system. In one popular choice, known as free evolution, a simulation is treated as a Cauchy problem, with the initial data constructed to satisfy the constraint equations. This initial data are then evolved via the evolution equations with no further enforcement of the constraint equations. The evolution, however, via the discretized evolution equations introduce constraint violating modes at the level of truncation error, and these constraint violating modes will behave in a formalism dependent way. This paper presents a generic method for incorporating the constraint equations into the evolution equations so that the off-constraint dynamics are biased toward the constraint satisfying solutions.