Fast algorithms for the computation of the pseudospectral abscissa and pseudospectral radius.

Series
Applied and Computational Mathematics Seminar
Time
Monday, November 9, 2009 - 1:00pm for 1 hour (actually 50 minutes)
Location
Skiles 255
Speaker
Nicola Guglielmi – Università di L'Aquila – guglielm@univaq.ithttp://univaq.it/~guglielm/
Organizer
Sung Ha Kang
This is a joint work with Michael Overton (Courant Institute, NYU). The epsilon-pseudospectral abscissa and radius of an n x n matrix are respectively the maximum real part and the maximal modulus of points in its epsilon-pseudospectrum. Existing techniques compute these quantities accurately but the cost is multiple SVDs of order n, which makesthe method suitable to middle size problems. We present a novel approach based on computing only the spectral abscissa or radius or a sequence of matrices, generating a monotonic sequence of lower bounds which, in many but not all cases, converges to the pseudospectral abscissa or radius.