Approximation of the Pseudospectral Abscissa via Eigenvalue Perturbation Theory

Abstract

Reliable and efficient computation of the pseudospectral abscissa in the large-scale setting is still not settled. Unlike the small-scale setting where there are globally convergent criss-cross algorithms, all algorithms in the large-scale setting proposed to date are at best locally convergent. We first describe how eigenvalue perturbation theory can be put in use to estimate the globally rightmost point in the ε-pseudospectrum if ε is small. Our treatment addresses both general nonlinear eigenvalue problems, and the standard eigenvalue problem as a special case. For small ε, the estimates by eigenvalue perturbation theory are quite accurate. In the standard eigenvalue case, we even derive a formula with an O(ε3) error. For larger ε, the estimates can be used to initialize the locally convergent algorithms. We also propose fixed-point iterations built on the the perturbation theory ideas for large ε that are suitable for the large-scale setting. The proposed fixed-point iterations initialized by using eigenvalue perturbation theory converge to the globally rightmost point in the pseudospectrum in a vast majority of the cases that we experiment with.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…