Probabilistic Numerical Methods for PDE-constrained Bayesian Inverse Problems

Abstract

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.

Cited by in corpus (16)

• Machine Learning of Linear Differential Equations using Gaussian Processes
• A Modern Retrospective on Probabilistic Numerics
• Strong convergence rates of probabilistic integrators for ordinary differential equations
• Random time step probabilistic methods for uncertainty quantification in chaotic and geometric numerical integration
• A probabilistic finite element method based on random meshes: Error estimators and Bayesian inverse problems
• Sparse Gaussian processes for solving nonlinear PDEs
• Parameter Estimation with Maximal Updated Densities
• Gaussian processes for Bayesian inverse problems associated with linear partial differential equations
• Randomised one-step time integration methods for deterministic operator differential equations
• Differentiable Likelihoods for Fast Inversion of 'Likelihood-Free' Dynamical Systems
• GPU-Accelerated Particle Methods for Evaluation of Sparse Observations for PDE-Constrained Inverse Problems
• Radiative Transfer as a Bayesian Linear Regression problem
• Images of Gaussian and other stochastic processes under closed, densely-defined, unbounded linear operators
• Optimality Criteria for Probabilistic Numerical Methods
• Rejoinder for "Probabilistic Integration: A Role in Statistical Computation?"
• Sampling Methods for Bayesian Inference Involving Convergent Noisy Approximations of Forward Maps