Computation of the exponential function of matrices by a formula without oscillatory integrals on infinite intervals
Fuente:
arXiv
Salvato in:
| Autori principali: | , |
|---|---|
| Natura: | Preprint |
| Pubblicazione: |
2024
|
| Soggetti: | |
| Accesso online: | |
| Tags: |
Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
|
| _version_ | 1866916498315935744 |
|---|---|
| author | Suzuki, Masato Tanaka, Ken'ichiro |
| author_facet | Suzuki, Masato Tanaka, Ken'ichiro |
| contents | We propose a quadrature-based formula for computing the exponential function of matrices with a non-oscillatory integral on an infinite interval and an oscillatory integral on a finite interval. In the literature, existing quadrature-based formulas are based on the inverse Laplace transform or the Fourier transform. We show these expressions are essentially equivalent in terms of complex integrals and choose the former as a starting point to reduce computational cost. By choosing a simple integral path, we derive an integral expression mentioned above. Then, we can easily apply the double-exponential formula and the Gauss-Legendre formula, which have rigorous error bounds. As numerical experiments show, the proposed formula outperforms the existing formulas when the imaginary parts of the eigenvalues of matrices have large absolute values. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2411_19086 |
| institution | arXiv |
| publishDate | 2024 |
| record_format | arxiv |
| spellingShingle | Computation of the exponential function of matrices by a formula without oscillatory integrals on infinite intervals Suzuki, Masato Tanaka, Ken'ichiro Numerical Analysis 65F60, 65D30 We propose a quadrature-based formula for computing the exponential function of matrices with a non-oscillatory integral on an infinite interval and an oscillatory integral on a finite interval. In the literature, existing quadrature-based formulas are based on the inverse Laplace transform or the Fourier transform. We show these expressions are essentially equivalent in terms of complex integrals and choose the former as a starting point to reduce computational cost. By choosing a simple integral path, we derive an integral expression mentioned above. Then, we can easily apply the double-exponential formula and the Gauss-Legendre formula, which have rigorous error bounds. As numerical experiments show, the proposed formula outperforms the existing formulas when the imaginary parts of the eigenvalues of matrices have large absolute values. |
| title | Computation of the exponential function of matrices by a formula without oscillatory integrals on infinite intervals |
| topic | Numerical Analysis 65F60, 65D30 |
| url | https://arxiv.org/abs/2411.19086 |