Related papers: Fast Multipole Method for Gravitational Lensing. A…
In this paper, a fast multipole method (FMM) is proposed to compute long-range interactions of wave sources embedded in 3-D layered media. The layered media Green's function for the Helmholtz equation, which satisfies the transmission…
Large matrix multiplication is a cornerstone of modern machine learning workloads, yet traditional approaches suffer from cubic computational complexity (e.g., $\mathcal{O}(n^3)$ for a matrix of size $n\times n$). We present Low-Rank GEMM,…
The line-of-sight structure of the Galactic magnetic field (GMF) can be studied using Faraday rotation measure (RM) grids. We analyze how the choice of interpolation kernel can affect the accuracy and reliability of reconstructed RM maps.…
Removing optical and atmospheric blur from galaxy images significantly improves galaxy shape measurements for weak gravitational lensing and galaxy evolution studies. This ill-posed linear inverse problem is usually solved with…
This paper is concerned with the modeling errors appeared in the numerical methods of inverse medium scattering problems (IMSP). Optimization based iterative methods are wildly employed to solve IMSP, which are computationally intensive due…
We propose a new, two-stage algorithm for inverting strong gravitational lenses. The key to the algorithm is decoupling the effects of lens magnification and intrinsic structure in the background source in the appearance of lensed arcs.…
The gravitational field of compact objects, such as primordial black holes, can create multiple images of background sources. For transients such as fast radio bursts (FRBs), these multiple images can be resolved in the time domain. Under…
We present the observational results of a survey for strong gravitational lens systems consisting of extended emission-line galaxies lensed by intervening early-type galaxies, conducted using integral field units (IFUs) of the Magellan…
Learning a Gaussian Mixture Model (GMM) is hard when the number of parameters is too large given the amount of available data. As a remedy, we propose restricting the GMM to a Gaussian Markov Random Field Mixture Model (GMRF-MM), as well as…
As research into magnetic thin films and spintronics devices is moving from single to multiple magnetic layers, there is a need for micromagnetics modelling tools specifically designed to efficiently handle magnetic multilayers. Here we…
We describe a new method for measuring galaxy magnification due to weak gravitational lensing. Our method makes use of a tight scaling relation between galaxy properties that are modified by gravitational lensing, such as apparent size, and…
We demonstrate a new computational illumination technique that achieves large space-bandwidth-time product, for quantitative phase imaging of unstained live samples in vitro. Microscope lenses can have either large field of view (FOV) or…
The kernel-independent fast multipole method (KIFMM) proposed in [1] is of almost linear complexity. In the original KIFMM the time-consuming M2L translations are accelerated by FFT. However, when more equivalent points are used to achieve…
We present an analytical formulation of gravitational lensing using familiar triaxial power-law mass distributions, where the 3-dimensional mass density is given by $\rho(X,Y,Z) = \rho_0 [1 + (\frac{X}{a})^2 + (\frac{Y}{b})^2 +…
Non-uniform fast Fourier Transform (NUFFT) and inverse NUFFT (INUFFT) algorithms, based on the Fast Multipole Method (FMM) are developed and tested. Our algorithms are based on a novel factorization of the FFT kernel, and are implemented…
In this work we introduced a new proposal to study the gravitational lensing theory by spherical lenses, starting from its surface mass density $\Sigma(x)$ written in terms of a decreasing function $f$ of a dimensionless coordinate $x$ on…
In this paper, we present a numerical algorithm for the accurate and efficient computation of the convolution of the frequency domain layered media Green's function with a given density function. Instead of compressing the convolution…
Given an observed gravitational lens mirage produced by a foreground deflector (cf. galaxy, quasar, cluster,...), it is possible via numerical lens inversion to retrieve the real source image, taking full advantage of the magnifying power…
Composite function minimization captures a wide spectrum of applications in both computer vision and machine learning. It includes bound constrained optimization and cardinality regularized optimization as special cases. This paper proposes…
A common problem in cosmology is to integrate the product of two or more spherical Bessel functions (sBFs) with different configuration-space arguments against the power spectrum or its square, weighted by powers of wavenumber. Naively…