A block bidiagonalization method for fixed-accuracy low-rank matrix approximation

From MaRDI portal



Abstract: We present randUBV, a randomized algorithm for matrix sketching based on the block Lanzcos bidiagonalization process. Given a matrix , it produces a low-rank approximation of the form , where and have orthonormal columns in exact arithmetic and is block bidiagonal. In finite precision, the columns of both and will be close to orthonormal. Our algorithm is closely related to the randQB algorithms of Yu, Gu, and Li (2018) in that the entries of are incrementally generated and the Frobenius norm approximation error may be efficiently estimated. Our algorithm is therefore suitable for the fixed-accuracy problem, and so is designed to terminate as soon as a user input error tolerance is reached. Numerical experiments suggest that the block Lanczos method is generally competitive with or superior to algorithms that use power iteration, even when has significant clusters of singular values.



Cites work



Describes a project that uses

Uses Software






This page was built for publication: A block bidiagonalization method for fixed-accuracy low-rank matrix approximation

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