Related papers: GRTresna: An open-source code to solve the initial…
Currently there is considerable interest in making use of many-core processor architectures, such as Nvidia and AMD graphics processing units (GPUs) for scientific computing. In this work we explore the use of the Open Computing Language…
We obtain the Schwarzschild solution based on teleparallel gravity (TG) theory formulated in a space-time with torsion only. The starting point is the Poincar\UNICODE{0xe9} gauge theory (PGT).The general structure of TG and its connection…
GMRES is a popular Krylov subspace method for solving linear systems of equations involving a general non-Hermitian coefficient matrix. The conventional bounds on GMRES convergence involve polynomial approximation problems in the complex…
Construction of astrophysically realistic initial data remains a central problem when modelling the merger and eventual coalescence of binary black holes in numerical relativity. The objective of this paper is to provide astrophysically…
We are entering an era where the numerical construction of generic spacetimes is becoming a reality. The use of computer simulations, in principle, allows us to solve Einstein equations in their full generality and unravel important…
Recent years have seen a significant progress in the development of general relativistic codes for the numerical solution of the equations of magnetohydrodynamics in spacetimes with high and dynamical curvature. These codes are valuable…
The progression of scientific computing resources has enabled the numerical approximation of mathematical models describing complex physical phenomena. A significant portion of researcher time is typically dedicated to the development of…
Detailed observations of phenomena involving black holes, be it via gravitational waves or more traditional electromagnetic means, can probe the strong field regime of the gravitational interaction. The prediction of features in such…
The standard paradigm of cosmology assumes General Relativity (GR) is a valid theory for gravity at scales in which it has not been properly tested. Developing novel tests of GR and its alternatives is crucial if we want to give strength to…
The construction of constraint-satisfying initial data is an essential element for the numerical exploration of the dynamics of compact-object binaries. While several codes have been developed over the years to compute generic…
Iharm3D is an open-source C code for simulating black hole accretion systems in arbitrary stationary spacetimes using ideal general-relativistic magnetohydrodynamics (GRMHD). It is an implementation of the HARM ("High Accuracy Relativistic…
I present the publicly available code GLADIS (GLobal Accretion Disk Instability Simulation) developed in my reserach group over the years 2002-2017. It can be freely downloaded and modified by the users via the link from the Astrophysics…
Modeling self-gravitating gas flows is essential to answering many fundamental questions in astrophysics. This spans many topics including planet-forming disks, star-forming clouds, galaxy formation, and the development of large-scale…
This paper shows that the OSGA algorithm -- which uses first-order information to solve convex optimization problems with optimal complexity -- can be used to efficiently solve arbitrary bound-constrained convex optimization problems. This…
We address the generation of initial conditions (ICs) for GRAMSES, a code for nonlinear general relativistic (GR) $N$-body cosmological simulations recently introduced in Ref. [1]. GRAMSES adopts a constant mean curvature slicing with a…
We introduce a new, open-source computational general relativity framework for the Wolfram Language called Gravitas, which boasts a number of novel and distinctive features as compared to the many pre-existing computational and numerical…
The need to smoothly cover a computational domain of interest generically requires the adoption of several grids. To solve the problem of interest under this grid-structure one must ensure the suitable transfer of information among the…
We develop an algorithm solving the 3x3 real symmetric eigenproblem. This is a common problem and in certain applications it must be solved many thousands of times, see for example \cite{tripref} where each element in a finite element grid…
We present approximate analytical solutions to the Hamiltonian and momentum constraint equations, corresponding to systems composed of two black holes with arbitrary linear and angular momentum. The analytical nature of these initial data…
We present the extension of GR-Athena++ to general-relativistic magnetohydrodynamics (GRMHD) for applications to neutron star spacetimes. The new solver couples the constrained transport implementation of Athena++ to the Z4c formulation of…