BAMPHI: Matrix-free and transpose-free action of linear combinations of φ-functions from exponential integrators

The time integration of stiff systems of Ordinary Differential Equations (ODEs), usually arising from the spatial discretization of Partial Differential Equations (PDEs), constitutes a hot topic in numerical analysis. In particular, exponential-type integrators have attracted much attention for thei...

Full description

Saved in:
Bibliographic Details
Published in:Journal of computational and applied mathematics Vol. 423; p. 114973
Main Authors: Caliari, Marco, Cassini, Fabio, Zivcovich, Franco
Format: Journal Article
Language:English
Published: Elsevier B.V 15.05.2023
Subjects:
Backward error analysis
Exponential integrators
Linear combinations of [formula omitted]-functions
Matrix exponential
Stiff ODEs
Backward error analysis
Linear combinations of φ-functions
Stiff ODEs
Exponential integrators
Matrix exponential
ISSN:0377-0427, 1879-1778
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:The time integration of stiff systems of Ordinary Differential Equations (ODEs), usually arising from the spatial discretization of Partial Differential Equations (PDEs), constitutes a hot topic in numerical analysis. In particular, exponential-type integrators have attracted much attention for their capacity to effectively handle the stiffness, allowing integration with large time steps. The efficiency of exponential-type integrators strongly relies on the fast computation of the action of the exponential of matrices which often change only slightly during the integration. The authors exploited this characteristic of exponential integrators to develop a backward stable algorithm, bamphi, which is designed to reuse the information gathered through the exponential integration steps, reaching unmatched levels of speed and accuracy on a variety of numerical experiments.
ISSN:0377-0427
1879-1778
DOI:10.1016/j.cam.2022.114973