Alphabeta Math
TheoremStatement: AI-adaptedProof: AI-generatedPipeline-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.

Sectional forward Euler converges globally with first-order error

Statement

Fix a finite grid and a final time T>0. Let y(t)=(n(t),L(t)) be the exact finite sectional solution from Finite sectional dynamics have a unique global positive solution, with nonnegative initial data and augmented mass M. Let N≥1 be an integer, set h=T/N, initialize y0=y(0), and apply Forward Euler for a finite sectional balance with a positivity CFL for N steps, with a uniform CFL bound hΛi(nk)≤1. Then there is a constant CT, independent of N, such that

max⁡0≤k≤N∥yk−y(tk)∥2≤CTh,tk=kh.

The constant may depend on the grid, rates, mass, and final time. This is convergence to the specified finite-dimensional sectional ODE; it does not include grid or cutoff error.

Facts & Assumptions

Given: The fixed finite-dimensional autonomous ODE, its nonnegative exact solution, the shared initial state y0=y(0), the fixed-time mesh h=T/N, and the uniform positivity CFL.

[F1]

Finite sectional dynamics have a unique global positive solution: the unique global nonnegative exact solution obeys 0≤ni(t)≤M/xi and 0≤L(t)≤M for initial augmented mass M.

[F2]

Forward Euler for a finite sectional balance with a positivity CFL: each successful update is forward Euler yk+1=yk+hF(yk).

[F3]

The sectional forward-Euler step is positive and conserves augmented mass: each accepted step stays in the same compact nonnegative mass simplex.

[F4]

Forward Euler has uniform local order one and a one-step Lipschitz bound: for a C1 field with the field and its first derivatives bounded on an open cylinder with convex state fibers, forward Euler has local defect O(h2) when the exact path and every exact-start Euler input/output pair lie in a common compact tube; the derivative bound gives the one-step factor 1+Lh for inputs in a common state fiber.

[F5]

Global error from uniform one-step defect and stability: with a common compact tube having positive margin in the open domain and containing the exact path, numerical states, and method stages, a uniform local defect Ch2 and one-step factor 1+Lh imply a fixed-time global error bound CTh.

Proof

1.1givenF1F3F4algebra

Use the finite event polynomials to extend F to all real augmented states, and use the Euler map Ψh(v)=v+hF(v) there. This map agrees with every successful algorithm step; the error analysis may evaluate it at other states without imposing the algorithm's acceptance test. Let BM={(n,L):0≤ni≤M/xi, 0≤L≤M}, a compact box since the grid is finite and each xi>0. The exact path lies in BM by [F1]; starting from y0=y(0), repeated application of [F3] puts every numerical state in the nonnegative augmented-mass simplex of mass M, which is contained in BM. Set RB=max⁡z∈BM∥z∥2 and BF=max⁡z∈BM∥F(z)∥2; these are finite because F is polynomial. Since h=T/N≤T, every exact-start Euler output obeys ∥y(tk)+hF(y(tk))∥2≤RB+TBF. Thus the exact path, all successful numerical states and Euler stage states, and the exact-start Euler inputs and outputs all lie in the compact tube K=[0,T]×B‾(0,RB+TBF+1). This tube lies with positive margin in the bounded open cylinder V=(−1,T+1)×B(0,RB+TBF+2), whose state fibers are convex. The time-independent finite polynomial field f(t,z)=F(z) and its first derivatives are bounded on V. Hence [F4] applies uniformly to every local defect, including the exact-start Euler output, and gives defect Ch2; its derivative bound also gives the one-step factor 1+Lh for inputs in K.

2.1givenF2F4F5step 1.1∎

Apply [F5] with initial error zero, the fixed mesh h=T/N, and the common compact tube from step 1.1. It gives max⁡k≤N∥yk−y(tk)∥2≤CTh, as claimed.

Depends on

Used by

Dependency tree · two levels

24 results within two dependency steps of this one, each drawn at its shortest distance from it. An arrow runs from a result to what uses it, so the chart reads left to right and ends at this result, which carries a heavier outline. Every node is a link to that result. Click elsewhere on the chart to enlarge it.

Sources