Subspace Methods For Computing The Pseudospectral Abscissa And The Stability Radius

The pseudospectral abscissa and the stability radius are well-established tools for quantifying the stability of a matrix under unstructured perturbations. Based on first-order eigenvalue expansions, Guglielmi and Overton [SIAM J. Matrix Anal. Appl., 32 (2011), pp. 1166-1192] recently proposed a linearly converging iterative method for computing the pseudospectral abscissa. In this paper, we propose to combine this method and its variants with subspace acceleration. Each extraction step computes the pseudospectral abscissa of a small rectangular matrix pencil, which is comparably cheap and guarantees monotonicity. We observe local quadratic convergence and prove local superlinear convergence of the resulting subspace methods. Moreover, these methods extend naturally to computing the stability radius. A number of numerical experiments demonstrate the robustness and efficiency of the subspace methods.

Published in:
SIAM Journal On Matrix Analysis And Applications, 35, 1, 292-313
Philadelphia, Siam Publications

 Record created 2014-05-02, last modified 2018-03-17

Rate this document:

Rate this document:
(Not yet reviewed)