Related papers: FAST-PT: a novel algorithm to calculate convolutio…
The polar decomposition of a matrix is a key element in the quantum linear algebra toolbox. We show that the problem of quantum polar decomposition, recently studied in Lloyd et al. [LBP+20], has a simple and concise implementation via the…
In this work we detail the application of a fast convolution algorithm computing high dimensional integrals to the context of multiplicative noise stochastic processes. The algorithm provides a numerical solution to the problem of…
We introduce a fast algorithm for computing volume potentials - that is, the convolution of a translation invariant, free-space Green's function with a compactly supported source distribution defined on a uniform grid. The algorithm relies…
We introduce a highly-parallelizable architecture for estimating parameters of compact binary coalescence using gravitational-wave data and waveform models. Using a spherical harmonic mode decomposition, the waveform is expressed as a sum…
We present a fast and accurate formulation for computing the nonlinear matter power spectrum at one-loop order based on Unified Lagrangian Perturbation Theory (ULPT). ULPT decomposes the density field into the Jacobian deviation, capturing…
In the previous paper [arXiv:2210.10435], the nonlinear perturbation theory of cosmological density field is generalized to include the tensor-valued bias of astronomical objects, such as spins and shapes of galaxies and any other tensors…
Over the last two decades, several fast, robust, and high-order accurate methods have been developed for solving the Poisson equation in complicated geometry using potential theory. In this approach, rather than discretizing the partial…
Standard cosmological perturbation theory (SPT) for the Large Scale Structure (LSS) of the Universe fails at small scales (UV) due to strong nonlinearities and to multistreaming effects. In Pietroni et al. 2011 a new framework was proposed…
We present quantum algorithms for the estimation of n-time correlation functions, the local and non-local density of states, and dynamical linear response functions. These algorithms are all based on block-encodings - a versatile technique…
Bayesian inference, while foundational to probabilistic reasoning, is often hampered by the computational intractability of posterior distributions, particularly through the challenging evidence integral. Conventional approaches like Markov…
We present a formulation of cosmological perturbation theory where the Boltzmann hierarchies that evolve the neutrino phase-space distributions are replaced by integrals that can be evaluated easily with fast Fourier transforms. The…
Image subtraction in astronomy is a tool for transient object discovery such as asteroids, extra-solar planets and supernovae. To match point spread functions (PSFs) between images of the same field taken at different times a convolution…
Calculating perturbation response properties of materials from first principles provides a vital link between theory and experiment, but is bottlenecked by the high computational cost. Here a general framework is proposed to perform density…
Modified gravity models with scale-dependent linear growth typically exhibit an enhancement in the power spectrum beyond a certain scale. The conventional methods for extracting cosmological information usually involve inferring modified…
A new approach to cosmological perturbation theory has been recently introduced by Bartelmann et al., relying on nonequilibrium statistical theory of classical particles, and treating the gravitational interaction perturbatively. They…
Nonuniform Fourier data are routinely collected in applications such as magnetic resonance imaging, synthetic aperture radar, and synthetic imaging in radio astronomy. To acquire a fast reconstruction that does not require an online inverse…
We propose two easy-to-implement fast algorithms based on moment-matching to compute the nonlocal potential $\varphi(\textbf{x})=(U\ast \rho)(\textbf{x})$ on bounded domain, where the kernel $U$ is singular at the origin and the density…
The nonlinear Fourier transform (NFT) has recently gained significant attention in fiber optic communications and other engineering fields. Although several numerical algorithms for computing the NFT have been published, the design of…
Optical imaging systems are inherently limited in their resolution due to the point spread function (PSF), which applies a static, yet spatially-varying, convolution to the image. This degradation can be addressed via Convolutional Neural…
The nonlinear mass is a characteristic scale in halo formation that has wide-ranging applications across cosmology. Naively, computing it requires repeated numerical integration to calculate the variance of the power spectrum on different…