{"title":"Akhiezer迭代","authors":"Cade Ballew, Thomas Trogdon","doi":"arxiv-2312.02384","DOIUrl":null,"url":null,"abstract":"We develop the Akhiezer iteration, a generalization of the classical\nChebyshev iteration, for the inner product-free, iterative solution of\nindefinite linear systems using orthogonal polynomials for measures supported\non multiple, disjoint intervals. The iteration applies to shifted linear solves\nand can then be used for efficient matrix function approximation. Using the\nasymptotics of orthogonal polynomials, error bounds are provided. A key\ncomponent in the efficiency of the method is the ability to compute the first\n$k$ orthogonal polynomial recurrence coefficients and the first $k$ weighted\nStieltjes transforms of these orthogonal polynomials in $\\mathrm{O}(k)$\ncomplexity using a numerical Riemann--Hilbert approach. For a special class of\northogonal polynomials, the Akhiezer polynomials, the method can be sped up\nsignificantly, with the greatest speedup occurring in the two interval case\nwhere important formulae of Akhiezer are employed and the Riemann--Hilbert\napproach is bypassed.","PeriodicalId":501061,"journal":{"name":"arXiv - CS - Numerical Analysis","volume":" 8","pages":""},"PeriodicalIF":0.0000,"publicationDate":"2023-12-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"The Akhiezer iteration\",\"authors\":\"Cade Ballew, Thomas Trogdon\",\"doi\":\"arxiv-2312.02384\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"We develop the Akhiezer iteration, a generalization of the classical\\nChebyshev iteration, for the inner product-free, iterative solution of\\nindefinite linear systems using orthogonal polynomials for measures supported\\non multiple, disjoint intervals. The iteration applies to shifted linear solves\\nand can then be used for efficient matrix function approximation. Using the\\nasymptotics of orthogonal polynomials, error bounds are provided. A key\\ncomponent in the efficiency of the method is the ability to compute the first\\n$k$ orthogonal polynomial recurrence coefficients and the first $k$ weighted\\nStieltjes transforms of these orthogonal polynomials in $\\\\mathrm{O}(k)$\\ncomplexity using a numerical Riemann--Hilbert approach. For a special class of\\northogonal polynomials, the Akhiezer polynomials, the method can be sped up\\nsignificantly, with the greatest speedup occurring in the two interval case\\nwhere important formulae of Akhiezer are employed and the Riemann--Hilbert\\napproach is bypassed.\",\"PeriodicalId\":501061,\"journal\":{\"name\":\"arXiv - CS - Numerical Analysis\",\"volume\":\" 8\",\"pages\":\"\"},\"PeriodicalIF\":0.0000,\"publicationDate\":\"2023-12-04\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"arXiv - CS - Numerical Analysis\",\"FirstCategoryId\":\"1085\",\"ListUrlMain\":\"https://doi.org/arxiv-2312.02384\",\"RegionNum\":0,\"RegionCategory\":null,\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"\",\"JCRName\":\"\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - CS - Numerical Analysis","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2312.02384","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
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.