Alphabeta Math
LemmaStatement: Literature-sourcedProof: AI-adaptedPipeline-generatedprecheck passjudge pass (gpt-6-luna)
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.

Fixed-memory L-Bfgs curvature and scale bounds imply uniform spectral bounds

Statement

Fix n,m≥1, constants 0<a≤b<∞, C>0, and 0<γmin⁡≤γmax⁡<∞. Suppose every stored pair in a limited-memory BFGS state satisfies a∥si∥22≤siTyi≤b∥si∥22,∥yi∥22siTyi≤C, and each base scale satisfies γmin⁡≤γ≤γmax⁡. Then the implicit SPD matrix Hk represented by the two-loop recursion of Limited-memory Bfgs two-loop recursion with curvature filtering has uniform bounds λ‾I⪯Hk⪯λ‾I,T:=nγmin⁡+mC,r:=min⁡{1,a/T},D:=γmax⁡−nrm,λ‾:=1T,λ‾:=Tn−1D. Here A⪯B means zTAz≤zTBz for all z. If the stored list is nonempty, the usual scale (slastTylast)/∥ylast∥22 lies in [1/C,b/a2].

Facts & Assumptions

Given: n,m,a,b,C,γmin⁡,γmax⁡, the stored-pair bounds, and the L-BFGS state in the Statement.

[F1]

The L-Bfgs two-loop recursion equals the stored inverse updates and costs linear memory work: the implicit matrix is obtained by at most m sequential inverse BFGS updates from H0=γI.

[F2]

Bfgs updates preserve the secant equations, symmetry, and positive definiteness: each positive-curvature inverse update preserves positive definiteness and is reciprocal to its Hessian update.

[F3]

Bfgs trace and determinant identities give a matrix potential inequality: for B=H−1, one update adds at most ∥y∥22/(sTy) to tr⁡(B) and multiplies det⁡(B) by (sTy)/(sTBs).

[F4]

In an orthonormal basis, self-adjoint means conjugate-transpose symmetry and normal means commuting with the conjugate transpose and Real spectral theorem: a self-adjoint endomorphism of a finite-dimensional real inner product space has an orthonormal eigenbasis: a real symmetric matrix represents a self-adjoint endomorphism of Euclidean Rn and has an orthonormal eigenbasis. If the matrix is SPD, each eigenvalue is positive by evaluating its quadratic form on a unit eigenvector.

[F5]

Limited-memory Bfgs two-loop recursion with curvature filtering: the scale from the newest pair is (slastTylast)/(ylastTylast).

Proof

Given: n,m,a,b,C,γmin⁡,γmax⁡, the stored-pair bounds, and the L-BFGS state in the Statement.

[F1] The L-Bfgs two-loop recursion equals the stored inverse updates and costs linear memory work: the implicit matrix is obtained by at most m sequential inverse BFGS updates from H0=γI.

[F2] Bfgs updates preserve the secant equations, symmetry, and positive definiteness: each positive-curvature inverse update preserves positive definiteness and is reciprocal to its Hessian update.

[F3] Bfgs trace and determinant identities give a matrix potential inequality: for B=H−1, one update adds at most ∥y∥22/(sTy) to tr⁡(B) and multiplies det⁡(B) by (sTy)/(sTBs).

[F4] In an orthonormal basis, self-adjoint means conjugate-transpose symmetry and normal means commuting with the conjugate transpose and Real spectral theorem: a self-adjoint endomorphism of a finite-dimensional real inner product space has an orthonormal eigenbasis: a real symmetric matrix represents a self-adjoint endomorphism of Euclidean Rn and has an orthonormal eigenbasis. If the matrix is SPD, each eigenvalue is positive by evaluating its quadratic form on a unit eigenvector.

[F5] Limited-memory Bfgs two-loop recursion with curvature filtering: the scale from the newest pair is (slastTylast)/(ylastTylast).

1.1givenF1F2F3algebra

Let Bk:=Hk−1. By [F1] and [F2], this is an SPD matrix obtained from B0=γ−1I by at most m direct BFGS updates. Its initial trace satisfies tr⁡(B0)=n/γ≤n/γmin⁡. Each update adds at most C to trace by [F3] and the pair bound. Hence tr⁡(Bk)≤T.

2.1givenF3F4F6step 1.1algebra

By [F4], choose an orthonormal eigenbasis for Bk and let its positive eigenvalues be λ1,…,λn. In this basis the characteristic polynomial factors as ∏i(t−λi), so [F6] gives tr⁡(Bk)=∑iλi and det⁡(Bk)=∏iλi. By step 1.1 each λi≤T, and expanding s in this orthonormal basis gives sTBks≤T∥s∥22. The determinant ratio in [F3] and the lower pair bound therefore give det⁡(Bk+1)det⁡(Bk)=sTysTBks≥aT≥r. Since 0<r≤1 and there are at most m updates, while det⁡(B0)=γ−n≥γmax⁡−n, it follows that det⁡(Bk)≥D.

3.1F4F6step 2.1algebra

Let λmin⁡(Bk) be the smallest eigenvalue. By step 2.1, det⁡(Bk) is the product of the n positive eigenvalues, and the other n−1 eigenvalues are at most T. Thus λmin⁡(Bk)≥D/Tn−1, while every eigenvalue is at most T. The eigenvalues of Hk=Bk−1 are the reciprocals, so they lie in [1/T,Tn−1/D]. The spectral theorem translates these bounds into the asserted quadratic-form inequalities.

4.1givenF5algebra∎

If the stored list is nonempty, take its newest pair. For the standard scale, [F5] and the upper curvature bound give γ=(sTy)/∥y∥22≤b∥s∥22/∥y∥22. By Cauchy–Schwarz and the lower curvature bound, ∥s∥2∥y∥2≥sTy≥a∥s∥22, so ∥y∥2≥a∥s∥2 and γ≤b/a2. The other pair bound gives ∥y∥22/(sTy)≤C, hence γ≥1/C.

Depends on

Used by

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.

Sources