The Langevin algorithms are frequently used to sample the posterior distributions in Bayesian inference. In many practical problems, however, the posterior distributions often consist of non-differentiable components, posing challenges for the standard Langevin algorithms, as they require to evaluate the gradient of the energy function in each iteration. To this end, a popular remedy is to utilize the proximity operator, and as a result one needs to solve a proximity subproblem in each iteration. The conventional practice is to solve the subproblems accurately, which can be exceedingly expensive, as the subproblem needs to be solved in each iteration. We propose an approximate primal-dual fixed-point algorithm for solving the subproblem, which only seeks an approximate solution of the subproblem and therefore reduces the computational cost considerably. We provide theoretical analysis of the proposed method and also demonstrate its performance with numerical examples.
翻译:Langevin 算法常用于贝叶斯推断中的后验分布采样。然而,在许多实际问题中,后验分布通常由不可微分的组分组成,这给标准 Langevin 算法带来了挑战,因为它们需要在每次迭代中计算能量函数的梯度。针对这个问题,一种常见的解决办法是使用近似算子,因此需要在每次迭代中解决一个近似子问题。传统的做法是准确地解决子问题,但这样做的代价非常昂贵,因为每次迭代都需要解决子问题。我们提出了使用近似原始-对偶不动点算法来解决子问题,该算法仅寻求近似解,因此大大降低了计算成本。我们对所提出的方法进行了理论分析,并在数值示例中展示了其性能。