△ Mechanism confirmed, baseline not beaten 2026

Compactified Burst Controller

Implementation & benchmark of arXiv:2607.14269 — Blow-up Parameter Landscapes for Polynomial Dynamical Systems

Usefulness7/10
Difficulty5/10
Novelty7/10

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

dr/dt = r^n a(u) + O(r^(n-1)), and du/dtau = G(u), with d tau = r^(n-1) dt.
G(u_star) = 0, alpha_star = a(u_star) > 0, and max Re eig(P D G(u_star) P restricted to the tangent space) > 0.
F_control(x) = F(x) - kappa(r)u, where u = x/r and kappa(r) = k softplus((r - r_c)/delta).
d||x||/dt is approximately r^n a(u) - kappa(r).

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

Mechanism confirmed, baseline not beaten

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
Limitations:

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
Limitations:

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

Implementation overview ⬇ Download all as ZIP 6 files · code, reports and structured results