Probabilistic Numerical Methods for PDE-constrained Bayesian Inverse Problems

This paper develops meshless methods for probabilistically describing discretisation error in the numerical solution of partial differential equations. This construction enables the solution of Bayesian inverse problems while accounting for the impact of the discretisation of the forward problem. In particular, this drives statistical inferences to be more conservative in the presence of significant solver error. Theoretical results are presented describing rates of convergence for the posteriors in both the forward and inverse problems. This method is tested on a challenging inverse problem with a nonlinear forward model.

[1]  Erich Novak,et al.  Essays on the Complexity of Continuous Problems , 2009 .

[2]  Michael A. Osborne,et al.  Probabilistic numerics and uncertainty in computations , 2015, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[3]  Qi Ye,et al.  Approximation of stochastic partial differential equations by a kernel-based collocation method , 2011, Int. J. Comput. Math..

[4]  A. Gelman.  Prior distributions for variance parameters in hierarchical models (comment on article by Browne and Draper) , 2004 .

[5]  Houman Owhadi,et al.  Multigrid with Rough Coefficients and Multiresolution Operator Decomposition from Hierarchical Information Games , 2015, SIAM Rev..

[6]  Andrew M. Stuart,et al.  Statistical analysis of differential equations: introducing probability measures on numerical solutions , 2016, Statistics and Computing.

[7]  Houman Owhadi,et al.  Bayesian Numerical Homogenization , 2014, Multiscale Model. Simul..

[8]  Mark A. Girolami,et al.  Probabilistic Meshless Methods for Partial Differential Equations and Bayesian Inverse Problems , 2016, ArXiv.

[9]  H. Owhadi,et al.  Bayesian Numerical Homogenization | Multiscale Modeling & Simulation | Vol. 13, No. 3 | Society for Industrial and Applied Mathematics , 2015 .

[10]  Gregory E. Fasshauer,et al.  Solving differential equations with radial basis functions: multilevel methods and smoothing , 1999, Adv. Comput. Math..

[11]  Simon W. Funke,et al.  Deflation Techniques for Finding Distinct Solutions of Nonlinear Partial Differential Equations , 2014, SIAM J. Sci. Comput..

[12]  Andrew M. Stuart,et al.  Inverse problems: A Bayesian perspective , 2010, Acta Numerica.