Efficient Bounds and Estimates for Canonical Angles in Randomized Subspace Approximations

Fuente: arXiv
Guardado en:
Detalles Bibliográficos
Autores principales: Dong, Yijun, Martinsson, Per-Gunnar, Nakatsukasa, Yuji
Formato: Preprint
Publicado: 2022
Materias:
Acceso en línea:
Etiquetas: Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
_version_ 1866913405797924864
author Dong, Yijun
Martinsson, Per-Gunnar
Nakatsukasa, Yuji
author_facet Dong, Yijun
Martinsson, Per-Gunnar
Nakatsukasa, Yuji
contents Randomized subspace approximation with "matrix sketching" is an effective approach for constructing approximate partial singular value decompositions (SVDs) of large matrices. The performance of such techniques has been extensively analyzed, and very precise estimates on the distribution of the residual errors have been derived. However, our understanding of the accuracy of the computed singular vectors (measured in terms of the canonical angles between the spaces spanned by the exact and the computed singular vectors, respectively) remains relatively limited. In this work, we present practical bounds and estimates for canonical angles of randomized subspace approximation that can be computed efficiently either a priori or a posteriori, without assuming prior knowledge of the true singular subspaces. Under moderate oversampling in the randomized SVD, our prior probabilistic bounds are asymptotically tight and can be computed efficiently, while bringing a clear insight into the balance between oversampling and power iterations given a fixed budget on the number of matrix-vector multiplications. The numerical experiments demonstrate the empirical effectiveness of these canonical angle bounds and estimates on different matrices under various algorithmic choices for the randomized SVD.
format Preprint
id arxiv_https___arxiv_org_abs_2211_04676
institution arXiv
publishDate 2022
record_format arxiv
spellingShingle Efficient Bounds and Estimates for Canonical Angles in Randomized Subspace Approximations
Dong, Yijun
Martinsson, Per-Gunnar
Nakatsukasa, Yuji
Numerical Analysis
15A18, 15A42, 65F15, 68W20
Randomized subspace approximation with "matrix sketching" is an effective approach for constructing approximate partial singular value decompositions (SVDs) of large matrices. The performance of such techniques has been extensively analyzed, and very precise estimates on the distribution of the residual errors have been derived. However, our understanding of the accuracy of the computed singular vectors (measured in terms of the canonical angles between the spaces spanned by the exact and the computed singular vectors, respectively) remains relatively limited. In this work, we present practical bounds and estimates for canonical angles of randomized subspace approximation that can be computed efficiently either a priori or a posteriori, without assuming prior knowledge of the true singular subspaces. Under moderate oversampling in the randomized SVD, our prior probabilistic bounds are asymptotically tight and can be computed efficiently, while bringing a clear insight into the balance between oversampling and power iterations given a fixed budget on the number of matrix-vector multiplications. The numerical experiments demonstrate the empirical effectiveness of these canonical angle bounds and estimates on different matrices under various algorithmic choices for the randomized SVD.
title Efficient Bounds and Estimates for Canonical Angles in Randomized Subspace Approximations
topic Numerical Analysis
15A18, 15A42, 65F15, 68W20
url https://arxiv.org/abs/2211.04676