Related papers: A new code for equilibriums and quasiequilibrium i…
In numerical relativity simulations with non-trivial matter configurations, one must solve the Hamiltonian and momentum constraints of the ADM formulation for the metric variables in the initial data. We introduce a new scheme based on the…
A new algorithms for computing discrete logarithms on elliptic curves defined over finite fields is suggested. It is based on a new method to find zeroes of summation polynomials. In binary elliptic curves one is to solve a cubic system of…
Based on a newly proposed spatial data model - spatial chromatic model (SCM), we developed a spatial coding scheme, called full-coded ordinary arranged chromatic diagram (full-OACD). Full-OACD is a type of spatial tessellation, where space…
I discuss the conformal approach to the numerical simulation of radiating isolated systems in general relativity. The method is based on conformal compactification and a reformulation of the Einstein equations in terms of rescaled…
Simulation of quasicircular compact binaries is a major goal in numerical relativity, as they are expected to constitute most gravitational wave observations. However, given that orbital eccentricity is not well-defined in general…
In this paper, we present a new version of the $i$QIST software package, which is capable of solving various quantum impurity models by using the hybridization expansion (or strong coupling expansion) continuous-time quantum Monte Carlo…
We present a numerical work aiming at the computation of excised initial data for black hole spacetimes in full general relativity, using Dirac gauge in the context of a constrained formalism for the Einstein equations. Introducing the…
We introduce a new Partition of Unity Method for the numerical homogenization of elliptic partial differential equations with arbitrarily rough coefficients. We do not restrict to a particular ansatz space or the existence of a finite…
We obtain an explicit solution of the momentum constraint for conformally flat, maximal slicing, initial data which gives an alternative to the purely longitudinal extrinsic curvature of Bowen and York. The new solution is related, in a…
Isolated horizon conditions enforce the time invariance of both the intrinsic and the extrinsic geometry of a (quasilocal) black hole horizon. Nonexpanding horizons, only requiring the invariance of the intrinsic geometry, have been…
We propose a numerical homogenization method for scalar linear partial differential equations with rough coefficients, that integrates classical coarse-scale solvers with quantum subroutines for fine-scale corrections. Inspired by the…
Partial differential equation (PDE)-constrained optimization, where an optimization problem is subject to PDE constraints, arises in various applications such as design, control, and inference. Solving such problems is computationally…
We propose a new approach, based on the puncture method, to construct black hole initial data in the so-called trumpet geometry, i.e. on slices that asymptote to a limiting surface of non-zero areal radius. Our approach is easy to implement…
We develop a numerical code to compute gravitational waves induced by a particle moving on eccentric inclined orbits around a Kerr black hole. For such systems, the black hole perturbation method is applicable. The gravitational waves can…
We present a characteristic algorithm for computing the perturbation of a Schwarzschild spacetime by means of solving the Teukolsky equation. We implement the algorithm as a characteristic evolution code and apply it to compute the advanced…
This article presents a new high-order accurate algorithm for finding a particular solution to a linear, constant-coefficient partial differential equation (PDE) by means of a convolution of the volumetric source function with the Green's…
A method is introduced for solving Einstein's equations using two distinct coordinate systems. The coordinate basis vectors associated with one system are used to project out components of the metric and other fields, in analogy with the…
We present a new numerical code developed for the evolution of binary black-hole spacetimes using different initial data and evolution techniques. The code is demonstrated to produce state-of-the-art simulations of orbiting and inspiralling…
CylindricalAlgebraicDecomposition.m2 is the first implementation of Cylindrical Algebraic Decomposition (CAD) in Macaulay2. CAD decomposes space into 'cells' where input polynomials are sign-invariant. This package computes an Open CAD…
Quantum error correction has given us a natural language for the emergence of spacetime, but the black hole interior poses a challenge for this framework: at late times the apparent number of interior degrees of freedom in effective field…