Abstract:Stiff systems of ordinary differential equations (ODEs) arise in a wide range of scientific and engineering disciplines and are traditionally solved using implicit integration methods due to their stability and efficiency. However, these methods are computationally expensive, particularly for applications requiring repeated integration, such as parameter estimation, Bayesian inference, neural ODEs, physics-informed neural networks, and MeshGraphNets. Explicit exponential integration methods have been proposed as a potential alternative, leveraging the matrix exponential to address stiffness without requiring nonlinear solvers. This study evaluates several state-of-the-art explicit single-step exponential schemes against classical implicit methods on benchmark stiff ODE problems, analyzing their accuracy, stability, and scalability with step size. Despite their initial appeal, our results reveal that explicit exponential methods significantly lag behind implicit schemes in accuracy and scalability for stiff ODEs. The backward Euler method consistently outperformed higher-order exponential methods in accuracy at small step sizes, with none surpassing the accuracy of the first-order integrating factor Euler method. Exponential methods fail to improve upon first-order accuracy, revealing the integrating factor Euler method as the only reliable choice for repeated, inexpensive integration in applications such as neural ODEs and parameter estimation. This study exposes the limitations of explicit exponential methods and calls for the development of improved algorithms.

Comparative Study of Matrix Exponential and Taylor Series Discretization Methods for Nonlinear Odes

Comparison Study of the Taylor Series Based Discretization Method for Nonlinear Input-Delay Systems

Time Discretization of Nonlinear Time-Delay System Using Matrix Exponential Method

Comparison Between First-Order Hold with Zero-Order Hold in Discretization of Input-Delay Nonlinear Systems

Time-Discretization of Delayed Multi-Input Nonlinear System Using A New Algorithm

Discrete the solving model of time-variant standard Sylvester-conjugate matrix equations using Euler-forward formula

High-Order Schemes of Exponential Time Differencing for Stiff Systems with Nondiagonal Linear Part

Numerical Discretization Methods for the Discounted Linear Quadratic Control Problem

Matrix analysis of discrete functionals in compact difference method for nonlinear problems with higher derivatives and program code (I: 1D problem)

A New Algorithm Applied in Control-Relevant Discretization of Nonlinear Systems

Exponential time differencing for matrix-valued dynamical systems

Comparison of Two Hybrid Functions for Numerical Solution of Nonlinear Mixed Partial Integro-Differential Equations

Effective equations for discrete systems: A time stepper based approach

Numerical Discretization Methods for Linear Quadratic Control Problems with Time Delays

Sampled-data Observer Design for a Class of Stochastic Nonlinear Systems Based on the Approximate Discretetime Models.

Numerical Discretization Methods for the Extended Linear Quadratic Control Problem

Performance Evaluation of Single-step Explicit Exponential Integration Methods on Stiff Ordinary Differential Equations

A Novel Fitted Method for a Class of Singularly Perturbed Differential-Difference Equations with Small Delay Exhibiting Twin Layer or Oscillatory Behaviour

Adaptive Time-Step Semi-Implicit One-Step Taylor Scheme for Stiff Ordinary Differential Equations

Black Holes with Weyl Charge and Non-Riemannian Waves

The Use of Adaptive Multiple Models for Predefined Exponential Stability of Control Systems of Uncertain Discrete‐time Nonlinear Objects