Gravitational waveforms with controlled accuracy

A partially first-order form of the characteristic formulation is introduced to control the accuracy in the computation of gravitational waveforms produced by highly distorted single black hole spacetimes. Our approach is to reduce the system of equations to first-order differential form on the angular derivatives, while retaining the proven radial and time integration schemes of the standard characteristic formulation. This results in significantly improved accuracy over the standard mixed-order approach in the extremely nonlinear post-merger regime of binary black hole collisions.