LeXInt: GPU-accelerated exponential integrators package

We present an open-source CUDA -based package, for the temporal integration of differential equations, that consists of a compilation of exponential integrators where the action of the matrix exponential or the φ l functions on a vector is approximated using the method of polynomial interpolation at Leja points. Using a couple of test examples on an NVIDIA A100 GPU, we show that one can achieve significant speedups using CUDA over the corresponding CPU code. LeXInt , written in a modular format, facilitates integration into existing software packages (written in C++ or CUDA ), for temporal integration of differential equations.