superB/NRPy: Scalable, Task-Based Numerical Relativity for 3G Gravitational Wave Science

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Jadoo, Nishita, Jacques, Terrence Pierre, Etienne, Zachariah B.
Format: Preprint
Published: 2025
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866913960131821568
author Jadoo, Nishita
Jacques, Terrence Pierre
Etienne, Zachariah B.
author_facet Jadoo, Nishita
Jacques, Terrence Pierre
Etienne, Zachariah B.
contents Modern gravitational-wave science demands increasingly accurate and computationally intensive numerical relativity (NR) simulations. The Python-based, open-source NRPy framework generates optimized C/C++ code for NR, including the complete NR code BlackHoles@Home (BH@H), which leverages curvilinear coordinates well-suited to many astrophysical scenarios. Historically, BH@H was limited to single-node OpenMP CPU parallelism. To address this, we introduce superB, an open-source extension to NRPy that enables automatic generation of scalable, task-based, distributed-memory Charm++ code from existing BH@H modules. The generated code partitions the structured grids used by NRPy/BH@H, managing communication between them. Its correctness is validated through bit-identical results with the standard OpenMP version on a single node and via a head-on binary black hole simulation in cylindrical-like coordinates, accurately reproducing quasi-normal modes (up to $\ell=8$). The superB/NRPy-generated code demonstrates excellent strong scaling, achieving an $\approx 45$x speedup on 64 nodes (7168 cores) compared to the original single-node OpenMP code for a large 3D vacuum test. This scalable infrastructure benefits demanding simulations and lays the groundwork for future multi-patch grid support, targeting long inspirals, extreme parameter studies, and rapid follow-ups. This infrastructure readily integrates with other NRPy/BH@H-based projects, enabling performant scaling for the general relativistic hydrodynamics code GRoovy, and facilitating future coupling with GPU acceleration via the NRPy-CUDA project.
format Preprint
id arxiv_https___arxiv_org_abs_2505_00097
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle superB/NRPy: Scalable, Task-Based Numerical Relativity for 3G Gravitational Wave Science
Jadoo, Nishita
Jacques, Terrence Pierre
Etienne, Zachariah B.
General Relativity and Quantum Cosmology
Instrumentation and Methods for Astrophysics
Computational Physics
83C57, 83C35
Modern gravitational-wave science demands increasingly accurate and computationally intensive numerical relativity (NR) simulations. The Python-based, open-source NRPy framework generates optimized C/C++ code for NR, including the complete NR code BlackHoles@Home (BH@H), which leverages curvilinear coordinates well-suited to many astrophysical scenarios. Historically, BH@H was limited to single-node OpenMP CPU parallelism. To address this, we introduce superB, an open-source extension to NRPy that enables automatic generation of scalable, task-based, distributed-memory Charm++ code from existing BH@H modules. The generated code partitions the structured grids used by NRPy/BH@H, managing communication between them. Its correctness is validated through bit-identical results with the standard OpenMP version on a single node and via a head-on binary black hole simulation in cylindrical-like coordinates, accurately reproducing quasi-normal modes (up to $\ell=8$). The superB/NRPy-generated code demonstrates excellent strong scaling, achieving an $\approx 45$x speedup on 64 nodes (7168 cores) compared to the original single-node OpenMP code for a large 3D vacuum test. This scalable infrastructure benefits demanding simulations and lays the groundwork for future multi-patch grid support, targeting long inspirals, extreme parameter studies, and rapid follow-ups. This infrastructure readily integrates with other NRPy/BH@H-based projects, enabling performant scaling for the general relativistic hydrodynamics code GRoovy, and facilitating future coupling with GPU acceleration via the NRPy-CUDA project.
title superB/NRPy: Scalable, Task-Based Numerical Relativity for 3G Gravitational Wave Science
topic General Relativity and Quantum Cosmology
Instrumentation and Methods for Astrophysics
Computational Physics
83C57, 83C35
url https://arxiv.org/abs/2505.00097