Abstract:We propose and analyze deterministic multilevel approximations for Bayesian inversion of operator equations with uncertain distributed parameters, subject to additive Gaussian measurement data. The algorithms use a multilevel (ML) approach based on deterministic, higher order quasi-Monte Carlo (HoQMC) quadrature for approximating the high-dimensional expectations, which arise in the Bayesian estimators, and a Petrov-Galerkin (PG) method for approximating the solution to the underlying partial differential equation (PDE). This extends the previous single-level approach from [J. Dick, R. N. Gantner, Q. T. Le Gia and Ch. Schwab, Higher order Quasi-Monte Carlo integration for Bayesian Estimation. Report 2016-13, Seminar for Applied Mathematics, ETH Zürich (in review)]. We obtain sufficient conditions which allow us to achieve arbitrarily high, algebraic convergence rates in terms of work, which are independent of the dimension of the parameter space. The convergence rates are limited only by the spatial regularity of the forward problem,the discretization order achieved by the Petrov Galerkin discretization, and by the sparsity of the uncertainty parametrization. We provide detailed numerical experiments for linear elliptic problems in two space dimensions, with $s=1024$ parameters characterizing the uncertain input, confirming the theory and showing that the ML HoQMC algorithms outperform, in terms of error vs.~computational work, both multilevel Monte Carlo (MLMC) methods and single-level (SL) HoQMC methods.

A Low-Rank Solver for Parameter Estimation and Uncertainty Quantification in Time-Dependent Systems of Partial Differential Equations

Low rank approximation method for perturbed linear systems with applications to elliptic type stochastic PDEs

A Low-Rank Approach for Nonlinear Parameter-Dependent Fluid-Structure Interaction Problems

Bayesian Parameter Inference for Partially Observed Diffusions using Multilevel Stochastic Runge-Kutta Methods

Dual parametric and state estimation for partial differential equations

A Low-rank solver for the Navier--Stokes equations with uncertain viscosity

Numerical Method for Parameter Inference of Systems of Nonlinear Ordinary Differential Equations with Partial Observations.

Efficient solution of ordinary differential equations with high-dimensional parametrized uncertainty

Parameter estimation in ODEs: assessing the potential of local and global solvers

CUR for Implicit Time Integration of Random Partial Differential Equations on Low-Rank Matrix Manifolds

Parameter Estimation for Partial Differential Equations by Collage-Based Numerical Approximation

A space-time adaptive low-rank method for high-dimensional parabolic partial differential equations

Multilevel higher order Quasi-Monte Carlo Bayesian Estimation

Parameter Estimation in Nonlinear Multivariate Stochastic Differential Equations Based on Splitting Schemes

Total Uncertainty Quantification in Inverse PDE Solutions Obtained with Reduced-Order Deep Learning Surrogate Models

High-order empirical interpolation methods for real time solution of parametrized nonlinear PDEs

Solving optimal control problems governed by random Navier-Stokes equations using low-rank methods

Parameter estimation for SPDEs based on discrete observations in time and space

Stochastic methods for solving high-dimensional partial differential equations

Joint Parameter and Parameterization Inference with Uncertainty Quantification through Differentiable Programming

Time-dependent stochastic basis adaptation for uncertainty quantification