Regge calculus as a fourth-order method in numerical relativity

The convergence properties of numerical Regge calculus as an approximation to continuum vacuum General Relativity is studied, both analytically and numerically. The Regge equations are evaluated on continuum spacetimes by assigning squared geodesic distances in the continuum manifold to the squared edge lengths in the simplicial manifold. It is found analytically that, individually, the Regge equations converge to zero as the second power of the lattice spacing, but that an average over local Regge equations converges to zero as (at the very least) the third power of the lattice spacing. Numerical studies using analytic solutions to the Einstein equations show that these averages actually converge to zero as the fourth power of the lattice spacing.