Certified Fold Map for Recurrent Fixed Points
Implementation & benchmark of arXiv:2608.07119 — Certified Detection of Bifurcation Candidates in Uncertain Nonlinear Systems using Interval Analysis
Source paper: Certified Detection of Bifurcation Candidates in Uncertain Nonlinear Systems using Interval Analysis arXiv:2608.07119 ⓘ · analyzed Aug 31, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Apply interval Krawczyk certification to the augmented equations for a recurrent-network fixed point and a singular state Jacobian. This produces a rigorous local certificate for the gain or feedback value at which two fixed points merge or disappear, allowing training or inference to avoid parameter boxes containing an uncertified fold.
Formulas
Mathematical statement
Let a recurrent update be x_(t+1)=Phi(x_t,mu), with state x in R^n and scalar parameter mu. Define F(x,mu)=Phi(x,mu)-x and J(x,mu)=D_x F(x,mu). A saddle-node candidate is a zero of g(z)=[F(x,mu); det(J(x,mu))], where z=(x,mu) is in R^(n+1). For a search box [z], let z_c=midpoint([z]), let C approximate the inverse of Dg(z_c), and let [Dg]([z]) enclose the Jacobian throughout the box. The Krawczyk set is K([z])=z_c-Cg(z_c)+(I-C[Dg]([z]))([z]-z_c). Strict inclusion K([z]) subset int([z]) certifies a unique zero; disjointness K([z]) intersect [z]=emptyset excludes the box. The contraction quantity q=norm(I-C[Dg]([z])) predicts whether certification succeeds.
Implementation notes
Integration point: use this as an offline certification layer for an RNN, implicit layer, or state-space model, with weights fixed and mu chosen as recurrent gain, residual scale, or feedback strength. For a small tanh RNN, set Phi(x,mu)=tanh(mu W x+b), evaluate F and J with interval arithmetic, and use automatic differentiation at box midpoints. For larger systems replace det(J)=0 with Jv=0 and c^T v-1=0, introducing a null vector v and fixed normalization vector c. Pseudocode: tile a box over x and mu; compute g(z_c), C=inverse(Dg(z_c)), and interval [Dg]([z]); form K; accept strict inclusion; discard disjoint boxes; otherwise bisect the widest coordinate and repeat. The paper contributes the Krawczyk construction, outward-rounded interval evaluation, and inclusion logic. The engineer estimates activation interval bounds and chooses search domains and stopping widths. First experiment: use a two-dimensional tanh RNN, compare certification against dense fixed-point continuation and Newton solves. Prediction: accepted boxes occur only when q=norm(I-C[Dg]) is below one and the Krawczyk image has positive interior margin; the certified mu interval must contain the continuation fold, while q decreases as box diameter shrinks near a regular augmented root.
Verification
Stage 1 · Toy mechanism gate: Failed ✗
Stage 2 · Mechanism transferred to benchmark: Not tested
Stage 2 · Practical benchmark result: Not run
Methodology: Toy-system gate first; the benchmark stage runs only after a pass. How verification works
Stage 1 — Mechanism check agent confidence 9/10
Built a 2D tanh-RNN augmented fold solver and Krawczyk sweep. The mechanism appeared numerically: q fell from 1.816 at box width 0.2 to 0.0407 at width 0.00625, strict inclusion began at width 0.1, and narrow boxes offset from the root were excluded. However, the implementation uses finite differences and corner sampling rather than outward-rounded interval arithmetic, so it is not a rigorous certificate; the continuation baseline also failed to localize the fold and reported mu=1.675 versus the augmented root mu=0.7434.
- Agent confidence
- 9/10
- Baseline
- Warm-start continuation reported closest_grid_mu=1.675 with absolute error 0.9316, and therefore did not reliably track the fold branch.
- Idea
- Augmented root [-0.57223, 0.15708, 0.74344], residual 1.18e-12; q(width 0.2 to 0.00625)=1.816 to 0.0407; first strict inclusion width=0.1; all tested offsets 0.02 to 0.40 were disjoint at radius 0.008.
Limitations: No outward-rounded interval library was used, so strict inclusion is only an empirical surrogate. No larger RNN, training/inference integration, certified box tiling/bisection, GPU experiment, or reliable branch-aware continuation comparison was tested.
How to run: python3 fold_map_experiment.py > toy_results.json
Verdict computed by deterministic test code from paired-seed statistics — not by the language model.
Artifacts
- fold_map_experiment.py 2.7 KB View Raw
- report.md 1.5 KB View
- toy_results.json 2.7 KB View Raw