A posteriori computation of the singular vectors in a preconditioned Jacobi SVD algorithm

Authors
Citation
Z. Drmac, A posteriori computation of the singular vectors in a preconditioned Jacobi SVD algorithm, IMA J NUM A, 19(2), 1999, pp. 191-213
Citations number
47
Categorie Soggetti
Mathematics
Journal title
IMA JOURNAL OF NUMERICAL ANALYSIS
ISSN journal
02724979 → ACNP
Volume
19
Issue
2
Year of publication
1999
Pages
191 - 213
Database
ISI
SICI code
0272-4979(199904)19:2<191:APCOTS>2.0.ZU;2-1
Abstract
This paper describes a novel way to implement the Jacobi algorithm for the singular value decomposition of full rank matrices. It is shown that the le ft and right singular vectors can be computed without explicit accumulation of Jacobi rotations. Instead, the accumulated product of Jacobi rotations is computed a posteriori as the solution of a certain well-conditioned matr ix equation. Theoretical analysis provides tools to estimate, check and, if necessary, to improve the accuracy of the computed decomposition. Experime ntal results show that the new technique performs very well in practice.