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

Full-memory Bfgs is locally superlinear when the Hessian is Lipschitz

Statement

Assume the hypotheses of Full-memory Bfgs converges globally and linearly on a uniformly strongly convex level set and consider an infinite, nonstationary exact-undamped BFGS run with successful strong Wolfe searches. Under these same hypotheses, Locally Lipschitz Hessians give the BFGS Dennis–Moré condition proves the Dennis–Moré condition used below. Suppose the Hessian is LG-Lipschitz on a neighborhood of its minimizer x∗: ∥∇2f(u)−∇2f(v)∥2≤LG∥u−v∥2. Then the BFGS matrices satisfy the Dennis–Moré condition ∥(Bk−G∗)sk∥2∥sk∥2⟶0,G∗:=∇2f(x∗), and the line search eventually accepts its initial unit trial. Consequently, ∥xk+1−x∗∥2∥xk−x∗∥2⟶0. Here a sequence is superlinear when this last error ratio tends to zero. The conclusion concerns the exact undamped branch; Powell-damped and rounded-direction variants are not covered by this theorem.

Facts & Assumptions

Given: The global BFGS hypotheses, the exact-undamped run, and the local Hessian-Lipschitz bound in the Statement.

[F1]

Full-memory Bfgs converges globally and linearly on a uniformly strongly convex level set: the iterates converge to x∗, the Hessian bounds hold on L0, and the initial inverse matrix is SPD.

[F2]

Locally Lipschitz Hessians give the BFGS Dennis–Moré condition: ∥(Bk−G∗)sk∥2/∥sk∥2→0; since sk=αkpk with αk>0, the same holds with pk.

[F3]

Full-memory Bfgs with Wolfe search, descent reset, and Powell damping: in exact mode pk=−Hkgk, Bk=Hk−1, and hence Bkpk=−gk.

[F4]

Uniform Hessian bounds control gradient changes and Taylor remainders: on a convex local neighborhood with Hessian bounds and Lipschitz Hessian, ∥g(x∗+z)−g∗−G∗z∥2≤LG∥z∥22 and ∥g(x∗+z)−g∗∥2≤Λ∥z∥2.

[F5]

Second-order Taylor expansion f(a+h)=f(a)+∇f(a)⋅h+12hTHf(a)h+o(∥h∥2): f(x∗+z)=f(x∗)+g∗Tz+12zTG∗z+o(∥z∥22) for a C2 scalar field.

[F6]

Finite-budget bracketing and zoom for a strong Wolfe step: the first expansion trial is α=1 and is returned if both strong Wolfe checks pass.

[F7]

Wolfe and strong Wolfe line-search conditions: strong Wolfe requires Armijo and ∣gk+1Tpk∣≤c2∣gkTpk∣.

[F10]

Clairaut--Schwarz theorem for continuous second partial derivatives: under the C2 hypothesis, the mixed second partials commute, so the Hessian quadratic form in [F5] uses the Hessian convention G∗=∇2f(x∗) of [F8].

Proof

Given: The hypotheses and notation in the Statement. Put ek:=xk−x∗, G∗:=G(x∗), and gk:=∇f(xk).

1.1givenF1F8algebra

The global minimizer x∗ is interior to Rn, so each coordinate restriction has a minimum at x∗ and g∗:=∇f(x∗)=0. By [F1], ek→0 and the eigenvalues of G∗ lie in [μ,Λ].

1.2givenF2F3algebra

The Dennis–Moré condition [F2] and sk=αkpk, αk>0, give ∥(Bk−G∗)pk∥2=o(∥pk∥2).

2.1givenF1F2F3F4F8F9step 1.2algebra

For all sufficiently large k, [F4] gives gk=G∗ek+rk with ∥rk∥2≤LG∥ek∥22. By [F3], Bkpk=−gk, so G∗pk=−gk+(G∗−Bk)pk. The smallest eigenvalue of G∗ is at least μ, hence μ∥pk∥2≤∥G∗pk∥2≤Λ∥ek∥2+o(∥pk∥2), using [F9]. For sufficiently large k, the o(∥pk∥2) term is at most μ∥pk∥2/2, so μ∥pk∥2/2≤Λ∥ek∥2 and ∥pk∥2≤2Λ∥ek∥2/μ. Substitution into G∗(ek+pk)=(G∗−Bk)pk−rk gives G∗(ek+pk)=o(∥ek∥2); the same eigenvalue lower bound yields ek+pk=o(∥ek∥2).

3.1givenF1F4F5F6F7F9F10step 2.1algebra

For sufficiently large k, ek and ek+pk lie in the local neighborhood. Apply [F4] and [F5] to z=ek and z=ek+pk, using [F10] to identify the Hessian convention. Since ek+pk=o(∥ek∥2), these expansions give gkTpk=−ekTG∗ek+o(∥ek∥22),f(xk+pk)−f(xk)=−12ekTG∗ek+o(∥ek∥22),g(xk+pk)=o(∥ek∥2). Because G∗⪰μI and c1<1/2, f(xk+pk)−f(xk)−c1gkTpk=−(1/2−c1)ekTG∗ek+o(∥ek∥22)<0 for all sufficiently large k, so Armijo holds. Also pk=O(∥ek∥2) by step 2.1; [F9] gives ∣g(xk+pk)Tpk∣=o(∥ek∥22), while ∣gkTpk∣≥(μ/2)∥ek∥22 eventually. Thus the strong curvature condition holds. The search first tries α=1 by [F6], and the two checks are exactly [F7], so the unit step is accepted.

4.1step 2.1step 3.1∎

Eventually xk+1=xk+pk, so ek+1=ek+pk=o(∥ek∥2) by step 2.1. This is the claimed superlinear error ratio.

Depends on

Used by

Nothing in the library uses this result yet.

Dependency tree · two levels

51 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