We present an open-source CUDA-based package that consists of a compilation of exponential integrators where the action of the matrix exponential or the $\varphi_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 easy integration into any existing software package, and can be used for temporal integration of any differential equation.
翻译:我们介绍一个开源的CUDA软件包,其中集成了多种指数积分器,通过使用Leja点上的多项式插值方法来近似矩阵指数或$\varphi_l$函数在向量上的作用。利用NVIDIA A100 GPU上的若干测试实例,我们展示了CUDA相较于对应CPU代码能够实现显著加速。LeXInt采用模块化设计,便于集成到任何现有软件包中,可用于任意微分方程的时间积分。