Related papers: Towards an exascale code for GRMHD on dynamical sp…

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 SphericalNR, a new framework for the publicly available Einstein Toolkit that numerically solves the Einstein field equations coupled to the equations of general relativistic magnetohydrodynamics (GRMHD) in a 3+1 split of…

Simulations of binary black hole systems using the Spectral Einstein Code (SpEC) are done on a computational domain that excises the regions inside the black holes. It is imperative that the excision boundaries are outflow boundaries with…

We present Horizon, a new graphics processing unit (GPU)-accelerated code to solve the equations of general relativistic magnetohydrodynamics in a given spacetime. We evaluate the code in several test cases, including magnetized Riemann…

(abridged) In the first two papers of this publication series, we present a comprehensive library of synthetic EHT observations and used this library to train and validate Bayesian neural networks for the parameter inference of accreting…

We present a sequence of optimizations to the performance-critical compute kernels of the high-order discontinuous Galerkin solver of the hyperbolic PDE engine ExaHyPE -- successively tackling bottlenecks due to SIMD operations, cache…

We have extended Cosmos++, a multi-dimensional unstructured adaptive mesh code for solving the covariant Newtonian and general relativistic radiation magnetohydrodynamic (MHD) equations, to accommodate both discrete finite volume and…

We extend a new finite element code, Einstein PHG (iPHG), to solve the evolution part of Einstein equations in first-order GH formalism. This paper is the third one of a systematic investigation of applying adaptive finite element method to…

The nonlinear behavior of higher dimensional black hole spacetimes is of interest in several contexts, ranging from an understanding of cosmic censorship to black hole production in high-energy collisions. However, nonlinear numerical…

This paper describes the development and testing of a general relativistic magnetohydrodynamic (GRMHD) code to study ideal MHD in the fixed background of a Kerr black hole. The code is a direct extension of the hydrodynamic code of Hawley,…

This paper presents a new computer code to solve the general relativistic magnetohydrodynamics (GRMHD) equations using distributed parallel adaptive mesh refinement (AMR). The fluid equations are solved using a finite difference Convex ENO…

We present the AGATE simulation code, a Python-based framework developed primarily for solving the magnetohydrodynamics (MHD) equations while maintaining adaptability to other equation sets. The code employs a modular, object-oriented…

In this paper we present a major update to the general relativistic magnetohydrodynamics (GRMHD) code cuHARM, which adds fully covariant treatment of radiation transport and the subsequent radiation backreaction on the dynamics of the…

The ExaNeSt and EuroExa H2020 EU-funded projects aim to design and develop an exascale ready computing platform prototype based on low-energy-consumption ARM64 cores and FPGA accelerators. We participate in the application-driven design of…

A new code and methodology are introduced for solving the general relativistic magnetohydrodynamic (GRMHD) equations in fixed background spacetimes using time-explicit, finite-volume discretization. The code has options for solving the…

The advent of high-fidelity imaging of supermassive black holes calls for efficient and robust data-analysis methods. In this work, we use $\texttt{Jipole}$, a differentiable, $\texttt{ipole}$-based radiative transfer code, to enable…

The Einstein and Maxwell equations are both systems of hyperbolic equations which need to satisfy a set of elliptic constraints throughout evolution. However, while electrodynamics (EM) and magnetohydrodynamics (MHD) have benefited from a…

We present a mixed-precision implementation of the high-order discontinuous Galerkin method with ADER time stepping (ADER-DG) for solving hyperbolic systems of partial differential equations (PDEs) in the hyperbolic PDE engine ExaHyPE. The…

Associated to every stationary extremal black hole is a unique near-horizon geometry, itself a solution of the field equations. These latter spacetimes are more tractable to analyze and most importantly, retain properties of the original…

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…