Abstract:Randomized algorithms for low-rank matrix approximation are investigated, with the emphasis on the fixed-precision problem and computational efficiency for handling large matrices. The algorithms are based on the so-called QB factorization, where Q is an orthonormal matrix. First, a mechanism for calculating the approximation error in the Frobenius norm is proposed, which enables efficient adaptive rank determination for a large and/or sparse matrix. It can be combined with any QB-form factorization algorithm in which B's rows are incrementally generated. Based on the blocked randQB algorithm by Martinsson and Voronin, this results in an algorithm called randQB_EI. Then, we further revise the algorithm to obtain a pass-efficient algorithm, randQB_FP, which is mathematically equivalent to the existing randQB algorithms and also suitable for the fixed-precision problem. Especially, randQB_FP can serve as a single-pass algorithm for calculating leading singular values, under a certain condition. With large and/or sparse test matrices, we have empirically validated the merits of the proposed techniques, which exhibit remarkable speedup and memory saving over the blocked randQB algorithm. We have also demonstrated that the single-pass algorithm derived by randQB_FP is much more accurate than an existing single-pass algorithm. And with data from a scenic image and an information retrieval application, we have shown the advantages of the proposed algorithms over the adaptive range finder algorithm for solving the fixed-precision problem.

High relative accuracy for rational q-Bernstein–Vandermonde matrices

Total Positivity and Accurate Computations Related to q-Abel Polynomials

Accurate computations with Said-Ball-Vandermonde matrices

Total positivity and accurate computations related to $q$-Abel polynomials

A fixed-point algorithm for matrix projections with applications in quantum information

Fast parallelizable methods for the Hermitian eigenvalue problem

Exact QR factorizations of rectangular matrices

Accurate Computations with Block Checkerboard Pattern Matrices

A Parallelizable Method for Computing the Singular Subspaces of a General Matrix

Efficient Randomized Algorithms for the Fixed-Precision Low-Rank Matrix Approximation

Diagonal transformations for the computation of the Perron root of nonnegative irreducible matrices

Towards a more robust algorithm for computing the restricted singular value decomposition

The Generation of All Regular Rational Orthogonal Matrices

Total positivity and high relative accuracy for several classes of Hankel matrices

Real dqds for the nonsymmetric tridiagonal eigenvalue problem

Exact Eigenvalues and Eigenvectors for Some n-Dimensional Matrices

Accurate computations of singular values and linear systems for Polynomial-Vandermonde-type matrices

A Cubically Convergent Parallelizable Method for the Hermitian Eigenvalue Problem

The Power of Bidiagonal Matrices

Qubit-Efficient Randomized Quantum Algorithms for Linear Algebra

A randomized algorithm for the QR decomposition-based approximate SVD