# Probabilistic Exponential Integrators

arXiv: https://arxiv.org/abs/2305.14978 · HTML: https://mlsanity.com/paper/2305.14978/

## Key facts

- Cited by 4 papers in the map (data through 14 Oct 2026).
- Builds on 20 of its 46 references that are in the map.
- First posted on arXiv 24 May 2023.

## Builds on (most cited first, 20 of 20)

- [Neural Ordinary Differential Equations](https://mlsanity.com/paper/1806.07366/) (2018)
- [Universal Differential Equations for Scientific Machine Learning](https://mlsanity.com/paper/2001.04385/) (2020)
- [Probabilistic Numerics and Uncertainty in Computations](https://mlsanity.com/paper/1506.01326/) (2015)
- [Bayesian Probabilistic Numerical Methods](https://mlsanity.com/paper/1702.03673/) (2017)
- [A probabilistic model for the numerical solution of initial value problems](https://mlsanity.com/paper/1610.05261/) (2016)
- [Probabilistic Solutions To Ordinary Differential Equations As Non-Linear Bayesian Filtering: A New Perspective](https://mlsanity.com/paper/1810.03440/) (2018)
- [Calibrated Adaptive Probabilistic ODE Solvers](https://mlsanity.com/paper/2012.08202/) (2020)
- [A Probabilistic State Space Model for Joint Inference from Differential Equations and Data](https://mlsanity.com/paper/2103.10153/) (2021)
- [Bayesian ODE Solvers: The Maximum A Posteriori Estimate](https://mlsanity.com/paper/2004.00623/) (2020)
- [Stable Implementation of Probabilistic ODE Solvers](https://mlsanity.com/paper/2012.10106/) (2020)
- [A Modern Retrospective on Probabilistic Numerics](https://mlsanity.com/paper/1901.04457/) (2019)
- [Convergence Rates of Gaussian ODE Filters](https://mlsanity.com/paper/1807.09737/) (2018)
- [Differentiable Likelihoods for Fast Inversion of 'Likelihood-Free' Dynamical Systems](https://mlsanity.com/paper/2002.09301/) (2020)
- [Probabilistic ODE Solutions in Millions of Dimensions](https://mlsanity.com/paper/2110.11812/) (2021)
- [Probabilistic Numerical Method of Lines for Time-Dependent Partial Differential Equations](https://mlsanity.com/paper/2110.11847/) (2021)
- [Fenrir: Physics-Enhanced Regression for Initial Value Problems](https://mlsanity.com/paper/2202.01287/) (2022)
- [Active Uncertainty Calibration in Bayesian ODE Solvers](https://mlsanity.com/paper/1605.03364/) (2016)
- [Pick-and-Mix Information Operators for Probabilistic ODE Solvers](https://mlsanity.com/paper/2110.10770/) (2021)
- [Bayesian Filtering for ODEs with Bounded Derivatives](https://mlsanity.com/paper/1709.08471/) (2017)
- [Adaptive high-order splitting schemes for large-scale differential Riccati equations](https://mlsanity.com/paper/1612.00677/) (2016)

## Built on by (most cited first, 4 of 4)

- [Parallel-in-Time Probabilistic Numerical ODE Solvers](https://mlsanity.com/paper/2310.01145/) (2023)
- [Adaptive Probabilistic ODE Solvers Without Adaptive Memory Requirements](https://mlsanity.com/paper/2410.10530/) (2024)
- [Propagating Model Uncertainty through Filtering-based Probabilistic Numerical ODE Solvers](https://mlsanity.com/paper/2503.04684/) (2025)
- [Stable and Scalable Probabilistic Numerical Solvers for Stiff and High-Dimensional ODEs](https://mlsanity.com/paper/2606.08203/) (2026)

## Abstract (arXiv)

Probabilistic solvers provide a flexible and efficient framework for simulation, uncertainty quantification, and inference in dynamical systems. However, like standard solvers, they suffer performance penalties for certain stiff systems, where small steps are required not for reasons of numerical accuracy but for the sake of stability. This issue is greatly alleviated in semi-linear problems by the probabilistic exponential integrators developed in this paper. By including the fast, linear dynamics in the prior, we arrive at a class of probabilistic integrators with favorable properties. Namely, they are proven to be L-stable, and in a certain case reduce to a classic exponential integrator -- with the added benefit of providing a probabilistic account of the numerical error. The method is also generalized to arbitrary non-linear systems by imposing piece-wise semi-linearity on the prior via Jacobians of the vector field at the previous estimates, resulting in probabilistic exponential Rosenbrock methods. We evaluate the proposed methods on multiple stiff differential equations and demonstrate their improved stability and efficiency over established probabilistic solvers. The present contribution thus expands the range of problems that can be effectively tackled within probabilistic numerics.