English

Posterior sampling with Adaptive Gaussian Processes in Bayesian parameter identification

Numerical Analysis 2024-11-28 v1 Numerical Analysis

Abstract

Posterior sampling by Monte Carlo methods provides a more comprehensive solution approach to inverse problems than computing point estimates such as the maximum posterior using optimization methods, at the expense of usually requiring many more evaluations of the forward model. Replacing computationally expensive forward models by fast surrogate models is an attractive option. However, computing the simulated training data for building a sufficiently accurate surrogate model can be computationally expensive in itself, leading to the design of computer experiments problem of finding evaluation points and accuracies such that the highest accuracy is obtained given a fixed computational budget. Here, we consider a fully adaptive greedy approach to this problem. Using Gaussian process regression as surrogate, samples are drawn from the available posterior approximation while designs are incrementally defined by solving a sequence of optimization problems for evaluation accuracy and positions. The selection of training designs is tailored towards representing the posterior to be sampled as good as possible, while the interleaved sampling steps discard old inaccurate samples in favor of new, more accurate ones. Numerical results show a significant reduction of the computational effort compared to just position-adaptive and static designs.

Keywords

Cite

@article{arxiv.2411.17858,
  title  = {Posterior sampling with Adaptive Gaussian Processes in Bayesian parameter identification},
  author = {Paolo Villani and Daniel Andrés-Arcones and Jörg F. Unger and Martin Weiser},
  journal= {arXiv preprint arXiv:2411.17858},
  year   = {2024}
}

Comments

23 pages, 10 figures

R2 v1 2026-06-28T20:13:47.335Z