Fast algorithms for the computation of the pseudospectral abscissa and pseudospectral radius.
- Series
- Applied and Computational Mathematics Seminar
- Time
- Monday, November 9, 2009 - 13:00 for 1 hour (actually 50 minutes)
- Location
- Skiles 255
- Speaker
- Nicola Guglielmi – Università di L'Aquila – guglielm@univaq.it
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.