import Mathlib
import Optlib


open Set InnerProductSpace Topology Filter InnerProduct

open scoped Pointwise

theorem extracted_formal_statement_164.{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] (n : ℕ) :
  let e' := e₁ (n + 1);
  let block :=
    -ey (n + 1) -
        ((1 - ADMM.τ E₁ E₂ F) * ADMM.ρ E₁ E₂ F) • ((OptProblem.A₁ E₂) (e₁ (n + 1)) + (OptProblem.A₂ E₁) (e₂ (n + 1))) -
      ADMM.ρ E₁ E₂ F • (OptProblem.A₂ E₁) (ADMM.x₂ E₁ F n - ADMM.x₂ E₁ F (n + 1));
  let u' :=
    -(ContinuousLinearMap.adjoint (OptProblem.A₁ E₂))
        (ADMM.y E₁ E₂ (n + 1) +
            ((1 - ADMM.τ E₁ E₂ F) * ADMM.ρ E₁ E₂ F) •
              ((OptProblem.A₁ E₂) (e₁ (n + 1)) + (OptProblem.A₂ E₁) (e₂ (n + 1))) +
          ADMM.ρ E₁ E₂ F • (OptProblem.A₂ E₁) (ADMM.x₂ E₁ F n - ADMM.x₂ E₁ F (n + 1)));
  let x_diff := ADMM.x₁ E₂ F (n + 1) - admm_kkt.x₁;
  let block₁ :=
    ADMM.y E₁ E₂ (n + 1) +
        ((1 - ADMM.τ E₁ E₂ F) * ADMM.ρ E₁ E₂ F) • ((OptProblem.A₁ E₂) (e₁ (n + 1)) + (OptProblem.A₂ E₁) (e₂ (n + 1))) +
      ADMM.ρ E₁ E₂ F • (OptProblem.A₂ E₁) (ADMM.x₂ E₁ F n - ADMM.x₂ E₁ F (n + 1));
  -block = block₁ - admm_kkt.y →
    -(ContinuousLinearMap.adjoint (OptProblem.A₁ E₂)) block₁ = u' →
      (ContinuousLinearMap.adjoint (OptProblem.A₁ E₂)) admm_kkt.y =
          (ContinuousLinearMap.adjoint (OptProblem.A₁ E₂)) admm_kkt.y →
        ⟪(ContinuousLinearMap.adjoint (OptProblem.A₁ E₂)) admm_kkt.y -
                (ContinuousLinearMap.adjoint (OptProblem.A₁ E₂)) block₁,
              e'⟫_ℝ =
            ⟪u' + (ContinuousLinearMap.adjoint (OptProblem.A₁ E₂)) admm_kkt.y, x_diff⟫_ℝ →
          ⟪(ContinuousLinearMap.adjoint (OptProblem.A₁ E₂)) admm_kkt.y +
                -(ContinuousLinearMap.adjoint (OptProblem.A₁ E₂)) block₁,
              e'⟫_ℝ =
            ⟪(ContinuousLinearMap.adjoint (OptProblem.A₁ E₂)) admm_kkt.y +
                -(ContinuousLinearMap.adjoint (OptProblem.A₁ E₂)) block₁,
              x_diff⟫_ℝ := 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
n : ℕ
Ae1 : F := (OptProblem.A₁ E₂) (e₁ (n + 1))
e' : E₁ := e₁ (n + 1)
block : F :=
  -ey (n + 1) -
      ((1 - ADMM.τ E₁ E₂ F) * ADMM.ρ E₁ E₂ F) • ((OptProblem.A₁ E₂) (e₁ (n + 1)) + (OptProblem.A₂ E₁) (e₂ (n + 1))) -
    ADMM.ρ E₁ E₂ F • (OptProblem.A₂ E₁) (ADMM.x₂ E₁ F n - ADMM.x₂ E₁ F (n + 1))
u' : E₁ :=
  -(ContinuousLinearMap.adjoint (OptProblem.A₁ E₂))
      (ADMM.y E₁ E₂ (n + 1) +
          ((1 - ADMM.τ E₁ E₂ F) * ADMM.ρ E₁ E₂ F) •
            ((OptProblem.A₁ E₂) (e₁ (n + 1)) + (OptProblem.A₂ E₁) (e₂ (n + 1))) +
        ADMM.ρ E₁ E₂ F • (OptProblem.A₂ E₁) (ADMM.x₂ E₁ F n - ADMM.x₂ E₁ F (n + 1)))
x_diff : E₁ := ADMM.x₁ E₂ F (n + 1) - admm_kkt.x₁
succ_n : ℕ+ := (n + 1).toPNat'
block₁ : F :=
  ADMM.y E₁ E₂ (n + 1) +
      ((1 - ADMM.τ E₁ E₂ F) * ADMM.ρ E₁ E₂ F) • ((OptProblem.A₁ E₂) (e₁ (n + 1)) + (OptProblem.A₂ E₁) (e₂ (n + 1))) +
    ADMM.ρ E₁ E₂ F • (OptProblem.A₂ E₁) (ADMM.x₂ E₁ F n - ADMM.x₂ E₁ F (n + 1))
split_block : -block = block₁ - admm_kkt.y
u'_eq : -(ContinuousLinearMap.adjoint (OptProblem.A₁ E₂)) block₁ = u'
Aty'_eq :
  (ContinuousLinearMap.adjoint (OptProblem.A₁ E₂)) admm_kkt.y =
    (ContinuousLinearMap.adjoint (OptProblem.A₁ E₂)) admm_kkt.y
⊢ ⟪(ContinuousLinearMap.adjoint (OptProblem.A₁ E₂)) admm_kkt.y -
        (ContinuousLinearMap.adjoint (OptProblem.A₁ E₂)) block₁,
      e'⟫_ℝ =
    ⟪u' + (ContinuousLinearMap.adjoint (OptProblem.A₁ E₂)) admm_kkt.y, x_diff⟫_ℝ

Tactic:
rw [← u'_eq, Aty'_eq, add_comm, sub_eq_add_neg]

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
 n : ℕ
 Ae1 : F := (OptProblem.A₁ E₂) (e₁ (n + 1))
 e' : E₁ := e₁ (n + 1)
 block : F :=
   -ey (n + 1) -
       ((1 - ADMM.τ E₁ E₂ F) * ADMM.ρ E₁ E₂ F) • ((OptProblem.A₁ E₂) (e₁ (n + 1)) + (OptProblem.A₂ E₁) (e₂ (n + 1))) -
     ADMM.ρ E₁ E₂ F • (OptProblem.A₂ E₁) (ADMM.x₂ E₁ F n - ADMM.x₂ E₁ F (n + 1))
 u' : E₁ :=
   -(ContinuousLinearMap.adjoint (OptProblem.A₁ E₂))
       (ADMM.y E₁ E₂ (n + 1) +
           ((1 - ADMM.τ E₁ E₂ F) * ADMM.ρ E₁ E₂ F) •
             ((OptProblem.A₁ E₂) (e₁ (n + 1)) + (OptProblem.A₂ E₁) (e₂ (n + 1))) +
         ADMM.ρ E₁ E₂ F • (OptProblem.A₂ E₁) (ADMM.x₂ E₁ F n - ADMM.x₂ E₁ F (n + 1)))
 x_diff : E₁ := ADMM.x₁ E₂ F (n + 1) - admm_kkt.x₁
 succ_n : ℕ+ := (n + 1).toPNat'
 block₁ : F :=
   ADMM.y E₁ E₂ (n + 1) +
       ((1 - ADMM.τ E₁ E₂ F) * ADMM.ρ E₁ E₂ F) • ((OptProblem.A₁ E₂) (e₁ (n + 1)) + (OptProblem.A₂ E₁) (e₂ (n + 1))) +
     ADMM.ρ E₁ E₂ F • (OptProblem.A₂ E₁) (ADMM.x₂ E₁ F n - ADMM.x₂ E₁ F (n + 1))
 split_block : -block = block₁ - admm_kkt.y
 u'_eq : -(ContinuousLinearMap.adjoint (OptProblem.A₁ E₂)) block₁ = u'
 Aty'_eq :
   (ContinuousLinearMap.adjoint (OptProblem.A₁ E₂)) admm_kkt.y =
     (ContinuousLinearMap.adjoint (OptProblem.A₁ E₂)) admm_kkt.y
 ⊢ ⟪(ContinuousLinearMap.adjoint (OptProblem.A₁ E₂)) admm_kkt.y +
         -(ContinuousLinearMap.adjoint (OptProblem.A₁ E₂)) block₁,
       e'⟫_ℝ =
     ⟪(ContinuousLinearMap.adjoint (OptProblem.A₁ E₂)) admm_kkt.y +
         -(ContinuousLinearMap.adjoint (OptProblem.A₁ E₂)) block₁,
       x_diff⟫_ℝ]