Documentation

ConvexAnalysisMonotoneOperators_BauschkeCombettes_2017.Chap05.Example_5_21

noncomputable def relaxedProjectionFamily {H : Type u} [NormedAddCommGroup H] [InnerProductSpace H] [CompleteSpace H] {I : Type v} (C : ISet H) (hC_closed : ∀ (j : I), IsClosed (C j)) (hC_convex : ∀ (j : I), Convex (C j)) (hC_inter_nonempty : (⋂ (j : I), C j).Nonempty) (β : I) :
IHH

The relaxed projector family Tᵢ = (1 - βᵢ) Id + βᵢ P_{Cᵢ} attached to a family of nonempty closed convex sets with nonempty common intersection.

Instances For
    @[simp]
    theorem relaxedProjectionFamily_apply {H : Type u} [NormedAddCommGroup H] [InnerProductSpace H] [CompleteSpace H] {I : Type v} (C : ISet H) (hC_closed : ∀ (j : I), IsClosed (C j)) (hC_convex : ∀ (j : I), Convex (C j)) (hC_inter_nonempty : (⋂ (j : I), C j).Nonempty) (β : I) (j : I) (x : H) :
    relaxedProjectionFamily C hC_closed hC_convex hC_inter_nonempty β j x = (1 - β j) x + β j P[C j, ] x

    Evaluating the relaxed projector family gives the textbook formula Tᵢ x = (1 - βᵢ) x + βᵢ P_{Cᵢ} x.

    @[reducible, inline]
    noncomputable abbrev stringBlockCompositionOperator {H : Type u} [NormedAddCommGroup H] [InnerProductSpace H] [CompleteSpace H] {I : Type v} {p : ℕ+} (m : Fin pℕ+) (i : (k : Fin p) → Fin (m k)I) (C : ISet H) (hC_closed : ∀ (j : I), IsClosed (C j)) (hC_convex : ∀ (j : I), Convex (C j)) (hC_inter_nonempty : (⋂ (j : I), C j).Nonempty) (β : I) (k : Fin p) :
    HH

    The k-th ordered block composition of the relaxed projector family.

    Instances For
      @[simp]
      theorem stringBlockCompositionOperator_eq_finiteComposition {H : Type u} [NormedAddCommGroup H] [InnerProductSpace H] [CompleteSpace H] {I : Type v} {p : ℕ+} (m : Fin pℕ+) (i : (k : Fin p) → Fin (m k)I) (C : ISet H) (hC_closed : ∀ (j : I), IsClosed (C j)) (hC_convex : ∀ (j : I), Convex (C j)) (hC_inter_nonempty : (⋂ (j : I), C j).Nonempty) (β : I) (k : Fin p) :
      stringBlockCompositionOperator m i C hC_closed hC_convex hC_inter_nonempty β k = finiteComposition fun (l : Fin (m k)) => relaxedProjectionFamily C hC_closed hC_convex hC_inter_nonempty β (i k l)

      The k-th string block operator is the ordered composition of the relaxed projectors indexed by that block.

      @[reducible, inline]
      noncomputable abbrev stringAveragedRelaxedProjectionOperator {H : Type u} [NormedAddCommGroup H] [InnerProductSpace H] [CompleteSpace H] {I : Type v} {p : ℕ+} (ω : Fin p) (m : Fin pℕ+) (i : (k : Fin p) → Fin (m k)I) (C : ISet H) (hC_closed : ∀ (j : I), IsClosed (C j)) (hC_convex : ∀ (j : I), Convex (C j)) (hC_inter_nonempty : (⋂ (j : I), C j).Nonempty) (β : I) :
      HH

      The string-averaged relaxed projection operator given by the weighted average of the ordered block compositions of the relaxed projectors.

      Instances For
        @[simp]
        theorem stringAveragedRelaxedProjectionOperator_apply {H : Type u} [NormedAddCommGroup H] [InnerProductSpace H] [CompleteSpace H] {I : Type v} {p : ℕ+} (ω : Fin p) (m : Fin pℕ+) (i : (k : Fin p) → Fin (m k)I) (C : ISet H) (hC_closed : ∀ (j : I), IsClosed (C j)) (hC_convex : ∀ (j : I), Convex (C j)) (hC_inter_nonempty : (⋂ (j : I), C j).Nonempty) (β : I) (x : H) :
        stringAveragedRelaxedProjectionOperator ω m i C hC_closed hC_convex hC_inter_nonempty β x = k : Fin p, ω k stringBlockCompositionOperator m i C hC_closed hC_convex hC_inter_nonempty β k x

        Evaluating the string-averaged relaxed projection operator gives the weighted sum of the block compositions of the relaxed projector family.

        noncomputable def stringAveragedRelaxedProjectionOrbit {H : Type u} [NormedAddCommGroup H] [InnerProductSpace H] [CompleteSpace H] {I : Type v} {p : ℕ+} (ω : Fin p) (m : Fin pℕ+) (i : (k : Fin p) → Fin (m k)I) (C : ISet H) (hC_closed : ∀ (j : I), IsClosed (C j)) (hC_convex : ∀ (j : I), Convex (C j)) (hC_inter_nonempty : (⋂ (j : I), C j).Nonempty) (β : I) (x₀ : H) :
        H

        The orbit generated by iterating the string-averaged relaxed projection operator from x₀.

        Instances For
          theorem stringAveragedRelaxedProjectionOrbit_succ {H : Type u} [NormedAddCommGroup H] [InnerProductSpace H] [CompleteSpace H] {I : Type v} {p : ℕ+} (ω : Fin p) (m : Fin pℕ+) (i : (k : Fin p) → Fin (m k)I) (C : ISet H) (hC_closed : ∀ (j : I), IsClosed (C j)) (hC_convex : ∀ (j : I), Convex (C j)) (hC_inter_nonempty : (⋂ (j : I), C j).Nonempty) (β : I) (x₀ : H) (n : ) :
          stringAveragedRelaxedProjectionOrbit ω m i C hC_closed hC_convex hC_inter_nonempty β x₀ (n + 1) = stringAveragedRelaxedProjectionOperator ω m i C hC_closed hC_convex hC_inter_nonempty β (stringAveragedRelaxedProjectionOrbit ω m i C hC_closed hC_convex hC_inter_nonempty β x₀ n)

          The string-averaged relaxed projection orbit satisfies the recursion xₙ₊₁ = S xₙ, where S is the associated string-averaged relaxed projection operator.

          theorem exists_tendsto_weakly_to_point_in_iInter_of_stringAveragedRelaxedProjectionOrbit {H : Type u} [NormedAddCommGroup H] [InnerProductSpace H] [CompleteSpace H] {I : Type v} {p : ℕ+} (m : Fin pℕ+) (i : (k : Fin p) → Fin (m k)I) (C : ISet H) (hC_closed : ∀ (j : I), IsClosed (C j)) (hC_convex : ∀ (j : I), Convex (C j)) (hC_inter_nonempty : (⋂ (j : I), C j).Nonempty) (β : I) ( : ∀ (j : I), β j Set.Ioo 0 2) (ω : Fin p) ( : ∀ (k : Fin p), ω k Set.Ioc 0 1) (hω_sum : k : Fin p, ω k = 1) (hcover : Function.Surjective fun (q : (k : Fin p) × Fin (m k)) => i q.fst q.snd) (x₀ : H) :
          z⋂ (j : I), C j, Filter.Tendsto (fun (n : ) => (toWeakSpace H) (stringAveragedRelaxedProjectionOrbit ω m i C hC_closed hC_convex hC_inter_nonempty β x₀ n)) Filter.atTop (nhds ((toWeakSpace H) z))

          Example 5.21: for a finite family of closed convex sets with nonempty intersection, the orbit generated by weighted string averages of the relaxed projectors Tᵢ = (1 - βᵢ) Id + βᵢ P_{Cᵢ} converges weakly to a point of the common intersection.