Alphabeta Math
ExampleConstruction: AI-generatedVerification: 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.

A scalar inexact Newton solve trades inner work for outer error

Example

Let f(x)=x2/2, so F(x)=f′(x)=x and F′(x)=1. At a nonzero iterate, solve the Newton equation approximately from p(0)=0 by

p(j+1)=p(j)−(1−ω)(p(j)+x),0<ω<1.

Stop at the first j with ωj≤η, where the forcing tolerance satisfies 0<η≤c2<1/2. With the exact strong Wolfe search from Finite-budget bracketing and zoom for a strong Wolfe step, the unit step is accepted. The number j of inner iterations grows as η is tightened, while the outer error contracts by a factor at most η; choosing ηk=∣xk∣ gives a quadratic outer error bound for 0<∣x0∣≤c2.

For each finite horizon N≥1, let jk be the least integer with ωjk≤ηk on the first N steps of the displayed recurrence, and choose the single inner budget MN≥max⁡0≤k<Njk. Use M=MN, ηˉ=c2, σ=1/2, J=0, E,Z≥1, γ=2, and λ0=λmax⁡=1. The deterministic solver on the scalar operator A=1 starts at zero and performs the displayed iteration until ∣p+x∣≤τ, where τ is its absolute requested residual tolerance, or until M iterations are used; it returns that scalar candidate and its iteration count. It returns no-candidate for other operators. For b=−x≠0 and any m≥1, Km(1,b)=R, so every returned candidate meets the Krylov input contract (and p(0)=0∈K0). This realizes every finite prefix. The infinite recurrence and its rate claim do not assert realization by one fixed finite inner budget.

Facts & Assumptions

Given: The scalar objective, the fixed ω∈(0,1), and Wolfe parameters 0<c1<c2<1/2.

[F1]

The Jacobian matrix of partial derivatives and the gradient in the scalar-valued case: for a scalar function the gradient is the vector of first partial derivatives.

[F2]

The Hessian matrix and critical points of a scalar field: the Hessian is the matrix of second partial derivatives.

[F3]

Finite-budget inexact Newton–Krylov with a descent safeguard: the unshifted solve is accepted only after recomputing H(x)p+g(x) and checking its forcing residual and descent threshold; the next iterate uses the returned Wolfe step.

[F4]

Finite-budget bracketing and zoom for a strong Wolfe step: the search tests t=1 first and returns it when both strong Wolfe inequalities hold.

[F5]

Forcing residuals determine local inexact-Newton rates: an unshifted full-step sequence with a relative residual bounded by ηk∥F(xk)∥2 converges quadratically when ηk=O(∥F(xk)∥2), under the stated local hypotheses.

[F7]

The Krylov subspace Km(A,b)=span⁡{b,Ab,…,Am−1b}: K0(A,b)={0} and Km(A,b)=span⁡{b,Ab,…,Am−1b} for m≥1.

Verification

Given: The objective and inner iteration of the Example, with x≠0 and 0<η≤c2.

1.1givenF1F2inductionalgebra

The exact scalar identities are F(x)=x, J(x)=1, and g(x)=x by [F1]–[F2]. Induction on j gives p(j)=−(1−ωj)x, so the linear residual is rj=J(x)p(j)+F(x)=ωjx.

1.2givenF6F7choosealgebra

Since 0<ω<1, [F6] gives ωj→0, so a least j with ωj≤η exists. As η<1, j≥1 and minimality gives ωj≤η<ωj−1. If the requested tolerance is reduced, the set of eligible j can only shrink, so the required inner iteration count cannot decrease. For any finite N, all jk on its first N nonzero iterates are finite, so their finite maximum defines MN. The specified solver therefore reaches its requested residual tolerance before its budget on that prefix, and [F7] verifies its scalar Krylov membership.

2.1givenF3step 1.1step 1.2algebra

At this j, ∣rj∣=ωj∣x∣≤η∣x∣, while g(x)p(j)=−(1−ωj)x2 and ∣g(x)∣ ∣p(j)∣=(1−ωj)x2. Because 0<σ<1, the descent safeguard passes.

3.1givenF3F4step 2.1algebra

Let q:=ωj≤η≤c2. The unit step gives x+=x+p(j)=qx. Its objective satisfies f(x+)=q2x2/2≤x2(1−2c1(1−q))/2=f(x)+c1g(x)p(j), because 1+q≥1>2c1. Also ∣ϕ′(1)∣=q(1−q)x2≤c2(1−q)x2=c2∣ϕ′(0)∣. Thus both strong Wolfe conditions hold, and [F4] accepts t=1.

4.1givenF5step 1.2step 3.1induction∎

For fixed η, q=ωj is a constant in (0,1) and xk+1=qxk, so the outer sequence is linearly convergent. If 0<∣x0∣≤c2 and ηk=∣xk∣, then ηk≤c2 at every step and ∣xk+1∣=ωjk∣xk∣≤∣xk∣2. Here ηk=∣F(xk)∣; the local tube in [F5] may contain [−c2,c2] because F(x)=x and J(x)=1 everywhere. Thus [F5] supplies the same local quadratic-rate conclusion.

Depends on

Used by

Nothing in the library uses this result yet.

Dependency tree · two levels

39 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.