Abstract:Rapid convergence of the shifted QR algorithm on symmetric matrices was shown more than fifty years ago. Since then, despite significant interest and its practical relevance, an understanding of the dynamics and convergence properties of the shifted QR algorithm on nonsymmetric matrices has remained elusive. We introduce a new family of shifting strategies for the Hessenberg shifted QR algorithm. We prove that when the input is a diagonalizable Hessenberg matrix $H$ of bounded eigenvector condition number $\kappa_V(H)$ -- defined as the minimum condition number of $V$ over all diagonalizations $VDV^{-1}$ of $H$ -- then the shifted QR algorithm with a certain strategy from our family is guaranteed to converge rapidly to a Hessenberg matrix with a zero subdiagonal entry, in exact arithmetic. Our convergence result is nonasymptotic, showing that the geometric mean of certain subdiagonal entries of $H$ decays by a fixed constant in every $QR$ iteration. The arithmetic cost of implementing each iteration of our strategy scales roughly logarithmically in the eigenvector condition number $\kappa_V(H)$, which is a measure of the nonnormality of $H$. The key ideas in the design and analysis of our strategy are: (1) We are able to precisely characterize when a certain shifting strategy based on Ritz values stagnates. We use this information to design certain ``exceptional shifts'' which are guaranteed to escape stagnation whenever it occurs. (2) We use higher degree shifts (of degree roughly $\log \kappa_V(H)$) to dampen transient effects due to nonnormality, allowing us to treat nonnormal matrices in a manner similar to normal matrices.

A Shifted Method for Solving Symmetric Eigenvalue Problems

Origin-shifted Algorithm for Matrix Eigenvalues

A Quadratically Convergent QR-like Method Without Shifts for the Hermitian Eigenvalue Problem

A Parallel Iterative Method For Solving Symmetrical Tridiagonal Extreme Eigenpair Problems

A Refined Shift-And-Invert Arnoldi Algorithm for Large Unsymmetric Generalized Eigenproblems

Global Convergence of Hessenberg Shifted QR I: Exact Arithmetic

Accelerated Subspace Iteration with Aggressive Shift

An Inverse Extreme Eigenpair Problem and Its Parallel Iterative Solution for Symmetric Tridiagonal Matrices

Shift-splitting Iteration Methods for a Class of Large Sparse Linear Matrix Equations

A parallel iterative method for solving symmetric tridiagonal extreme eigenpair problems

Algorithm 1019: A Task-based Multi-shift QR/QZ Algorithm with Aggressive Early Deflation

A New Shift of the QL Algorithm for Irreducible Symmetric Tridiagonal Matrices

Neural networks based approach for computing eigenvectors and eigenvalues of symmetric matrix

Solving Inverse Eigenvalue Problems Via Householder and Rank-One Matrices

A Symmetric Structure-Preserving ΓQR Algorithm for Linear Response Eigenvalue Problems

A Communication-Avoiding Parallel Algorithm for the Symmetric Eigenvalue Problem

A greedy algorithm for computing eigenvalues of a symmetric matrix

A Shift Selection Strategy for Parallel Shift-Invert Spectrum Slicing in Symmetric Self-Consistent Eigenvalue Computation

An Improved Diagonal Transformation Algorithm for the Maximum Eigenvalue of Zero Symmetric Nonnegative Matrices

Fast parallelizable methods for the Hermitian eigenvalue problem

An Efficient Parallel Krylov-Schur Method for Eigen-Analysis of Large-Scale Power Systems