The Parallel Compact Object CALculator: An Efficient General Relativistic Initial Data Solver for Compact Objects

Abstract

Every numerical general relativistic investigation starts from the solution of the initial value equations at a given time. Astrophysically relevant initial values for different systems lead to distinct set of equations that obey specific assumptions tied to the particular problem. Therefore a robust and efficient solver for a variety of strongly gravitating sources is needed. In this work we present the OpenMP version of the Compact Object CALculator (COCAL) on shared memory processors. We performed extensive profiling of the core COCAL modules in order to identify bottlenecks in efficiency which we addressed. Using modest resources, the new parallel code achieves speedups approximately one order of magnitude relative to the original serial COCAL code, which is crucial for parameter studies of computationally expensive systems, such as magnetized neutron stars, as well as its further development towards more realistic scenarios. As a novel example of our new code we compute a binary quark system where each companion has a dimensionless spin of $0.43$ aligned with the orbital angular momentum.

26 pages, 6 figures

References in corpus (15)

• GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral
• Multi-messenger Observations of a Binary Neutron Star Merger
• Gravitational Waves and Gamma-rays from a Binary Neutron Star Merger: GW170817 and GRB 170817A
• Estimating the Contribution of Dynamical Ejecta in the Kilonova Associated with GW170817
• A constrained scheme for Einstein equations based on Dirac gauge and spherical coordinates
• Waveless Approximation Theories of Gravity
• Initial data for binary neutron stars with arbitrary spins
• Modeling differential rotations of compact stars in equilibriums
• Equilibrium solutions of relativistic rotating stars with mixed poloidal and toroidal magnetic fields
• Numerical method for binary black hole/neutron star initial data: Code test
• Quasi-equilibrium models for triaxially deformed rotating compact stars
• Evolution of equal mass binary bare quark stars in full general relativity: could a supramassive merger remnant experience prompt collapse?
• Evolution of bare quark stars in full general relativity: Single star case
• NRPyElliptic: A Fast Hyperbolic Relaxation Elliptic Solver for Numerical Relativity, I: Conformally Flat, Binary Puncture Initial Data
• Equilibriums of extremely magnetized compact stars with force-free magnetotunnels