Related papers: GRTresna: An open-source code to solve the initial…
We present Lenstronomy, a multi-purpose open-source gravitational lens modeling python package. Lenstronomy is able to reconstruct the lens mass and surface brightness distributions of strong lensing systems using forward modelling.…
Real-space grids are a powerful alternative for the simulation of electronic systems. One of the main advantages of the approach is the flexibility and simplicity of working directly in real space where the different fields are discretized…
We construct approximate analytical solutions to the constraint equations of general relativity for binary black holes of arbitrary mass ratio in quasicircular orbit. We adopt the puncture method to solve the constraint equations in the…
Black holes and other compact objects are powerful tools to observationally test Einsteins theory of General Relativity. We develop raytracing code to create visual images of compact objects that are solutions of Einsteins field equations.…
The rapid development of gravitational wave astronomy provides the unique opportunity of exploring the dynamics of the Universe using clustering properties of coalescing binary black hole mergers. Gravitational wave data, along with…
This thesis focuses on the application of numerical relativity methods to the solutions of problems in strong gravity. Our goal is the study of mergers of compact objects in the strong field regime where non-linear dynamics manifest and…
We present the new open-source spherically-symmetric general-relativistic (GR) hydrodynamics code GR1D. It is based on the Eulerian formulation of GR hydrodynamics (GRHD) put forth by Romero-Ibanez-Gourgoulhon and employs radial-gauge,…
We present a new open source C library \texttt{msolve} dedicated to solving multivariate polynomial systems of dimension zero through computer algebra methods. The core algorithmic framework of \texttt{msolve} relies on Gr\''obner bases and…
We use rigorous techniques from numerical analysis of hyperbolic equations in bounded domains to construct stable finite-difference schemes for Numerical Relativity, in particular for their use in black hole excision. As an application, we…
Low-luminosity Active Galactic Nuclei (AGN) are believed to be surrounded by a collisionless, highly magnetized accretion flow. As a result, Particle-in-Cell simulations are the best tools to study the immediate vicinity of the event…
We report on a new 3D numerical code designed to solve the Einstein equations for general vacuum spacetimes. This code is based on the standard 3+1 approach using cartesian coordinates. We discuss the numerical techniques used in developing…
We introduce a new code, ECOSMOG, to run N-body simulations for a wide class of modified gravity and dynamical dark energy theories. These theories generally have one or more new dynamical degrees of freedom, the dynamics of which are…
A benchmark problem for numerical relativity has been the head-on collision of two black holes starting from the ``Misner initial data,'' a closed form momentarily stationary solution to the constraint equations with an adjustable closeness…
We describe the Blacklight code, intended for post-processing general-relativistic magnetohydrodynamic simulation data. Beyond polarized ray tracing of synchrotron radiation, it can produce a number of outputs that aid in analyzing data…
Many applications of generalised linear models (GLMs) can be improved by applying constraints that impose assumptions on the associations or improve consistency of the estimators. Yet, there are still barriers to the implementation and…
We present a new open-source axisymmetric general relativistic hydrodynamics code Gmunu (General-relativistic multigrid numerical solver) which uses a multigrid method to solve the elliptic metric equations in the conformally flat condition…
The solution of large, sparse constrained least-squares problems is a staple in scientific and engineering applications. However, currently available codes for such problems are proprietary or based on MATLAB. We announce a freely available…
We give a pedagogical introduction of the essential features of General Theory of Relativity (GTR) in the format of an undergraduate (UG) project. A set of simple MATHEMATICA code is developed which enables the UG students to calculate the…
The Einstein constraint equations describe the space of initial data for the evolution equations, dictating how space should curve within spacetime. Under certain assumptions, the constraints reduce to a scalar quasilinear parabolic…
This paper presents a fast, economical particle-multiple-mesh N-body code optimized for large-N modelling of collisionless dynamical processes, such as black-hole wandering or bar-halo interactions, occurring within isolated galaxies. The…