def uniformNetwork (J : ℕ) (ρ : ℝ) : ComplementarityNetwork J where
w := fun j k => if j = k then 0 else (1 - ρ)
symmetric := fun j k => by
by_cases h : j = k
· simp [h]
· simp [h, Ne.symm h]
zero_diag := fun j => by simpCES on Networks: Heterogeneous Pairwise Complementarity