  In Brodbeck et al 1999 it has been shown that the linearised time evolution
equations of general relativity can be extended to a system whose solutions
asymptotically approach solutions of the constraints. In this paper we extend
the non-linear equations in similar ways and investigate the effect of various
possibilities by numerical means. Although we were not able to make the
constraint submanifold an attractor for all solutions of the extended system,
we were able to significantly reduce the growth of the numerical violation of
the constraints. Contrary to our expectations this improvement did not imply a
numerical solution closer to the exact solution, and therefore did not improve
the quality of the numerical solution.
