Computation of eigenvalues and eigenvectors for the discrete ordinate and matrix operator methods in radiative transfer

Knut Stamnes, Si Chee Tsay, Teruyuki Nakajima

Research output: Contribution to journalArticlepeer-review

23 Scopus citations

Abstract

Nakajima and Tanaka showed that the algebraic eigenvalue problem occurring in the discrete ordinate and matrix operator methods can be reduced to finding eigenvalues and eigenvectors of the product of two symmetric matrices, one of which is positive definite. Here, we show that the Cholesky decomposition of this positive definite matrix can be used to convert the eigenvalue problem into one involving a symmetric matrix. The Cholesky decomposition is extremely stable and is expected to improve the speed of the eigenvalue/eigenvector computation. After a careful comparison of the Nakajima and Tanaka procedure, our new Cholesky decomposition method and the original procedure suggested by Stamnes and Swanson, we find (contrary to our expectations) that the Stamnes and Swanson prescription is still the most accurate because it avoids round-off errors due to matrix multiplications needed to symmetrize the matrix in the two other procedures. We also find that, when the QR algorithm (used to solve the asymmetric eigenvalue problem in the Stamnes and Swanson procedure) is changed to avoid complex arithmetic, the speed becomes comparable to that of the two other procedures based on reduction to symmetric matrices.

Original languageEnglish
Pages (from-to)415-419
Number of pages5
JournalJournal of Quantitative Spectroscopy and Radiative Transfer
Volume39
Issue number5
DOIs
StatePublished - May 1988

Fingerprint

Dive into the research topics of 'Computation of eigenvalues and eigenvectors for the discrete ordinate and matrix operator methods in radiative transfer'. Together they form a unique fingerprint.

Cite this