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 X, Y ∈Rn× m. The algorithm aims to find two matrices Q, P ∈Rn× m such that range(X) = range(Q), range(Y) = range(P) and ( Q)T P = I, where ∈Rs × n is a sketching matrix satisfying an oblivious subspace -embedding property; in other words, the biorthogonality condition on the columns of Q and P 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 Q, P 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.
Turn this paper into a full lesson
ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.