2007.229: Simultaneous tridiagonalization of two symmetric matrices

2007.229: Seamus D. Garvey, Françoise Tisseur, Michael I. Friswell, John E. T. Penny and Uwe Prells (2003) Simultaneous tridiagonalization of two symmetric matrices. International Journal for Numerical Methods in Engineering, 57 (12). pp. 1643-1660. ISSN 0029-5981

Full text available as:

PDF - Archive staff only - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader 168 Kb

DOI: 10.1002/nme.733

Abstract

We show how to simultaneously reduce a pair of symmetric matrices to tridiagonal form by congruence transformations. No assumptions are made on the non-singularity or definiteness of the two matrices. The reduction follows a strategy similar to the one used for the tridiagonalization of a single symmetric matrix via Householder reflectors. Two algorithms are proposed, one using non-orthogonal rank-one modifications of the identity matrix and the other, more costly but more stable, using a combination of Householder reflectors and non-orthogonal rank-one modifications of the identity matrix with minimal condition numbers. Each of these tridiagonalization processes requires O(n3) arithmetic operations and respects the symmetry of the problem. We illustrate and compare the two algorithms with some numerical experiments.

Download Statistics: last 4 weeks
Repository Staff Only: edit this item