After this procedure is done, one obtains a system which is symmetric hyperbolic for most choices of given lapse-shift functions, once they are suitably re-scaled. Subsequent arguments go very much on adding equations for the lapse-shift vector in order to make the whole system well posed, and presumably useful for some application. It is instructive to think of these modifications of the evolution equations from the point of view of the initial value formulation. There one starts by solving the constraint equations, the time-time and time-space components of the Einstein tensor, at the initial surface. With the initial data thus obtained, one finds the solution to the evolution equations which are taken to be the space-space components of the Einstein tensor. Since that evolution preserves the constraints, (The vector field generating the flow in phase-space is tangent to the constraint sub-manifold.), one can forget about the constraint equations and think of the evolution equations as providing an evolution for the whole phase-space. In this sense, the modification one is making affects the evolution vector outside the constraint sub-manifold, leaving the vector intact at it. Uniqueness of solutions, which follows from the well posedness of the system, then implies that the solutions stay on the sub-manifold. Nevertheless, and we shall return to this point, as shown in , there is no guarantee that the sub-manifold of constraint solutions is stable with respect to the evolution vector field as extended on the whole phase-space. This is an important point for numerical simulations.
© Max Planck Society and the author(s)