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

Residual bounds for exact and approximate-state adjoint hypergradients

Statement

Let p,d≥1, let Θ⊆Rp and W⊆Rd be open, and let f:Θ×W→R be C2 and g:W×Θ→R be C1, where C1 means total differentiability with continuous derivative and C2 means the derivative map is C1. Let w∗:P→W be a C1 map on an open neighborhood P of θ0, with w∗(θ0)=w0 and F(θ,w∗(θ))=0 for F(θ,w):=∇wf(θ,w), and define J(θ):=g(w∗(θ),θ), H(θ,w):=DwF(θ,w), and B(θ,w):=DθF(θ,w). Assume H(θ0,w0) is invertible. Write c(w):=∂θg(w,θ0) and d(w):=∇wg(w,θ0), and use the induced Euclidean operator norm for matrices. If the adjoint algorithm at w0 returns a vector v~, define ρ:=d(w0)−H(θ0,w0)Tv~ and h^:=c(w0)−B(θ0,w0)Tv~. Suppose the evaluated Hessian is symmetric and has a valid lower-bound certificate xTH(θ0,w0)x≥μ∥x∥22 for every x, for some μ>0. Then

∥h^−∇θJ(θ0)∥2≤∥B(θ0,w∗(θ0))∥ ∥ρ∥2μ.

There is also a separate approximate-state bound under stronger hypotheses. Put R(w,θ):=∇wf(θ,w) and A(w):=DwR(w,θ0)=H(θ0,w). Suppose R is C2 on an open neighborhood of (w0,θ0). Let w~ lie in a C1 local inverse chart: for r:=R(w~,θ0), there is a local inverse Ψ of w↦R(w,θ0) with Ψ(0)=w0, defined on a convex open residual neighborhood containing 0 and r, with w~=Ψ(r) and ∥DΨ(y)∥≤κ along the segment from 0 to r. Suppose the segment from w0 to w~ lies in a convex state neighborhood where ∥A(w)−1∥≤M, ∥B(θ0,w)∥≤B0, and A(w),B(θ0,w) are Lipschitz with constants LA,LB. Suppose also that c,d are Lipschitz there with constants Lc,Ld and ∥d(w0)∥2≤D∗. If the adjoint algorithm at w~ returns v~, set ρ:=d(w~)−A(w~)Tv~ and h^:=c(w~)−B(θ0,w~)Tv~. Assume the certificate A(w~)=A(w~)T and xTA(w~)x≥μ∥x∥22 for all x, with the same μ>0. Then

∥h^−∇θJ(θ0)∥2≤κ(Lc+D∗(M2LAB0+MLB)+MB0Ld)∥r∥2+B0μ∥ρ∥2.

The first term bounds the implicit formula evaluated at an approximate state; it does not identify that formula with the derivative of a stopped iteration. The second term bounds the separate linear-solve error.

Facts & Assumptions

Given: The exact stationary branch, its exact and approximate states, the local inverse chart and Lipschitz bounds, and a solver output vector with the stated actual residual.

[F1]

Smooth bilevel stationary problems and their hypergradients: R=∇wf, A=H=DwR, B=DθR, and J(θ)=g(w∗(θ),θ).

[F2]

Hypergradient of a smooth local stationary branch: at an exact branch state, ∇θJ=c−BTH−Td.

[F3]

An implicit sensitivity evaluated at an approximate state has a local residual bound: under a common invertible residual chart and the stated inverse and Lipschitz bounds, ∥S(w~)−S(w∗)∥≤κ(M2LAB0+MLB)∥r∥2, where S(w)=−A(w)−1B(w).

[F4]

The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a): the derivative of a composition of finite-dimensional totally differentiable maps is the composition of their derivatives.

[F5]

On a convex open set, a uniform bound ∥Df(z)v∥2≤M∥v∥2 implies ∥f(y)−f(x)∥2≤M∥y−x∥2: on a convex open set, a uniform bound ∥Df(z)v∥2≤M∥v∥2 implies ∥f(y)−f(x)∥2≤M∥y−x∥2.

[F6]

Adjoint hypergradient with a measured linear-solve residual: for any returned vector v~, the algorithm computes the actual residual ρ=d−HTv~ and estimate h^=c−BTv~.

Proof

1.1F1F3given

Put S(w):=−A(w)−1B(w) on the stated state neighborhood. The chart and Lipschitz hypotheses with [F3] give ∥S(w~)−S(w∗)∥≤κ(M2LAB0+MLB)∥r∥2.

1.2F4F5given

For the state bound, set γ(t):=Ψ(tr) for 0<t<1. By [F4], γ′(t)=DΨ(tr)r, so ∥γ′(t)∥2≤κ∥r∥2. For 0<ε<1/2, apply [F5] to γ on the convex open interval (ε/2,1−ε/2) at ε and 1−ε; it gives ∥Ψ((1−ε)r)−Ψ(εr)∥2≤κ(1−2ε)∥r∥2. Letting ε↓0 and using continuity of Ψ, Ψ(0)=w∗ and Ψ(r)=w~ give ∥w~−w∗∥2≤κ∥r∥2.

1.3givenalgebra

The inverse and operator-norm hypotheses give ∥S(w~)∥≤∥A(w~)−1∥ ∥B(w~)∥≤MB0.

1.4F1F2givenalgebra

Define the implicit formula evaluated at a state by hˉ(w):=c(w)+S(w)Td(w). By [F1] and [F2], hˉ(w∗)=∇θJ(θ0). Expansion gives hˉ(w~)−hˉ(w∗)=(c(w~)−c(w∗))+(S(w~)−S(w∗))Td(w∗)+S(w~)T(d(w~)−d(w∗)).

1.5F6F7givenalgebra

At the evaluated state, write H=A(w~) and B=B(θ0,w~). Let v∗=H−Td(w~). The residual identity gives v∗−v~=H−Tρ. For any vector x, the certificate and [F7] imply μ∥x∥22≤xTHx≤∥x∥2∥Hx∥2, hence ∥H−1∥≤1/μ; symmetry gives the same bound for H−T. The Euclidean dual formula and [F7] give ∥BTy∥2≤∥B∥∥y∥2, so ∥h^−hˉ(w~)∥2=∥BTH−Tρ∥2≤(B0/μ)∥ρ∥2.

1.6F2F6F7givenalgebra

For the exact-state claim, put H:=H(θ0,w0), B:=B(θ0,w0), c:=c(w0), and d:=d(w0), and let v∗:=H−Td. Since w0=w∗(θ0), [F2] gives ∇θJ(θ0)=c−BTv∗. The residual identity gives v∗−v~=H−Tρ. For every x, the stated symmetry and lower-bound certificate together with [F7] imply μ∥x∥22≤xTHx≤∥x∥2∥Hx∥2, so ∥H−1∥≤1/μ; symmetry gives ∥H−T∥≤1/μ. The Euclidean dual formula and [F7] give ∥BTy∥2≤∥B∥∥y∥2. Therefore ∥h^−∇θJ(θ0)∥2=∥BTH−Tρ∥2≤∥B∥ ∥ρ∥2μ=∥B(θ0,w∗(θ0))∥ ∥ρ∥2μ. This proves the exact-state claim using only its stated hypotheses.

2.1F1F2F7step 1.1step 1.3step 1.4givenalgebra

For any matrix M and vector x, the Euclidean dual formula ∥MTx∥2=sup⁡∥y∥2=1∣yTMTx∣ and [F7] give ∥MTx∥2≤∥M∥∥x∥2. Put δ:=∥w~−w∗∥2. Apply this norm bound to the three terms in step 1.4. The first and third terms have norms at most Lcδ and MB0Ldδ, respectively, by the stated Lipschitz bounds and step 1.3. The middle term is at most κD∗(M2LAB0+MLB)∥r∥2 by step 1.1 and ∥d(w∗)∥2≤D∗. Therefore ∥hˉ(w~)−∇θJ(θ0)∥2≤(Lc+MB0Ld)δ+κD∗(M2LAB0+MLB)∥r∥2.

3.1step 1.2step 2.1algebra

Substitute the state-residual estimate from step 1.2 into step 2.1. This gives the first term in the asserted approximate-state bound for the implicit formula evaluation.

4.1step 3.1step 1.5givenalgebra∎

Combine the approximate-state formula bound from step 3.1 and the solve-error bound from step 1.5 by the triangle inequality. This gives the approximate-state bound under its separate chart and Lipschitz hypotheses.

Depends on

Used by

Nothing in the library uses this result yet.

Dependency tree · two levels

37 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