DeepAI AI Chat
Log In Sign Up

Calibrated Adaptive Probabilistic ODE Solvers

by   Nathanael Bosch, et al.

Probabilistic solvers for ordinary differential equations (ODEs) assign a posterior measure to the solution of an initial value problem. The joint covariance of this distribution provides an estimate of the (global) approximation error. The contraction rate of this error estimate as a function of the solver's step size identifies it as a well-calibrated worst-case error. But its explicit numerical value for a certain step size, which depends on certain parameters of this class of solvers, is not automatically a good estimate of the explicit error. Addressing this issue, we introduce, discuss, and assess several probabilistically motivated ways to calibrate the uncertainty estimate. Numerical experiments demonstrate that these calibration methods interact efficiently with adaptive step-size selection, resulting in descriptive, and efficiently computable posteriors. We demonstrate the efficiency of the methodology by benchmarking against the classic, widely used Dormand-Prince 4/5 Runge-Kutta method.


page 7

page 13

page 14

page 15

page 16


Stability of step size control based on a posteriori error estimates

A posteriori error estimates based on residuals can be used for reliable...

A probabilistic model for the numerical solution of initial value problems

Like many numerical methods, solvers for initial value problems (IVPs) o...

Efficient adaptive step size control for exponential integrators

Traditional step size controllers make the tacit assumption that the cos...

Understanding the impact of numerical solvers on inference for differential equation models

Most ordinary differential equation (ODE) models used to describe biolog...

A Probabilistic Least-Mean-Squares Filter

We introduce a probabilistic approach to the LMS filter. By means of an ...

Extraction of a computer-certified ODE solver

Reliably determining system trajectories is essential in many analysis a...

Stable Implementation of Probabilistic ODE Solvers

Probabilistic solvers for ordinary differential equations (ODEs) provide...