  Combining incoming and outgoing characteristic formulations can provide
numerical relativists with a natural implementation of Einstein's equations
that better exploits the causal properties of the spacetime and gives access to
both null infinity and the interior region simultaneously (assuming the
foliation is free of caustics and crossovers). We discuss how this combination
can be performed and illustrate its behavior in the Einstein-Klein-Gordon field
in 1D.
