Alphabeta Math
Pipeline-generated
How statement and proof provenance work

The first chip identifies the source of the statement or construction; the second identifies the source of its local proof or verification.

  • Literature-sourced: the exact statement appears in a cited source; only wording and notation differ.
  • AI-adapted: a semantically identical restatement of literature-sourced material, modulo indexing, notation, and boundary cases adopted by the library.
  • AI-generated: a genuinely novel statement formulated by AI, with no source for the claim itself.

These labels describe origin, not correctness: citations and verification chips remain separate evidence.

✓ 3 results · all verified · 3 also independently AI-judged
Every result on this page is machine-checked by a proof checker and read in full and owner-audited; the judge is an additional, independent cross-model AI review of the proofs; all 3 also cleared it.

Adaptive ODE solvers and error estimation — Examples

1 · Prerequisites

2 · Summary

Summary

These examples isolate three limits of an embedded controller: the leading estimator term can vanish, a fast-decay step can be rejected down to the minimum-step safeguard, and a rounded implementation can accept while its exact endpoint error exceeds the requested tolerance. They illustrate the stated method and do not supply global adaptive accuracy or a general stiffness diagnosis.

3 · Logical flowchart

4 · Definitions, theorems and proofs

None yet.

5 · Examples, counterexamples and false statements

ExampleConstruction: AI-generatedVerification: AI-generatedprecheck passjudge pass (gpt-6-luna)Open item page →

A tolerance below roundoff can pass with a zero embedded difference

Example

Use the toy binary format whose nonzero representable numbers are ±m2e with integer e and m∈{4,5,6,7} (three significant bits), rounded to nearest with ties to an even significand. The next larger representable value after 8=4⋅2 is 10=5⋅2, so their midpoint 9 rounds to 8 because 4 is even and 5 is odd. Apply a rounded Euler–Heun implementation to y′=1, y(0)=8 on [0,1], with one step h=1, hmin⁡=hmax⁡=hinit=1, attempt budget N=1, atol=1/4, rtol=0, safety 0.9, and scale bounds 0.2 and 5. Both computed endpoints round to 8; the stored difference is zero and the step is accepted, although the exact endpoint is 9 and the actual error is 1.

This example concerns the rounded realization described in Finite-budget adaptive Runge–Kutta controller. It does not contradict its exact-arithmetic specification and shows why that acceptance test is not a certificate.

Facts & Assumptions

Given: The toy format and rounding rule in the Example, field f(y)=1, initial state 8, and controller tolerance parameters.

[F1]

Euler–Heun embedded Runge–Kutta pair: both pair members use K1=f(yn) and K2=f(yn+hK1), with Euler and Heun weighted endpoints.

[F2]

Finite-budget adaptive Runge–Kutta controller: a rounded realization uses stored high and low outputs in the difference and accepts when the scaled difference Q≤1.

[F3]

Local orders and leading estimator term for the Euler–Heun pair: for this pair the low member is order p=1, which is the controller exponent input.

Verification

1.1givenF1F2F3algebra

The two stages evaluate to K1=K2=1; each exact endpoint formula is 8+1=9. By the specified tie rule, the Euler endpoint rounds to 8 and the Heun endpoint also rounds to 8. Hence the stored difference is 0 and the scaled error is Q=0/(1/4)=0≤1, so the rounded controller accepts.

2.1givenstep 1.1algebra∎

The exact solution is y(t)=8+t, so y(1)=9. The accepted stored endpoint is 8, giving absolute error ∣9−8∣=1>1/4=atol. Thus the embedded difference can be zero while the actual error exceeds the requested tolerance.

ExampleConstruction: AI-generatedVerification: AI-generatedprecheck passjudge pass (gpt-6-luna)Open item page →

A fast-decay test problem is rejected at the minimum step

Example

Run the Euler–Heun controller on the scalar fast-decay equation y′=−1000y, y(0)=1, with T=0.004, hmin⁡=0.002, hmax⁡=h0=0.004, atol=10−4, rtol=0, and attempt budget N=2, p=1, safety σ=0.9, and scale factors γmin⁡=0.2, γmax⁡=5. The computed pair difference rejects both attempted steps, so the procedure stops at its minimum-step safeguard. This scalar fast-decay equation is used as a stiffness illustration: the calculation verifies repeated rejection of this explicit trial only, and gives no general stiffness diagnosis.

Facts & Assumptions

Given: The scalar field f(y)=−1000y, starting value 1, trial sizes 0.004 and 0.002, and controller parameters in the Example.

[F1]

Euler–Heun embedded Runge–Kutta pair: K1=f(yn), K2=f(yn+hK1), YE=yn+hK1, and YH=yn+h2(K1+K2).

[F2]

Finite-budget adaptive Runge–Kutta controller: reject if Q>1, scale by max⁡(γmin⁡,min⁡(γmax⁡,σQ−1/(p+1))), clamp the next proposal to the step bounds, and return minimum-step-rejected for a rejection at hmin⁡.

[F3]

The bounded adaptive Runge–Kutta controller terminates: this finite-budget controller reaches a terminal status after finitely many attempts.

[F4]

Local orders and leading estimator term for the Euler–Heun pair: the Euler–Heun pair has low order p=1.

Verification

1.1givenF1F2F4algebra

At h=0.004, K1=−1000, the Euler stage value is −3, and K2=3000. Thus YE=−3, YH=5, the difference magnitude is 8, and Q=8/10−4=80000>1. The proposed scale is clamped to 0.2, so the next step is max⁡(0.002,0.2⋅0.004)=0.002.

2.1F1F2F3F4step 1.1algebra∎

At h=0.002, K1=−1000, the Euler stage value is −1, and K2=1000. Thus YE=−1, YH=1, the difference magnitude is 2, and Q=20000>1. This rejection occurs at hmin⁡, so [F2] returns minimum-step-rejected; finite termination is also covered by [F3].

ExampleConstruction: AI-generatedVerification: AI-generatedprecheck passjudge pass (gpt-6-luna)Open item page →

The leading Euler–Heun estimator coefficient can vanish

Example

Apply the Euler–Heun pair to y′=t2, y(0)=0, on one step of size h>0. The low-order defect begins at order h3, not h2, so the embedded difference is not asymptotically equal to the low-order defect in ratio.

Facts & Assumptions

Given: The field f(t,y)=t2, initial time and value t0=y0=0, and h>0.

[F1]

Euler–Heun embedded Runge–Kutta pair: K1=f(0,0), K2=f(h,hK1), YE=hK1, and YH=h2(K1+K2).

[F2]

Local orders and leading estimator term for the Euler–Heun pair: the local error expansions compare each one-step value with the exact endpoint.

[F3]

A nonzero leading defect makes the embedded difference asymptotically match the low-order defect: a ratio-one conclusion requires a nonzero leading low-defect coefficient.

Verification

1.1givenF1F2algebra

Here K1=0 and K2=h2, so YE=0 and YH=h3/2. Since y(h)=∫0ht2dt=h3/3, the Euler defect is −h3/3, the Heun defect is h3/6, and the embedded estimator is ∣YH−YE∣=h3/2.

2.1F1F3step 1.1algebra∎

The Euler defect has zero coefficient at order h2 and the estimator-to-Euler-defect ratio is (h3/2)/(h3/3)=3/2 for every h>0. This concrete case shows why [F3] requires a nonzero leading coefficient before asserting ratio one.

Sources