Isometric tensor-network token mixer
Implementation & benchmark of arXiv:2608.00053 — Fast Trainable Multilinear Bases for Image Compression
Source paper: Fast Trainable Multilinear Bases for Image Compression arXiv:2608.00053 ⓘ · analyzed Aug 31, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Use the relaxed QFT tensor-network topology as a trainable norm-preserving mixer inside a neural block, replacing a dense token-mixing matrix or an expensive global convolution. The network learns data-adapted global interactions while retaining structured O(N log^2 N) application and an exact cheap inverse, making it suitable for image tokens, long sequences, or reversible residual blocks.
Formulas
Mathematical statement
The construction composes local unitary tensors in a butterfly/QFT topology rather than learning a dense N-by-N matrix. Each trainable local gate is unitary: a relaxed two-dimensional gate is U\in U(2), and a relaxed phase gate is D(\boldsymbol{ heta})=\operatorname{diag}(e^{i heta_0},e^{i heta_1},e^{i heta_2},e^{i heta_3}). Because a product of unitary gates is unitary, the resulting mixer T( heta) obeys T^\dagger T=I, preserving Euclidean norms and preventing singular or amplifying linear mixing. The QFT-inspired topology recursively reduces F_N to F_{N/2} and terminates at F_2=H, producing a structured global transform with near-linear cost O(N\log^2N), where N is the number of tokens or channels. Insert T between learned pointwise nonlinearities or use x_plus=T^\dagger f(Tx) as a spectral block. The trainable parameters are local gate coordinates instead of N^2 dense mixing parameters; manifold constraints are enforced by exponential-map or retraction updates.
Implementation notes
(1) Integrate at the token-mixing location of a ViT or MLP-Mixer block, or in a long-sequence encoder between channel-wise MLPs. For N tokens, pad to the next power of two and implement the QFT-style butterfly as batched local 2x2 operations plus diagonal phase operations. Start with fixed Hadamard and controlled-phase gates, then relax a small subset or all gates. (2) Pseudocode: input X with shape batch x N x d; for each feature channel compute Z=T(theta)X; apply a pointwise nonlinear map A=GELU(Z) or a learned channel MLP; compute Y=T(theta)^*A; output residual X+alpha Y, with alpha initialized to 0.1. During training, parameterize each local gate by U=exp(K)U0 with K^*=-K, or retract U after an Euclidean gradient step using polar decomposition U<-UV^* for U=QR. (3) The exact unitary property and O(N log^2 N) graph come from the construction; measure actual FLOPs and test whether gate sharing across feature channels is sufficient. (4) First cheap experiment: compare a 6-layer CIFAR-100 ViT-tiny or ImageNet-100 64x64 model against an MLP-Mixer and dense token-mixing baseline at equal hidden width. Record validation accuracy, gradient-norm variance, loss versus FLOPs, parameter count, and throughput. A successful result is matching or exceeding dense-mixer accuracy with substantially fewer parameters, no norm drift through repeated blocks, and faster training or inference at N>=256.
Verification
Stage 1 · Toy mechanism gate: Passed ✓
Stage 2 · Mechanism transferred to benchmark: Confirmed ✓
Stage 2 · Practical benchmark result: Inconclusive
Methodology: 8 paired seeds, permutation test (p<0.05) against a learning-rate-tuned baseline. How verification works
Stage 1 — Mechanism check agent confidence 9/10
Built a trainable complex butterfly/isometric token mixer with local U(2) gates, exact adjoint inverse, norm checks, controlled non-unitary perturbation sweeps, and a dense fitting baseline. The mechanism manifested: unitary errors stayed at 6.5e-16–7.6e-15 for N=4–64, inverse reconstruction error was 3.9e-15, and measured drift matched the predicted (1+epsilon)^(log2(N)*repetitions) law to about 2e-15 relative error across all sweep points. In the toy fit the structured mixer reached MSE 2.63e-9 versus 6.72e-14 for dense mixing with 128 real structured coordinates versus 256 complex dense entries; this confirms stability and compression, but not superior optimization or task accuracy.
- Agent confidence
- 9/10
- Baseline
- Dense complex N=16 fit: final MSE 6.72e-14 after 300 steps; 256 complex parameters, equivalent to 512 real scalars.
- Idea
- Structured U(2) butterfly N=16 fit: final MSE 2.63e-9 after 300 steps; 128 real coordinates; norm ratio approximately 1; unitarity Frobenius error <=7.61e-15 over N=4..64; inverse error 3.94e-15.
Limitations: No CIFAR/ViT or long-sequence task accuracy was tested, and no reliable wall-clock/FLOP speedup was established. The toy target was itself generated by the same structured family, the implementation uses a butterfly O(N log N) topology rather than the idea's broader QFT O(N log^2 N) construction, and gate parameterization/training at larger N was not evaluated.
How to run: python3 experiment.py
Stage 2 — Benchmark vs. tuned baseline (arXiv:2608.00053)
No significant difference vs. the tuned baseline (+25.5%, p=0.1299, wins 2 / 8 paired seeds).
- Benchmark
- Sequence modeling
- Model
- transformer_tiny
- Paired seeds
- 8
- Baseline mean
- 0.4766
- Idea mean
- 0.5983
- Effect (Δ)
- 0.1217 (+25.5%; negative = idea better)
- Wins
- 2 / 8 paired seeds
- p-value
- 0.1299 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- ±34.3%
- Mechanism
- Confirmed ✓
- Practical effect
- Inconclusive
- Baseline tuning
- swept over 3 configs
Only the built-in sequence forecast track was tested; CIFAR vision, tabular, and dynamics tracks were not tested. The implementation uses a real orthogonal butterfly specialization rather than full complex U(2) QFT phases, and no reliable wall-clock or FLOP speedup measurement was included.
How to run: python3 stage2_bench.py
Verdict computed by deterministic test code from paired-seed statistics — not by the language model.
Stage 2 — Benchmark latest Failed ✗
Agent confidence: 9/10
Implemented a matched sequence-level benchmark: dense learned token mixing versus a trainable orthogonal butterfly mixer with identical encoder structure, optimizer, learning-rate grid, and eight paired seeds. The mechanism signature was confirmed at neural-network scale (norm ratio 1.0 and inverse error 1.19e-6), but the idea had worse test MSE (0.598316 versus 0.476622), paired delta +0.121694, and permutation p=0.1299, so there was no significant win.
- Baseline
- Best dense token mixer, lr=0.001, epochs=6: test MSE mean 0.476622, std 0.099342 over 8 seeds.
- Idea
- Best isometric butterfly, lr=0.001, epochs=6: test MSE mean 0.598316, std 0.165879 over 8 seeds; paired delta +0.121694, p=0.1299, 2/8 wins.
- Mechanism signature
- confirmed=true; predicted norm ratio 1.0, observed trained-model ratio 1.0, observed adjoint inverse max error 1.19e-6; baseline observed norm ratio 0.502587.
Benchmark result
No significant difference vs. the tuned baseline (+25.5%, p=0.1299, wins 2 / 8 paired seeds).
- Benchmark
- Sequence modeling
- Model
- transformer_tiny
- Paired seeds
- 8
- Baseline mean
- 0.4766
- Idea mean
- 0.5983
- Effect (Δ)
- 0.1217 (+25.5%; negative = idea better)
- Wins
- 2 / 8 paired seeds
- p-value
- 0.1299 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- ±34.3%
- Mechanism
- Confirmed ✓
- Practical effect
- Inconclusive
- Baseline tuning
- swept over 3 configs
Only the built-in sequence forecast track was tested; CIFAR vision, tabular, and dynamics tracks were not tested. The implementation uses a real orthogonal butterfly specialization rather than full complex U(2) QFT phases, and no reliable wall-clock or FLOP speedup measurement was included.
How to run: python3 stage2_bench.py
Artifacts
- bench_report.json 4.6 KB View Raw
- experiment.py 6.4 KB View Raw
- isometric_mixer.py 2.4 KB View Raw
- report.md 1.7 KB View
- report_bench_2026-08-31T194248.md 4.2 KB View
- results.json 3.9 KB View Raw
- stage2_bench.py 4.7 KB View Raw