Open AccessComputer SciencePhysicsMathematics

Thejs Brinckmann, Julien Lesgourgues

2018.4.19Physics of the Dark Universe

DOI: 10.1016/j.dark.2018.100260

tlooto Summary

Two new ingredients both contributing to improve the performance of Metropolis-Hastings sampling are explained: an adaptation algorithm for the jumping factor, and a calculation of the inverse Fisher matrix, which can be used as a proposal density.

Abstract

MontePython is a parameter inference package for cosmology. We present the latest development of the code over the past couple of years. We explain, in particular, two new ingredients both contributing to improve the performance of Metropolis-Hastings sampling: an adaptation algorithm for the jumping factor, and a calculation of the inverse Fisher matrix, which can be used as a proposal density. We present several examples to show that these features speed up convergence and can save many hundreds of CPU-hours in the case of difficult runs, with a poor prior knowledge of the covariance matrix. We also summarise all the functionalities of MontePython in the current release, including new likelihoods and plotting options.

Citation format

BRINCKMANN, Thejs; LESGOURGUES, Julien. Montepython 3: Boosted MCMC sampler and other features [preprint]. arXiv, 2018. arXiv:1804.07261.