HMC: Reducing the number of rejections by not using leapfrog and some results on the acceptance rate
HMC: Reducing the number of rejections by not using leapfrog and some results on the acceptance rate
复制标题
DOI:
10.1016/j.jcp.2021.110333
复制
发表时间:
2021-04-13
影响因子:
4.1
通讯作者:
Sanz-Serna, J. M.
中科院分区:
文献类型:
--
作者:
Calvo, M. P.;Sanz-Alonso, D.;Sanz-Serna, J. M.
The leapfrog integrator is routinely used within the Hamiltonian Monte Carlo method and its variants. We give strong numerical evidence that alternative, easy to implement algorithms yield fewer rejections with a given computational effort. When the dimensionality of the target distribution is high, the number of accepted proposals may be multiplied by a factor of three or more. This increase in the number of accepted proposals is not achieved by impairing any positive features of the sampling. We also establish new non-asymptotic and asymptotic results on the monotonic relationship between the expected acceptance rate and the expected energy error. These results further validate the derivation of one of the integrators we consider and are of independent interest. (C) 2021 Elsevier Inc. All rights reserved.