Numerical Approach for Corvino-Type Gluing of Brill-Lindquist Initial Data.

Building on the work of Giulini and Holzegel (2005), a new numerical approach is developed for computing Cauchy data for Einstein's equations by gluing a Schwarzschild end to a Brill-Lindquist metric via a Corvino-type construction. In contrast to, and in extension of, the numerical strategy of Doulis and Rinne (2016), the overdetermined Poisson problem resulting from the Brill wave ansatz is decomposed to obtain two uniquely solvable problems. A pseudospectral method and Newton-Krylov root finder are utilized to perform the gluing. The convergence analysis strongly indicates that the numerical strategy developed here is able to produce highly accurate results. It is observed that Schwarzschild ends of various ADM masses can be glued to the same interior configuration using the same gluing radius.