Related papers: ExaGRyPE: Numerical General Relativity Solvers Bas…
We present the first results on the construction of an exascale hyperbolic PDE engine (ExaHyPE), a code for the next generation of supercomputers with the objective to evolve dynamical spacetimes of black holes, neutron stars and binaries.…
ExaHyPE ("An Exascale Hyperbolic PDE Engine") is a software engine for solving systems of first-order hyperbolic partial differential equations (PDEs). Hyperbolic PDEs are typically derived from the conservation laws of physics and are…
We introduce NRPyElliptic, an elliptic solver for numerical relativity (NR) built within the NRPy+ framework. As its first application, NRPyElliptic sets up conformally flat, binary black hole (BBH) puncture initial data (ID) on a single…
Next-generation gravitational wave detectors such as Cosmic Explorer, the Einstein Telescope, and LISA, demand highly accurate and extensive gravitational wave (GW) catalogs to faithfully extract physical parameters from observed signals.…
The ISO C++17 standard introduces \emph{parallel algorithms}, a parallel programming model promising portability across a wide variety of parallel hardware including multi-core CPUs, GPUs, and FPGAs. Since 2019, the NVIDIA HPC SDK compiler…
We extend a recently developed numerical code to obtain stationary, axisymmetric solutions that describe rotating black hole spacetimes in a wide class of modified theories of gravity. The code utilizes a relaxed Newton-Raphson method to…
Wave equations help us to understand phenomena ranging from earthquakes to tsunamis. These phenomena materialise over very large scales. It would be computationally infeasible to track them over a regular mesh. Yet, since the phenomena are…
The development of a high performance PDE solver requires the combined expertise of interdisciplinary teams with respect to application domain, numerical scheme and low-level optimization. In this paper, we present how the ExaHyPE engine…
Novel applications of Numerical Relativity demand for more flexible algorithms and tools. In this paper, I develop and test a multigrid solver, based on the infrastructure provided by the Einstein Toolkit, for elliptic partial differential…
In the Exa-Dune project we have developed, implemented and optimised numerical algorithms and software for the scalable solution of partial differential equations (PDEs) on future exascale systems exhibiting a heterogeneous massively…
Elliptic partial differential equations must be solved numerically for many problems in numerical relativity, such as initial data for every simulation of merging black holes and neutron stars. Existing elliptic solvers can take multiple…
In an era where we can not afford to checkpoint frequently, replication is a generic way forward to construct numerical simulations that can continue to run even if hardware parts fail. Yet, replication often is not employed on larger…
This work presents a novel methodology for deriving stationary and axially symmetric solutions to Einstein field equations using the 1+3 tetrad formalism. This approach reformulates the Einstein equations into first order scalar equations,…
We study the performance behaviour of a seismic simulation using the ExaHyPE engine with a specific focus on memory characteristics and energy needs. ExaHyPE combines dynamically adaptive mesh refinement (AMR) with ADER-DG. It is…
Memory bound applications such as solvers for large sparse systems of equations remain a challenge for GPUs. Fast solvers should be based on numerically efficient algorithms and implemented such that global memory access is minimised. To…
We present a new numerical code designed to solve the Einstein field equations for axisymmetric spacetimes. The long term goal of this project is to construct a code that will be capable of studying many problems of interest in axisymmetry,…
We describe a high-order ADER-DG solver for the compressible Euler equations within the ExaHyPE framework. The implementation combines a high-order ADER-DG polynomial representation, a local space-time DG predictor, adaptive mesh…
Excision techniques are used in order to deal with black holes in numerical simulations of Einstein equations and consist in removing a topological sphere containing the physical singularity from the numerical domain, applying instead…
We present the numerical relativity module within AthenaK, an open source performance-portable astrophysics code designed for exascale computing applications. This module employs the Z4c formulation to solve the Einstein equations. We…
The Poisson-Boltzmann equation (PBE) is a relevant partial differential equation commonly used in biophysical applications to estimate the electrostatic energy of biomolecular systems immersed in electrolytic solutions. A conventional mean…