Lemma 3.3 —
ProvedAndersonAccel.Safe.norm_windowB_leanderson-accelerationoperator-normp2o-batch-p200ap2o-gran-per-chapterp2o-plan-paperp2o-v1rank-one-update
Let and let be the max-memory. Let and be vectors in with , and let be the Powell-regularized matrix of (3.5). Assume
- every update is well-defined: ;
- ;
- ,
for . Then, in the induced operator norm,
The bound is uniform in the iteration and depends only on , and the max-memory ; together with Lemma 3.2 it controls the inverse matrices used by the algorithm.
Formalization Note The paper assumes " is chosen by rule (3.7)" and uses it, in the first line of its proof, only through and ; these two consequences are the hypotheses here. The connection with rule (3.7) along the algorithm is Corollary 3.5. The exponent is the max-memory , not .
Preamble
import Mathlib import Definitions.Def_AndersonAccel_Safe_Basic import Definitions.Def_AndersonAccel_Safe_windowB
Formal statement
namespace AndersonAccel.Safe
/-- Lemma 3.3 (p. 3178). Under the hypotheses of Lemma 3.2, if `‖y_i‖ ≤ 2‖s_i‖` on the window and
the window satisfies what rule (3.7) guarantees (`‖ŝ_i‖ ≥ τ‖s_i‖` and `m_k ≤ m`), then
`‖B‖ ≤ 3((1 + θ̄ + τ)/τ)^m - 2` in the operator 2-norm. -/
theorem norm_windowB_le {n : ℕ} (θbar τ : ℝ) (m : ℕ) (hθ0 : 0 < θbar) (hθ1 : θbar < 1)
(hτ0 : 0 < τ) (hτ1 : τ < 1) (hm : 1 ≤ m)
(s y : ℕ → EuclideanSpace ℝ (Fin n)) (mk : ℕ) (hmk : mk ≤ m)
(hwd : ∀ i < mk, inner ℝ (windowShat s i) (s i) ≠ 0)
(hy : ∀ i < mk, ‖y i‖ ≤ 2 * ‖s i‖)
(hτs : ∀ i < mk, τ * ‖s i‖ ≤ ‖windowShat s i‖) :
‖windowB θbar s y mk‖ ≤ 3 * ((1 + θbar + τ) / τ) ^ m - 2 := by sorry
end AndersonAccel.Safe
Source
Zhang, O'Donoghue, Boyd, SIAM J. Optim. 30 (2020), p. 3178, Lemma 3.3 (rule (3.7), p. 3177)
Human review
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.