A randomized preconditioned Cholesky-QR algorithm

Fuente: arXiv
Enregistré dans:
Détails bibliographiques
Auteurs principaux: Garrison, James E., Ipsen, Ilse C. F.
Format: Preprint
Publié: 2024
Sujets:
Accès en ligne:
Tags: Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
_version_ 1866911945293037568
author Garrison, James E.
Ipsen, Ilse C. F.
author_facet Garrison, James E.
Ipsen, Ilse C. F.
contents We a present and analyze rpCholesky-QR, a randomized preconditioned Cholesky-QR algorithm for computing the thin QR factorization of real mxn matrices with rank n. rpCholesky-QR has a low orthogonalization error, a residual on the order of machine precision, and does not break down for highly singular matrices. We derive rigorous and interpretable two-norm perturbation bounds for rpCholesky-QR that require a minimum of assumptions. Numerical experiments corroborate the accuracy of rpCholesky-QR for preconditioners sampled from as few as 3n rows, and illustrate that the two-norm deviation from orthonormality increases with only the condition number of the preconditioned matrix, rather than its square -- even if the original matrix is numerically singular.
format Preprint
id arxiv_https___arxiv_org_abs_2406_11751
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle A randomized preconditioned Cholesky-QR algorithm
Garrison, James E.
Ipsen, Ilse C. F.
Numerical Analysis
65F35, 68W20, 60B20, 15A12, 15A18, 15A42, 15B10
We a present and analyze rpCholesky-QR, a randomized preconditioned Cholesky-QR algorithm for computing the thin QR factorization of real mxn matrices with rank n. rpCholesky-QR has a low orthogonalization error, a residual on the order of machine precision, and does not break down for highly singular matrices. We derive rigorous and interpretable two-norm perturbation bounds for rpCholesky-QR that require a minimum of assumptions. Numerical experiments corroborate the accuracy of rpCholesky-QR for preconditioners sampled from as few as 3n rows, and illustrate that the two-norm deviation from orthonormality increases with only the condition number of the preconditioned matrix, rather than its square -- even if the original matrix is numerically singular.
title A randomized preconditioned Cholesky-QR algorithm
topic Numerical Analysis
65F35, 68W20, 60B20, 15A12, 15A18, 15A42, 15B10
url https://arxiv.org/abs/2406.11751