LeXInt: GPU-accelerated exponential integrators package

被引：1

作者：

Deka, Pranab J. ^{[1
,2
]}

Moriggl, Alexander ^{[1
]}

Einkemmer, Lukas ^{[1
]}

机构：

[1] Univ Innsbruck, Dept Math, A-6020 Innsbruck, Austria

[2] Katholieke Univ Leuven, Leuven, Belgium

来源：

SOFTWAREX | 2025年 / 29卷

基金：

奥地利科学基金会;

关键词：

Accelerator; GPGPU (General Purpose computation on Graphics Processing Unit); NVIDIA; CUDA; Numerical methods; Time integration; Exponential integrators; Polynomial interpolation; STIFF SYSTEMS; EFFICIENT; INTERPOLATION; PERFORMANCE; SOFTWARE; IMPLICIT; SCHEMES;

D O I：

10.1016/j.softx.2024.101949

中图分类号：

TP31 [计算机软件];

学科分类号：

081202 ; 0835 ;

摘要：

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 phi(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.

引用

页数：14

共 38 条

[1] THE PRINCIPLE OF MINIMIZED ITERATIONS IN THE SOLUTION OF THE MATRIX EIGENVALUE PROBLEM [J].

ARNOLDI, WE .

QUARTERLY OF APPLIED MATHEMATICS, 1951, 9 (01) :17-29

[2]

Baglama J., 1998, Electronic Transactions on Numerical Analysis, V7, P124

[3]

Bergamaschi L, 2006, LECT NOTES COMPUT SC, V3994, P685

[4] A new class of time discretization schemes for the solution of nonlinear PDEs [J].

Beylkin, G ;

Keiser, JM ;

Vozovoi, L .

JOURNAL OF COMPUTATIONAL PHYSICS, 1998, 147 (02) :362-387

[5]

Boyd J., 2001, Chebyshev and Fourier Spectral Methods, V2nd

[6] Accurate evaluation of divided differences for polynomial interpolation of exponential propagators [J].

Caliari, M. .

COMPUTING, 2007, 80 (02) :189-201

[7] Interpolating discrete advection-diffusion propagators at Leja sequences [J].

Caliari, M ;

Vianello, M ;

Bergamaschi, L .

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS, 2004, 172 (01) :79-99

[8] The LEM exponential integrator for advection-diffusion-reaction equations [J].

Caliari, Marco ;

Vianello, Marco ;

Bergamaschi, Luca .

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS, 2007, 210 (1-2) :56-63

[9] Comparison of software for computing the action of the matrix exponential [J].

Caliari, Marco ;

Kandolf, Peter ;

Ostermann, Alexander ;

Rainer, Stefan .

BIT NUMERICAL MATHEMATICS, 2014, 54 (01) :113-128

[10] Implementation of exponential Rosenbrock-type integrators [J].

Caliari, Marco ;

Ostermann, Alexander .

APPLIED NUMERICAL MATHEMATICS, 2009, 59 (3-4) :568-581

← 1 2 3 4 →