The Akhiezer iteration
Fuente:
arXiv
Saved in:
| Main Authors: | , |
|---|---|
| Format: | Preprint |
| Published: |
2023
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866909075496763392 |
|---|---|
| author | Ballew, Cade Trogdon, Thomas |
| author_facet | Ballew, Cade Trogdon, Thomas |
| contents | We develop the Akhiezer iteration, a generalization of the classical Chebyshev iteration, for the inner product-free, iterative solution of indefinite linear systems using orthogonal polynomials for measures supported on multiple, disjoint intervals. The iteration applies to shifted linear solves and can then be used for efficient matrix function approximation. Using the asymptotics of orthogonal polynomials, error bounds are provided. A key component in the efficiency of the method is the ability to compute the first $k$ orthogonal polynomial recurrence coefficients and the first $k$ weighted Stieltjes transforms of these orthogonal polynomials in $\mathrm{O}(k)$ complexity using a numerical Riemann--Hilbert approach. For a special class of orthogonal polynomials, the Akhiezer polynomials, the method can be sped up significantly, with the greatest speedup occurring in the two interval case where important formulae of Akhiezer are employed and the Riemann--Hilbert approach is bypassed. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2312_02384 |
| institution | arXiv |
| publishDate | 2023 |
| record_format | arxiv |
| spellingShingle | The Akhiezer iteration Ballew, Cade Trogdon, Thomas Numerical Analysis Complex Variables 42C05, 65E05, 33C47, 65F10 We develop the Akhiezer iteration, a generalization of the classical Chebyshev iteration, for the inner product-free, iterative solution of indefinite linear systems using orthogonal polynomials for measures supported on multiple, disjoint intervals. The iteration applies to shifted linear solves and can then be used for efficient matrix function approximation. Using the asymptotics of orthogonal polynomials, error bounds are provided. A key component in the efficiency of the method is the ability to compute the first $k$ orthogonal polynomial recurrence coefficients and the first $k$ weighted Stieltjes transforms of these orthogonal polynomials in $\mathrm{O}(k)$ complexity using a numerical Riemann--Hilbert approach. For a special class of orthogonal polynomials, the Akhiezer polynomials, the method can be sped up significantly, with the greatest speedup occurring in the two interval case where important formulae of Akhiezer are employed and the Riemann--Hilbert approach is bypassed. |
| title | The Akhiezer iteration |
| topic | Numerical Analysis Complex Variables 42C05, 65E05, 33C47, 65F10 |
| url | https://arxiv.org/abs/2312.02384 |