A Sampling Strategy for High-Dimensional Spaces Applied to Free-Form Gravitational Lensing
Instrumentation and Methods for Astrophysics
2015-06-05 v1
Abstract
We present a novel proposal strategy for the Metropolis-Hastings algorithm designed to efficiently sample general convex polytopes in 100 or more dimensions. This improves upon previous sampling strategies used for free-form reconstruction of gravitational lenses, but is general enough to be applied to other fields. We have written a parallel implementation within the lens modeling framework GLASS. Testing shows that we are able to produce uniform uncorrelated random samples which are necessary for exploring the degeneracies inherent in lens reconstruction.
Cite
@article{arxiv.1207.1722,
title = {A Sampling Strategy for High-Dimensional Spaces Applied to Free-Form Gravitational Lensing},
author = {Mario Lubini and Jonathan Coles},
journal= {arXiv preprint arXiv:1207.1722},
year = {2015}
}
Comments
10 pages, 9 figures. Accepted for publication in MNRAS