Computing accurate eigensystems of scaled diagonally dominant matrices

Jesse Barlow, James Demmel

Research output: Contribution to journalArticle

98 Scopus citations

Abstract

When eigenvalues of symmetric matrices and singular values of general matrices are computed in finite-precision arithmetic, it is generally expected that they will be computed with an error bound proportional to the product of machine precision and the norm of the matrix. In particular, tiny eigenvalues and singular values are usually not computed to high relative accuracy. There are some important classes of matrices that provide a higher level of precision, including bidiagonal matrices, scaled diagonally dominant matrices, and scaled diagonally dominant definite pencils. These classes include many graded matrices and all symmetric positive-definite matrices that can be consistently ordered (i.e., all symmetric positive-definite tridiagonal matrices). In particular, singular values and eigenvalues are determined to high relative precision, independent of their magnitudes, and there are algorithms to compute them accurately. The eigenvectors are also determined more accurately than for general matrices, and may be computed more accurately as well. This work extends results of Kahan and Demmel on bidiagonal and tridiagonal matrices.

Original languageEnglish (US)
Pages (from-to)762-791
Number of pages30
JournalSIAM Journal on Numerical Analysis
Volume27
Issue number3
DOIs
StatePublished - Jan 1 1990

All Science Journal Classification (ASJC) codes

  • Numerical Analysis
  • Computational Mathematics
  • Applied Mathematics

Fingerprint Dive into the research topics of 'Computing accurate eigensystems of scaled diagonally dominant matrices'. Together they form a unique fingerprint.

  • Cite this