---
title: 'RandRAND: Randomized Range-Deflation Preconditioners'
url: https://www.emergentmind.com/topics/randrand
type: topic
---

# RandRAND: Randomized Range-Deflation Preconditioners

Searching arXiv for papers on “RandRAND” and closely related work.
RandRAND denotes a class of randomized deflation preconditioners for large linear systems. Introduced in “Preconditioning via Randomized Range Deflation (RandRAND),” it deflates the spectrum via efficient orthogonal projections onto random subspaces, without computing eigenpairs or low-rank approximations. The method targets systems of the form \((A+\mu I)x=b\), uses randomized range sketches of the shifted operator, and yields projector-based preconditioners whose condition number bounds depend only weakly on the problem size and reduce to a small constant when the dimension of the deflated subspace is comparable to the effective spectral dimension [2509.19747].

## 1. Problem setting and randomized range deflation

The basic object is the shifted matrix
\[
A_\mu := A+\mu I,\qquad A\in\mathbb{R}^{n\times n},\ \mu\in\mathbb{R},
\]
with the goal of accelerating Krylov solvers for
\[
(A+\mu I)x=b.
\]
The regime of interest is one in which \(A\) or \(A_\mu\) has a small number of large singular or eigenvalues compared to the rest. In that setting, classical deflation would build a basis of dominant eigenvectors and project against it. RandRAND replaces this explicit spectral computation by randomized range approximation [2509.19747].

The construction starts from a sketch matrix \(\Omega\in\mathbb{R}^{n\times l}\), typically generated from an oblivious subspace embedding. For symmetric \(A\),
\[
\Omega := A^q X^\top,
\]
and for non-symmetric \(A\),
\[
\Omega := (A^\top A)^q X^\top,
\]
where \(q\in\{0,1,\dots\}\) is a small power and \(X\in\mathbb{R}^{l\times n}\) is an OSE. The randomized deflation space is then
\[
V := A_\mu \Omega.
\]
If the dominant eigendirections or singular directions are well represented in \(\mathrm{range}(V)\), projecting onto that range suppresses the large-scale spectral components that drive poor conditioning.

This formulation is projection-based rather than low-rank-approximation-based. A plausible implication is that the method is naturally suited to settings where matvecs with \(A\) are available but explicit eigenspaces, factorizations, or matrix entries are not.

## 2. Orthogonal projectors and the three RandRAND preconditioners

Let \(\Pi=\Pi_V\) be the orthogonal projector onto \(\mathrm{range}(V)\), so that
\[
\Pi^2=\Pi,\qquad \Pi^\top=\Pi,\qquad \mathrm{range}(\Pi)=\mathrm{range}(V).
\]
With an explicit QR factorization
\[
A_\mu\Omega = QR,\quad Q^\top Q=I_l,
\]
the projector is
\[
\Pi u = QQ^\top u.
\]
A central identity is
\[
A_\mu^{-1}\Pi = (A_\mu\Omega)^\dagger,
\]
and, when \(A_\mu\Omega=QR\),
\[
(A_\mu\Omega)^\dagger u = R^{-\top}Q^\top u.
\]
These identities let RandRAND act on the deflated range without forming \(A_\mu^{-1}\) explicitly [2509.19747].

Three preconditioner families are introduced.

The **R-RandRAND** preconditioner for SPD \(A_\mu\) is
\[
P := A_\mu^{-1}\left((I-\Pi)A_\mu(I-\Pi)+\tau\Pi\right),
\]
with preconditioned operator
\[
B := A_\mu P = (I-\Pi)A_\mu(I-\Pi)+\tau\Pi.
\]
On \(\mathrm{range}(\Pi)\), one has \(A_\mu Pv=\tau v\). Thus the large eigenvalues captured by \(\Pi\) are flattened to \(\tau\), while the orthogonal complement is governed by the projected operator \((I-\Pi)A_\mu(I-\Pi)\).

The **C-RandRAND** preconditioner for SPD \(A_\mu\) is
\[
P := (I-\Pi)+\tau A_\mu^{-1}\Pi = (I-\Pi)+\tau (A_\mu\Omega)^\dagger.
\]
This is a correction-type construction: it leaves the complement essentially untouched and corrects the action of \(A_\mu\) on the projected range.

The **G-RandRAND** preconditioner for indefinite or non-symmetric systems is
\[
P := (I-\Pi)+\tau\left((A_\mu^{-1}\Pi)^\top A_\mu^{-1}\Pi\right)^{1/2}.
\]
This form is SPD and is used as a left preconditioner for symmetric indefinite systems with MINRES and for general non-symmetric systems in left-preconditioned GMRES/BiCG-type methods.

The distinction among the three forms is structural. R-RandRAND modifies the operator directly on the projected range; C-RandRAND introduces an additive correction; G-RandRAND transfers the same idea to singular-value control and indefinite or non-normal regimes.

## 3. Spectral effect, effective spectral dimension, and condition-number theory

The phrase “deflating the spectrum” is used in a specific sense. RandRAND does not zero out eigenvalues. Instead, it flattens the action of \(A_\mu\) on \(\mathrm{range}(\Pi)\) to a controlled scale \(\tau\), while the complement is governed by the projected operator. For R-RandRAND,
\[
B=(I-\Pi)A_\mu(I-\Pi)+\tau\Pi,
\]
so the large eigenvalues represented in \(\mathrm{range}(\Pi)\) are replaced by \(\tau\) [2509.19747].

Let
\[
E := (I-\Pi)A_\mu(I-\Pi).
\]
For SPD \(A_\mu\), if
\[
\rho^{-1}\lambda_{\min}(A_\mu)\le \tau\le \rho\|E\|
\]
for some \(\rho\ge 1\), then
\[
\mathrm{cond}(A_\mu P)\le \rho\, \frac{\|E\|}{\lambda_{\min}(A_\mu)}.
\]
This shows that the preconditioned condition number is controlled by the norm of the projected residual operator \(E\), rather than by the full spectral spread of \(A_\mu\).

The randomized theory then links \(\|E\|\) to the quality of the sketch. For Gaussian \(X\in\mathbb{R}^{l\times n}\), \(l=2ck+4\), \(1<c\le 2\), \(k\ge 1\), and \(\Omega=A^qX^\top\), with probability at least \(1-6^{-k}\),
\[
\mathrm{cond}(A_\mu P)\le 1+\left(C_3\frac{n}{(c-1)k}\right)^{1/2}\mathrm{cond}_k(A_\mu^2)^{1/2}
\]
for \(q=0\), where \(C_3=2\cdot 18^2\), and \(4e^2\) in expectation. For \(q\ge 1\), the paper gives a corresponding bound involving \((A_\mu^2A^{2q})^{1/(2q+2)}\) and \(\mathrm{cond}_{k,\mu}\). The dependence on \(n\) appears only through a mild factor of the form \((n/k)^{1/(2q+2)}\).

The analysis is organized around the **effective spectral dimension**. Informally, this is the number of eigenvalues above the regularization scale \(\mu\), or above a comparable threshold relative to the smallest shifted eigenvalue. If the sketch dimension is chosen comparable to that dimension, then the projected operator norm \(\|(I-\Pi)A_\mu\|\) becomes close to the trailing spectral scale, and the condition number can reduce to a small constant.

A further convergence statement is given for R-RandRAND. If
\[
T := \left(\frac{\|(I-\Pi)A_\mu(I-\Pi)\|}{\lambda_n+\mu}\right)^{1/2},
\]
then CG or MINRES on the preconditioned operator satisfies an iteration estimate of the form
\[
t \lesssim \frac{1}{2}T(\log 2+\log(1/u)).
\]
For sketch size \(l\ge 2.4d+4\), the paper states
\[
\mathbb{E}(T)\le 1+\left(148\frac{n}{d}\right)^{1/(4q+4)},
\]
and, with probability at least \(1-6^{-d}\),
\[
T \le 1+\left(3240\frac{n}{d}\right)^{1/(4q+4)}.
\]
This suggests that the iteration count depends only weakly on \(n/d\) rather than on \(\mathrm{cond}(A_\mu)\).

## 4. Basis-explicit and basis-less algorithms

An explicit implementation computes \(A_\mu\Omega=QR\), stores \(Q\), and applies
\[
\Pi u = QQ^\top u,\qquad (A_\mu\Omega)^\dagger u = R^{-\top}Q^\top u.
\]
For symmetric \(A\), the paper gives construction cost
\[
O(T_{\rm mv}\cdot l(q+1)+nl^2+l^3)\ \text{flops},\qquad O(nl)\ \text{memory},
\]
and application cost \(O(T_{\rm mv}+nl)\) per preconditioned operator call. This is viable when \(l\) is small and explicit orthogonalization is affordable [2509.19747].

A major contribution is the **basis-less** implementation. Instead of storing \(Q\), RandRAND retains only Q-less information and uses
\[
Qu = A_\mu\Omega R^{-1}u,\qquad Q^\top u = R^{-\top}\Omega^\top A_\mu^\top u.
\]
This makes it possible to exploit fast randomized transforms and bypass explicit basis orthogonalization.

The first basis-less route is **fast Q-less Cholesky QR**. It forms
\[
G=\Omega^\top A_\mu^\top A_\mu\Omega,
\]
then sets \(R=\mathrm{chol}(G)\). With structured sketches such as SRHT or sparse embeddings, the Gram matrix can be built at cost
\[
O(2(q+1)T_{\rm mv}l + nl\log n).
\]
However, the paper notes that for ill-conditioned \(A_\mu\Omega\), the orthogonality error of Cholesky QR scales like \(u\,\mathrm{cond}(A_\mu\Omega)^2\).

The stabilization strategy is **preconditioned Cholesky QR**. A second-level sketch \(\Omega_{\rm sk}\) is used to form
\[
S=\Omega_{\rm sk}^\top A_\mu\Omega,\qquad S=\tilde Q R_{\rm sk},
\]
followed by
\[
G = R_{\rm sk}^{-T}\Omega^\top A_\mu^\top A_\mu\Omega R_{\rm sk}^{-1},\qquad
R_{\rm chol}=\mathrm{chol}(G),\qquad
R=R_{\rm sk}R_{\rm chol}.
\]
If \(\Omega_{\rm sk}\) is an \(\varepsilon\)-embedding with \(\varepsilon<1\), then
\[
\mathrm{cond}(A_\mu\Omega R^{-1}) \le \sqrt{\frac{1+\varepsilon}{1-\varepsilon}},
\]
so the orthogonality error becomes \(O(u\,\mathrm{cond}(A_\mu))\).

A third route is **factorization-free iterative projection**. Instead of factoring \(\Omega^\top A_\mu^2\Omega\), the projection step is recast as a least-squares problem:
\[
Q^\dagger u = \arg\min_v \|A_\mu\Omega R^{-1}v-u\|.
\]
With \(R=I\), this becomes
\[
Q^\dagger u = \arg\min_v \|A_\mu\Omega v-u\|,
\]
and can be solved by iterative sketch-based least-squares methods, such as Newton-Sketch-type refinement. This route is particularly attractive in GPU or streaming settings.

## 5. Numerical stability, related methods, and empirical behavior

The numerical-stability analysis centers on projector accuracy, loss of orthogonality, and preservation of symmetry or SPD structure. The paper proposes several remedies: preconditioned Cholesky QR, repeated application of \((I-\Pi)\) in the projected operator, and refined projector approximations such as
\[
\Pi u \approx Q(2I-Q^\top Q)Q^\top u,
\]
which reduce projector error from first order to second order in \(\|I-Q^\top Q\|\) [2509.19747].

RandRAND is closely related to two existing lines of work. Relative to **classical deflation**, it serves the same spectral role but does not compute eigenvectors. Relative to **Nyström-based randomized preconditioning**, it uses the same randomized information in a different way. The paper rewrites the Nyström pseudo-inverse as
\[
\hat A_\mu^\dagger = A_\mu^{-1}\Pi,
\]
and observes that
\[
P = (I-\Pi)+\tau A_\mu^{-1}\Pi
\]
may be viewed as “Nyström pseudo-inverse plus a projector.” The theoretical distinction is that Nyström analysis is driven by \(\|A-\hat A\|\), whereas RandRAND is driven by projection errors such as \(\|(I-\Pi)A_\mu\|\). The paper states that, for heavy-tailed spectra, Nyström’s bound may be worse by a factor
\[
O\Bigl(\bigl(n/d\bigr)^{(2q+2)/(2q+1)}\Bigr).
\]

The empirical study covers synthetic SPD and indefinite matrices, ridge regression with random Fourier features, kernel ridge regression, portfolio optimization, PDE-constrained inverse problems, and limited-memory kernel ridge regression. The reported pattern is consistent across these settings. For synthetic SPD systems, RandRAND reduces condition numbers by 2–4 orders of magnitude for moderate \(l\). For portfolio optimization, RandRAND cuts iterations by approximately \(3\)–\(5\times\), and R-RandRAND gives approximately \(15\)–\(20\%\) fewer iterations than Nyström and C-RandRAND. In basis-less PDE inverse and limited-memory kernel settings, basis-less RandRAND remains robust, while basis-less Nyström is reported as numerically unstable or prone to breakdown for large \(l\). This suggests that the projection-based design is especially advantageous when memory is limited or when explicit orthogonalization is infeasible.

## 6. Terminology and other uses of the label

In the supplied arXiv literature, “RandRAND” is not uniformly used for a single concept. The formal named method is the randomized range-deflation preconditioner described above [2509.19747]. The label is also used informally or hypothetically in several distinct contexts.

In “Batched Ranged Random Integer Generation,” the paper itself states that it does not name the algorithm “RandRAND,” but the label is used there to denote efficient unbiased conversion of 64-bit random words into multiple independent bounded integers by full-width multiplication and a single rejection criterion [2408.06213]. In “The RAM equivalent of P vs. RP,” the core primitive is \(\mathrm{RAND}(y)\), which returns a uniformly distributed random integer in the range \([0,y)\); the supplied discussion explicitly treats “RandRAND” as shorthand for the equivalence between \(\mathrm{RAND}(y)\) and \(\mathrm{RAND}(2^k)\) in suitably equipped RAM models [1306.0400]. Two random-number-generation papers use the term only as a hypothetical design label: one for a chaos-based PRNG built from the Rössler attractor and a Maiorana-type bent Boolean function [1702.07502], and one for an enhanced-NRBG architecture that combines a quantum entropy source with SHA-256 Hash-DRBG and Toeplitz hashing [2108.04331].

This suggests that, in current arXiv usage, “RandRAND” is a polysemous label. In the strict bibliographic sense, it refers to randomized range deflation for preconditioning large linear systems. In adjacent discussions, it also serves as an informal name for ranged random integer generation, RAND-based RAM stochasticity, or hypothetical RNG architectures.

Source: https://www.emergentmind.com/topics/randrand