import Mathlib
import Optlib


open Set Bornology Filter BigOperators Topology BCD

open BCD

theorem extracted_formal_statement_8.{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} [inst_6 : ProperSpace E] [inst_7 : ProperSpace F] (γ : ℝ)
  (hγ : γ > 1) (ck : ∀ (k : ℕ), alg.c k = 1 / (γ * ↑l)) (dk : ∀ (k : ℕ), alg.d k = 1 / (γ * ↑l))
  (bd : Bornology.IsBounded (z '' univ)) (lbdψ : BddBelow (ψ '' univ)) (a b : Set (WithLp 2 (E × F)))
  (closea : IsClosed a) (closeb : IsClosed b) (sub_ab : limit_set z ⊆ a ∪ b) (nez_a : (limit_set z ∩ a).Nonempty)
  (nez_b : (limit_set z ∩ b).Nonempty) (nz_ab : limit_set z ∩ (a ∩ b) = ∅) :
  let A := limit_set z ∩ a;
  let B := limit_set z ∩ b;
  A ∪ B = limit_set z →
    A ∩ B = ∅ →
      ∀ (z : WithLp 2 (E × F)),
        (EMetric.infEdist z A).toReal + (EMetric.infEdist z B).toReal = 0 → z ∈ A → IsClosed A := 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
inst✝¹ : ProperSpace E
inst✝ : ProperSpace F
γ : ℝ
hγ : γ > 1
ck : ∀ (k : ℕ), alg.c k = 1 / (γ * ↑l)
dk : ∀ (k : ℕ), alg.d k = 1 / (γ * ↑l)
bd : Bornology.IsBounded (BCD.z '' univ)
lbdψ : BddBelow (ψ '' univ)
a b : Set (WithLp 2 (E × F))
closea : IsClosed a
closeb : IsClosed b
sub_ab : limit_set BCD.z ⊆ a ∪ b
nez_a : (limit_set BCD.z ∩ a).Nonempty
nez_b : (limit_set BCD.z ∩ b).Nonempty
nz_ab : limit_set BCD.z ∩ (a ∩ b) = ∅
A : Set (WithLp 2 (E × F)) := limit_set BCD.z ∩ a
B : Set (WithLp 2 (E × F)) := limit_set BCD.z ∩ b
eqAB : A ∪ B = limit_set BCD.z
disjoint_AB : A ∩ B = ∅
ω : WithLp 2 (E × F) → ℝ :=
  fun z => (EMetric.infEdist z A).toReal / ((EMetric.infEdist z A).toReal + (EMetric.infEdist z B).toReal)
z : WithLp 2 (E × F)
eq0 : (EMetric.infEdist z A).toReal + (EMetric.infEdist z B).toReal = 0
⊢ z ∈ A

Tactic:
apply EMetric.mem_closure_iff_infEdist_zero.mpr
    have : (EMetric.infEdist z A).toReal = 0 := by
  linarith [eq0, @ENNReal.toReal_nonneg (EMetric.infEdist z A), @ENNReal.toReal_nonneg (EMetric.infEdist z B)]
    exact
  (((fun {x y} hx hy ↦ (ENNReal.toReal_eq_toReal_iff' hx hy).mp) ENNReal.top_ne_zero.symm (Metric.infEdist_ne_top nez_a)
      (id (Eq.symm this)))).symm
    simp
    constructor

Tactic state after:
[case a.left
 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
 inst✝¹ : ProperSpace E
 inst✝ : ProperSpace F
 γ : ℝ
 hγ : γ > 1
 ck : ∀ (k : ℕ), alg.c k = 1 / (γ * ↑l)
 dk : ∀ (k : ℕ), alg.d k = 1 / (γ * ↑l)
 bd : Bornology.IsBounded (BCD.z '' univ)
 lbdψ : BddBelow (ψ '' univ)
 a b : Set (WithLp 2 (E × F))
 closea : IsClosed a
 closeb : IsClosed b
 sub_ab : limit_set BCD.z ⊆ a ∪ b
 nez_a : (limit_set BCD.z ∩ a).Nonempty
 nez_b : (limit_set BCD.z ∩ b).Nonempty
 nz_ab : limit_set BCD.z ∩ (a ∩ b) = ∅
 A : Set (WithLp 2 (E × F)) := limit_set BCD.z ∩ a
 B : Set (WithLp 2 (E × F)) := limit_set BCD.z ∩ b
 eqAB : A ∪ B = limit_set BCD.z
 disjoint_AB : A ∩ B = ∅
 ω : WithLp 2 (E × F) → ℝ :=
   fun z => (EMetric.infEdist z A).toReal / ((EMetric.infEdist z A).toReal + (EMetric.infEdist z B).toReal)
 z : WithLp 2 (E × F)
 eq0 : (EMetric.infEdist z A).toReal + (EMetric.infEdist z B).toReal = 0
 ⊢ IsClosed A,
 case a.right
 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
 inst✝¹ : ProperSpace E
 inst✝ : ProperSpace F
 γ : ℝ
 hγ : γ > 1
 ck : ∀ (k : ℕ), alg.c k = 1 / (γ * ↑l)
 dk : ∀ (k : ℕ), alg.d k = 1 / (γ * ↑l)
 bd : Bornology.IsBounded (BCD.z '' univ)
 lbdψ : BddBelow (ψ '' univ)
 a b : Set (WithLp 2 (E × F))
 closea : IsClosed a
 closeb : IsClosed b
 sub_ab : limit_set BCD.z ⊆ a ∪ b
 nez_a : (limit_set BCD.z ∩ a).Nonempty
 nez_b : (limit_set BCD.z ∩ b).Nonempty
 nz_ab : limit_set BCD.z ∩ (a ∩ b) = ∅
 A : Set (WithLp 2 (E × F)) := limit_set BCD.z ∩ a
 B : Set (WithLp 2 (E × F)) := limit_set BCD.z ∩ b
 eqAB : A ∪ B = limit_set BCD.z
 disjoint_AB : A ∩ B = ∅
 ω : WithLp 2 (E × F) → ℝ :=
   fun z => (EMetric.infEdist z A).toReal / ((EMetric.infEdist z A).toReal + (EMetric.infEdist z B).toReal)
 z : WithLp 2 (E × F)
 eq0 : (EMetric.infEdist z A).toReal + (EMetric.infEdist z B).toReal = 0
 ⊢ limit_set BCD.z ∩ a ⊆ A]