---
title: Pareto-Guided SVD Compression for LLMs and VLMs
url: https://www.emergentmind.com/topics/pareto-guided-singular-value-decomposition-pgsvd
type: topic
---

# Pareto-Guided SVD Compression for LLMs and VLMs

Searching arXiv for the specified paper and closely related work on PGSVD and activation-aware low-rank compression.
Pareto-Guided Singular Value Decomposition (PGSVD) is a zero-shot, activation-informed low-rank compression framework for large language models (LLMs) and vision–language models (VLMs). It combines a loss-sensitivity upper bound that ties layer-wise activation-based compression to network-level loss, a bi-objective formulation of compression as a trade-off between model size and loss change, and a Pareto-guided rule that selects heterogeneous ranks through a single uniform tolerance. In the formulation introduced in "Activation-Informed Pareto-Guided Low-Rank Compression for Efficient LLM/VLM" [2510.05544], the resulting pipeline uses sampled activations to guide factor optimization and an alternating least-squares (ALS) solver to improve stability and accuracy over truncation-based SVD.

## 1. Definition and problem setting

PGSVD addresses the deployment problem posed by foundation models that deliver strong performance but impose severe memory and compute costs during inference. The framework targets zero-shot compression—without retraining—in order to reduce both memory and FLOPs while preserving model quality. Its basic mechanism is low-rank factorization, in which a weight matrix \(W\) is replaced by factors \(A B\), with \(A \in \mathbb{R}^{N\times r}\) and \(B \in \mathbb{R}^{r\times M}\). This reduces parameters from \(N M\) to \(r(N+M)\) and lowers per-forward FLOPs from \(O(NM)\) to \(O(r(N+M))\) [2510.05544].

The framework is motivated by a limitation of direct weight-space compression. Many pre-trained layers do not exhibit low-rank weights, so direct truncation of \(W\) can degrade accuracy substantially. By contrast, activations \(X\) often exhibit stronger low-rank structure. PGSVD therefore treats activation-informed distortion, rather than weight-space distortion alone, as the relevant compression target. The conventional truncated SVD objective is
\[
\min_{A\in \mathbb{R}^{N\times r},\; B\in \mathbb{R}^{r\times M}}
\;\bigl\| W - A B \bigr\|_{F}^{2},
\]
whereas the activation-aware objective is
\[
\min_{A\in \mathbb{R}^{N\times r},\; B\in \mathbb{R}^{r\times M}}
\;\bigl\| W X - A B X \bigr\|_{F}^{2}.
\]

The distinction is consequential. The activation-informed distortion measures how the compressed layer perturbs the layer outputs under observed inputs. This suggests that PGSVD is not merely a matrix approximation procedure, but a deployment-oriented compression scheme in which the optimization target is aligned with downstream model behavior.

## 2. Loss sensitivity and surrogate compression theory

The theoretical core of PGSVD is a loss-sensitivity upper bound linking activation-based layer compression errors to the network-level change in loss. Consider a feedforward layer \(l\) with pre-activations \(z_l = W_l x_l\), an elementwise nonlinearity \(\sigma\), and batch activations
\[
X_l = \bigl[x_l^{(1)}~\cdots~x_l^{(B)}\bigr].
\]
Let the compressed weight be \(\hat{W}_l = W_l + \Delta W_l\), and assume \(\sigma\) has bounded slope \(\sup_{t\in\mathbb R}|\sigma'(t)|\le c<\infty\). For a differentiable scalar loss \(\mathcal{L}\), the stated theorem gives
\[
\bigl|\Delta\mathcal{L}\bigr|
~\le~
G \sum_{l=1}^{L}
\Biggl(\prod_{m=l+1}^{L} \mathcal{K}_m\Biggr)
c\,\|\Delta W_l X_l\|_F,
\]
where \(G := \bigl\|\nabla_{Y}\mathcal{L}\bigr\|_{F}\), \(Y=X_{L+1}\), and
\[
\mathcal{K}_l:=\sup_{1\le i\le B}\|J_l^{(i)}\|,\qquad
J_l^{(i)}
=\operatorname{diag}\!\bigl(\sigma'(W_l x_l^{(i)})\bigr)\,W_l.
\]

The proof sketch is explicitly first-order: a first-order perturbation analysis bounds the propagated activation perturbation by layer-dependent Jacobian norms and the compression-induced distortion \(\|\Delta W_l X_l\|_F\), and a first-order Taylor bound converts activation perturbation into loss change via \(\|\nabla_Y \mathcal{L}\|_F\). The theoretical implication stated for the method is that, for fixed data and model, the coefficients are constants, so reducing \(\sum_l \|\Delta W_l X_l\|_F\) provably reduces a surrogate upper bound on the loss change.

PGSVD then formulates compression as a bi-objective optimization. For a rank-\(r_l\) approximation \(\hat{W}_l^{(r_l)}\), define the relative weight-space error
\[
e_l(r_l)=\frac{\bigl\|\hat{W}_l^{(r_l)}-W_l\bigr\|_F}{\|W_l\|_F},\qquad
\Gamma_l:=\{r_l\in\mathbb{Z}\mid e_l(r_l)\le \varepsilon_l\},
\]
where \(\varepsilon_l\in[0,1]\) is a tolerance. Let \(r=[r_1,\ldots,r_L]\) collect ranks and let
\[
S(r)=\sum_{l=1}^L P_l(r_l)
\]
denote the network parameter count. The bi-objective problem is
\[
\min_{r \in \prod_{l=1}^L \Gamma_l}
\;\;\bigl( S(r), \;\; |\Delta {\cal L}(r)| \bigr).
\tag{B}
\]

Using \(\|\Delta W_l X_l\|_F\le \|\Delta W_l\|_F\,\|X_l\|_F\) and the loss-sensitivity theorem yields the scalar surrogate
\[
|\Delta \mathcal{L}(r)| \le \sum_{l=1}^{L}\alpha_l\, e_l(r_l),
\]
with
\[
\alpha_l=\|\nabla_{Y} \mathcal{L}\|_F
\Bigg( \prod_{m = l+1}^L \mathcal{K}_m \Bigg)c\, \|X_l\|_F\,\|W_l\|_F.
\]
Under a compression budget \(b\), the surrogate rank-allocation program becomes
\[
\min_{r\,\in\prod_{l=1}^{L}\Gamma_l} \;\; |\Delta \mathcal{L}(r)| \le \sum_{l=1}^{L}\alpha_l e_l(r_l)
\quad\text{s.t.}\quad \sum_{l=1}^{L}P_{l}(r_l)\le b.
\tag{P}
\]

This formulation places PGSVD within a class of compression methods that optimize a tractable surrogate rather than the exact post-compression loss. A common misconception is that the framework’s theory directly optimizes the true loss; the stated result is instead an upper-bound argument, and the paper explicitly notes that the loss bound is first-order and conservative.

## 3. Pareto-guided rank allocation via a uniform tolerance

A distinctive aspect of PGSVD is its conversion of rank allocation into tolerance allocation. For any \(W\in\mathbb{R}^{N\times M}\) and \(\varepsilon \in [0,1]\), the framework defines a unique minimal rank \(r^\star(\varepsilon)\) such that
\[
\frac{\|W-\hat{W}^{(r^\star)}\|_F}{\|W\|_F} \;\le\; \varepsilon,
\]
with parameter count \(P(r)=r(M+N)\). The associated map is
\[
h:[0,1]\to\mathbb{Z}_{\ge 0},\qquad h(\varepsilon):=P\bigl(r^\star(\varepsilon)\bigr).
\]
The paper states that the rank-allocation problem and the \(\varepsilon\)-allocation problem
\[
\min_{\,0\le \varepsilon_1,\dots,\varepsilon_L\le 1}
\ \sum_{l=1}^L \alpha_l\,\varepsilon_l \quad \text{s.t.}\quad \sum_{l=1}^L h_l(\varepsilon_l)\le b
\tag{E}
\]
have the same optimal value [2510.05544].

The central Pareto-guided claim is then made under a robust surrogate in which \(\alpha_l\equiv\alpha\) for all \(l\), together with bounded, nonincreasing convex envelopes for each layer’s parameter function \(h_l\). Under those assumptions, there exists a uniform solution \(\varepsilon_1=\cdots=\varepsilon_L=\varepsilon^\star\), and every uniform tolerance \(\varepsilon\in[0,1]\) corresponds to a point on the surrogate Pareto frontier of the bi-objective problem. The per-layer ranks are nevertheless heterogeneous because each layer’s spectrum determines
\[
r_l(\varepsilon) = \min\{r\in\mathbb{Z}\mid e_l(r)\le\varepsilon\}.
\]

This directly addresses another common misunderstanding: a single tolerance does not imply a single rank. The tolerance is uniform, but the ranks are heterogeneous because layers have different SVD profiles \(e_l(r)\). In the PGSVD formulation, the “single-knob” control variable is the tolerance \(\varepsilon\), while the rank pattern is induced by layer-specific spectra. For VLMs, the same idea is extended to separate tolerances \(\varepsilon_v\) and \(\varepsilon_t\) for the vision and text towers, respectively.

The budget-to-tolerance mapping is also explicit. Given a compression budget \(b\), one solves for \(\varepsilon\) using the upper envelope \(\overline h\) via
\[
L\,\overline h(\varepsilon) = b.
\]
The layer-wise ranks then follow by the minimal-rank condition \(e_l(r_l)\le\varepsilon\). This suggests that PGSVD replaces ad hoc layer-wise rank search with a structured Pareto-guided rule, although the theoretical guarantee is stated only for the homogeneous-sensitivity surrogate and bounded-profile assumptions.

## 4. Algorithmic pipeline and ALS factorization

PGSVD is presented as a zero-shot pipeline with four stages: activation sampling, tolerance selection, rank allocation, and activation-weighted ALS factorization. The procedure requires sampled activations but no retraining. The step-by-step pipeline is as follows.

First, representative inputs are run through the model to collect per-layer activations \(X_l\), stacked by columns. The empirical activation covariance is then formed as
\[
M_l=X_l X_l^\top.
\]

Second, a single uniform tolerance \(\varepsilon\) is chosen, or per-modality tolerances \(\varepsilon_v, \varepsilon_t\) in VLMs. If a parameter or FLOP budget \(b\) is given, the mapping uses the convex envelope relation \(L\,\overline h(\varepsilon)=b\).

Third, each layer receives the minimal rank
\[
r_l = \min\{r\in \mathbb{Z} \mid e_l(r) \le \varepsilon\},
\]
where \(e_l(r) = \| \hat{W}_l^{(r)} - W_l\|_F / \|W_l\|_F\) is the weight-space relative truncation error.

Fourth, low-rank factors are initialized by SVD on \(W_l\). With
\[
W_l \approx U_r \Sigma_r V_r^\top,
\]
the initialization is
\[
A_l=U_r\Sigma_r^{1/2},\qquad B_l=\Sigma_r^{1/2}V_r^\top.
\]
These factors are then refined by ALS to minimize
\[
\min_{A,B}\ \bigl\|W_l X_l - A B X_l\bigr\|_F^2.
\]

The ALS updates use the empirical activation covariance \(M=X X^\top\):
\[
A = W M B^{\top} \left( B M B^{\top}\right)^{\dagger},\qquad
B = \left(A^{\top}A \right)^{\dagger}A^{\top}W.
\]
Only \(r\times r\) pseudo-inverses are required. The distinction from plain truncated SVD is explicit: truncated SVD optimizes \(\|W-AB\|_F\), whereas ALS optimizes \(\|W X - A B X\|_F\), thereby weighting errors by the empirical activation covariance \(M\). The paper further states that this avoids the numerical fragility, including Cholesky failures, of whitening-based methods and is more efficient than EVD-based implementations.

The algorithmic outline takes as input weights \(\{W_l\}_{l=1}^{L}\), activation covariances \(\{M_l\}_{l=1}^L\), tolerance \(\varepsilon\) or \(\varepsilon_v/\varepsilon_t\), and ALS iterations \(\tau\). For each layer, it computes \(r_l\), initializes \(A_l,B_l\) by rank-\(r_l\) SVD, performs \(t=1\) to \(\tau\) ALS updates, and returns the factorized model \(\{A_l, B_l\}_{l=1}^L\). The paper emphasizes that activation sampling affects factor optimization but not the rank allocation when using uniform \(\varepsilon\) [2510.05544].

## 5. Scope, computational characteristics, and implementation practice

The reported compression scope covers all linear layers in self-attention—query, key, value, and output projections—and all linear projections in MLP blocks. In CLIP-like VLMs, both the image encoder, such as ViT, and the text encoder transformers are compressed, typically with separate tolerances per modality.

The parameter and FLOP accounting are straightforward. For layer \(l\), parameters drop from \(N_l M_l\) to \(r_l(N_l+M_l)\), so the total parameter count is
\[
S(r)=\sum_{l} r_l(N_l+M_l).
\]
For inference on a single vector \(x\), FLOPs reduce from \(O(N_l M_l)\) to \(O(r_l M_l + r_l N_l)\) through
\[
y \approx A_l(B_l x).
\]

The main hyperparameter is the uniform tolerance \(\varepsilon\), or \(\varepsilon_v,\varepsilon_t\) in multimodal settings. It may be chosen via a target budget \(b\) using \(L\,\overline h(\varepsilon)=b\), or by sweeping \(\varepsilon\) to trace the surrogate Pareto front. Activation sampling requires only a small batch of representative data per modality to estimate \(M_l=X_l X_l^\top\). The reported implementation guidance is to use full precision for inference during activation collection and double precision during factorization for numerical stability. The recommended ALS setting is \(\tau \approx 10\), with empirical improvements plateauing between 5–10 iterations; even a single iteration yields acceptable perplexity. Stopping may be based on fixed \(\tau\) or on stabilization of \(\|W_l X_l - A_l B_l X_l\|_F\). No regularization is required in the reported experiments.

These implementation choices place PGSVD between purely analytic compression and post-training adaptation. It is zero-shot and deployment-oriented, yet still activation-informed. A plausible implication is that its operational cost is concentrated in one pass for activation collection and a relatively small number of ALS refinements, rather than in any fine-tuning loop.

## 6. Empirical behavior, comparative position, and limitations

The reported evaluation covers LLMs—LLaMA-2-7B, LLaMA-2-13B, and Mistral-7B—and VLMs based on CLIP [2510.05544]. For LLMs, metrics include perplexity on WikiText-2 and zero-shot accuracy on ARC-E, CSQA, LAMBADA, PIQA, and WinoGrande, with additional tasks ARC-C, BoolQ, HellaSwag, MathQA, MMLU, and RACE. On LLaMA-2-7B at 20% compression, PGSVD reports PPL \(= 7.38\) versus SVD-LLM \(7.70\) and SVD-ALS \(7.72\), with average accuracy \(56.08\%\) versus \(53.93\%\) and \(54.42\%\). At 40% compression, it reports PPL \(= 13.46\) versus \(14.95\) and \(15.03\), with average accuracy \(43.66\%\) versus \(41.09\%\) and \(41.35\%\). On LLaMA-2-13B, the reported values are PPL \(= 5.96\) and average accuracy \(60.38\%\) at 20%, and PPL \(= 9.55\) with average accuracy \(47.57\%\) at 40%, outperforming both SVD-LLM and SVD-ALS. On Mistral-7B, the reported values are PPL \(= 6.71\) and average accuracy \(57.14\%\) at 20%, and PPL \(= 14.43\) with average accuracy \(40.98\%\) at 40%.

Throughput is reported to be similar to uniform-rank SVD-ALS under naïve Python implementations, with substantial gains over base models across sequence lengths. The ALS ablation states that compression times are lower than EVD-based SVD variants and that perplexity improvements plateau after approximately 5–10 ALS iterations.

For VLMs, the evaluation uses zero-shot Top-1 and Top-5 over Caltech101, Food101, OxfordPets, StanfordCars, EuroSAT, and DTD. At 20% compression, PGSVD reports average Top-1/Top-5 \(= 62.60\%/83.95\%\), compared with SVD-ALS \(= 45.53\%/74.05\%\) and naive SVD \(\approx 3.09\%/16.22\%\). At 40% compression, the reported values are \(50.89\%/76.18\%\), compared with \(37.08\%/65.94\%\) for SVD-ALS and \(\approx 3.28\%/14.55\%\) for SVD. The paper further reports that modality-aware tolerances \(\varepsilon_v\) and \(\varepsilon_t\) yield better trade-offs than a single global \(\varepsilon\), acknowledging cross-modality disparities.

The comparison with pruning baselines is restricted in the paper to WikiText-2 perplexity. For LLaMA-2-7B at 10%/30%/50% compression, PGSVD yields \(6.52/9.20/27.46\), outperforming LLM-Pruner and ShortGPT, and is competitive with SliceGPT at lower ratios. For LLaMA-2-13B, it reports \(5.36/7.09/18.04\), consistently better than LLM-Pruner and ShortGPT and competitive with SliceGPT at high compression. On reasoning accuracy at 20% compression for LLaMA-2-7B, PGSVD is reported to maintain balanced performance, with PIQA \(71.27\) and ARC-E \(70.75\), while SliceGPT drops sharply.

In comparative positioning, PGSVD is distinguished from truncated SVD/PCA, LoRA, and prior activation-aware methods such as Weighted SVD, ASVD, and SVD-LLM. The stated differences are threefold: a network loss bound linking \(\|\Delta W_l X_l\|_F\) to \(|\Delta \mathcal{L}|\); a bi-objective formulation with a Pareto-guided uniform tolerance that induces heterogeneous ranks with theoretical guarantees; and an ALS solver that avoids the numerical fragility of whitening-based methods and is more efficient than EVD-based implementations.

The paper also states the principal limitations. The Pareto guidance assumes homogeneous sensitivity and bounded, nonincreasing convex envelopes for layer SVD profiles. The loss bound is first-order and conservative, so extreme compressions may exceed its validity. The uniform \(\varepsilon\) rule is surrogate-optimal, whereas data-driven per-layer \(\varepsilon\) choices could further improve accuracy at the cost of robustness and complexity. Future directions explicitly identified include learning tolerance allocations from data, combining PGSVD with quantization and lightweight fine-tuning, refining the loss bounds through tighter higher-order analysis, and extending the approach to larger and more diverse VLM architectures.

Source: https://www.emergentmind.com/topics/pareto-guided-singular-value-decomposition-pgsvd