More Numerically Accurate Algorithm for Stiff Matrix Exponential

Teddy Lazebnik () and Svetlana Bunimovich-Mendrazitsky
Additional contact information
Teddy Lazebnik: Department of Mathematics, Ariel University, Ariel 4070000, Israel
Svetlana Bunimovich-Mendrazitsky: Department of Mathematics, Ariel University, Ariel 4070000, Israel

Mathematics, 2024, vol. 12, issue 8, 1-13

Abstract: In this paper, we propose a novel, highly accurate numerical algorithm for matrix exponentials (MEs). The algorithm is based on approximating Putzer’s algorithm by analytically solving the ordinary differential equation (ODE)-based coefficients and approximating them. We show that the algorithm outperforms other ME algorithms for stiff matrices for several matrix sizes while keeping the computation and memory consumption asymptotically similar to these algorithms. In addition, we propose a numerical-error- and complexity-optimized decision tree model for efficient ME computation based on machine learning and genetic programming methods. We show that, while there is not one ME algorithm that outperforms the others, one can find a good algorithm for any given matrix according to its properties.

Keywords: decision tree for a numerical algorithm; stiff matrix exponential; Putzer approximation (search for similar items in EconPapers)
JEL-codes: C (search for similar items in EconPapers)
Date: 2024
References: View complete reference list from CitEc
Citations:

Downloads: (external link)
https://www.mdpi.com/2227-7390/12/8/1151/pdf (application/pdf)
https://www.mdpi.com/2227-7390/12/8/1151/ (text/html)

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:gam:jmathe:v:12:y:2024:i:8:p:1151-:d:1373993

Access Statistics for this article

Mathematics is currently edited by Ms. Emma He

More articles in Mathematics from MDPI
Bibliographic data for series maintained by MDPI Indexing Manager ().