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.

Newton, Cauchy, and truncated-CG steps on one indefinite model

Example

Let f(z)=gTz+12zTBz on R2, where B=diag⁡(−1,1), g=(−1,−3)/2, and the trust-region radius is Δ=1. The raw Newton direction is outside the ball but satisfies the exact inexact-Newton descent safeguard and a strong Wolfe line search accepts its full step. The Cauchy and first truncated-CG boundary steps coincide and have smaller predicted reduction. For the inexact-Newton call take M≥2, N≥1, ηk=ηˉ=0, σ=1/2, J=0, E,Z≥1, γ=2, and λ0=λmax⁡=1; use the deterministic two-dimensional projected solver specified in step 1.1 with iteration count m=2, returning no-candidate if its two-dimensional basis or projected matrix is singular or if M<2.

Facts & Assumptions

Given: The quadratic trust-region model, the displayed B,g,Δ, the Cauchy and truncated-CG procedures, and the supplied inexact-Newton and strong-Wolfe line-search algorithms.

[F1]

Quadratic trust-region models and the actual-to-predicted reduction ratio: m(p)=f(0)+gTp+12pTBp and feasible steps satisfy ∥p∥2≤1.

[F2]

Global trust-region minimizers satisfy the shifted-Hessian conditions, including the hard case: every feasible step satisfying (B+λI)p=−g, B+λI⪰0, and λ(1−∥p∥2)=0 is a global trust-region minimizer.

[F3]

The Cauchy step for a quadratic trust-region model: if gTBg>0, then pC=−min⁡{∥g∥22/(gTBg),1/∥g∥2}g.

[F4]

Truncated conjugate gradients with boundary and negative-curvature exits: start with z0=0, r0=g, d0=−g; if d0TBd0>0, set α0=∥r0∥22/(d0TBd0) and, when ∥z0+α0d0∥2≥Δ, return the positive boundary root on z0+td0 as candidate before comparing with the Cauchy step. The unique positive root and feasibility of the selected output are established by Truncated-CG iterates are feasible and lie in generated Krylov subspaces.

[F5]

Finite-budget inexact Newton–Krylov with a descent safeguard: with zero shift, a candidate p is accepted as a linear solve when ∥Bp+g∥2≤ηk∥g∥2 and gTp<−σ∥g∥2∥p∥2; the descent fallback is p=−g.

[F6]

Wolfe and strong Wolfe line-search conditions: a strong Wolfe step obeys ϕ(t)≤ϕ(0)+c1tϕ′(0) and ∣ϕ′(t)∣≤c2∣ϕ′(0)∣, with 0<c1<c2<1/2.

[F7]

Finite-budget bracketing and zoom for a strong Wolfe step: from a descent direction, the first expansion trial is t=1 and it returns that point when both strong Wolfe inequalities hold.

[F8]

The Krylov subspace Km(A,b)=span⁡{b,Ab,…,Am−1b}: K2(B,−g)=span⁡{−g,B(−g)}.

Proof

Given: B=diag⁡(−1,1), g=(−1,−3)/2, Δ=1, and f(z)=gTz+12zTBz.

1.1givenF5F8algebra

The matrix B is invertible and the raw Newton solution is pN=−B−1g=(−1,3)/2. Direct calculation gives BpN+g=0, ∥pN∥2=5, gTpN=−4, ∥g∥2∥pN∥2=5, and the unrestricted model decrease at pN is −gTpN−12pNTBpN=2. Also −g=(1,3)/2 and B(−g)=(−1,3)/2 are linearly independent, so [F8] gives K2(B,−g)=R2. Put V=[−g B(−g)]; its determinant is 3, so the projected matrix VTBV is invertible because its determinant is det⁡(V)2det⁡(B)=−9. Solving VTBVy=VT(−g) gives p=Vy∈K2 with Bp+g orthogonal to K2=R2, hence Bp+g=0 and p=pN. Thus a deterministic two-dimensional Galerkin Krylov solver may return pN. With η0=0 and descent parameter σ=1/2, [F5]'s residual and strict descent tests pass because 0≤0 and −4<−5/2.

1.2givenF1F3F4algebra

We have ∥g∥22=5 and gTBg=4. By [F3], pC=(1,3)/10 and ∥pC∥2=1. Its predicted reduction is −gTpC−12pCTBpC=5−25. For truncated CG with residual tolerance zero and inner cap at least one, d0=−g, d0TBd0=4, and α0=5/4. The unconstrained point z^=−(5/4)g has norm 55/4>1, so the boundary exit is z0+t0d0 with t0=1/5. It equals pC, and the final comparison in [F4] returns a step with the same predicted reduction 5−2/5.

2.1givenF6F7step 1.1algebra

For ϕ(t):=f(tpN), direct substitution gives ϕ(t)=−4t+2t2 and ϕ′(t)=−4+4t. Set c1=1/10 and c2=1/4. At t=1, ϕ(1)=−2≤−2/5=c1ϕ′(0) and ∣ϕ′(1)∣=0≤1=c2∣ϕ′(0)∣. Thus the full step satisfies [F6]; with expansion budget at least one, [F7] returns t=1.

3.1givenF1F2step 2.1step 1.2algebra∎

The point p∗=(1,1)/2 has norm one and (B+2I)p∗=(1,3)/2=−g, while B+2I=diag⁡(1,3)⪰0. By [F2], it is a global trust- region minimizer. Its predicted reduction is −gTp∗−12p∗TBp∗=2. The Cauchy decrease is strictly smaller because 5<12/5, which follows from 5<144/25. Thus the raw Newton line-search step is accepted at full length, whereas the feasible Cauchy and truncated-CG steps have the strictly smaller reduction 5−2/5<2.

Depends on

Used by

Nothing in the library uses this result yet.

Dependency tree · two levels

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