  The dual-frame formalism leads to an approach to extend numerical relativity
simulations in generalized harmonic gauge (GHG) all the way to null infinity. A
major setback is that without care, even simple choices of initial data give
rise to logarithmically divergent terms that would result in irregular
variables and equations on the compactified domain, which would in turn prevent
accurate numerical approximation. It has been shown, however, that a suitable
choice of gauge and constraint addition can be used to prevent their
appearance. Presently we give a first order symmetric hyperbolic reduction of
general relativity in GHG on compactified hyperboloidal slices that exploits
this knowledge and eradicates these log-terms at leading orders. Because of
their effect on the asymptotic solution space, specific formally singular terms
are systematically chosen to remain. Such formally singular terms have been
successfully treated numerically in toy models and result in a formulation with
the desirable property that unphysical radiation content near infinity is
suppressed.
