constant_config_total_defect

formal source

 127  rw [Finset.sum_const, Finset.card_univ, Fintype.card_fin, nsmul_eq_mul]
 128
 129/-- The constant configuration has total defect `N * Jlog μ`. -/
 130private theorem constant_config_total_defect {N : ℕ} (μ : ℝ) :
 131    total_defect (constant_config μ : Configuration N) = (N : ℝ) * Jlog μ := by
 132  unfold total_defect constant_config
 133  simp only [Finset.sum_const, Finset.card_univ, Fintype.card_fin, nsmul_eq_mul]
 134  rfl
 135
 136/-- Weighted average of the logs equals `log_charge / N`. -/
 137private theorem weighted_log_average {N : ℕ} (hN : 0 < N) (c : Configuration N) :
 138    (∑ i ∈ (Finset.univ : Finset (Fin N)), (1 / (N : ℝ)) * Real.log (c.entries i)) =
 139      log_charge c / N := by
 140  unfold log_charge
 141  rw [← Finset.mul_sum]
 142  ring
 143
 144/-- Weighted average of `Jlog(log x_i)` equals `total_defect / N`. -/
 145private theorem weighted_Jlog_average {N : ℕ} (c : Configuration N) :
 146    (∑ i ∈ (Finset.univ : Finset (Fin N)), (1 / (N : ℝ)) * Jlog (Real.log (c.entries i))) =
 147      (1 / (N : ℝ)) * total_defect c := by
 148  calc
 149    (∑ i ∈ (Finset.univ : Finset (Fin N)), (1 / (N : ℝ)) * Jlog (Real.log (c.entries i)))
 150        = ∑ i ∈ (Finset.univ : Finset (Fin N)), (1 / (N : ℝ)) * defect (c.entries i) := by
 151            apply Finset.sum_congr rfl
 152            intro i _
 153            unfold Jlog defect J Jcost
 154            rw [Real.exp_log (c.entries_pos i)]
 155    _ = (1 / (N : ℝ)) * total_defect c := by
 156          unfold total_defect
 157          rw [← Finset.mul_sum]
 158
 159/-- Jensen lower bound: fixed log-charge implies a defect lower bound. -/
 160private theorem total_defect_lower_bound {N : ℕ} (hN : 0 < N) (c : Configuration N) :