# Self-Correcting Euler Horizon Rule

- ID: 904
- Canonical URL: https://synthcore.org/idea/904/self-correcting-euler-horizon-rule
- API JSON: https://synthcore.org/api/idea/904.json
- API Markdown: https://synthcore.org/api/idea/904.md
- Verification status: mech_ok_no_baseline
- Source: [arXiv:2607.14291](https://arxiv.org/abs/2607.14291)
- Category: dynamics
- Solves: stability, speedup, accuracy
- ML areas: diffusion-sampling, inference-speedup, ssm, world-model
- Math tags: dynamical-systems, numerical-analysis, control-theory, spectral-theory
- Ratings: usefulness 7/10; difficulty 4/10; novelty 6/10

## Idea description

Use contraction-aware integration rather than assuming that Euler discretization error grows monotonically with sampling time. For a contracting neural ODE, permit a transient error peak but choose the step size and terminal horizon using the predicted peak time and subsequent exponential decay.

## Mathematical statement

Consider the exact contracting flow \(\dot{x}=v(x)\) and Euler iterates \(x_{k+1}=x_k+h v(x_k)\), where \(h\) is the step size and \(\lambda>0\) is a contraction rate. Local truncation errors injected at early steps are propagated by the contraction, so the global error has a forcing-versus-decay profile rather than monotone accumulation. In the linear surrogate \(\dot{x}=-\lambda x\), the exact multiplier is \(e^{-\lambda h}\) and Euler's multiplier is \(1-\lambda h\), requiring \(0<h<2/\lambda\) for discrete stability. For nonlinear contracting fields, the paper predicts a universal transient peak whose time is controlled by \(\lambda\), followed by exponential decay; this creates a testable criterion for safe coarse integration and for distinguishing numerical error from genuine model instability.

## Key formulas

- $$x_{k+1}=x_k+h v(x_k),\qquad \dot{x}=v(x),$$
- $$\dot{x}=-\lambda x:\quad x(t)=e^{-\lambda t}x_0,\qquad x_k=(1-\lambda h)^kx_0,$$
- $$0<h<\frac{2}{\lambda}\quad\Longleftrightarrow\quad |1-\lambda h|<1,$$
- $$e(t)\ \text{is predicted to grow initially, peak at }t_{\mathrm{peak}}=O(\lambda^{-1}),\text{ then decay as }e(t)\propto e^{-\lambda(t-t_{\mathrm{peak}})}.$

## Implementation notes

Add a contraction-aware solver controller to neural ODE, continuous normalizing flow, or diffusion-ODE sampling. Estimate a conservative local contraction rate \(\hat\lambda_t\) from minibatch secant pairs or from the largest eigenvalue of the symmetric Jacobian part, \(\lambda_t=-\lambda_{\max}((J_v+J_v^\top)/2)\). Before each Euler or Heun step, enforce \(h_t\leq\alpha/\max(\hat\lambda_t,\epsilon)\) with \(\alpha<2\), and log the discrepancy between one full step and two half steps as an error proxy. Do not reject a trajectory merely because this proxy increases: maintain a running error curve, identify its peak near \(t_{\mathrm{peak}}\approx c/\hat\lambda\), and permit larger later steps when the proxy begins decaying at the predicted rate. Pseudocode: estimate \(\hat\lambda_t\); clip \(h_t\); compute full and half-step solutions; record \(e_t\); update the peak detector; after the peak, increase \(h_t\) only while \(e_{t+\Delta}/e_t\) is below \(\exp(-\hat\lambda_t\Delta)\) up to tolerance. The paper supplies the self-correction mechanism and contraction dependence; Jacobian rates and constants must be estimated empirically for the neural field. First cheap experiment: use a 2D contracting MLP neural ODE and a CIFAR-10 latent diffusion ODE with fixed Euler versus this controller. Plot error against a high-accuracy RK reference. The prediction is a non-monotone error curve with a peak at a time proportional to \(1/\lambda\), followed by a log-linear decay slope near \(-\lambda\); violating \(h<2/\lambda\) should remove this decay and produce persistent or divergent error.

## Verification

- Status: mech_ok_no_baseline
- Mechanism evidence: yes
- Mechanism confirmed: yes
- Practical verdict: no_effect
- Verdict: Built a reproducible scalar contracting-ODE MVP with Euler error sweeps, stability-boundary verification, peak-time scaling, decay-rate fitting, and a contraction-aware controller. The three mechanism predictions manifested: the stability transition was observed at lambda*h=2.00, peak time scaled as 1/lambda with 9.7% deviation from the analytic estimate, and post-peak log-error slopes matched -lambda within 0.91%. The controller used 33 versus 40 steps but had worse terminal error (3.26e-4 versus 2.03e-4), so this is evidence for the mechanism, not an end-to-end accuracy or speed win.

### Mechanism check

- Verdict: Built a reproducible scalar contracting-ODE MVP with Euler error sweeps, stability-boundary verification, peak-time scaling, decay-rate fitting, and a contraction-aware controller. The three mechanism predictions manifested: the stability transition was observed at lambda*h=2.00, peak time scaled as 1/lambda with 9.7% deviation from the analytic estimate, and post-peak log-error slopes matched -lambda within 0.91%. The controller used 33 versus 40 steps but had worse terminal error (3.26e-4 versus 2.03e-4), so this is evidence for the mechanism, not an end-to-end accuracy or speed win.
- Confidence: 9/10
- Limitations: Only the scalar linear contracting surrogate was tested; no nonlinear MLP neural ODE, CIFAR-10 latent diffusion, Jacobian/secant-rate estimation, wall-clock benchmark, adaptive horizon extension, or GPU experiment was performed.

### Practical benchmark

- Paired seeds: 8
- Baseline mean: 0
- Idea mean: 0
- p-value: 1
- Paired wins: 0/8
- Benchmark verdict: no measurable effect

## Artifacts

- [bench_experiment.py](https://synthcore.org/code/376/bench_experiment.py)
- [bench_report.json](https://synthcore.org/code/376/bench_report.json)
- [experiment.py](https://synthcore.org/code/376/experiment.py)
- [report.md](https://synthcore.org/code/376/report.md)
- [report_bench_2026-08-31T101317.md](https://synthcore.org/code/376/report_bench_2026-08-31T101317.md)
- [results.json](https://synthcore.org/code/376/results.json)
- [Download all files as ZIP](https://synthcore.org/download/376)

## Disclaimer

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