A scaling and recovering algorithm for the matrix $φ$-functions
Fuente:
arXiv
Gespeichert in:
| Hauptverfasser: | , |
|---|---|
| Format: | Preprint |
| Veröffentlicht: |
2025
|
| Schlagworte: | |
| Online-Zugang: | |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| _version_ | 1866911170739306496 |
|---|---|
| author | Al-Mohy, Awad H. Liu, Xiaobo |
| author_facet | Al-Mohy, Awad H. Liu, Xiaobo |
| contents | A new scaling and recovering algorithm is proposed for simultaneously computing the matrix $φ$-functions that arise in exponential integrator methods for the numerical solution of certain first-order systems of ordinary differential equations. The algorithm initially scales the input matrix down by a nonnegative integer power of two, and then evaluates the $[m/m]$ diagonal Padé approximant to $φ_p$, where $p$ is the largest index of interest. The remaining $[m+p{-}j/m]$ Padé approximants to $φ_j$, $0 \le j < p$, are obtained implicitly via a recurrence relation. The effect of scaling is subsequently recovered using the double-argument formula. A rigorous backward error analysis, based on the $[m+p/m]$ Padé approximant to the exponential, enables sharp bounds on the relative backward errors. These bounds are expressed in terms of the sequence $\|A^k\|^{1/k}$, which can be much smaller than $\|A\|$ for nonnormal matrices. The scaling parameter and the degrees of the Padé approximants are selected to minimize the overall computational cost, which benefits from the sharp bounds and the optimal evaluation schemes for diagonal Padé approximants. Furthermore, if the input matrix is (quasi-)triangular, the algorithm exploits its structure in the recovering phase. Numerical experiments demonstrate the superiority of the proposed algorithm over existing alternatives in both accuracy and efficiency. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2506_01193 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | A scaling and recovering algorithm for the matrix $φ$-functions Al-Mohy, Awad H. Liu, Xiaobo Numerical Analysis 15A16, 65F60, 65L05 A new scaling and recovering algorithm is proposed for simultaneously computing the matrix $φ$-functions that arise in exponential integrator methods for the numerical solution of certain first-order systems of ordinary differential equations. The algorithm initially scales the input matrix down by a nonnegative integer power of two, and then evaluates the $[m/m]$ diagonal Padé approximant to $φ_p$, where $p$ is the largest index of interest. The remaining $[m+p{-}j/m]$ Padé approximants to $φ_j$, $0 \le j < p$, are obtained implicitly via a recurrence relation. The effect of scaling is subsequently recovered using the double-argument formula. A rigorous backward error analysis, based on the $[m+p/m]$ Padé approximant to the exponential, enables sharp bounds on the relative backward errors. These bounds are expressed in terms of the sequence $\|A^k\|^{1/k}$, which can be much smaller than $\|A\|$ for nonnormal matrices. The scaling parameter and the degrees of the Padé approximants are selected to minimize the overall computational cost, which benefits from the sharp bounds and the optimal evaluation schemes for diagonal Padé approximants. Furthermore, if the input matrix is (quasi-)triangular, the algorithm exploits its structure in the recovering phase. Numerical experiments demonstrate the superiority of the proposed algorithm over existing alternatives in both accuracy and efficiency. |
| title | A scaling and recovering algorithm for the matrix $φ$-functions |
| topic | Numerical Analysis 15A16, 65F60, 65L05 |
| url | https://arxiv.org/abs/2506.01193 |