Abstract:ExaGRyPE describes a suite of solvers for numerical relativity, based upon ExaHyPE 2, the second generation of our Exascale Hyperbolic PDE Engine. The presented generation of ExaGRyPE solves the Einstein equations in the CCZ4 formulation under a 3+1 foliation and focuses on black hole spacetimes. It employs a block-structured Cartesian grid carrying a higher-order order Finite Difference scheme with adaptive mesh refinement, it facilitates massive parallelism combining message passing, domain decomposition and task parallelism, and it supports the injection of particles into the grid as data probes or tracers. We introduce the ExaGRyPE-specific building blocks within ExaHyPE 2, and discuss its software architecture and compute-n-feel: For this, we formalize the creation of any specific simulation with ExaGRyPE as a sequence of lowering operations, where abstract logical tasks are successively broken into finer tasks until we obtain an abstraction level that can be mapped onto a C++ executable. The overall program logic is fully specified via a domain-specific Python interface, we map this logic onto a more detailed set of numerical tasks, subsequently lower this representation onto technical tasks that the underlying ExaHyPE engine uses to parallelize the application, before eventually the technical tasks are mapped onto task graphs including the actual PDE term evaluations, initial conditions, boundary conditions, and so forth. We end up with a rigorous separation of concerns which shields ExaGRyPE users from technical details and hence simplifies the development of novel physical models. We present the simulations and data for the gauge wave, static single black holes and rotating binary black hole systems, demonstrating that the code base is mature and usable. However, we also uncover domain-specific numerical challenges that need further study by the community in future work.

Graphical Processing of Geodesic Propagation in Python

GYOTO 2.0: a polarized relativistic ray-tracing code

Algorithms and Programs for Strong Gravitational Lensing in Kerr Space-time Including Polarization

Interactive visualization of a thin disc around a Schwarzschild black hole

Adaptive Analytical Ray Tracing of Black Hole Photon Rings

KerrGeoPy: A Python Package for Computing Timelike Geodesics in Kerr Spacetime

The Role of Adaptive Ray Tracing in Analyzing Black Hole Structure

Mahakala: a Python-based Modular Ray-tracing and Radiative Transfer Algorithm for Curved Space-times

Gravitational lensing and ghost images in the regular Bardeen no-horizon spacetimes

The Evolution of Distorted Rotating Black Holes I: Methods and Tests

Kerr-fully Diving into the Abyss: Analytic Solutions to Plunging Geodesics in Kerr

Mechanics of geodesics in Information geometry and Black Hole Thermodynamics

GPU-Accelerated Simulations of Isolated Black Holes

Comparison of Polarized Radiative Transfer Codes used by the EHT Collaboration

Propagation of light in the presence of gravity generated by static and spherically symmetric curved space-times using Maxwell equations

A Tetrad-First Approach to Robust Numerical Algorithms in General Relativity

A CODE TO COMPUTE THE EMISSION OF THIN ACCRETION DISKS IN NON-KERR SPACETIMES AND TEST THE NATURE OF BLACK HOLE CANDIDATES

Observing supermassive black holes in virtual reality

ExaGRyPE: Numerical General Relativity Solvers Based upon the Hyperbolic PDEs Solver Engine ExaHyPE

Initial Data and Eccentricity Reduction Toolkit for Binary Black Hole Numerical Relativity Waveforms

From Newton to Einstein –N-Body Dynamics in Galactic Nuclei and SPH Using New Special Hardware and Astrogrid-D