Applied and Computational Mathematics Seminar: Fast algorithms for the computation of the pseudospectral abscissa and pseudospectral radius.
Mon, 11/09/2009 - 1:00pm, Skiles 255
Nicola Guglielmi, Università di L'Aquila , email 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.