Unverified 2026

Inverse-Eigenvector Tight-Frame Codebook

Implementation & benchmark of arXiv:2607.26048 — Isotropic Decompositions via Inverse Eigenvectors

Usefulness5/10
Difficulty6/10
Novelty6/10

Source paper: Isotropic Decompositions via Inverse Eigenvectors arXiv:2607.26048 ⓘ · analyzed Aug 31, 2026

AI-generated research hypothesis, automatically tested. Not peer-reviewed.

Idea description

Construct a finite neural prototype dictionary from solutions of Mα = α⁻¹, where the inverse is coordinatewise, and assign positive weights so the dictionary obeys the isotropy identity Σᵢ cᵢαᵢαᵢᵀ = I. Use the resulting frame as the initialization or fixed geometry for embedding prototypes, attention directions, or MoE router experts instead of initializing those vectors independently. The isotropy guarantee should reduce directional collapse and make early optimization gradients more uniform.

Formulas

$$M\alpha=\alpha^{-1},\qquad (\alpha^{-1})_j=\frac{1}{\alpha_j},\qquad \alpha_j(M\alpha)_j=1\quad (j=1,\ldots,n).$$
$$I_n=\sum_{i=1}^{N}c_i(\alpha_i\otimes\alpha_i)=\sum_{i=1}^{N}c_i\alpha_i\alpha_i^{\top},\qquad c_i>0,\qquad \sum_{i=1}^{N}c_i=1.$$
$$z=\sum_{i=1}^{N}c_i(\alpha_i^{\top}z)\alpha_i,\qquad \|z\|_2^2=\sum_{i=1}^{N}c_i(\alpha_i^{\top}z)^2.$$
$$\mathcal{L}_{\mathrm{frame}}=\left\|I_n-\sum_{i=1}^{N}c_i\alpha_i\alpha_i^{\top}\right\|_F^2+\lambda\sum_{i=1}^{N}\left\|M\alpha_i-\alpha_i^{-1}\right\|_2^2,\qquad c_i=\frac{e^{u_i}}{\sum_k e^{u_k}}.$$

Mathematical statement

An inverse eigenvector is a vector α ∈ Rⁿ with all coordinates nonzero satisfying Mα = α⁻¹, equivalently αⱼ(Mα)ⱼ = 1 for every coordinate j, where M ∈ Rⁿˣⁿ is a real correlation matrix and α⁻¹ = (1/α₁, …, 1/αₙ). The paper's central construction gives the identity decomposition Iₙ = Σᵢ₌₁ᴺ cᵢ(αᵢ ⊗ αᵢ) = Σᵢ₌₁ᴺ cᵢαᵢαᵢᵀ, with positive coefficients cᵢ > 0 satisfying Σᵢcᵢ = 1 for real correlation matrices whose off-diagonal entries have modulus strictly below one. Here Iₙ is the n × n identity, N is the number of selected inverse eigenvectors, and αᵢ are vectors in Rⁿ. The tensor is bilinear: α ⊗ α = ααᵀ. Consequently, for every z ∈ Rⁿ, Σᵢcᵢ(αᵢᵀz)αᵢ = z and Σᵢcᵢ(αᵢᵀz)² = ||z||₂². For approximate roots and weights, monitor the root residuals rᵢ = Mαᵢ − αᵢ⁻¹ and the frame residual R = Iₙ − Σᵢcᵢαᵢαᵢᵀ.

Implementation notes

(1) Exact integration point: use the construction as an offline-generated codebook for an embedding layer, an MoE expert/router direction matrix, or the key and value direction matrices of a small attention head. Let the feature width be n and choose N at least n. Estimate M as the empirical correlation matrix of normalized calibration embeddings, then regularize it with M ← (1 − ε)M + εI to avoid singularity. Store the αᵢ vectors and normalized weights cᵢ as an initialization; optionally allow a small trainable correction to each αᵢ. (2) Pseudocode: normalize calibration samples h; compute M; initialize every αᵢ with nonzero coordinates; repeat, gᵢ = Mαᵢ − αᵢ⁻¹, update αᵢ ← αᵢ − η times a damped Adam or Newton step on gᵢ, and clip |αᵢⱼ| below by δ; then optimize unconstrained logits uᵢ with c = softmax(u) using the frame loss. Alternate root-residual and frame-residual updates. For a router, compute logits ℓᵢ = hᵀαᵢ divided by square root of n, then apply softmax or top-k routing and add β times the frame loss to the neural training objective. (3) Computed from the mathematics: inverse-eigenvector residuals, simplex-constrained positive weights, and the Frobenius tight-frame residual. Estimated empirically: M, ε, damping, clipping δ, and the regularization coefficient β. (4) First experiment: train a width-64 MLP-Mixer or four-expert top-2 MoE on CIFAR-10, comparing this codebook with Gaussian and orthogonal initialization. Measure loss versus optimizer steps, feature-coordinate gradient variance, frame residual, router expert-load entropy, and final accuracy. Success means faster early loss descent at equal FLOPs, lower gradient anisotropy, and fewer collapsed MoE experts without sacrificing final accuracy.

Verification

This idea has not been verified yet.

Verification happens in two stages: Stage 1 — a mechanism check on a toy system confirms the claimed mathematical phenomenon reproduces; Stage 2 — a benchmark implements the idea on a real (small) neural network task and compares it against a tuned baseline over 8 paired seeds with a permutation test.

Artifacts

Artifacts unavailable.