In a recent paper, Rostami [SIAM J. Sci. Comput, 37(2015), pp. S447-S471] has presented an interesting algorithm for the computation of the real pseudospectral abscissa and the real stability radius (aka the distance to instability) of a square matrix A ∈ ℝn, n in the spectral norm. The method is particularly well suited for large sparse matrices. The algorithm to compute the stability radius is based on the alternation of an iterative method in the rank-2 manifold of real matrices and a Lyapunov inverse iteration. The algorithm shows quadratic convergence in all experiments, but this peculiar property is only conjectured. In this paper, we provide a rigorous proof and propose a modification which replaces the Lyapunov inverse iteration and speeds up the method. Moreover, we interpret the algorithm to compute the real ϵ-pseudospectral abscissa as a nonstandard discretization of the ODE-based method given in Guglielmi and Lubich [SIAM J. Matrix Anal. Appl., 34(2013), pp. 40-66]. Several numerical illustrations comparing the modified algorithm with that proposed by Rostami conclude the paper.

On the method by Rostami for computing the real stability radius of large and sparse matrices

GUGLIELMI, NICOLA
2016-01-01

Abstract

In a recent paper, Rostami [SIAM J. Sci. Comput, 37(2015), pp. S447-S471] has presented an interesting algorithm for the computation of the real pseudospectral abscissa and the real stability radius (aka the distance to instability) of a square matrix A ∈ ℝn, n in the spectral norm. The method is particularly well suited for large sparse matrices. The algorithm to compute the stability radius is based on the alternation of an iterative method in the rank-2 manifold of real matrices and a Lyapunov inverse iteration. The algorithm shows quadratic convergence in all experiments, but this peculiar property is only conjectured. In this paper, we provide a rigorous proof and propose a modification which replaces the Lyapunov inverse iteration and speeds up the method. Moreover, we interpret the algorithm to compute the real ϵ-pseudospectral abscissa as a nonstandard discretization of the ODE-based method given in Guglielmi and Lubich [SIAM J. Matrix Anal. Appl., 34(2013), pp. 40-66]. Several numerical illustrations comparing the modified algorithm with that proposed by Rostami conclude the paper.
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11697/108392
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 3
  • ???jsp.display-item.citation.isi??? 3
social impact