IQC-Certified Training Dynamics
Source paper: Generalization as a robust performance property of learning-enabled dynamical systems arXiv:2608.30431 ⓘ · analyzed Sep 1, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Represent a learned optimizer or recurrent training controller as a discrete-time feedback system and certify its sensitivity to one-sample dataset replacement using an IQC dissipativity inequality. Penalize the smallest certified disturbance-to-state gain during meta-training or use it as a post-training acceptance test, favoring update dynamics that do not amplify microscopic data perturbations over many iterations.
Formulas
Mathematical statement
Let x_k in R^n be the optimizer or recurrent-module state, y_k = C x_k + D u_k its feedback signal, and u_k = Phi(y_k,S) the data-dependent neural update. For neighboring datasets S and S', decompose the operator difference as u_k - u'_k = w_k + e_k, where w_k is the incremental response to the same feedback signal and e_k is the one-sample replacement disturbance. The sensitivity dynamics are delta_x_(k+1) = A delta_x_k + B w_k + B e_k and delta_y_k = C delta_x_k + D w_k + D e_k. Assume the incremental operator satisfies the IQC q_k^T Q q_k >= 0, where q_k = [delta_y_k; w_k], Q is a symmetric matrix, and delta_y_k = y_k - y'_k. A sufficient dissipativity certificate is an LMI in a storage matrix P, decay factor rho, disturbance gain gamma, and IQC multiplier lambda. It implies V_(k+1) - rho V_k + ||delta_x_k||^2 - gamma^2 ||e_k||^2 + lambda q_k^T Q q_k <= 0, with V_k = delta_x_k^T P delta_x_k. Summation gives a finite-horizon sensitivity bound, so if the one-sample operator sensitivity is at most epsilon_N, cumulative algorithmic stability scales as O(gamma epsilon_N).
Implementation notes
First, integrate the mechanism into a recurrent optimizer or learned training controller whose state contains parameters, momentum, adaptive statistics, or a compact controller state. A cheap baseline is heavy-ball training with state x_k = [theta_k; theta_(k-1)] and a data-dependent gradient operator Phi; the same construction can later be applied to an RNN or state-space model used for long-horizon prediction. Second, use the following MVP procedure: sample a task and two datasets differing by one example; run the controller on both; estimate e_k = Phi(y_k,S) - Phi(y_k,S') after aligning the feedback state; estimate an incremental slope bound L from pairs of outputs and updates; choose an IQC such as Q = diag(L^2 I_p, -I_m) for the condition ||w_k|| <= L ||delta_y_k||; construct the lifted matrices T and H; solve an SDP for P, lambda, rho, and the smallest gamma; and add alpha log(gamma + 1e-8) or alpha max(0, lambda_max(M))^2 to the controller objective. Initially stop gradients through the SDP and alternate controller updates with certificate updates. The matrices A, B, C, and D come from the optimizer realization; L and one-sample disturbance magnitudes are estimated empirically, while the SDP is the formal certificate conditional on those estimates. Test on a d = 20 finite-sum quadratic regression problem with N = 100, comparing vanilla gradient descent, heavy-ball, Nesterov, and a learned two-state optimizer. Sweep learning rate and momentum, measure leave-one-sample-out sensitivity, loss divergence, and certified gamma. The prediction is that the empirical divergence or sensitivity boundary appears close to the first hyperparameter value where the LMI becomes infeasible, and cumulative sensitivity grows approximately linearly with gamma, with empirical sensitivity within a factor of two of the certificate.
Verification
Mechanism evidence: Not confirmed
Practical benchmark: Did not beat baseline
Stage 1 — Mechanism check agent confidence 8/10
Built a scalar IQC-certified feedback-dynamics MVP with lifted LMI construction, numerical certification, and reproducible parameter sweeps. The mechanism manifested: the predicted stability boundary L*=1-A=0.8 was observed between L=0.79 (stable) and L=0.81 (unstable), while disturbance sensitivity scaled exactly linearly with epsilon. The certificate was feasible on stable grid points and infeasible on unstable points, although conservative near the boundary.
- Agent confidence
- 8/10
- Baseline
- No-feedback L=0: empirical finite-horizon state gain 1.0206.
- Idea
- Feedback L=0.6: empirical gain 1.6667 and IQC certificate gamma 5.1285, or 1.633x the baseline gain. Predicted boundary L*=0.8; observed transition [0.79,0.81]. At L=0.4, epsilon values [0,.25,.5,1,2] produced state norms [0,.3125,.625,1.25,2.5], with fitted slope 1.25 exactly matching the analytic gain.
Limitations: This is a scalar toy realization, not a learned optimizer or d=20 finite-sum regression meta-training experiment. The certificate uses fixed rho=0.99 and a bounded SciPy SLSQP fallback because CVXPY was unavailable; its gamma is conservative and should be independently checked with a production SDP solver. No generalization-error or wall-clock comparison was tested.
How to run: python3 iqc_experiment.py
Stage 2 — Benchmark vs. tuned baseline
Ahead of the baseline (−20.3%) but not statistically significant (p=0.0522, wins 6 / 8 paired seeds).
- Benchmark
- Dynamics forecasting
- Model
- rnn_small
- Paired seeds
- 8
- Baseline mean
- 0.0004
- Idea mean
- 0.0003
- Effect (Δ)
- -0.0001 (−20.3%; negative = idea better)
- Wins
- 6 / 8 paired seeds
- p-value
- 0.0522 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- ±19.7%
- Mechanism
- Not confirmed ✗
- Practical effect
- Inconclusive
- Baseline tuning
- swept over 3 configs
The gain is an empirical IQC-inspired proxy, not a production SDP certificate for the nonlinear GRU. Only the built-in dynamics track was tested; no cross-domain transfer was tested. The benchmark effect narrowly missed the required significance threshold.
How to run: python3 bench_iqc.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 IQC-inspired gradient clipping as the training-loop intervention on the structurally matched dynamics/rnn_small benchmark. The idea improved mean test MSE from 0.0003780504 to 0.0003013260 and won 6/8 paired seeds, but the permutation p-value was 0.0522, just above the required 0.05 threshold; therefore the required significant-win verdict was not achieved. The trained-model mechanism signature was marked confirmed by its quantitative proxy check, but this does not override the task-metric significance criterion.
- Baseline
- Adam best lr=0.01: mean test MSE 0.0003780504303, std 0.0000731425047.
- Idea
- IQC-inspired clipped Adam, lr=0.01 and clip=1.0: mean test MSE 0.0003013259811, std 0.0000658611425.
- Idea wins
- 6/8
- Paired delta
- -0.0000767244492
- Permutation p value
- 0.0522
Benchmark result
Ahead of the baseline (−20.3%) but not statistically significant (p=0.0522, wins 6 / 8 paired seeds).
- Benchmark
- Dynamics forecasting
- Model
- rnn_small
- Paired seeds
- 8
- Baseline mean
- 0.0004
- Idea mean
- 0.0003
- Effect (Δ)
- -0.0001 (−20.3%; negative = idea better)
- Wins
- 6 / 8 paired seeds
- p-value
- 0.0522 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- ±19.7%
- Mechanism
- Not confirmed ✗
- Practical effect
- Inconclusive
- Baseline tuning
- swept over 3 configs
The gain is an empirical IQC-inspired proxy, not a production SDP certificate for the nonlinear GRU. Only the built-in dynamics track was tested; no cross-domain transfer was tested. The benchmark effect narrowly missed the required significance threshold.
How to run: python3 bench_iqc.py
Artifacts
- bench_iqc.py 5.2 KB View Raw JSON
- bench_report.json 3.3 KB View Raw JSON
- iqc_experiment.py 5.5 KB View Raw JSON
- report.md 1.5 KB View
- report_bench_2026-09-01T234656.md 5.1 KB View
- results.json 3.7 KB View Raw JSON