AIMech AI

Frequent-Directions Deterministic Streaming Matrix Sketching

Algorithms AIMech:FrequentDirectionsDeterministicStreamingMatrixSketching DRAFT
Aliases: Frequent-Directions Matrix Sketching, Liberty Deterministic Covariance Matrix Sketching, Spectral Misra-Gries / Space-Saving Matrix Sketch
Frequent-Directions Deterministic Streaming Matrix Sketching is a deterministic one-pass row-streaming matrix sketching algorithm introduced by Liberty (2013) that maintains an l x m sketch matrix B of an n x m input matrix A by generalizing the Misra-Gries / Space-Saving frequent-items counter decrement rule to singular values and right singular vectors. In the basic algorithm (Algorithm 1), each streaming row a_i in R^m is inserted into an all-zero row of B in R^{l x m}; whenever B has no zero rows, the algorithm computes the SVD B = U Sigma V^T, sets delta_i = sigma_l^2, shrinks all l squared singular values uniformly via bar{Sigma} = sqrt(Sigma^2 - I_l delta_i), and updates B <- bar{Sigma} V^T so that the last row of B is guaranteed to be zero. In the amortized fast variant (Algorithm 2), B is zeroed out in its bottom l/2 rows by shrinking at threshold delta_i = sigma_{l/2}^2 only once every l/2 row insertions, reducing the total runtime from O(n m l^2) to O(n m l) at space O(m l). Claims 1 and 2 and Theorems 4.1 and 4.2 prove that for every unit vector x in R^m (||x||_2 = 1), 0 <= ||A x||_2^2 - ||B x||_2^2 <= 2 (||A||_F^2 - ||B||_F^2) / l <= 2 ||A||_F^2 / l (with factor 1 for Algorithm 1 and factor 2 for Algorithm 2), guaranteeing A^T A >= B^T B and ||A^T A - B^T B||_2 <= epsilon ||A||_F^2 for l = ceil(2 / epsilon).
Identity Notes & Concept Scope Distinct from batch Principal Component Analysis (`AIO:PrincipalComponentAnalysis`), which requires storing A or A^T A in O(n m) or O(m^2) memory, and from randomized oblivious subspace embeddings (`AIMech:SubsampledRandomizedHadamardTransformSubspaceEmbedding`) or randomized column sampling (`AIMech:NystromGramMatrixLowRankApproximationWithDiagonalSampling`). Specifically captures Liberty's (2013) deterministic SVD singular-value shrinkage rule bar{Sigma} = sqrt(Sigma^2 - I_l delta_i) (Algorithms 1 and 2), its exact reduction to the Misra-Gries / Space-Saving heavy-hitters counter algorithm on orthogonal basis streams (Section 3, Theorem 3.1), and the deterministic Loewner positive semidefinite covariance bound 0 <= A^T A - B^T B <= (2 / l)(||A||_F^2 - ||B||_F^2) I_m (Claims 1-2, Theorems 4.1-4.2).

Claims & Grounded Evidence

3 Claims
Mechanism #claim-01
In each step i of Frequent-Directions (Algorithm 1), shrinking the squared singular values of the l x m sketch by delta_i = sigma_l^2 via bar{Sigma} = sqrt(Sigma^2 - I_l delta_i) zeroes out the l-th row and satisfies ||C_i||_F^2 - ||B_i||_F^2 = l delta_i, which telescopes over all n streaming rows to give the exact identity sum_{i=1}^n delta_i = (||A||_F^2 - ||B||_F^2) / l <= ||A||_F^2 / l (Claim 1), while for every unit vector x in R^m, Claim 2 proves 0 <= ||A x||_2^2 - ||B x||_2^2 <= sum_{i=1}^n delta_i, implying both A^T A >= B^T B and ||A^T A - B^T B||_2 <= (||A||_F^2 - ||B||_F^2) / l (Theorem 4.1).
Scope #claim-02
By doubling the sketch capacity and shrinking by delta_i = sigma_{l/2}^2 only when all l/2 zero rows are filled (Algorithm 2), Theorem 4.2 proves that the SVD cost O(m l^2) is amortized over l/2 row insertions to achieve total streaming runtime O(n m l) = O(n m / epsilon) for l = ceil(2 / epsilon) while preserving the deterministic guarantee ||A^T A - B^T B||_2 <= 2 (||A||_F^2 - ||B||_F^2) / l <= epsilon ||A||_F^2; moreover, when input rows are indicator vectors e_{j_i}, Frequent-Directions reduces exactly to the Misra-Gries / Space-Saving item frequency estimation algorithm (Section 3, Theorem 3.1).
Limitation #claim-03
Frequent-Directions does not preserve input sparsity (the sketch B = bar{Sigma} V^T is dense even when the streaming rows a_i are sparse, requiring O(n m l) time rather than input-sparsity time), and its worst-case spectral error bound scales with ||A||_F^2 / l rather than the tail residual ||A - A_k||_F^2 / (l - k) in the original Liberty (2013) analysis when the spectrum of A is completely flat.

⚠️ Curated Limitations & Negative Results

Artificial Intelligence Ontology (AIO) Grounding

Mapped AIO Classes
AIO Target Class Mapping Relation Mapping Rationale
AIO:DimensionalityReduction
Dimensionality Reduction
narrow_match Compresses an n x m streaming data matrix into an l x m deterministic row sketch B approximating A^T A.
AIO:PrincipalComponentAnalysis
Principal Component Analysis
related_match Maintains a streaming SVD sketch whose right singular vectors approximate the dominant principal components of A.

Evidence Sources & Primary Literature

3 Sources
Source ID Publication / Document Type / Access Reference & Checksum
r330-primary-pdf
Edo Liberty
primary paper
full_text
ID: https://api.crossref.org/works/10.1145/2487575.2487623
sha256:11d3293e2dd5…
r330-arxiv-abs
Edo Liberty
primary paper
abstract_only
ID: https://arxiv.org/abs/1206.0594v4
sha256:b5a8ffc45ea2…
r330-metadata-cr
Edo Liberty
primary paper
abstract_only
ID: https://api.crossref.org/works/10.1145/2487575.2487623
sha256:e243a474b5c3…

Ontology & Grounding Topology

Open in interactive network graph →
CONCEPT KIND Algorithms AIMech:FrequentDirectionsDeterministicStreamingMatrixSketching Frequent-Directions... 3 claims · 3 sources AIO GROUNDING AIO:DimensionalityReduction (2)

Curation Provenance

Curator: Antigravity
Method: ai_assisted
Model: Gemini
Created: 2026-10-19T06:52:00Z
Updated: 2026-10-19T06:52:00Z
Source File: data/curated/FrequentDirectionsDeterministicStreamingMatrixSketching.yaml · Schema: curated_concept.yaml

Explore across X-Mech

Related computational & biological Mechs: Embedding Map Network Graph ProteinTraitsMech (ESM-2 pLM) DUFMech (Domain pLM) TraitMech CultureMech All 14 X-Mech Repositories