Strong convergence rates of probabilistic integrators for ordinary differential equations

Abstract

Probabilistic integration of a continuous dynamical system is a way of systematically introducing model error, at scales no larger than errors introduced by standard numerical discretisation, in order to enable thorough exploration of possible responses of the system to inputs. It is thus a potentially useful approach in a number of applications such as forward uncertainty quantification, inverse problems, and data assimilation. We extend the convergence analysis of probabilistic integrators for deterministic ordinary differential equations, as proposed by Conrad et al.\ (\textit{Stat.\ Comput.}, 2017), to establish mean-square convergence in the uniform norm on discrete- or continuous-time solutions under relaxed regularity assumptions on the driving vector fields and their induced flows. Specifically, we show that randomised high-order integrators for globally Lipschitz flows and randomised Euler integrators for dissipative vector fields with polynomially-bounded local Lipschitz constants all have the same mean-square convergence rate as their deterministic counterparts, provided that the variance of the integration noise is not of higher order than the corresponding deterministic integrator. These and similar results are proven for probabilistic integrators where the random perturbations may be state-dependent, non-Gaussian, or non-centred random variables.

25 pages

References in corpus (5)

• Probabilistic Numerical Methods for PDE-constrained Bayesian Inverse Problems
• Error analysis of randomized Runge-Kutta methods for differential equations with time-irregular coefficients
• Random forward models and log-likelihoods in Bayesian inverse problems
• Implicit Probabilistic Integrators for ODEs
• On the Bayesian Solution of Differential Equations

Cited by in corpus (9)

• A Modern Retrospective on Probabilistic Numerics
• 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
• Estimation of ordinary differential equation models with discretization error quantification
• Randomised one-step time integration methods for deterministic operator differential equations
• Differentiable Likelihoods for Fast Inversion of 'Likelihood-Free' Dynamical Systems
• Calibrated Adaptive Probabilistic ODE Solvers
• Piecewise monotone estimation in one-parameter exponential family
• Sampling Methods for Bayesian Inference Involving Convergent Noisy Approximations of Forward Maps