Probability Measures for Numerical Solutions of Differential Equations

Abstract

In this paper, we present a formal quantification of epistemic uncertainty induced by numerical solutions of ordinary and partial differential equation models. Numerical solutions of differential equations contain inherent uncertainties due to the finite dimensional approximation of an unknown and implicitly defined function. When statistically analysing models based on differential equations describing physical, or other naturally occurring, phenomena, it is therefore important to explicitly account for the uncertainty introduced by the numerical method. This enables objective determination of its importance relative to other uncertainties, such as those caused by data contaminated with noise or model error induced by missing physical or inadequate descriptors. To this end we show that a wide variety of existing solvers can be randomised, inducing a probability measure over the solutions of such differential equations. These measures exhibit contraction to a Dirac measure around the true unknown solution, where the rates of convergence are consistent with the underlying deterministic numerical method. Ordinary differential equations and elliptic partial differential equations are used to illustrate the approach to quantifying uncertainty in both the statistical analysis of the forward and inverse problems.

Cited by in corpus (9)

• Gamblets for opening the complexity-bottleneck of implicit schemes for hyperbolic and parabolic ODEs/PDEs with rough coefficients
• Frank-Wolfe Bayesian Quadrature: Probabilistic Integration with Theoretical Guarantees
• Active Uncertainty Calibration in Bayesian ODE Solvers
• A hierarchical Bayesian setting for an inverse problem in linear parabolic PDEs with noisy boundary conditions
• Quantitative Kinetic Models from Intravital Microcopy: A Case Study Using Hepatic Transport
• Probabilistic solvers for partial differential equations
• Consistency and convergence of simulation schemes in Information field dynamics
• A probabilistic model for the numerical solution of initial value problems
• Confidence Intervals for Finite Difference Solutions