  We present our numerical comparisons between the BSSN formulation widely usedin numerical relativity today and its adjusted versions using constraints. Weperformed three testbeds: gauge-wave, linear wave, and Gowdy-wave tests,proposed by the Mexico workshop on the formulation problem of the Einsteinequations. We tried three kinds of adjustments, which were previously proposedfrom the analysis of the constraint propagation equations, and investigated howthey improve the accuracy and stability of evolutions. We observed that thesignature of the proposed Lagrange multipliers are always right and theadjustments improve the convergence and stability of the simulations. When theoriginal BSSN system already shows satisfactory good evolutions (e.g., linearwave test), the adjusted versions also coincide with those evolutions; while insome cases (e.g., gauge-wave or Gowdy-wave tests) the simulations using theadjusted systems last 10 times as long as those using the original BSSNequations. Our demonstrations imply a potential to construct a robust evolutionsystem against constraint violations even in highly dynamical situations.
We performed three testbeds: gauge-wave, linear wave, and Gowdy-wave tests, proposed by the Mexico workshop on the formulation problem of the Einstein equations. We tried three kinds of adjustments, which were previously proposed from the analysis of the constraint propagation equations, and investigated how they improve the accuracy and stability of evolutions.