import Mathlib
import Optlib


open Set InnerProductSpace Topology Filter Bornology Metric Real

theorem extracted_formal_statement_21.{u_3, u_2, u_1} {E₁ : Type u_1} {E₂ : Type u_2} {F : Type u_3}
  [inst : NormedAddCommGroup E₁] [inst_1 : InnerProductSpace ℝ E₁] [inst_2 : FiniteDimensional ℝ E₁]
  [inst_3 : NormedAddCommGroup E₂] [inst_4 : InnerProductSpace ℝ E₂] [inst_5 : FiniteDimensional ℝ E₂]
  [inst_6 : NormedAddCommGroup F] [inst_7 : InnerProductSpace ℝ F] [inst_8 : FiniteDimensional ℝ F]
  {admm : ADMM E₁ E₂ F} {admm_kkt : Existance_of_kkt admm} [inst_9 : Setting E₁ E₂ F admm admm_kkt] (r_Φ : ℝ)
  (hΦ : range Φ ⊆ ball 0 r_Φ)
  (h1 :
    ∀ (k : ℕ+),
      (1 ⊓ (1 + 1 / ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F)) * ADMM.ρ E₁ E₂ F *
          ‖(OptProblem.A₁ E₂) (e₁ (↑k + 1)) + (OptProblem.A₂ E₁) (e₂ (↑k + 1))‖ ^ 2 ≤
        (ADMM.τ E₁ E₂ F ⊓ (1 + ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F ^ 2)) * ADMM.ρ E₁ E₂ F *
            ‖(OptProblem.A₂ E₁) (ADMM.x₂ E₁ F ↑k - ADMM.x₂ E₁ F (↑k + 1))‖ ^ 2 +
          (1 ⊓ (1 + 1 / ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F)) * ADMM.ρ E₁ E₂ F *
            ‖(OptProblem.A₁ E₂) (e₁ (↑k + 1)) + (OptProblem.A₂ E₁) (e₂ (↑k + 1))‖ ^ 2)
  (h2 :
    ∀ (k : ℕ+),
      (ADMM.τ E₁ E₂ F ⊓ (1 + ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F ^ 2)) * ADMM.ρ E₁ E₂ F *
            ‖(OptProblem.A₂ E₁) (ADMM.x₂ E₁ F ↑k - ADMM.x₂ E₁ F (↑k + 1))‖ ^ 2 +
          (1 ⊓ (1 + 1 / ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F)) * ADMM.ρ E₁ E₂ F *
            ‖(OptProblem.A₁ E₂) (e₁ (↑k + 1)) + (OptProblem.A₂ E₁) (e₂ (↑k + 1))‖ ^ 2 ≤
        Φ ↑k - Φ (↑k + 1))
  (h3 : ∀ (k : ℕ+), Φ ↑k - Φ (↑k + 1) < 2 * r_Φ)
  (h4 :
    ∀ (k : ℕ+),
      (1 ⊓ (1 + 1 / ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F)) * ADMM.ρ E₁ E₂ F *
          ‖(OptProblem.A₁ E₂) (e₁ (↑k + 1)) + (OptProblem.A₂ E₁) (e₂ (↑k + 1))‖ ^ 2 <
        2 * r_Φ)
  (h5 :
    ∀ (k : ℕ+),
      ‖(OptProblem.A₁ E₂) (e₁ (↑k + 1)) + (OptProblem.A₂ E₁) (e₂ (↑k + 1))‖ <
        √(2 * r_Φ / ((1 ⊓ (1 + 1 / ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F)) * ADMM.ρ E₁ E₂ F)))
  (n : ℕ) (hn : 1 < n) (h_pos : 0 < n - 1) :
  let k := (n - 1).toPNat h_pos;
  ‖(OptProblem.A₁ E₂) (e₁ (↑k + 1)) + (OptProblem.A₂ E₁) (e₂ (↑k + 1))‖ <
      √(2 * r_Φ / ((1 ⊓ (1 + 1 / ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F)) * ADMM.ρ E₁ E₂ F)) →
    ↑k = n - 1 →
      ‖(OptProblem.A₁ E₂) (e₁ n) + (OptProblem.A₂ E₁) (e₂ n)‖ <
          √(2 * r_Φ / ((1 ⊓ (1 + 1 / ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F)) * ADMM.ρ E₁ E₂ F)) →
        1 ≤ n := sorry

Tactic state before:
E₁ : Type u_1
E₂ : Type u_2
F : Type u_3
inst✝⁹ : NormedAddCommGroup E₁
inst✝⁸ : InnerProductSpace ℝ E₁
inst✝⁷ : FiniteDimensional ℝ E₁
inst✝⁶ : NormedAddCommGroup E₂
inst✝⁵ : InnerProductSpace ℝ E₂
inst✝⁴ : FiniteDimensional ℝ E₂
inst✝³ : NormedAddCommGroup F
inst✝² : InnerProductSpace ℝ F
inst✝¹ : FiniteDimensional ℝ F
admm : ADMM E₁ E₂ F
admm_kkt : Existance_of_kkt admm
inst✝ : Setting E₁ E₂ F admm admm_kkt
r_Φ : ℝ
hΦ : range Φ ⊆ ball 0 r_Φ
h1 :
  ∀ (k : ℕ+),
    (1 ⊓ (1 + 1 / ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F)) * ADMM.ρ E₁ E₂ F *
        ‖(OptProblem.A₁ E₂) (e₁ (↑k + 1)) + (OptProblem.A₂ E₁) (e₂ (↑k + 1))‖ ^ 2 ≤
      (ADMM.τ E₁ E₂ F ⊓ (1 + ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F ^ 2)) * ADMM.ρ E₁ E₂ F *
          ‖(OptProblem.A₂ E₁) (ADMM.x₂ E₁ F ↑k - ADMM.x₂ E₁ F (↑k + 1))‖ ^ 2 +
        (1 ⊓ (1 + 1 / ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F)) * ADMM.ρ E₁ E₂ F *
          ‖(OptProblem.A₁ E₂) (e₁ (↑k + 1)) + (OptProblem.A₂ E₁) (e₂ (↑k + 1))‖ ^ 2
h2 :
  ∀ (k : ℕ+),
    (ADMM.τ E₁ E₂ F ⊓ (1 + ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F ^ 2)) * ADMM.ρ E₁ E₂ F *
          ‖(OptProblem.A₂ E₁) (ADMM.x₂ E₁ F ↑k - ADMM.x₂ E₁ F (↑k + 1))‖ ^ 2 +
        (1 ⊓ (1 + 1 / ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F)) * ADMM.ρ E₁ E₂ F *
          ‖(OptProblem.A₁ E₂) (e₁ (↑k + 1)) + (OptProblem.A₂ E₁) (e₂ (↑k + 1))‖ ^ 2 ≤
      Φ ↑k - Φ (↑k + 1)
h3 : ∀ (k : ℕ+), Φ ↑k - Φ (↑k + 1) < 2 * r_Φ
h4 :
  ∀ (k : ℕ+),
    (1 ⊓ (1 + 1 / ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F)) * ADMM.ρ E₁ E₂ F *
        ‖(OptProblem.A₁ E₂) (e₁ (↑k + 1)) + (OptProblem.A₂ E₁) (e₂ (↑k + 1))‖ ^ 2 <
      2 * r_Φ
h5 :
  ∀ (k : ℕ+),
    ‖(OptProblem.A₁ E₂) (e₁ (↑k + 1)) + (OptProblem.A₂ E₁) (e₂ (↑k + 1))‖ <
      √(2 * r_Φ / ((1 ⊓ (1 + 1 / ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F)) * ADMM.ρ E₁ E₂ F))
n : ℕ
hn : 1 < n
h_pos : 0 < n - 1
k : ℕ+ := (n - 1).toPNat h_pos
h_k :
  ‖(OptProblem.A₁ E₂) (e₁ (↑k + 1)) + (OptProblem.A₂ E₁) (e₂ (↑k + 1))‖ <
    √(2 * r_Φ / ((1 ⊓ (1 + 1 / ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F)) * ADMM.ρ E₁ E₂ F))
k_to_n : ↑k = n - 1
⊢ ‖(OptProblem.A₁ E₂) (e₁ n) + (OptProblem.A₂ E₁) (e₂ n)‖ <
    √(2 * r_Φ / ((1 ⊓ (1 + 1 / ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F)) * ADMM.ρ E₁ E₂ F))

Tactic:
rw [k_to_n] at h_k
    rw [Nat.sub_add_cancel] at h_k
    exact h_k

Tactic state after:
[E₁ : Type u_1
 E₂ : Type u_2
 F : Type u_3
 inst✝⁹ : NormedAddCommGroup E₁
 inst✝⁸ : InnerProductSpace ℝ E₁
 inst✝⁷ : FiniteDimensional ℝ E₁
 inst✝⁶ : NormedAddCommGroup E₂
 inst✝⁵ : InnerProductSpace ℝ E₂
 inst✝⁴ : FiniteDimensional ℝ E₂
 inst✝³ : NormedAddCommGroup F
 inst✝² : InnerProductSpace ℝ F
 inst✝¹ : FiniteDimensional ℝ F
 admm : ADMM E₁ E₂ F
 admm_kkt : Existance_of_kkt admm
 inst✝ : Setting E₁ E₂ F admm admm_kkt
 r_Φ : ℝ
 hΦ : range Φ ⊆ ball 0 r_Φ
 h1 :
   ∀ (k : ℕ+),
     (1 ⊓ (1 + 1 / ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F)) * ADMM.ρ E₁ E₂ F *
         ‖(OptProblem.A₁ E₂) (e₁ (↑k + 1)) + (OptProblem.A₂ E₁) (e₂ (↑k + 1))‖ ^ 2 ≤
       (ADMM.τ E₁ E₂ F ⊓ (1 + ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F ^ 2)) * ADMM.ρ E₁ E₂ F *
           ‖(OptProblem.A₂ E₁) (ADMM.x₂ E₁ F ↑k - ADMM.x₂ E₁ F (↑k + 1))‖ ^ 2 +
         (1 ⊓ (1 + 1 / ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F)) * ADMM.ρ E₁ E₂ F *
           ‖(OptProblem.A₁ E₂) (e₁ (↑k + 1)) + (OptProblem.A₂ E₁) (e₂ (↑k + 1))‖ ^ 2
 h2 :
   ∀ (k : ℕ+),
     (ADMM.τ E₁ E₂ F ⊓ (1 + ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F ^ 2)) * ADMM.ρ E₁ E₂ F *
           ‖(OptProblem.A₂ E₁) (ADMM.x₂ E₁ F ↑k - ADMM.x₂ E₁ F (↑k + 1))‖ ^ 2 +
         (1 ⊓ (1 + 1 / ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F)) * ADMM.ρ E₁ E₂ F *
           ‖(OptProblem.A₁ E₂) (e₁ (↑k + 1)) + (OptProblem.A₂ E₁) (e₂ (↑k + 1))‖ ^ 2 ≤
       Φ ↑k - Φ (↑k + 1)
 h3 : ∀ (k : ℕ+), Φ ↑k - Φ (↑k + 1) < 2 * r_Φ
 h4 :
   ∀ (k : ℕ+),
     (1 ⊓ (1 + 1 / ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F)) * ADMM.ρ E₁ E₂ F *
         ‖(OptProblem.A₁ E₂) (e₁ (↑k + 1)) + (OptProblem.A₂ E₁) (e₂ (↑k + 1))‖ ^ 2 <
       2 * r_Φ
 h5 :
   ∀ (k : ℕ+),
     ‖(OptProblem.A₁ E₂) (e₁ (↑k + 1)) + (OptProblem.A₂ E₁) (e₂ (↑k + 1))‖ <
       √(2 * r_Φ / ((1 ⊓ (1 + 1 / ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F)) * ADMM.ρ E₁ E₂ F))
 n : ℕ
 hn : 1 < n
 h_pos : 0 < n - 1
 k : ℕ+ := (n - 1).toPNat h_pos
 h_k :
   ‖(OptProblem.A₁ E₂) (e₁ (n - 1 + 1)) + (OptProblem.A₂ E₁) (e₂ (n - 1 + 1))‖ <
     √(2 * r_Φ / ((1 ⊓ (1 + 1 / ADMM.τ E₁ E₂ F - ADMM.τ E₁ E₂ F)) * ADMM.ρ E₁ E₂ F))
 k_to_n : ↑k = n - 1
 ⊢ 1 ≤ n]