Related papers: Fast and Slow solutions in General Relativity: The…
Approximations are commonly employed to find approximate solutions to the Einstein equations. These solutions, however, are usually only valid in some specific spacetime region. A global solution can be constructed by gluing approximate…
An optimization program is used to re-adjust the initial conditions, in order to reproduce as closely as possible the predictions of a complete ephemeris by using simplified equations of motion in the numerical integration. The adjustment…
We present a model reduction approach for the real-time solution of time-dependent nonlinear partial differential equations (PDEs) with parametric dependencies. The approach integrates several ingredients to develop efficient and accurate…
Strongly coupled gravitational systems describe Einstein gravity and matter in the limit that Newton's constant G is assumed to be very large. The nonlinear evolution of these systems may be solved analytically in the classical and…
Spacetime is foliated by spatial hypersurfaces in the 3+1 split of General Relativity. The initial value problem then consists of specifying initial data for all relevant fields on one such a spatial hypersurface. These fields are the…
Successive quadratic approximations, or second-order proximal methods, are useful for minimizing functions that are a sum of a smooth part and a convex, possibly nonsmooth part that promotes regularization. Most analyses of iteration…
A delayed term in a differential equation reflects the fact that information takes significant time to travel from one place to another within a process being studied. Despite de apparent similarity with ordinary differential equations,…
We consider parameterized variational inverse problems that are constrained by partial differential equations (PDEs). We seek to efficiently compute the solution of the inverse problem when auxiliary model parameters, which appear in the…
We consider the approximation to the solution of the initial boundary value problem for the heat equation with right hand side and initial condition that merely belong to $L^1$. Due to the low integrability of the data, to guarantee…
We introduce a generic scheme for accelerating first-order optimization methods in the sense of Nesterov, which builds upon a new analysis of the accelerated proximal point algorithm. Our approach consists of minimizing a convex objective…
In Einstein theory of gravity the initial configuration of metric field and its time derivative are related to matter configuration by four equations called constraints. We use the method of conformal metrics (York Method) to solve…
We exploit analogies between first-order algorithms for constrained optimization and non-smooth dynamical systems to design a new class of accelerated first-order algorithms for constrained optimization. Unlike Frank-Wolfe or projected…
There is a tendency to write the equations of general relativity as a first order symmetric system of time dependent partial differential equations. However, for numerical reasons, it might be advantageous to use a second order formulation…
An efficient approximate version of implicit Taylor methods for initial-value problems of systems of ordinary differential equations (ODEs) is introduced. The approach, based on an approximate formulation of Taylor methods, produces a…
There are many complementary approaches to the construction of solutions to the field equations of general relativity. Among these, numerical approximation offers the only possibility to compute a variety of dynamical spacetimes, and so has…
We present here the explicit parametric solutions of second order differential equations invariant under time translation and rescaling and third order differential equations invariant under time translation and the two homogeneity…
The purpose of this paper is to demonstrate a new method of generating exact solutions to the Einstein's equations obtained by the Hamiltonian reduction. The key element to the successful Hamiltonian reduction is finding the privileged…
Nonlinear systems of partial differential equations (PDEs) may permit several distinct solutions. The typical current approach to finding distinct solutions is to start Newton's method with many different initial guesses, hoping to find…
The Einstein equations have proven surprisingly difficult to solve numerically. A standard diagnostic of the problems which plague the field is the failure of computational schemes to satisfy the constraints, which are known to be…
The conformal method for constructing initial data for Einstein's equations is presented in both the Hamiltonian and Lagrangian picture (extrinsic curvature decomposition and conformal thin sandwich formalism, respectively), and advantages…