Related papers: A Sampling Strategy for High-Dimensional Spaces Ap…
We introduce a deep generative framework for high-dimensional Bayesian inference that enables efficient posterior sampling. As telescopes and simulations rapidly expand the volume and resolution of astrophysical data, fast simulation-based…
We study a class of Metropolis-Hastings algorithms for target measures that are absolutely continuous with respect to a large class of non-Gaussian prior measures on Banach spaces. The algorithm is shown to have a spectral gap in a…
An easy-to-implement form of the Metropolis Algorithm is described which, unlike most standard techniques, is well suited to sampling from multi-modal distributions on spaces with moderate numbers of dimensions (order ten) in environments…
We study a generalized framework for structured sparsity. It extends the well-known methods of Lasso and Group Lasso by incorporating additional constraints on the variables as part of a convex optimization problem. This framework provides…
Uncertainty quantification for large-scale inverse problems remains a challenging task. For linear inverse problems with additive Gaussian noise and Gaussian priors, the posterior is Gaussian but sampling can be challenging, especially for…
We propose a new random method to minimize deterministic continuous functions over subsets $\mathcal{S}$ of high-dimensional space $\mathbb{R}^K$ without assuming convexity. Our procedure alternates between a Global Search (GS) regime to…
Coulomb and log-gases are exchangeable singular Boltzmann-Gibbs measures appearing in mathematical physics at many places, in particular in random matrix theory. We explore experimentally an efficient numerical method for simulating such…
We present Lenstronomy, a multi-purpose open-source gravitational lens modeling python package. Lenstronomy is able to reconstruct the lens mass and surface brightness distributions of strong lensing systems using forward modelling.…
Learning to sample from complex unnormalized distributions is a fundamental challenge in computational physics and machine learning. While score-based and variational methods have achieved success in continuous domains, extending them to…
Most surveys for multiply-imaged gravitational lenses, outside of rich galaxy clusters, are based on sifting through large samples of distant sources to identify the rare examples of lensing. An alternative strategy, based on the selection…
We revisit the issue of non-parametric gravitational lens reconstruction and present a new method to obtain the cluster mass distribution using strong lensing data without using any prior information on the underlying mass. The method…
Gravitational lensing and stellar dynamics are two independent methods, based solely on gravity, to study the mass distributions of galaxies. Both methods suffer from degeneracies, however, that are difficult to break. In this paper, we…
The maximum-entropy method is applied to the problem of reconstructing the projected mass density of a galaxy cluster using its gravitational lensing effects on background galaxies. We demonstrate the method by reconstructing the mass…
Markov Chain Monte Carlo (MCMC) techniques have long been studied in computational geometry subjects whereabouts the problems to be studied are complex geometric objects which by their nature require optimized techniques to be deployed or…
In this paper, we introduce a new and efficient data augmentation approach to the posterior inference of the models with shape parameters when the reciprocal gamma function appears in full conditional densities. Our approach is to…
In this work, we examine sampling problems with non-smooth potentials. We propose a novel Markov chain Monte Carlo algorithm for sampling from non-smooth potentials. We provide a non-asymptotical analysis of our algorithm and establish a…
Traditional methods for unsupervised learning of finite mixture models require to evaluate the likelihood of all components of the mixture. This becomes computationally prohibitive when the number of components is large, as it is, for…
Context. Existing samples of strong lenses have been assembled by giving priority to sample size, at the cost of having a complex selection function. With the advent of the next generation of wide-field photometric surveys, however, it…
In this paper, we study the development of efficient multiscale methods for flows in heterogeneous media. Our approach uses the Generalized Multiscale Finite Element (GMsFEM) framework. The main idea of GMsFEM is to approximate the solution…
Sparse regression based on global-local shrinkage priors are increasingly used for Bayesian modeling of modern high-dimensional data, but scaling up the Gibbs sampler for posterior inference remains a challenge. While much effort has gone…