ebook img

Implicit Jacobi Algorithm for the SVD with High Relative Accuracy PDF

111 Pages·2011·0.65 MB·English
by  
Most books are stored in the elastic cloud where traffic is expensive. For this reason, we have a limit on daily download.

Preview Implicit Jacobi Algorithm for the SVD with High Relative Accuracy

Implicit Jacobi Algorithm for the SVD with High Relative Accuracy Johan A. Ceballos, Froil´an M. Dopico and Juan Manuel Molera IWASEP 8 Berlin, June 28-July 1 JuanM.Molera(U.CarlosIIIdeMadrid) IWASEP82010-Berlin, June28-July1,2010 1/32 Outline 1 Objectives of this work 2 Jacobi methods for the SVD 3 Implicit Jacobi SVD 4 Error Analysis 5 Numerical Experiments 6 Conclusions JuanM.Molera(U.CarlosIIIdeMadrid) IWASEP82010-Berlin, June28-July1,2010 2/32 Outline 1 Objectives of this work 2 Jacobi methods for the SVD 3 Implicit Jacobi SVD 4 Error Analysis 5 Numerical Experiments 6 Conclusions JuanM.Molera(U.CarlosIIIdeMadrid) IWASEP82010-Berlin, June28-July1,2010 3/32 2 Both start from a Rank Revealing Decomposition (RRD) A=XDYT. 3 Alg. 3.1, QR-decomposition+one application of one-sided Jacobi. The error bound depends of the condition number (usually small) of an intermediate factor: κ(R(cid:48)). 4 Alg. 3.2, Two applications of one-sided Jacobi. It does not have the previous factor. • The implicit Jacobi SVD algorithm gets hra, starting from an RRD, with error bounds that reflect the sensitivity of the problem, with only one run of Jacobi. Objectives of this work • Demmel, Gu, Eisenstat, Slapni˘car, Ves´elic and Drma˘c [1999] presented 1 Two algorithms that compute SVD with high relative accuracy (hra). JuanM.Molera(U.CarlosIIIdeMadrid) IWASEP82010-Berlin, June28-July1,2010 4/32 3 Alg. 3.1, QR-decomposition+one application of one-sided Jacobi. The error bound depends of the condition number (usually small) of an intermediate factor: κ(R(cid:48)). 4 Alg. 3.2, Two applications of one-sided Jacobi. It does not have the previous factor. • The implicit Jacobi SVD algorithm gets hra, starting from an RRD, with error bounds that reflect the sensitivity of the problem, with only one run of Jacobi. Objectives of this work • Demmel, Gu, Eisenstat, Slapni˘car, Ves´elic and Drma˘c [1999] presented 1 Two algorithms that compute SVD with high relative accuracy (hra). 2 Both start from a Rank Revealing Decomposition (RRD) A=XDYT. JuanM.Molera(U.CarlosIIIdeMadrid) IWASEP82010-Berlin, June28-July1,2010 4/32 4 Alg. 3.2, Two applications of one-sided Jacobi. It does not have the previous factor. • The implicit Jacobi SVD algorithm gets hra, starting from an RRD, with error bounds that reflect the sensitivity of the problem, with only one run of Jacobi. Objectives of this work • Demmel, Gu, Eisenstat, Slapni˘car, Ves´elic and Drma˘c [1999] presented 1 Two algorithms that compute SVD with high relative accuracy (hra). 2 Both start from a Rank Revealing Decomposition (RRD) A=XDYT. 3 Alg. 3.1, QR-decomposition+one application of one-sided Jacobi. The error bound depends of the condition number (usually small) of an intermediate factor: κ(R(cid:48)). JuanM.Molera(U.CarlosIIIdeMadrid) IWASEP82010-Berlin, June28-July1,2010 4/32 • The implicit Jacobi SVD algorithm gets hra, starting from an RRD, with error bounds that reflect the sensitivity of the problem, with only one run of Jacobi. Objectives of this work • Demmel, Gu, Eisenstat, Slapni˘car, Ves´elic and Drma˘c [1999] presented 1 Two algorithms that compute SVD with high relative accuracy (hra). 2 Both start from a Rank Revealing Decomposition (RRD) A=XDYT. 3 Alg. 3.1, QR-decomposition+one application of one-sided Jacobi. The error bound depends of the condition number (usually small) of an intermediate factor: κ(R(cid:48)). 4 Alg. 3.2, Two applications of one-sided Jacobi. It does not have the previous factor. JuanM.Molera(U.CarlosIIIdeMadrid) IWASEP82010-Berlin, June28-July1,2010 4/32 Objectives of this work • Demmel, Gu, Eisenstat, Slapni˘car, Ves´elic and Drma˘c [1999] presented 1 Two algorithms that compute SVD with high relative accuracy (hra). 2 Both start from a Rank Revealing Decomposition (RRD) A=XDYT. 3 Alg. 3.1, QR-decomposition+one application of one-sided Jacobi. The error bound depends of the condition number (usually small) of an intermediate factor: κ(R(cid:48)). 4 Alg. 3.2, Two applications of one-sided Jacobi. It does not have the previous factor. • The implicit Jacobi SVD algorithm gets hra, starting from an RRD, with error bounds that reflect the sensitivity of the problem, with only one run of Jacobi. JuanM.Molera(U.CarlosIIIdeMadrid) IWASEP82010-Berlin, June28-July1,2010 4/32 1 SVD [Demmel et al; Ceballos, Dopico, M] 2 Eigenvalues and Eigenvectors [Dopico, M, Moro; Dopico, Koev, M] 3 Solutions of Linear Systems (for almost any rhs) [Dopico, M] • With this algorithm we have now algorithms that useRRD to compute with hra Objectives of this work • This work is a natural extension of the implicit Jacobi method for the symmetric eigenvalue problem [Dopico, Koev, M, 2009] JuanM.Molera(U.CarlosIIIdeMadrid) IWASEP82010-Berlin, June28-July1,2010 5/32 1 SVD [Demmel et al; Ceballos, Dopico, M] 2 Eigenvalues and Eigenvectors [Dopico, M, Moro; Dopico, Koev, M] 3 Solutions of Linear Systems (for almost any rhs) [Dopico, M] Objectives of this work • This work is a natural extension of the implicit Jacobi method for the symmetric eigenvalue problem [Dopico, Koev, M, 2009] • With this algorithm we have now algorithms that useRRD to compute with hra JuanM.Molera(U.CarlosIIIdeMadrid) IWASEP82010-Berlin, June28-July1,2010 5/32

Description:
Demmel, Gu, Eisenstat, Slapni˘car, Vesélic and Drma˘c [1999] presented. 1 Two algorithms that compute SVD with high relative accuracy (hra).
See more

The list of books you might like

book image

Haunting Adeline

H. D. Carlton
·2021
·3.65 MB

book image

$100m Offers

Alex Hormozi
·205 Pages
·2021
·3.18 MB

book image

The Mountain Is You

Brianna Wiest
·2020
·0.34 MB

book image

An Introduction to Fully Integrated Mixed Methods Research

Elizabeth G. Creamer
·404 Pages
·2017
·12.451 MB

book image

Design of Observational Studies

Paul R. Rosenbaum (auth.)
·382 Pages
·2010
·2.027 MB

book image

Greek Government Gazette: Part 7, 2006 no. 892

The Government of the Hellenic Republic
·2006
·0.16 MB

book image

Centenary Magazine

Centenary College of Louisiana
·2006
·345.5 MB

book image

O Evangelho segundo o Espiritismo

Allan Kardec
·562 Pages
·2007
·1.323 MB

book image

The Yellow Frigate by James Grant

240 Pages
·2021
·1.72 MB

book image

Capitalism from Below

Opper, Sonja; Nee, Victor
·2012
·7.264 MB

book image

The Western United States

125 Pages
·2021
·0.68 MB

book image

The Works of Lord Byron Volume 4

339 Pages
·2021
·2.65 MB

book image

Thomism: the philosophy of Thomas Aquinas

Étienne Gilson, Armand Augustine Maurer, Laurence K. Shook, Pontifical Institute of Mediaeval Studies
·468 Pages
·2002
·45.994 MB