Random time step probabilistic methods for uncertainty quantification in chaotic and geometric numerical integration

Abstract

A novel probabilistic numerical method for quantifying the uncertainty induced by the time integration of ordinary differential equations (ODEs) is introduced. Departing from the classical strategy to randomize ODE solvers by adding a random forcing term, we show that a probability measure over the numerical solution of ODEs can be obtained by introducing suitable random time-steps in a classical time integrator. This intrinsic randomization allows for the conservation of geometric properties of the underlying deterministic integrator such as mass conservation, symplecticity or conservation of first integrals. Weak and mean-square convergence analysis are derived. We also analyse the convergence of the Monte Carlo estimator for the proposed random time step method and show that the measure obtained with repeated sampling converges in the mean-square sense independently of the number of samples. Numerical examples including chaotic Hamiltonian systems, chemical reactions and Bayesian inferential problems illustrate the accuracy, robustness and versatility of our probabilistic numerical method.

References in corpus (7)

• The pseudo-marginal approach for efficient Monte Carlo computations
• Bayesian Probabilistic Numerical Methods
• Probabilistic Numerical Methods for PDE-constrained Bayesian Inverse Problems
• A Modern Retrospective on Probabilistic Numerics
• Well-posed Bayesian inverse problems and heavy-tailed stable quasi-Banach space priors
• Strong convergence rates of probabilistic integrators for ordinary differential equations
• Random forward models and log-likelihoods in Bayesian inverse problems

Cited by in corpus (5)

• 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
• Choosing observation operators to mitigate model error in Bayesian inverse problems
• Gradient flows and randomised thresholding: sparse inversion and classification