Randomized Cholesky QR factorizations

by   Oleg Balabanov, et al.

This article proposes and analyzes several variants of the randomized Cholesky QR factorization of a matrix X. Instead of computing the R factor from X^T X, as is done by standard methods, we obtain it from a small, efficiently computable random sketch of X, thus saving computational cost and improving numerical stability. The proposed direct variant of the randomized Cholesky QR requires only half the flops and the same communication cost as the classical Cholesky QR. At the same time, it is more robust since it is guaranteed to be stable whenever the input matrix is numerically full-rank. The rank-revealing randomized Cholesky QR variant has the ability to sort out the linearly dependent columns of X, which allows to have an unconditional numerical stability and reduce the computational cost when X is rank-deficient. We also depict a column-oriented randomized Cholesky QR that establishes the connection with the randomized Gram-Schmidt process, and a reduced variant that outputs a low-dimensional projection of the Q factor rather than the full factor and therefore yields drastic computational savings. It is shown that performing minor operations in higher precision in the proposed algorithms can allow stability with working unit roundoff independent of the dominant matrix dimension. This feature may be of particular interest for a QR factorization of tall-and-skinny matrices on low-precision architectures.


page 1

page 2

page 3

page 4


Randomized Gram-Schmidt process with application to GMRES

A randomized Gram-Schmidt algorithm is developed for orthonormalization ...

Analysis of Randomized Householder-Cholesky QR Factorization with Multisketching

CholeskyQR2 and shifted CholeskyQR3 are two state-of-the-art algorithms ...

Randomized block Gram-Schmidt process for solution of linear systems and eigenvalue problems

We propose a block version of the randomized Gram-Schmidt process for co...

A Generalized Randomized Rank-Revealing Factorization

We introduce a Generalized Randomized QR-decomposition that may be appli...

Fast and stable randomized low-rank matrix approximation

Randomized SVD has become an extremely successful approach for efficient...

New robust ScaLAPACK routine for computing the QR factorization with column pivoting

In this note we describe two modifications of the ScaLAPACK subroutines ...

Matrix-free construction of HSS representation using adaptive randomized sampling

We present new algorithms for the randomized construction of hierarchica...

Please sign up or login with your details

Forgot password? Click here to reset