Randomized biorthogonalization through a two-sided Gram-Schmidt process
Abstract
We propose and analyze a randomized two-sided Gram-Schmidt process for the biorthogonalization of two given matrices . The algorithm aims to find two matrices such that , and , where is a sketching matrix satisfying an oblivious subspace -embedding property; in other words, the biorthogonality condition on the columns of and is replaced by an equivalent condition on their sketches. This randomized approach is computationally less expensive than the classical two-sided Gram-Schmidt process, has better numerical stability, and the condition number of the computed bases is often smaller than in the deterministic case. Several different implementations of the randomized algorithm are analyzed and compared numerically. The randomized two-sided Gram-Schmidt process is applied to the nonsymmetric Lancozs algorithm for the approximation of eigenvalues and both left and right eigenvectors.
Cite
@article{arxiv.2509.04386,
title = {Randomized biorthogonalization through a two-sided Gram-Schmidt process},
author = {Laura Grigori and Lorenzo Piccinini and Igor Simunec},
journal= {arXiv preprint arXiv:2509.04386},
year = {2025}
}
Comments
26 pages, 5 figures, 3 tables