Randomized methods for computing joint eigenvalues, with applications to multiparameter eigenvalue problems and root finding
Fuente:
arXiv
Enregistré dans:
| Auteurs principaux: | , , |
|---|---|
| Format: | Preprint |
| Publié: |
2024
|
| Sujets: | |
| Accès en ligne: | |
| Tags: |
Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
|
| _version_ | 1866917826513600512 |
|---|---|
| author | He, Haoze Kressner, Daniel Plestenjak, Bor |
| author_facet | He, Haoze Kressner, Daniel Plestenjak, Bor |
| contents | It is well known that a family of $n\times n$ commuting matrices can be simultaneously triangularized by a unitary similarity transformation. The diagonal entries of the triangular matrices define the $n$ joint eigenvalues of the family. In this work, we consider the task of numerically computing approximations to such joint eigenvalues for a family of (nearly) commuting matrices. This task arises, for example, in solvers for multiparameter eigenvalue problems and systems of multivariate polynomials, which are our main motivations. We propose and analyze a simple approach that computes eigenvalues as one-sided or two-sided Rayleigh quotients from eigenvectors of a random linear combination of the matrices in the family. We provide some analysis and numerous numerical examples, showing that such randomized approaches can compute semisimple joint eigenvalues accurately and lead to improved performance of existing solvers. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2409_00500 |
| institution | arXiv |
| publishDate | 2024 |
| record_format | arxiv |
| spellingShingle | Randomized methods for computing joint eigenvalues, with applications to multiparameter eigenvalue problems and root finding He, Haoze Kressner, Daniel Plestenjak, Bor Numerical Analysis 65F15, 15A27, 68W20, 15A69, 65H04 It is well known that a family of $n\times n$ commuting matrices can be simultaneously triangularized by a unitary similarity transformation. The diagonal entries of the triangular matrices define the $n$ joint eigenvalues of the family. In this work, we consider the task of numerically computing approximations to such joint eigenvalues for a family of (nearly) commuting matrices. This task arises, for example, in solvers for multiparameter eigenvalue problems and systems of multivariate polynomials, which are our main motivations. We propose and analyze a simple approach that computes eigenvalues as one-sided or two-sided Rayleigh quotients from eigenvectors of a random linear combination of the matrices in the family. We provide some analysis and numerous numerical examples, showing that such randomized approaches can compute semisimple joint eigenvalues accurately and lead to improved performance of existing solvers. |
| title | Randomized methods for computing joint eigenvalues, with applications to multiparameter eigenvalue problems and root finding |
| topic | Numerical Analysis 65F15, 15A27, 68W20, 15A69, 65H04 |
| url | https://arxiv.org/abs/2409.00500 |