Compactified Burst Controller
Implementation & benchmark of arXiv:2607.14269 — Blow-up Parameter Landscapes for Polynomial Dynamical Systems
Source paper: Blow-up Parameter Landscapes for Polynomial Dynamical Systems arXiv:2607.14269 ⓘ · analyzed Aug 30, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Use the paper's distinction between radial attraction and tangential instability at infinity to detect impending hidden-state bursts before they cause numerical failure. When the state approaches a radially growing directional equilibrium, temporarily add radial damping or switch to a bounded fallback update, then restore the original dynamics after angular ejection.
Formulas
Mathematical statement
For x = r u and a homogeneous leading field F_n, define a(u) = u transpose F_n(u) and G(u) = F_n(u) - a(u)u. A candidate equilibrium at infinity satisfies G(u_star) = 0. Positive a(u_star) causes radial growth because dr/dt is approximately r^n a(u_star). Tangential instability is diagnosed from the projected Jacobian J_T = P D G(u_star) P on the tangent space, where P = I - u_star u_star transpose. If the largest real part of an eigenvalue of J_T is positive, the direction is angularly unstable: the state can grow while nearly aligned with the direction and later be ejected, producing a large but transient burst. A radial controller subtracts kappa(r)u from the vector field, reducing radial growth without directly suppressing angular motion.
Implementation notes
1. Integration point: insert the controller inside every neural ODE function evaluation or recurrent hidden-state update, before advancing the solver. It is intended for long-horizon RNN and SSM rollouts, DEQ iterations, and polynomial residual blocks where rare norm bursts cause NaNs or excessive adaptive solver steps. 2. Pseudocode: compute r = norm(h) plus epsilon and u = h/r; evaluate the estimated leading field v = F_n(u); compute a = dot(u,v) and q = norm(v - a u). Maintain an exponential moving average of a. If r exceeds r_c, the averaged radial growth exceeds a margin, and q is below an angular-equilibrium threshold, replace dh/dt = F(h) by dh/dt = F(h) - kappa(r)u. Estimate angular instability every K steps using finite differences or automatic differentiation of G. 3. The paper supplies the radial and tangential infinity mechanism; r_c, the growth margin, q threshold, damping gain k, and transition width delta are calibrated on non-divergent training rollouts. 4. First cheap experiment: compare vanilla updates, global gradient clipping, and this controller on a cubic RNN trained on the copy and addition tasks, evaluating at ten times the training horizon. Sweep the cubic gain lambda. Prediction: uncontrolled bursts begin when averaged a becomes positive while q is initially small; controller interventions occur before the norm maximum, reduce solver failures, and preserve trajectories with negative a. The measured critical lambda for burst onset should agree with the infinity-atlas prediction within 20 percent.
Verification
Stage 1 · Toy mechanism gate: Passed ✓
Stage 2 · Mechanism transferred to benchmark: Confirmed ✓
Stage 2 · Practical benchmark result: Significantly worse ✗
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 readable 2D homogeneous cubic burst simulator with infinity-atlas diagnostics and gated softplus radial damping. The core math manifested: finite-difference a(u*) and tangent eigenvalue matched predictions to 1e-12, radial threshold-time scaling matched the cubic blow-up law within 8e-5 relative error, and angular slopes matched beta within 1.25%, including the stable negative-beta case. The predicted finite-horizon onset boundary alpha=0.8203 was observed between 0.8 and 0.9; in burst regimes the controller prevented numerical failure and held max radius near 2.2-3.0 versus vanilla failure or large excursions.
- Agent confidence
- 9/10
- Baseline
- Vanilla cubic field: alpha=1.2, 1.6, 2.0 all failed before T=0.8 with recorded max radii 2.96e16, 1.84e18, and 5.69e9; alpha=0.8 remained bounded at max r=1.881.
- Idea
- Controlled field: alpha=1.2, 1.6, 2.0 did not fail, with max radii 2.209, 2.411, and 3.009; interventions began at t=0.547, 0.410, and 0.328 respectively. At alpha=0.8 the controller did not intervene and matched vanilla exactly.
Limitations: This is a deterministic 2D toy cubic vector field, not a trained RNN, copy/addition task, adaptive ODE solver, or FLOP/speed benchmark. The controller uses an instantaneous angular-equilibrium gate rather than implementing EMA growth estimates or periodic Jacobian estimation, and the burst prevention result is not evidence of end-to-end neural-network training gains.
How to run: python3 burst_controller_experiment.py
Stage 2 — Benchmark vs. tuned baseline (arXiv:2607.14269)
Significantly WORSE than the tuned baseline (+176.9%, p=0.0147, wins 1 / 8 paired seeds).
- Benchmark
- Dynamics forecasting
- Model
- rnn_small
- Paired seeds
- 8
- Baseline mean
- 0.0015
- Idea mean
- 0.0042
- Effect (Δ)
- 0.0027 (+176.9%; negative = idea better)
- Wins
- 1 / 8 paired seeds
- p-value
- 0.0147 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- ±116.3%
- Mechanism
- Confirmed ✓
- Practical effect
- Hurts
- Baseline tuning
- swept over 6 configs
Only the registered dynamics track was tested. Copy/addition tasks, longer-horizon rollout, adaptive ODE solver failures, inference speed or FLOP savings, EMA growth estimation, periodic Jacobian estimation, and broader controller hyperparameter sweeps were not tested.
How to run: python3 compactified_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
The Compactified Burst Controller was evaluated on the registered dynamics track with the matched rnn_small task and eight paired seeds. Its trained-model mechanism signature was confirmed: mean maximum hidden norm decreased from 2.679 to 2.058 with 1919 interventions. Nevertheless, test MSE significantly worsened from 0.001530 to 0.004237, with paired delta +0.002706 and permutation p=0.01475; therefore the system did not work on this benchmark.
- Baseline
- Best config lr=0.009, epochs=12, weight_decay=0.0001; full 8-seed MSE mean 0.0015301889.
- Idea
- Best tested config lr=0.009, epochs=12, weight_decay=0.0001; full 8-seed MSE mean 0.0042366631; paired delta +0.0027064742, p=0.01475.
- Mechanism
- Mean maximum hidden norm reduced by 0.6204301, from 2.6786648 to 2.0582347; 1919 controller interventions; confirmed=true.
Benchmark result
Significantly WORSE than the tuned baseline (+176.9%, p=0.0147, wins 1 / 8 paired seeds).
- Benchmark
- Dynamics forecasting
- Model
- rnn_small
- Paired seeds
- 8
- Baseline mean
- 0.0015
- Idea mean
- 0.0042
- Effect (Δ)
- 0.0027 (+176.9%; negative = idea better)
- Wins
- 1 / 8 paired seeds
- p-value
- 0.0147 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- ±116.3%
- Mechanism
- Confirmed ✓
- Practical effect
- Hurts
- Baseline tuning
- swept over 6 configs
Only the registered dynamics track was tested. Copy/addition tasks, longer-horizon rollout, adaptive ODE solver failures, inference speed or FLOP savings, EMA growth estimation, periodic Jacobian estimation, and broader controller hyperparameter sweeps were not tested.
How to run: python3 compactified_bench.py
Artifacts
- bench_report.json 4.4 KB View Raw
- burst_controller_experiment.py 5.0 KB View Raw
- compactified_bench.py 5.3 KB View Raw
- report.md 1.7 KB View
- report_bench_2026-08-31T095644.md 4.7 KB View
- results.json 3.4 KB View Raw