Theorem 2.1 — the hyperplane projection method converges to a solution of
ProvedSolodovSvaiterVI.Alg21.theorem2_1_converges_to_solutionLet be a closed convex subset of and continuous. Suppose that the solution set of the variational inequality is nonempty and that condition (1.2) holds:
Let . Then every sequence generated by Algorithm 2.1 converges to a solution of : there is with
The theorem requires neither Lipschitz continuity nor monotonicity of ; condition (1.2) holds for pseudomonotone maps and more generally. Each iteration of the method needs only one projection onto and one onto .
Formalization Note is EuclideanSpace ℝ (Fin n); projections are projOnto, which is the true projection on the nonempty closed convex sets the algorithm uses. A run that stops at a solution () is continued as a constant sequence (the stall encoding of IsAlg21Run), so finitely terminating runs are included and converge trivially. The whole sequence is required to converge to a single solution.
import Mathlib import Definitions.Def_SolodovSvaiterVI_Alg21_viSol import Definitions.Def_SolodovSvaiterVI_Alg21_SatisfiesCond12 import Definitions.Def_SolodovSvaiterVI_Alg21_IsAlg21Run
namespace SolodovSvaiterVI.Alg21
/-- Theorem 2.1 of Solodov–Svaiter (p. 769): let `C ⊆ ℝⁿ` be closed and convex and `F`
continuous; suppose the solution set `S` of `VI(F, C)` is nonempty and condition (1.2) holds.
Then any sequence `{xⁱ}` generated by Algorithm 2.1 (parameters `γ, σ ∈ (0, 1)`) converges to
a solution of `VI(F, C)`.
**Formalization Note.** A run that stops at a solution is continued as a constant sequence
(stall encoding, see `IsAlg21Run`); `P_K` is `projOnto K`. -/
theorem theorem2_1_converges_to_solution {n : ℕ}
(F : EuclideanSpace ℝ (Fin n) → EuclideanSpace ℝ (Fin n))
(C : Set (EuclideanSpace ℝ (Fin n))) (hCc : IsClosed C) (hCcv : Convex ℝ C)
(hF : Continuous F) (hS : (viSol F C).Nonempty) (h12 : SatisfiesCond12 F C)
(gamma sigma : ℝ) (hγ0 : 0 < gamma) (hγ1 : gamma < 1) (hσ0 : 0 < sigma) (hσ1 : sigma < 1)
(x : ℕ → EuclideanSpace ℝ (Fin n)) (k : ℕ → ℕ) (hrun : IsAlg21Run F C gamma sigma x k) :
∃ xs ∈ viSol F C, Filter.Tendsto x Filter.atTop (nhds xs) := by sorry
end SolodovSvaiterVI.Alg21
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.