Schur-Riesz Greedy Adapter Expansion
Implementation & benchmark of arXiv:2608.11764 — Schur-Riesz Variational Enrichment: A Generalized Refinement Framework for Finite Elements
Source paper: Schur-Riesz Variational Enrichment: A Generalized Refinement Framework for Finite Elements arXiv:2608.11764 ⓘ · analyzed Sep 1, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Replace blind width or rank expansion with a residualized candidate-selection procedure for low-rank adapters or expert banks. Each candidate block is measured after removing the component already represented by the incumbent network, filtered for numerical conditioning, and accepted only when its projected response produces sufficient loss-relevant gain. This should produce smaller adaptive networks than adding candidates by fixed order or magnitude heuristics.
Formulas
Mathematical statement
The paper defines the residual-space innovations Q_hp=(I-P_{Z_m}^Y)(T\circ C_hp) and Q_op=(I-P_{Z_m}^Y)(T\circ C_op), where V_m is the current approximation space, Z_m=T(V_m) is its image under the variational operator T, C_hp and C_op are candidate blocks, P_{Z_m}^Y is projection onto the incumbent image under the Y-inner product, and I-P_{Z_m}^Y removes already represented content. For a neural adaptation, let B_m in R^{n x r_m} contain incumbent responses on n calibration examples, let C_j in R^{n x r_j} contain responses of candidate block j, let T be either the identity on logits or features or a minibatch Jacobian mapping parameter perturbations to logits, and let Y be a positive diagonal example or output weighting matrix. Compute P_{B_m}^Y=B_m(B_m^T Y B_m)^dagger B_m^T Y and Q_j=(I-P_{B_m}^Y)T C_j. Rank by projected gain g_j=||Q_j||_Y^2=tr(Q_j^T Y Q_j). A practical Riesz-style stability test retains j only if the nonzero singular values of Y^{1/2}Q_j lie in [alpha,beta], with alpha>0 preventing nearly null or redundant coordinates and beta detecting badly scaled candidates. The exact paper-specific Schur gain and finite-element certificates are not present in the extraction; g_j is an explicit neural surrogate whose success is falsifiable.
Implementation notes
(1) Integration point: apply the procedure during progressive construction of LoRA adapters, feature banks, or MoE experts, while keeping the base Transformer frozen in the first experiment. For a selected layer, candidate j is a low-rank update Delta W_j=A_j B_j, an extra MLP feature group, or an expert. Collect its response C_j on a fixed calibration batch. Use T=I on logits for a cheap version, or use Jacobian-vector products T C_j to measure the candidate's first-order effect on logits. The incumbent matrix B_m contains responses from accepted blocks. (2) Pseudocode: initialize B as an empty matrix. At each growth round, generate K candidates; evaluate R_j=T C_j; compute the Y-orthogonal projection using P=B(B^T Y B)^dagger B^T Y; form Q_j=(I-P)R_j; calculate g_j=trace(Q_j^T Y Q_j) and the singular values of Y^{1/2}Q_j; discard candidates with sigma_min_plus below alpha or sigma_max above beta; choose the largest remaining g_j; train only that block for s steps; accept it if held-out loss decreases by at least threshold tau; append its response columns to B; recompute the projection and continue. Stop when no candidate passes tau or a parameter budget is reached. (3) Compute P, Q_j, and g_j from the residualization formula. Estimate alpha, beta, Y, s, and tau empirically. Use QR or randomized SVD rather than explicitly constructing I-P. (4) First cheap experiment: compare against fixed-rank LoRA and random or magnitude-based adapter expansion on a 125M-parameter Transformer fine-tuned on a 10k-example instruction dataset or GLUE task. Plot validation loss versus trainable parameters and FLOPs. Success means at least 20% fewer adapter parameters at matched validation loss, or faster loss descent at equal parameter count, with accepted blocks having substantially larger projected response energy than rejected blocks.
Verification
Stage 1 · Toy mechanism gate: Passed ✓
Stage 2 · Mechanism transferred to benchmark: Not reproduced ✗
Stage 2 · Practical benchmark result: No measurable effect
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 NumPy MVP of Schur-Riesz greedy adapter expansion with Y-weighted projection, residual candidate responses, gain ranking, singular-value filtering, and fixed-order comparison. The math checks manifested quantitatively: projection orthogonality was 2.85e-14, idempotence error was 1.53e-15, residual gain was approximately zero at zero novel amplitude and followed exact quadratic scaling (ratios 0.0625, 0.25, 1, 4 for amplitudes 0.25, 0.5, 1, 2), while predicted and observed least-squares loss drops matched to 1.83e-16 relative error. In the controlled toy regression, greedy selected all three useful blocks and reached loss 9.51 versus 483.09 for fixed-order expansion, a 98.0% lower final loss at the same three added scalar blocks; this is a clear mechanism signal, not evidence of production-scale neural-adapter gains.
- Agent confidence
- 9/10
- Baseline
- Fixed-order expansion selected blocks [0,1,2], all labeled redundant; loss 501.51 -> 483.09 after 3 additions.
- Idea
- Projected-gain greedy selected blocks [7,5,11], all labeled useful; loss 501.51 -> 9.51 after 3 additions; final loss ratio idea/baseline = 0.01968. Math: orthogonality 2.85e-14, idempotence 1.53e-15, gain amplitude ratios exactly [0.0625,0.25,1,4], loss-drop relative prediction error 1.83e-16.
Limitations: This is a deterministic synthetic weighted least-squares experiment, not a Transformer, LoRA, MoE, or held-out neural fine-tuning benchmark. It does not measure GPU runtime, FLOPs, parameter-memory savings, train-after-selection acceptance thresholds, minibatch Jacobian operators, or robustness across datasets and random seeds.
How to run: python3 schur_riesz_mvp.py > results.json
Stage 2 — Benchmark vs. tuned baseline (arXiv:2608.11764)
No significant difference vs. the tuned baseline (+0.4%, p=0.185, wins 3 / 8 paired seeds).
- Benchmark
- Tabular regression (Friedman)
- Model
- candidate_mlp_shared
- Paired seeds
- 8
- Baseline mean
- 21.9049
- Idea mean
- 21.9884
- Effect (Δ)
- 0.0835 (+0.4%; negative = idea better)
- Wins
- 3 / 8 paired seeds
- p-value
- 0.185 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- ±0.6%
- Mechanism
- Not confirmed ✗
- Practical effect
- No effect
- Baseline tuning
- swept over 3 configs
Only the tabular track and a small candidate-bank MLP were tested; no LoRA, Transformer, MoE, vision, sequence, or dynamics implementation was evaluated. The candidate responses were direct output-block responses rather than minibatch Jacobian-mapped logits, and parameter-count/runtime savings were not measured.
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 Schur-Riesz greedy candidate expansion as a trained candidate-bank MLP and ran the full 8-seed paired protocol with an equal three-point learning-rate sweep. The idea had test MSE 21.9884 versus baseline 21.9049, paired delta +0.0835 and permutation p=0.18505, so it was not a significant win. The behavior-based mechanism signature was not confirmed because predicted projected gain correlated negatively with observed held-out selected-response energy (correlation -0.4155).
- Baseline
- Best lr=0.01; full 8-seed test MSE mean 21.9048938751, std 1.9880925148.
- Comparison
- delta_mean=0.0834848881, p_value=0.18505, idea_wins=3/8, verdict=no measurable effect.
- Idea
- Best lr=0.01; full 8-seed test MSE mean 21.9883787632, std 1.8636409833.
- Mechanism signature
- predicted_projected_gain_mean=26.9075886458, observed_heldout_selected_response_mse_mean=65.3522763252, correlation=-0.4154841401, confirmed=false
Benchmark result
No significant difference vs. the tuned baseline (+0.4%, p=0.185, wins 3 / 8 paired seeds).
- Benchmark
- Tabular regression (Friedman)
- Model
- candidate_mlp_shared
- Paired seeds
- 8
- Baseline mean
- 21.9049
- Idea mean
- 21.9884
- Effect (Δ)
- 0.0835 (+0.4%; negative = idea better)
- Wins
- 3 / 8 paired seeds
- p-value
- 0.185 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- ±0.6%
- Mechanism
- Not confirmed ✗
- Practical effect
- No effect
- Baseline tuning
- swept over 3 configs
Only the tabular track and a small candidate-bank MLP were tested; no LoRA, Transformer, MoE, vision, sequence, or dynamics implementation was evaluated. The candidate responses were direct output-block responses rather than minibatch Jacobian-mapped logits, and parameter-count/runtime savings were not measured.
How to run: python3 stage2_bench.py
Artifacts
- bench_report.json 3.9 KB View Raw
- report.md 1.9 KB View
- report_bench_2026-09-01T124330.md 4.3 KB View
- results.json 1.7 KB View Raw
- schur_riesz_mvp.py 5.6 KB View Raw
- stage2_bench.py 5.3 KB View Raw