---
title: Computable Distortion Index
url: https://www.emergentmind.com/topics/computable-distortion-index
type: topic
---

# Computable Distortion Index

Searching arXiv for the cited papers and closely related work on computable distortion indices.
A computable distortion index is a distortion-related quantity that can be evaluated algorithmically from a formal specification of a source, representation, estimator, or signal pair. In the cited literature, the term is not confined to a single construction. It appears as the rate–distortion value \(R(D)\) or its finite-blocklength feed-forward analogue \(R_n(D)\), as a distortion–perception tradeoff function \(D(P)\), as the Class Angular Distortion Index (CADI) for dimensionality reduction, as the Feature Distortion Index (FDI) for model-extraction detection, and as the Gammachirp Envelope Distortion Index (GEDI) for speech intelligibility prediction. A recurring distinction is that scalar performance values can be computable even when the exact optimizing object that attains them is not [2511.09412].

## 1. Conceptual scope and recurring forms

In the cited work, a computable distortion index is always an explicit map from structured inputs to either a scalar or a curve. The inputs vary by domain: a source law and distortion matrix in rate–distortion theory, a joint law \(p_{X,Y}\) and a perceptual metric in distortion–perception analysis, a labeled high-dimensional dataset and its embedding in dimensionality reduction, a query and anchor features in MLaaS security, or a clean/enhanced speech pair in auditory modeling. What is common is that the index is specified by an algorithmic pipeline and is intended to be reproducible from the available mathematical or signal representation.

| Domain | Index | Computable object |
|---|---|---|
| Rate–distortion theory | \(R(D)\) | Minimum achievable rate at distortion \(D\) |
| Feed-forward source coding | \(R_n(D)\) | Finite-\(n\) directed-information rate |
| Distortion–perception theory | \(D(P)\) | Minimum distortion under perception constraint |
| Dimensionality reduction | CADI | Mean squared cosine-angle distortion over between-class triplets |
| MLaaS security | FDI | Layerwise feature-distance vector to class anchors |
| Speech intelligibility | GEDI / mr-GEDI | Envelope-domain signal-to-distortion ratio |

This range of usages indicates that computability can refer either to a value function produced by optimization or to a directly evaluated residual/error functional. In some settings the output is a single scalar, as with CADI or GEDI. In others it is a parametric tradeoff curve, as with \(R(D)\), \(R_n(D)\), or \(D(P)\). The cited papers also show that computability is domain-specific: in finite-alphabet information theory it is formulated through computable analysis, while in applied settings it typically means that the index can be evaluated with explicit finite procedures, sampling schemes, or convex programs [2107.02555].

## 2. Information-theoretic computability and the optimizer gap

For a discrete memoryless source over finite alphabets, with computable source law \(P_X\), computable single-letter distortion \(d\), and computable distortion level \(D\), the rate–distortion function is
\[
R(D) = \inf_{P_{Y|X}:\,\mathbb{E}[d(X,Y)]\le D} I(X;Y).
\]
The corresponding optimizing test channel is any
\[
P_{Y|X}^* \in \arg\inf_{P_{Y|X}:\,\mathbb{E}[d(X,Y)]\le D} I(X;Y).
\]
The paper "Computability of the Optimizer for Rate Distortion Functions" establishes a precise asymmetry: \(R(D)\) is computable, but the optimizer is in general non-computable, even for simple finite distortion measures including Hamming-type settings [2511.09412].

The computability framework used there is computable analysis in the sense of Pour-El and Richards and Weihrauch. A computable real is represented by an effectively convergent rational sequence; Banach–Mazur computability and uniform effective continuity give the stronger function notion on compact domains. Within that framework, a classic minimum-versus-argmin dichotomy is recalled: the minimum of a computable function on a compact set is computable, but an argmin need not be. The rate–distortion problem inherits exactly that pattern. Mutual information is continuous and computable in \(P_X\) and \(P_{Y|X}\), the feasible set \(\{P_{Y|X}:\mathbb{E}[d(X,Y)]\le D\}\) is compact and effectively specified, and therefore the value \(R(D)\) is computable. By contrast, there is no universal Banach–Mazur, hence no universal Turing-computable, selector
\[
F_{\mathrm{opt}}:(P_X,d,D)\mapsto P_{Y|X}^*
\]
for all computable instances with finite alphabets \(|\mathcal X|\ge 2\), \(|\mathcal Y|\ge 2\), and nontrivial \(P_X\) [2511.09412].

The negative result is strengthened in two directions. First, even if \(P_X\) is fixed and only the distortion measure varies, a computable optimizer map need not exist. Second, even if the distortion measure is fixed—including binary weighted Hamming distortion and more general normal distortion measures—there is no general computable map from \((P_X,D)\) to an exact optimal test channel. The constructions encode a recursively enumerable but non-recursive set into the structure of the optimizer, using explicit distortion matrices and source perturbations. The paper also states that there is not even a uniformly computable approximator that recovers the optimizer to arbitrary precision in total variation distance [2511.09412].

A related but operationally distinct construction appears in source coding with feed-forward. For stationary, ergodic, discrete-time sources, the finite-blocklength quantity
\[
R_n(D)=\frac{1}{n}\min I(\hat{X}^n\rightarrow X^n)
\]
is achievable for any \(n\), and \(\lim_{n\to\infty}R_n(D)\) is the feed-forward rate–distortion function [1106.0895]. Here the distortion index is not the optimizer but the sequence of computable rates themselves. The paper provides an alternating minimization procedure, derived as an extended Blahut–Arimoto algorithm, and also a dual formulation that becomes a geometric programming problem. This yields a sequence of achievable and computable rates converging to the operational limit [1106.0895].

Taken together, these results separate value computability from structure computability. This suggests that, in information-theoretic usage, a computable distortion index is safest when defined through optimal values such as \(R(D)\), \(R_n(D)\), or directly computable thresholds such as
\[
D_{\max}=\min_{y\in\mathcal Y}\mathbb E[d(X,y)],
\]
rather than through exact optimal test channels or induced codebook distributions [2511.09412].

## 3. Distortion–perception functions as computable indices

A second major line of work defines a computable distortion index as a distortion–perception tradeoff function. For random vectors \(X\) and \(Y\), with estimator \(\hat X\) satisfying the Markov condition \(X\to Y\to \hat X\), the general distortion–perception function is
\[
D(P)=\min_{p_{\hat X|Y}}\left\{\mathbb E[d(X,\hat X)] : d_p(p_X,p_{\hat X})\le P\right\}.
\]
In the MSE–Wasserstein-2 setting, with
\[
d(x,\hat x)=\|x-\hat x\|^2,\qquad d_p(p_X,p_{\hat X})=W_2(p_X,p_{\hat X}),
\]
the paper "A Theory of the Distortion-Perception Tradeoff in Wasserstein Space" proves that the DP function is always quadratic:
\[
D(P)=D^*+\left[(P^*-P)_+\right]^2,
\]
where \(D^*=\mathbb E[\|X-\mathbb E[X|Y]\|^2]\) and \(P^*=W_2(p_X,p_{X^*})\) with \(X^*=\mathbb E[X|Y]\) [2107.02555].

That result is geometric rather than merely variational. The minimizers on the DP curve form a Wasserstein-2 geodesic between the distribution of the signal and that of the minimum-MSE estimator. If \(\hat X_0\) is an optimal perfect-perception estimator satisfying \(p_{\hat X_0}=p_X\), then every optimal estimator for \(P\in[0,P^*]\) can be written as
\[
\hat X_P=\left(1-\frac{P}{P^*}\right)\hat X_0+\frac{P}{P^*}X^*.
\]
In the Gaussian case, the paper gives closed-form optimal estimators and a computable lower bound via the Gelbrich distance
\[
D(P)\ge D^*+\bigl[(G^*-P)_+\bigr]^2,
\]
with equality for jointly Gaussian \((X,X^*)\) [2107.02555].

For finite-alphabet channels, the same idea becomes a linear-programming object. The paper "Characterization of the Distortion-Perception Tradeoff for Finite Channels with Arbitrary Metrics" studies arbitrary distortion matrices and Wasserstein-\(1\) perception induced by an arbitrary metric on a finite source alphabet. In that setting, computing both the DP function and optimal reconstructions is equivalent to solving linear programs. The distortion criterion is affine in the stochastic decoder \(q_{\hat X|Y}\), and the Wasserstein-\(1\) perception constraint is encoded through an explicit coupling \(\Pi\). The resulting DP function is non-increasing, convex, and piecewise linear in the perception index \(P\) [2402.02265].

The finite-alphabet paper also derives a dual characterization in which
\[
D(P)=\max_i (p_0^i+p_1^i P),
\]
so the curve is the upper envelope of finitely many affine functions. In the binary case, it provides closed-form breakpoints, slopes, and optimal deterministic threshold reconstructions, with randomized interpolation between adjacent breakpoints. This makes the distortion index fully computable in the sense of explicit LP evaluation or explicit closed form, depending on the alphabet size [2402.02265].

These two DP frameworks reveal two distinct computable geometries. In Wasserstein-2 with MSE, the tradeoff is quadratic and geodesic. In finite alphabets with Wasserstein-1 and arbitrary metrics, it is polyhedral and piecewise linear. In both cases, the index is the optimal value function, not a heuristic proxy.

## 4. Angular distortion and class geometry in dimensionality reduction

In dimensionality reduction, the computable distortion index can be defined directly on the geometry of labeled point sets rather than through an optimization over encoders or decoders. The Class Angular Distortion Index (CADI) is defined for high-dimensional data
\[
X=\{x_1,\dots,x_n\}\subset\mathbb R^d,
\]
a low-dimensional embedding
\[
Y=\{y_1,\dots,y_n\}\subset\mathbb R^t,\quad t\in\{2,3\},
\]
and a class partition \(C=\{C_1,\dots,C_m\}\). For between-class triplets with reference point \(i\in C_a\) and unordered pair \(j,k\in C_b\neq C_a\), CADI is
\[
\mathrm{CADI}(X,Y,C)
=
\frac{1}{T}
\sum
\left(
\cos\theta_X(j,i,k)-\cos\theta_Y(j,i,k)
\right)^2.
\]
The cosines are computed from inner products and norms of the corresponding edge vectors in the original and embedded spaces [2605.00637].

CADI measures distortion at the level of class-to-class geometry. Because it uses between-class triplets, it is sensitive to how one class is seen from another: nested shells, rings, interlinked manifolds, anisotropic spreads, and non-globular cluster shapes. The paper positions it against metrics such as Silhouette Score, Davies–Bouldin Index, NMI, ARI, centroid-based Cluster Distance Score, Label Trustworthiness/Continuity, and Steadiness & Cohesiveness. Its stated role is not separability alone, but the preservation of inter-class manifold geometry [2605.00637].

The naïve computation is at least \(\mathcal O(n^3)\), since the number of between-class triplets is on that order. To make the index computably practical, the paper uniformly samples
\[
k=\mathcal O(n)
\]
between-class triplets. It reports that \(k=10n\) already gives stable estimates to approximately three decimal digits, and uses \(k=40n\) in its evaluations. With fixed dimensionalities, this makes runtime linear in \(n\), and the paper reports evaluation times around a second for its largest datasets on a consumer laptop [2605.00637].

CADI is also differentiable because each triplet loss
\[
\ell_{ijk}(Y)=\left(\cos\theta_X(j,i,k)-\cos\theta_Y(j,i,k)\right)^2
\]
depends smoothly on dot products and norms of the embedded coordinates, except at degenerate zero-norm cases, which are handled explicitly. This enables AngleEmbedding, a supervised dimensionality reduction method that minimizes CADI with a feedforward MLP and Adam at learning rate \(10^{-3}\), using sampled triplets during training [2605.00637].

Several limitations are explicit. CADI is not a separability metric; it can assign low distortion to embeddings with overlap if the class geometry is preserved. It requires class labels or a clustering. Larger classes contribute more triplets by design, because the normalization is by sampled triplet count rather than by class size. It is invariant to translation, rotation, and uniform scaling, but sensitive to non-uniform feature scaling and nonlinear transforms [2605.00637].

## 5. Feature-space distortion indices in MLaaS security

In MLaaS model-extraction defense, the Feature Distortion Index (FDI) instantiates a computable distortion index as a layerwise deviation between the internal feature representation of an incoming query and anchor features derived from the victim model’s training set. The setting is a private DNN \(F_V\) deployed behind a prediction API. For each class \(c\), the defender selects \(K\) high-confidence anchor samples from the training set and precomputes their features at \(L\) layers. For a query \(x\), the victim’s predicted class \(c=F_V(x)\) determines which anchor set is used [2306.11338].

The per-layer, per-anchor index is
\[
\mathcal I_j^l
=
d\!\left(F_V^l(x)-F_V^l(x_{c,j})\right),
\]
with \(d\) chosen as the \(\ell_2\) norm, so concretely
\[
\mathcal I_j^l=\left\|F_V^l(x)-F_V^l(x_{c,j})\right\|_2.
\]
The full FDI vector is the concatenation over all anchors \(j=1,\dots,K\) and considered layers \(l=1,\dots,L\):
\[
\mathcal I(x,c;\mathcal D_{\mathrm{anc}},F_V)\in\mathbb R^{K\times L}.
\]
For VGG19 on CIFAR-10, the paper uses \(K=100\) and \(L=5\), giving a 500-dimensional FDI vector [2306.11338].

This index is not a divergence between explicitly estimated continuous distributions. Instead, the distributional claim is operationalized through distances to class-specific anchor features. The anchor set approximates the center of each class cluster in feature space, and the FDI vector records how much an incoming query deviates from that benign feature distribution. FDINet uses these vectors to train a binary detector with architecture Linear\([K\times L,256]\), Linear\([256,32]\), Linear\([32,2]\), followed by Softmax, optimized by cross-entropy on benign training data and an auxiliary out-of-distribution positive set [2306.11338].

Detection is performed over batches rather than single queries. For a batch of size \(bs\), each FDI vector is passed through the classifier, converted to a binary malicious/benign decision by \(\arg\max\), and averaged:
\[
ac=\frac{1}{bs}\sum_{i=1}^{bs}\arg\max \mathcal C\!\left(\mathcal I(x_i,c_i;\mathcal D_{\mathrm{anc}},F_V)\right).
\]
A threshold \(\tau_1=\mu_{ac}+3\sigma_{ac}\), estimated from benign validation batches, is used to flag malicious clients. The same FDI representation is also used for colluding-adversary detection via two-sample \(t\)-tests on FDI statistics across clients, with a significance threshold \(\tau_2\) such as \(0.05\) [2306.11338].

The empirical claims are specific. FDINet achieves 100% detection accuracy on DFME and DaST, uses 50 queries to raise an extraction alarm with an average confidence of 96.08% for GTSRB, identifies colluding adversaries with accuracy exceeding 91%, and detects two adaptive attacks, Dummy Query and Feature Correction, though with some degradation on attacks whose queries are closer to benign data such as Knockoff and ActiveThief [2306.11338].

As a distortion index, FDI is fully computable online because anchor features are precomputed and the per-query computation reduces to forward feature extraction and \(O(KL)\) Euclidean distances. Its limitations are also domain-specific: it depends on label-conditioned anchors, it was evaluated in computer vision rather than in other modalities, and similarity-based colluding detection is weaker when distinct attacks induce very similar feature distributions [2306.11338].

## 6. Auditory-envelope distortion indices and cross-domain synthesis

GEDI and mr-GEDI define a computable distortion index for speech intelligibility prediction by operating in the temporal envelope after auditory and modulation filtering. The pipeline begins with a dynamic compressive gammachirp filterbank (dcGC-FB) with \(I=100\) channels, equally spaced in ERB\(_\mathrm N\)-number from 100 to 6000 Hz. For each auditory channel, the clean and enhanced signals are transformed into Hilbert envelopes and low-pass filtered at 150 Hz. Envelope distortion is then defined, with \(p=2\), as
\[
e_{D,i}(n)=\Bigl(\bigl|\{e_{S,i}(n)\}^p-\{e_{\hat S,i}(n)\}^p\bigr|\Bigr)^{1/p}.
\]
This produces an explicit distortion signal in the auditory-envelope domain [1904.02096].

GEDI then analyzes modulation power. For clean and distortion envelopes, modulation-filter outputs are summarized as powers \(P_{\mathrm{env},S,i,j}\) and \(P_{\mathrm{env},D,i,j}\), with an ERB\(_\mathrm N\)-based weighting
\[
W_i=\frac{\mathrm{ERB_N}(1000)}{\mathrm{ERB_N}(f_i)}.
\]
For each modulation channel \(j\), the weighted envelope-domain signal-to-distortion ratio is
\[
\mathrm{SDR}_{\mathrm{env},j}^W
=
\frac{\sum_{i=1}^{100} W_i P_{\mathrm{env},S,i,j}}
{\sum_{i=1}^{100} W_i P_{\mathrm{env},D,i,j}},
\]
and the global GEDI core quantity is
\[
\mathrm{SDR}_{\mathrm{env}}
=
\sqrt{\sum_{j=1}^{J}\left(\mathrm{SDR}_{\mathrm{env},j}^W\right)^2}.
\]
This scalar is then mapped, via an ideal-observer model, to predicted intelligibility [1904.02096].

Mr-GEDI extends the same construction to non-stationary noise by using multi-resolution temporal analysis. Instead of a single FFT-based modulation analysis, it employs one third-order low-pass and eight second-order band-pass IIR modulation filters, with frame lengths inversely matched to modulation center frequency. It computes framewise powers and framewise
\[
\mathrm{SDR}_{\mathrm{env},j,t}^W
=
\frac{\sum_{i=1}^{100} W_i P_{\mathrm{env},S,i,j,t}}
{\sum_{i=1}^{100} W_i P_{\mathrm{env},D,i,j,t}},
\]
averages them over frames, and then pools across modulation bands. The paper reports that mr-GEDI predicts intelligibility curves better than STOI, ESTOI, and HASPI under pink-noise conditions, and better than HASPI under babble-noise conditions. It further states that mr-GEDI does not present an overestimation tendency and is considered a more conservative approach than STOI and ESTOI [1904.02096].

Across the cited literature, several misconceptions are explicitly corrected. First, computable does not imply that the exact optimizer is computable; in finite-alphabet rate–distortion theory, the opposite is shown for \(P_{Y|X}^*\) [2511.09412]. Second, a low-distortion geometric index need not mean strong class separability; CADI is designed to preserve class geometry, not to optimize blob-like separation [2605.00637]. Third, an index can be practically computable without being a universal exact solver; the feed-forward \(R_n(D)\) sequence, LP-based finite-channel DP functions, sampled CADI, and FDI all rely on explicit finite procedures rather than inaccessible analytic formulas [1106.0895].

A plausible synthesis is that computable distortion indices fall into two broad categories. The first consists of optimal-value functions, such as \(R(D)\), \(R_n(D)\), and \(D(P)\), where computability attaches to the optimum value of a convex or linear program. The second consists of direct residual measures, such as CADI, FDI, and GEDI, where computability attaches to an explicitly evaluated discrepancy in a chosen representation. The cited papers jointly show that the design choice is consequential: computable values are often robustly definable, whereas exact optimizing structures may be non-computable, non-unique, or only practically approximable.

Source: https://www.emergentmind.com/topics/computable-distortion-index