[{"data":1,"prerenderedAt":-1},["ShallowReactive",2],{"doc-detail-85207-en":3,"doc-seo-85207-105":30,"detail-sidebar-cat-0-en-105":91},{"code":4,"msg":5,"data":6},0,"success",{"doc_id":7,"user_id":8,"nickname":9,"user_avatar":10,"doc_module":4,"category_id":11,"category_name":12,"doc_title":13,"doc_description":14,"doc_content":15,"file_id":16,"file_url":17,"file_type":18,"file_size":19,"view_count":20,"is_deleted":4,"is_public":21,"is_downloadable":21,"audit_status":21,"page_count":22,"language":23,"language_code":24,"site_id":25,"html_lang":24,"table_of_contents":26,"faqs":27,"seo_title":13,"seo_description":14,"update_tm":28,"read_time":29},85207,962075114101,"Seraphina","https://ap-avatar.wpscdn.com/avatar/e000253a75eb197efd?x-image-process=image/resize,m_fixed,w_180,h_180&k=1780044092746381165",8,"Research & Report","Convergence Analysis of a Nonlinear Eigensolver Based on Rational Approximation of the Resolvent","Given a holomorphic matrix-valued function, eigenvalues emerge as the poles of a sketched resolvent. After constructing a high-quality rational approximation to this sketched resolvent, the poles of the rational approximant typically cluster near the true eigenvalues, enabling a flexible route for linear and nonlinear eigenvalue problems. The paper studies convergence of this pole-based eigensolver and introduces block probing and zooming-in to improve eigenvalue accuracy, with proven backward/forward stability and supporting numerical experiments.","arXiv :2607 . 10377v1 [math .NA] 11 Jul 2026  \nCONVERGENCE ANALYSIS OF A NONLINEAR EIGENSOLVER BASED ON RATIONAL APPROXIMATION OF THE RESOLVENT  \nNIAN SHAO˚ AND YUJI NAKATSUKASA:  \nAbstract. Given a holomorphic matrix-valued function, the poles of its sketched resolvent are generically its eigenvalues. Once a good rational approximation of the sketched resolvent is obtained, the poles of this rational approximation typically lie close to those eigenvalues, thus providing a flexible framework for solving both linear and nonlinear eigenvalue problems. However, the accuracy of the computed eigenvalues is limited and remains poorly understood. This paper analyzes the convergence of this approach and demonstrates the effectiveness of two techniques to improve accuracy: block probing and zooming in. We also establish the backward and forward stability of polefinding for a barycentric rational form via a generalized eigenproblem. Numerical experiments demonstrate the sharpness of our theoretical results.  \nKey words. Nonlinear eigenvalue problem, rational approximation, pole/root finding.  \nAMS subject classifications. 65H17, 41A20, 65F15  \n1. Introduction. This work is concerned with nonlinear eigenvalue problems (NEPs) [19, 28] of the form  \n(1.1) Tpλqv “ 0, λ P D and v P Cnzt0u.  \nThroughout this paper, we assume that D Ă C is a nonempty domain enclosed by a Jordan curve, and T is an n ˆ n holomorphic matrix-valued function on some domain that contains the closure of D. We also assume that the NEP (1.1) is regular, that is, that det Tpz0 q ‰ 0 for some z0 P D. In this paper, we consider an approach that computes eigenvalues through the poles of rational approximations of the resolvent T ´1 . To our knowledge, such an approach was first proposed in [1] for symmetric eigenvalue problems, and later developed in [7] for the evaluation of resonances. A general discussion can be found in [33] .  \nLet ΩL and ΩR be two n ˆ b probing matrices, such as Gaussian random matrices. Our approach starts with the following surrogate function for the resolvent:  \nFpzq :“ ΩHLTpzq´1ΩR , z P D.  \nFor generic ΩL and ΩR , the poles of the surrogate function F are exact eigenvalues of T. Let R be a good rational approximation of F in D. Then we expect that the eigenvalues of T are approximated by the poles of R. Once we obtain an approximate eigenvalue λ, the approximate eigenvector can be computed by SVD of Tpλq or residual inverse iteration [34] .  \nThe main ingredient of our approach is computing a rational approximation. The AAA (adaptive Antoulas-Anderson) algorithm [32] provides an efficient tool for rational approximation, and has been employed in a number of applications [33] . The AAA algorithm basically computes a barycentric representation of the rational approximation by a greedy strategy. For a matrix-valued function, the set-valued AAA method [27] computes a matrix-valued rational approximation in a similar approach. Unlike existing contour integral-based approaches [5, 14 , 24 , 35 , 38], which interpret  \n˚ Institute of Mathematics, EPF Lausanne, 1015 Lausanne, Switzerland ([nian.shao@epfl.ch](nian.shao@epfl.ch))  \n: Mathematical Institute, University of Oxford, Oxford, OX2 6GG, UK (nakat[sukasa@maths.ox.ac.uk](sukasa@maths.ox.ac.uk)).  \n2 NIAN SHAO AND YUJI NAKATSUKASA  \nthe quadrature rule as a rational approximation of the indicator function of D, our method directly approximates the resolvent itself and extracts the eigenvalues from its poles. Moreover, our method differs from the nonlinear eigensolvers based on rational approximation proposed in [17, 18 , 20 , 37], which approximate surrogate functions of T instead of its resolvent T ´1 .  \nEven though the (set-valued) AAA can compute a near-best rational approximation effectively, the accuracy of the (nonlinear) eigensolver is not always satisfactory. Consider a toy example where Tpzq “ diagt0 .1 ,..., 0.9u ´ zI9 and ΩL and ΩR are two complex Gaussian random vectors. We uniformly sa","cbCaibHTSa9p6Uup","https://ap.wps.com/l/cbCaibHTSa9p6Uup","pdf",3283197,3,1,24,"English","en",105,"# Introduction\n## Nonlinear eigenvalue problems and resolvent-pole approach\n## Rational approximation via AAA\n## Limitations of naive polefinding\n# Convergence strategy improvements\n## Zooming-in\n## Block probing\n## Stability of polefinding","[{\"question\":\"How are eigenvalues obtained in this method?\",\"answer\":\"Eigenvalues are extracted from the poles of rational approximations of the resolvent sketched from a holomorphic matrix-valued function.\"},{\"question\":\"Why is eigenvalue accuracy limited in the naive approach?\",\"answer\":\"Even with well-separated eigenvalues, numerical examples show that the poles computed from a basic rational approximation can have incorrect digits, especially for interior eigenvalues.\"},{\"question\":\"What techniques does the paper propose to improve accuracy?\",\"answer\":\"The paper demonstrates two strategies: zooming in (recursively subdividing the domain) and block probing (using probing matrices instead of vectors).\"}]",1784201755,60,{"code":4,"msg":31,"data":32},"ok",{"site_id":25,"language":24,"slug":33,"title":13,"keywords":34,"description":14,"schema_data":35,"social_meta":86,"head_meta":88,"extra_data":90,"updated_unix":28},"convergence-analysis-of-a-nonlinear-eigensolver-based-on-rational-approximation-of-the-resolvent","",{"@graph":36,"@context":85},[37,53,68],{"@type":38,"itemListElement":39},"BreadcrumbList",[40,44,48,50],{"item":41,"name":42,"@type":43,"position":21},"https://docshare.wps.com","Home","ListItem",{"item":45,"name":46,"@type":43,"position":47},"https://docshare.wps.com/document/","Document",2,{"item":49,"name":12,"@type":43,"position":20},"https://docshare.wps.com/document/research-report/",{"item":51,"name":13,"@type":43,"position":52},"https://docshare.wps.com/document/convergence-analysis-of-a-nonlinear-eigensolver-based-on-rational-approximation-of-the-resolvent/85207/",4,{"url":51,"name":13,"@type":54,"author":55,"headline":13,"publisher":57,"fileFormat":60,"inLanguage":24,"description":14,"dateModified":61,"datePublished":62,"encodingFormat":60,"isAccessibleForFree":63,"interactionStatistic":64},"DigitalDocument",{"name":9,"@type":56},"Person",{"url":41,"name":58,"@type":59},"DocShare","Organization","application/pdf","2026-07-24","2026-07-16",true,{"@type":65,"interactionType":66,"userInteractionCount":20},"InteractionCounter",{"@type":67},"ViewAction",{"@type":69,"mainEntity":70},"FAQPage",[71,77,81],{"name":72,"@type":73,"acceptedAnswer":74},"How are eigenvalues obtained in this method?","Question",{"text":75,"@type":76},"Eigenvalues are extracted from the poles of rational approximations of the resolvent sketched from a holomorphic matrix-valued function.","Answer",{"name":78,"@type":73,"acceptedAnswer":79},"Why is eigenvalue accuracy limited in the naive approach?",{"text":80,"@type":76},"Even with well-separated eigenvalues, numerical examples show that the poles computed from a basic rational approximation can have incorrect digits, especially for interior eigenvalues.",{"name":82,"@type":73,"acceptedAnswer":83},"What techniques does the paper propose to improve accuracy?",{"text":84,"@type":76},"The paper demonstrates two strategies: zooming in (recursively subdividing the domain) and block probing (using probing matrices instead of vectors).","https://schema.org",{"og:url":51,"og:type":87,"og:title":13,"og:site_name":58,"og:description":14},"article",{"robots":89,"canonical":51},"index,follow",{"doc_id":7,"site_id":25},{"code":4,"msg":5,"data":92},[93,97,101,105,109,114,119,122,127,130,134],{"id":21,"doc_module":4,"doc_module_name":46,"category_name":94,"show_sort_weight":95,"slug":96},"Story & Novel",90,"story-novel",{"id":47,"doc_module":4,"doc_module_name":46,"category_name":98,"show_sort_weight":99,"slug":100},"Literature",80,"literature",{"id":52,"doc_module":4,"doc_module_name":46,"category_name":102,"show_sort_weight":103,"slug":104},"Exam",70,"exam",{"id":106,"doc_module":4,"doc_module_name":46,"category_name":107,"show_sort_weight":29,"slug":108},5,"Comic","comic",{"id":110,"doc_module":4,"doc_module_name":46,"category_name":111,"show_sort_weight":112,"slug":113},6,"Technology",50,"technology",{"id":115,"doc_module":4,"doc_module_name":46,"category_name":116,"show_sort_weight":117,"slug":118},7,"Healthcare",40,"healthcare",{"id":11,"doc_module":4,"doc_module_name":46,"category_name":12,"show_sort_weight":120,"slug":121},30,"research-report",{"id":123,"doc_module":4,"doc_module_name":46,"category_name":124,"show_sort_weight":125,"slug":126},9,"Religion & Spirituality",20,"religion-spirituality",{"id":125,"doc_module":4,"doc_module_name":46,"category_name":128,"show_sort_weight":125,"slug":129},"World Cup","world-cup",{"id":131,"doc_module":4,"doc_module_name":46,"category_name":132,"show_sort_weight":131,"slug":133},10,"Lifestyle","lifestyle",{"id":135,"doc_module":4,"doc_module_name":46,"category_name":136,"show_sort_weight":106,"slug":137},19,"General","general"]