import Mathlib
import Optlib


open Set Bornology Filter BigOperators Topology BCD

open BCD

theorem extracted_formal_statement_87.{u_2, u_1} {E : Type u_1} {F : Type u_2} [inst : NormedAddCommGroup E]
  [inst_1 : InnerProductSpace ℝ E] [inst_2 : CompleteSpace E] [inst_3 : NormedAddCommGroup F]
  [inst_4 : InnerProductSpace ℝ F] [inst_5 : CompleteSpace F] {f : E → ℝ} {g : F → ℝ} {H : WithLp 2 (E × F) → ℝ}
  {x0 : E} {y0 : F} {l : NNReal} {alg : BCD f g H l x0 y0} (γ : ℝ) (hγ : γ > 1) (ck : ∀ (k : ℕ), alg.c k = 1 / (γ * ↑l))
  (dk : ∀ (k : ℕ), alg.d k = 1 / (γ * ↑l)) (k : ℕ)
  (lip : ∀ (x y : WithLp 2 (E × F)), ‖gradient H x - gradient H y‖ ≤ ↑l * ‖x - y‖) (cpos' : (alg.c k)⁻¹ ≥ 0)
  (dpos' : (alg.d k)⁻¹ ≥ 0) (h1 : ‖(subdiff k).1‖ ≤ ↑l * (γ + 1) * ‖z (k + 1) - z k‖) :
  let a := (alg.d k)⁻¹ • (alg.y k - alg.y (k + 1));
  ‖a + (gradient H (alg.x (k + 1), alg.y (k + 1))).2 - (gradient H (alg.x (k + 1), alg.y k)).2‖ =
      ‖a + (gradient H (alg.x (k + 1), alg.y (k + 1))).2 - grad_snd H (alg.x (k + 1)) (alg.y k)‖ →
    ‖a +
            ((WithLp.equiv 2 (E × F)).symm
                (grad_fst H (alg.x (k + 1), alg.y (k + 1)).2 (alg.x (k + 1), alg.y (k + 1)).1,
                  grad_snd H (alg.x (k + 1), alg.y (k + 1)).1 (alg.x (k + 1), alg.y (k + 1)).2)).2 -
          ((WithLp.equiv 2 (E × F)).symm
              (grad_fst H (alg.x (k + 1), alg.y k).2 (alg.x (k + 1), alg.y k).1,
                grad_snd H (alg.x (k + 1), alg.y k).1 (alg.x (k + 1), alg.y k).2)).2‖ =
      ‖a +
            ((WithLp.equiv 2 (E × F)).symm
                (grad_fst H (alg.x (k + 1), alg.y (k + 1)).2 (alg.x (k + 1), alg.y (k + 1)).1,
                  grad_snd H (alg.x (k + 1), alg.y (k + 1)).1 (alg.x (k + 1), alg.y (k + 1)).2)).2 -
          grad_snd H (alg.x (k + 1)) (alg.y k)‖ := sorry

Tactic state before:
E : Type u_1
F : Type u_2
inst✝⁵ : NormedAddCommGroup E
inst✝⁴ : InnerProductSpace ℝ E
inst✝³ : CompleteSpace E
inst✝² : NormedAddCommGroup F
inst✝¹ : InnerProductSpace ℝ F
inst✝ : CompleteSpace F
f : E → ℝ
g : F → ℝ
H : WithLp 2 (E × F) → ℝ
x0 : E
y0 : F
l : NNReal
alg : BCD f g H l x0 y0
γ : ℝ
hγ : γ > 1
ck : ∀ (k : ℕ), alg.c k = 1 / (γ * ↑l)
dk : ∀ (k : ℕ), alg.d k = 1 / (γ * ↑l)
k : ℕ
lip : ∀ (x y : WithLp 2 (E × F)), ‖gradient H x - gradient H y‖ ≤ ↑l * ‖x - y‖
cpos' : (alg.c k)⁻¹ ≥ 0
dpos' : (alg.d k)⁻¹ ≥ 0
h1 : ‖(subdiff k).1‖ ≤ ↑l * (γ + 1) * ‖z (k + 1) - z k‖
a : F := (alg.d k)⁻¹ • (alg.y k - alg.y (k + 1))
⊢ ‖a + (gradient H (alg.x (k + 1), alg.y (k + 1))).2 - (gradient H (alg.x (k + 1), alg.y k)).2‖ =
    ‖a + (gradient H (alg.x (k + 1), alg.y (k + 1))).2 - grad_snd H (alg.x (k + 1)) (alg.y k)‖

Tactic:
rw [grad_eq_block_grad, grad_fun_comp, grad_comp, grad_fun_comp, grad_comp]

Tactic state after:
[E : Type u_1
 F : Type u_2
 inst✝⁵ : NormedAddCommGroup E
 inst✝⁴ : InnerProductSpace ℝ E
 inst✝³ : CompleteSpace E
 inst✝² : NormedAddCommGroup F
 inst✝¹ : InnerProductSpace ℝ F
 inst✝ : CompleteSpace F
 f : E → ℝ
 g : F → ℝ
 H : WithLp 2 (E × F) → ℝ
 x0 : E
 y0 : F
 l : NNReal
 alg : BCD f g H l x0 y0
 γ : ℝ
 hγ : γ > 1
 ck : ∀ (k : ℕ), alg.c k = 1 / (γ * ↑l)
 dk : ∀ (k : ℕ), alg.d k = 1 / (γ * ↑l)
 k : ℕ
 lip : ∀ (x y : WithLp 2 (E × F)), ‖gradient H x - gradient H y‖ ≤ ↑l * ‖x - y‖
 cpos' : (alg.c k)⁻¹ ≥ 0
 dpos' : (alg.d k)⁻¹ ≥ 0
 h1 : ‖(subdiff k).1‖ ≤ ↑l * (γ + 1) * ‖z (k + 1) - z k‖
 a : F := (alg.d k)⁻¹ • (alg.y k - alg.y (k + 1))
 ⊢ ‖a +
           ((WithLp.equiv 2 (E × F)).symm
               (grad_fst H (alg.x (k + 1), alg.y (k + 1)).2 (alg.x (k + 1), alg.y (k + 1)).1,
                 grad_snd H (alg.x (k + 1), alg.y (k + 1)).1 (alg.x (k + 1), alg.y (k + 1)).2)).2 -
         ((WithLp.equiv 2 (E × F)).symm
             (grad_fst H (alg.x (k + 1), alg.y k).2 (alg.x (k + 1), alg.y k).1,
               grad_snd H (alg.x (k + 1), alg.y k).1 (alg.x (k + 1), alg.y k).2)).2‖ =
     ‖a +
           ((WithLp.equiv 2 (E × F)).symm
               (grad_fst H (alg.x (k + 1), alg.y (k + 1)).2 (alg.x (k + 1), alg.y (k + 1)).1,
                 grad_snd H (alg.x (k + 1), alg.y (k + 1)).1 (alg.x (k + 1), alg.y (k + 1)).2)).2 -
         grad_snd H (alg.x (k + 1)) (alg.y k)‖,
 E : Type u_1
 F : Type u_2
 inst✝⁵ : NormedAddCommGroup E
 inst✝⁴ : InnerProductSpace ℝ E
 inst✝³ : CompleteSpace E
 inst✝² : NormedAddCommGroup F
 inst✝¹ : InnerProductSpace ℝ F
 inst✝ : CompleteSpace F
 f : E → ℝ
 g : F → ℝ
 H : WithLp 2 (E × F) → ℝ
 x0 : E
 y0 : F
 l : NNReal
 alg : BCD f g H l x0 y0
 γ : ℝ
 hγ : γ > 1
 ck : ∀ (k : ℕ), alg.c k = 1 / (γ * ↑l)
 dk : ∀ (k : ℕ), alg.d k = 1 / (γ * ↑l)
 k : ℕ
 lip : ∀ (x y : WithLp 2 (E × F)), ‖gradient H x - gradient H y‖ ≤ ↑l * ‖x - y‖
 cpos' : (alg.c k)⁻¹ ≥ 0
 dpos' : (alg.d k)⁻¹ ≥ 0
 h1 : ‖(subdiff k).1‖ ≤ ↑l * (γ + 1) * ‖z (k + 1) - z k‖
 a : F := (alg.d k)⁻¹ • (alg.y k - alg.y (k + 1))
 ⊢ Differentiable ℝ H]