A generalization of the randomized singular value decomposition

MaRDI QIDQ6368690zbMATH OpenFDO

Authors Nicolas Boullé, Alex Townsend

Publication date 27 May 2021

Abstract: The randomized singular value decomposition (SVD) is a popular and effective algorithm for computing a near-best rank $k$ approximation of a matrix $A$ using matrix-vector products with standard Gaussian vectors. Here, we generalize the randomized SVD to multivariate Gaussian vectors, allowing one to incorporate prior knowledge of $A$ into the algorithm. This enables us to explore the continuous analogue of the randomized SVD for Hilbert--Schmidt (HS) operators using operator-function products with functions drawn from a Gaussian process (GP). We then construct a new covariance kernel for GPs, based on weighted Jacobi polynomials, which allows us to rapidly sample the GP and control the smoothness of the randomly generated functions. Numerical examples on matrices and HS operators demonstrate the applicability of the algorithm.

This page was built for publication: A generalization of the randomized singular value decomposition

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6368690)