△ Mechanism confirmed, baseline not beaten 2026

Topological Reachable-Set Coverage Scheduler

Implementation & benchmark of arXiv:2608.10332 — Topological Feasibility Guarantees for Differentiable Predictive Control

Usefulness7/10
Difficulty6/10
Novelty7/10

Source paper: Topological Feasibility Guarantees for Differentiable Predictive Control arXiv:2608.10332 ⓘ · analyzed Sep 1, 2026

AI-generated research hypothesis, automatically tested. Not peer-reviewed.

Idea description

Use the reachable-safe-set viewpoint to make training data generation adaptive: maintain an approximation of the states reached by the current neural policy, identify boundary regions with weak barrier margin, and sample there until the set is sufficiently covered. This replaces random rollout expansion with a measurable coverage condition that can support finite-sample safety claims.

Formulas

R_theta = {x(t; x_0, pi_theta) : x_0 in X_0 and 0 <= t <= T}, C = {x : h(x) >= 0}
delta(S, R_theta) = sup over x in R_theta of min over x_i in S of ||x - x_i||
min over x_i in S of r_theta(x_i) >= epsilon and r_theta is L_r-Lipschitz implies inf over x in R_theta of r_theta(x) >= epsilon - L_r delta(S, R_theta)
delta(S, R_theta) < epsilon / L_r implies r_theta(x) > 0 for every x in R_theta

Mathematical statement

Let R_theta denote the set of states reachable under policy pi_theta from the initial-state set over the training horizon. Let r_theta(x) be the CBF residual and suppose it is L_r-Lipschitz on R_theta. A finite sample set S is a delta-net of R_theta if every x in R_theta has a sample x_i in S with distance at most delta. If all samples have residual at least epsilon, then every reachable state has residual at least epsilon - L_r delta. Define the empirical covering radius as delta_hat = max over validation states x of the distance from x to its nearest training sample. Training or certification stops only when delta_hat is below epsilon divided by the estimated L_r.

Implementation notes

1. Integration point: attach the scheduler to offline policy training or differentiable model-predictive rollouts. Store state-action trajectories, barrier residuals, and a k-d tree over training states. 2. Pseudocode: initialize a seed set; train the policy for several epochs; generate a large validation cloud of reachable states; for each validation state compute nearest-training-sample distance and barrier residual; select states with high distance, low residual, or both; add them to the training set; retrain; stop only when delta_hat is below epsilon divided by estimated L_r. A practical priority score is q(x) = distance_to_nearest_sample + gamma times max(0, epsilon - r_theta(x)). 3. The theory supplies the margin relation and threshold epsilon/L_r. Estimate the reachable set, coverage radius, and L_r empirically, and report conservative upper bounds. For high-dimensional states, use a latent metric only when a decoder reconstruction-error bound is available. 4. First cheap experiment: use a two-dimensional pendulum or planar robot with a known safe region. Compare random sampling, uniform initial-state sampling, and adaptive sampling at equal simulator-call budgets. Measure maximum dense-grid violation, estimated coverage radius, minimum sampled residual, and training-set size. The predicted signature is monotonic reduction in violation as delta_hat decreases, with a transition near delta_hat = epsilon / L_r. Adaptive sampling should reach zero violation with fewer trajectories than random sampling.

Verification

Mechanism confirmed, baseline not beaten

Stage 1 · Toy mechanism gate: Passed ✓

Stage 2 · Mechanism transferred to benchmark: Confirmed ✓

Stage 2 · Practical benchmark result: Inconclusive

Methodology: 8 paired seeds, permutation test (p<0.05) against a learning-rate-tuned baseline. How verification works

Stage 1 — Mechanism check agent confidence 8/10

Built a readable 2D reachable-set coverage scheduler with random and adaptive/farthest sampling, exact KD-tree coverage measurement, and reproducible sweeps. The Lipschitz safety bound held in 80/80 checks, residual-gap scaling was linear in L with constant gap/L=0.141421, and adaptive sampling reduced median coverage radius and reached certification substantially sooner for feasible thresholds. This validates the mechanism in the controlled toy setting, but not yet in a learned neural-policy or physical rollout setting.

Agent confidence
8/10
Baseline
Random sampling: median coverage radius at 12/25/50/80 samples = 0.43013/0.32473/0.22326/0.17309; certification median samples = 93 at L=0.25 and 266.5 at L=0.5 among successful runs.
Idea
Adaptive sampling: median coverage radius at 12/25/50/80 samples = 0.32998/0.18856/0.10541/0.10000; certification median samples = 36 at L=0.25 and 88 at L=0.5. Bound validity was 100% (80/80), and median residual gap/L stayed 0.141421 across L=0.25,0.5,1,2,4.

Limitations: The experiment uses a known static 2D unit-square reachable set and an analytically Lipschitz residual on a finite grid. It does not test a learned policy, simulator-generated trajectories, estimated reachable-set error, noisy residuals, high-dimensional metrics, latent reconstruction bounds, or confidence intervals; the low-margin score was inactive because all residuals exceeded epsilon.

How to run: python3 coverage_scheduler_experiment.py

Stage 2 — Benchmark vs. tuned baseline (arXiv:2608.10332)

Ahead of the baseline (−13.5%) but not statistically significant (p=0.5151, wins 4 / 8 paired seeds).

Benchmark
Dynamics forecasting
Model
rnn_small
Paired seeds
8
Baseline mean
0.0022
Idea mean
0.0019
Effect (Δ)
-0.0003 (−13.5%; negative = idea better)
Wins
4 / 8 paired seeds
p-value
0.5151 (permutation test, 20 000 shuffles)
Smallest detectable effect
±46.5%
Mechanism
Confirmed ✓
Practical effect
Inconclusive
Baseline tuning
swept over 3 configs
Limitations:

Only the built-in dynamics task was tested; the scheduler used an empirical state-space pool and simple residual proxy, without learned reachable-set confidence bounds, physical safety certification, larger models, longer horizons, or more than 8 paired seeds.

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 the Topological Reachable-Set Coverage Scheduler on the structurally matched dynamics track with the standard rnn_small model and canonical train_model loop. Adaptive sampling reduced empirical reachable-state coverage radius by 54.70% at equal simulator budget, confirming the mechanism signature, and the Lipschitz sanity check held. However, the 8-seed test-MSE improvement was not statistically significant (paired delta=-0.00030233, permutation p=0.51515), so the idea did not demonstrate a significant benchmark win.

Baseline
Best lr=0.006; 8-seed test MSE mean=0.00224437, std=0.00068261.
Idea
Adaptive scheduler, lr=0.006; 8-seed test MSE mean=0.00194204, std=0.00059465; paired delta=-0.00030233, p=0.51515; coverage radius 0.34452 versus 0.76052 baseline, 54.70% lower.

Benchmark result

Ahead of the baseline (−13.5%) but not statistically significant (p=0.5151, wins 4 / 8 paired seeds).

Benchmark
Dynamics forecasting
Model
rnn_small
Paired seeds
8
Baseline mean
0.0022
Idea mean
0.0019
Effect (Δ)
-0.0003 (−13.5%; negative = idea better)
Wins
4 / 8 paired seeds
p-value
0.5151 (permutation test, 20 000 shuffles)
Smallest detectable effect
±46.5%
Mechanism
Confirmed ✓
Practical effect
Inconclusive
Baseline tuning
swept over 3 configs
Limitations:

Only the built-in dynamics task was tested; the scheduler used an empirical state-space pool and simple residual proxy, without learned reachable-set confidence bounds, physical safety certification, larger models, longer horizons, or more than 8 paired seeds.

How to run: python3 stage2_bench.py

Artifacts

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