Oscar Reula (F.A.M.A.F-U.N.C. y Mirta Iriondo (F.A.M.A.F-U.N.C.)
The evolution equations can be writen
in a very simple form and are symmetric hyperbolic in Eddington-Finkelstein
coordinates.
The numerical evolution is performed
with a standard method of lines
fourth order in space and time.
The time algorithm is Runge-Kuta while the
space discrete derivative is symmetric
(non-dissipative). The constraints are
preserved under evolution (within
numerical errors) and we are able to
reproduce several known results.