A GPU-accelerated mixed-precision WENO method for extremal black hole and gravitational wave physics computations

Abstract

We develop and use a novel mixed-precision weighted essentially non-oscillatory (WENO) method for solving the Teukolsky equation, which arises when modeling perturbations of Kerr black holes. We show that WENO methods outperform higher-order finite-difference methods, standard in the discretization of the Teukolsky equation, due to the need to add dissipation for stability purposes in the latter. In particular, as the WENO scheme uses no additional dissipation it is well-suited for scenarios requiring long-time evolution such as the study of Price tails and gravitational wave emission from extreme mass ratio binaries. In the mixed-precision approach, the expensive computation of the WENO weights is performed in reduced floating-point precision that results in a significant speedup factor of 3.3. In addition, we use state-of-the-art Nvidia general-purpose graphics processing units and cluster parallelism to further accelerate the WENO computations. Our optimized WENO solver can be used to quickly generate accurate results of significance in the field of black hole and gravitational wave physics. We apply our solver to study the behavior of the Aretakis charge -- a conserved quantity, that if detected by a gravitational wave observatory like LIGO/Virgo would prove the existence of extremal black holes.

19 pages, 6 plots

References in corpus (12)

• First M87 Event Horizon Telescope Results. IV. Imaging the Central Supermassive Black Hole
• Hyperboloidal foliations and scri-fixing
• Towards adiabatic waveforms for inspiral into Kerr black holes: I. A new model of the source for the time domain perturbation equation
• A geometric framework for black hole perturbations
• Towards adiabatic waveforms for inspiral into Kerr black holes: II. Dynamical sources and generic orbits
• A new gravitational wave generation algorithm for particle perturbations of the Kerr spacetime
• Binary black hole merger gravitational waves and recoil in the large mass ratio limit
• Hyperboloidal evolution with the Einstein equations
• Spacelike matching to null infinity
• Gravitational perturbations of Schwarzschild spacetime at null infinity and the hyperboloidal initial value problem
• A hyperboloidal study of tail decay rates for scalar and Yang-Mills fields
• Scalar and gravitational hair for extreme Kerr black holes

Cited by in corpus (4)

• Surrogate model for gravitational wave signals from non-spinning, comparable- to large-mass-ratio black hole binaries built on black hole perturbation theory waveforms calibrated to numerical relativity
• Measuring quasi-normal mode amplitudes with misaligned binary black hole ringdowns
• Scalar and Gravitational Transient "Hair" for Near-Extremal Black Holes
• Discontinuous Galerkin method for linear wave equations involving derivatives of the Dirac delta distribution