Related papers: Improved Aberth-Ehrlich root-finding algorithm and…
We introduce a new adaptive and fully Bayesian grid-based method to model strong gravitational lenses with extended images. The primary goal of this method is to quantify the level of luminous and dark-mass substructure in massive galaxies,…
This study examines the characterization of binary star systems using Spectral Energy Distributions (SEDs), a technique increasingly essential with the rise of large-scale astronomical surveys. Binaries can emit flux at different regions of…
Evolutionary algorithms (EAs) are general-purpose problem solvers that usually perform an unbiased search. This is reasonable and desirable in a black-box scenario. For combinatorial optimization problems, often more knowledge about the…
We show that the higher dimensional Weierstrass and Ehrlich-Aberth methods for finding roots of polynomials have infinite orbits that diverge to infinity. This is possible for the Jacobi update scheme (all coordinates are updated in…
We develop a reduction multigrid based on approximate ideal restriction (AIR) for use with asymmetric linear systems. We use fixed-order GMRES polynomials to approximate $A_\textrm{ff}^{-1}$ and we use these polynomials to build grid…
Gaudi & Gould (1997) showed that close companions of remote binary systems can be efficiently detected by using gravitational microlensing via the deviations in the lensing light curves induced by the existence of the lens companions. In…
In modern astronomy, the quantity of data collected has vastly exceeded the capacity for manual analysis, necessitating the use of advanced artificial intelligence (AI) techniques to assist scientists with the most labor-intensive tasks. AI…
A numerical method ADER-DG with a local DG predictor for solving a DAE system has been developed, which was based on the formulation of ADER-DG methods using a local DG predictor for solving ODE and PDE systems. The basis functions were…
We present the first light curve solutions of the four binary systems observed and catalogued by the All Sky Automated Survey program. The light curves are analysed by using PHOEBE software, which is based on the Wilson-Devinney method. The…
Massive black hole binaries are expected to be observable as periodic AGN in time-domain photometric surveys. Periodicities may originate from different physical processes, including the intermittent gas feeding of the black holes caused by…
The phase retrieval problem, where one aims to recover a complex-valued image from far-field intensity measurements, is a classic problem encountered in a range of imaging applications. Modern phase retrieval approaches usually rely on…
We revisit the problem of model-based object recognition for intensity images and attempt to address some of the shortcomings of existing Bayesian methods, such as unsuitable priors and the treatment of residuals with a non-robust error…
We present a new spectral-method-based algorithm for finding apparent horizons in three-dimensional space-like hypersurfaces without symmetries. While there are already a wide variety of algorithms for finding apparent horizons, our new…
Equilibria of binary neutron stars in close circular orbits are computed numerically in a waveless formulation: The full Einstein-relativistic-Euler system is solved on an initial hypersurface to obtain an asymptotically flat form of the…
An extra-solar planet can be detected by microlensing because the planet can perturb the smooth lensing light curve created by the primary lens. However, it was shown by Gaudi that a subset of binary-source events can produce light curves…
The Bayesian gravitational shear estimation algorithm developed by Bernstein and Armstrong (2014) can potentially be used to overcome multiplicative noise bias and recover shear using very low signal-to-noise ratio (S/N) galaxy images. In…
Path planning is an active area of research essential for many applications in robotics. Popular techniques include graph-based searches and sampling-based planners. These approaches are powerful but have limitations. This paper continues…
Mass measurements of gravitational microlenses require one to determine the microlens parallax $\pie$, but precise $\pie$ measurement, in many cases, is hampered due to the subtlety of the microlens-parallax signal combined with the…
Integration operational matrix methods based on Zernike polynomials are used to determine approximate solutions of a class of non-homogeneous partial differential equations (PDEs) of first and second order. Due to the nature of the Zernike…
In this paper we propose a new class of iterative regularization methods for solving ill-posed linear operator equations. The prototype of these iterative regularization methods is in the form of second order evolution equation with a…