Related papers: Pseudo-spectral apparent horizon finders: an effic…
The search for continuous gravitational waves demands computationally efficient algorithms that can handle highly non-linear parameter spaces. Loosely coherent algorithms establish upper limits and detect signals by analyzing families of…
Grover's unstructured search algorithm is one of the best examples to date for the superiority of quantum algorithms over classical ones. Its applicability, however, has been questioned by many due to its oracular nature. We propose a…
We present a new method of extracting gravitational radiation from three-dimensional numerical relativity codes and providing outer boundary conditions. Our approach matches the solution of a Cauchy evolution of Einstein's equations to a…
We present a detailed analytical study of spherically symmetric self-similar solutions in the SU(2) sigma model coupled to gravity. Using a shooting argument we prove that there is a countable family of solutions which are analytic inside…
We introduce a new computational technique for searching for faint moving sources in astronomical images. Starting from a maximum likelihood estimate for the probability of the detection of a source within a series of images, we develop a…
In this article, we consider the numerical approximation of far-field statistics for acoustic scattering problems in the case of random obstacles. In particular, we consider the computation of the expected far-field pattern and the expected…
A method for detecting and approximating fault lines or surfaces, respectively, or decision curves in two and three dimensions with guaranteed accuracy is presented. Reformulated as a classification problem, our method starts from a set of…
Motivated by nonconvex, inconsistent feasibility problems in imaging, the relaxed alternating averaged reflections algorithm, or relaxed Douglas-Rachford algorithm (DR$\lambda$), was first proposed over a decade ago. Convergence results for…
The approximate solution of the Cauchy problem for second-order evolution equations is performed, first of all, using three-level time approximations. Such approximations are easily constructed and relatively uncomplicated to investigate…
We give efficient data-oblivious algorithms for several fundamental geometric problems that are relevant to geographic information systems, including planar convex hulls and all-nearest neighbors. Our methods are "data-oblivious" in that…
We present two graph drawing algorithms based on the recently defined "grand-Schnyder woods", which are a far-reaching generalization of the classical Schnyder woods. The first is a straight-line drawing algorithm for plane graphs with…
This paper investigates the global dynamics of the apparent horizon. We present an approach to establish its existence and its long-term behaviors. Our apparent horizon is constructed by solving the marginally outer trapped surface (MOTS)…
We advocate for a new paradigm of cosmological likelihood-based inference, leveraging recent developments in machine learning and its underlying technology, to accelerate Bayesian inference in high-dimensional settings. Specifically, we…
We develop and analyze several different second-order algorithms for computing a near-optimal solution path of a convex parametric optimization problem with smooth Hessian. Our algorithms are inspired by a differential equation perspective…
The isolated horizon formalism recently introduced by Ashtekar et al. aims at providing a quasi-local concept of a black hole in equilibrium in an otherwise possibly dynamical spacetime. In this formalism, a hierarchy of geometrical…
We prove that the surface gravity of a compact non-degenerate Cauchy horizon in a smooth vacuum spacetime, can be normalized to a non-zero constant. This result, combined with a recent result by Oliver Petersen and Istv\'an R\'acz, end up…
We propose an efficient semi-Lagrangian method for solving the two-dimensional incompressible Euler equations with high precision on a coarse grid. The new approach evolves the flow map using the gradient-augmented level set method (GALSM).…
Structured sparse optimization is an important and challenging problem for analyzing high-dimensional data in a variety of applications such as bioinformatics, medical imaging, social networks, and astronomy. Although a number of structured…
High-dimensional linear regression under heavy-tailed noise or outlier corruption is challenging, both computationally and statistically. Convex approaches have been proven statistically optimal but suffer from high computational costs,…
Detection of interesting (e.g., coherent or anomalous) clusters has been studied extensively on plain or univariate networks, with various applications. Recently, algorithms have been extended to networks with multiple attributes for each…