Dynamics of the symmetric eigenvalue problem with shift strategies
Abstract
A common algorithm for the computation of eigenvalues of real symmetric tridiagonal matrices is the iteration of certain special maps Fσ called shifted QR steps. Such maps preserve spectrum and a natural common domain is TΛ, the manifold of real symmetric tridiagonal matrices conjugate to the diagonal matrix Λ. More precisely, a (generic) shift s ∈ defines a map Fs: TΛ TΛ. A strategy σ: TΛ specifies the shift to be applied at T so that Fσ(T) = Fσ(T)(T). Good shift strategies should lead to fast deflation: some off-diagonal coordinate tends to zero, allowing for reducing of the problem to submatrices. For topological reasons, continuous shift strategies do not obtain fast deflation; many standard strategies are indeed discontinuous. Practical implementation only gives rise systematically to bottom deflation, convergence to zero of the lowest off-diagonal entry b(T). For most shift strategies, convergence to zero of b(T) is cubic, |b(Fσ(T))| = Θ(|b(T)|k) for k = 3. The existence of arithmetic progressions in the spectrum of T sometimes implies instead quadratic convergence, k = 2. The complete integrability of the Toda lattice and the dynamics at non-smooth points are central to our discussion. The text does not assume knowledge of numerical linear algebra.
Turn this paper into a lesson
ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.