  Full relativistic simulations in three dimensions are known to develop
runaway modes that grow exponentially and are accompanied by violations of the
Hamiltonian and momentum constraints. We present here a method that controls
the violation of these constraints and is tested with simulations of binary
neutron stars in circular orbits. We show that this technique improves the
overall quality of the simulations.
