Alphabeta Math
TheoremStatement: Literature-sourcedProof: AI-adaptedPipeline-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.

Gauss-Newton converges quadratically at a full-rank zero-residual root

Statement

Let m,n≥1, and let r:U→Rm be C1 on an open U⊆Rn, and let x∗∈U satisfy r(x∗)=0. Let B={x:∥x−x∗∥2≤R}⊆U be convex. Suppose J=Dr is L-Lipschitz on B and, for some σ>0, ∥J(x)v∥2≥σ∥v∥2(x∈B, v∈Rn). Thus every J(x) has full column rank. For the exact Gauss-Newton iteration xk+1=xk+pk, where pk is the unique minimizer of ∥r(xk)+J(xk)p∥22. When the iteration uses the exact branch of Gauss-Newton linearized least-squares step with QR/SVD rank safeguards, assume each invocation is supplied with an exact reduced QR factorization of J(xk). Set C=L/(2σ). Then every xk∈B satisfies ∥xk+1−x∗∥2≤C∥xk−x∗∥22. If C>0 and ∥x0−x∗∥2≤min⁡{R,1/(2C)}, all iterates stay in the ball, converge to x∗, and contract at least by 1/2 per iteration after the quadratic estimate. If C=0 and x0∈B, one exact step reaches x∗. This result concerns exact solves; it asserts no finite-precision convergence rate.

Facts & Assumptions

Given: The residual, root, convex ball, Lipschitz Jacobian, uniform lower bound σ, and exact Gauss-Newton steps in the Statement.

[F1]

Rectangular residual linearization error under a Lipschitz Jacobian: on a convex segment, ∥r(x)−r(x∗)−J(x)(x−x∗)∥2≤(L/2)∥x−x∗∥22.

[F2]

Gauss-Newton linearized least-squares step with QR/SVD rank safeguards: an exact full-rank branch returns the unique minimizer of the linearized residual least-squares problem.

[F3]

The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn: the Euclidean inner product is bilinear and symmetric, and ∥v∥22=vTv.

[A1]

r(x∗)=0.

[A2]

∥J(x)v∥2≥σ∥v∥2 with σ>0 for every x∈B and every v.

Proof

1.1givenA1F1construct

Fix x∈B, put e=x−x∗, and write E=r(x)−r(x∗)−J(x)e. By [F1] and [A1], r(x)=J(x)e+E with ∥E∥2≤(L/2)∥e∥22.

2.1givenA2F2F3step 1.1algebra

Let p be the exact Gauss-Newton step and set h:=p+e. For every v∈Rn, minimality of p gives 12∥r(x)+J(x)(p+tv)∥22≥12∥r(x)+J(x)p∥22 for all real t. Expanding in t shows the linear coefficient vTJ(x)T(r(x)+J(x)p) must be zero (otherwise a sufficiently small t of the opposite sign decreases the quadratic). Since this holds for every v, J(x)T(J(x)h+E)=0 by step 1.1.

3.1step 1.1step 2.1A2F3algebra

Set q:=−h. Step 2.1 gives J(x)T(E−J(x)q)=0, so E−J(x)q is orthogonal to J(x)q; hence ∥E∥22=∥J(x)q∥22+∥E−J(x)q∥22≥∥J(x)q∥22. By [A2], ∥J(x)q∥2≥σ∥q∥2. Therefore ∥x+p−x∗∥2=∥h∥2=∥q∥2≤∥E∥2/σ≤C∥e∥22.

4.1step 3.1algebra∎

If C=0 and x0∈B, step 3.1 gives x1=x∗. If C>0 and ∥x0−x∗∥2≤1/(2C), the stated hypothesis also puts x0∈B, and step 3.1 gives ∥x1−x∗∥2≤12∥x0−x∗∥2, so x1∈B; induction repeats the same bound for every k, and ∥xk−x∗∥2≤2−k∥x0−x∗∥2→0.

Depends on

Used by

Nothing in the library uses this result yet.

Dependency tree · two levels

36 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