Documentation

LeanPool.HighContrastHomogenization.Entry.Multiscale.OneGrid.ExponentAndLossArithmetic

High-contrast homogenization: Entry.Multiscale.OneGrid.ExponentAndLossArithmetic #

Imported from the Apache-2.0 HighContrastHomogenization development at commit 7a13dbcd8d6609264a713373f5c69ceeac870472.

Exponent And Loss Arithmetic #

This module proves the exponent arithmetic of the one-grid propagation and the algebra of the log-determinant loss it is applied to. It establishes the numeric facts about Q = 𝒬(d,Ξ³) and the contraction, span and averaging coefficients, the elementary exponential inequalities and the geometric profile weights; and it proves additivity of the log-determinant loss, its comparison with the synchronized loss, and the synchronized multiplicity bound. It serves p.fixed.geometry.one.grid.propagation.

Printed exponent and coefficient arithmetic #

Group A of the printed proof (p.fixed.geometry.one.grid.propagation): the numeric facts about Q = 𝒬(d,Ξ³), the contraction, span and averaging coefficients, the elementary exponential inequalities and the geometric profile weights.

Part of the proof of the statement in HCPoly/Entry/Statements/OneGridPropagation.lean. Conventions of the group: q = 𝒬(π”ͺ) is Geometry.explicitRoundedGrid jStar metric at every loss, history, profile and drift; metric (not m) names the positive matrix, because m, n, mβ‚€ are generations; the source lower scale e.source.lower.scale is carried exactly by the module HCPoly/Entry/OneGridPropagation.lean; this group's stronger internal algebraic and history lemmas omit an unused threshold, and this group's generic integration helpers take finiteness, integrability or measurability inputs that are proved at their actual use sites, not extra premises of the printed proposition.

Group A. The printed exponent and coefficient arithmetic #

(p.fixed.geometry.one.grid.propagation)

theorem HCPolySupport.HighContrast.Multiscale.bigQ_twelve_le (d : β„•) (hd : 2 ≀ d) (Ξ³ : ℝ) (hΞ³ : Ξ³ ∈ Set.Ico 0 1) :
12 ≀ bigQ d Ξ³

Q β‰₯ 12 for d β‰₯ 2 and Ξ³ ∈ [0,1) (p.fixed.geometry.one.grid.propagation).

theorem HCPolySupport.HighContrast.Multiscale.one_sub_gamma_mul_span_ge (d : β„•) (hd : 2 ≀ d) (Ξ³ : ℝ) (hΞ³ : Ξ³ ∈ Set.Ico 0 1) (h : β„•) (hh : 2 * bigQ d Ξ³ ≀ h) :
8 * (↑d + 1) ≀ (1 - Ξ³) * ↑h

(1-Ξ³)h β‰₯ 2(1-Ξ³)Q β‰₯ 8(d+1) for every natural h β‰₯ 2Q (p.fixed.geometry.one.grid.propagation).

theorem HCPolySupport.HighContrast.Multiscale.three_rpow_neg_eighth_span_lt (d : β„•) (hd : 2 ≀ d) (Ξ³ : ℝ) (hΞ³ : Ξ³ ∈ Set.Ico 0 1) (h : β„•) (hh : 2 * bigQ d Ξ³ ≀ h) :
3 ^ (-((1 - Ξ³) / 8) * ↑h) < 1 / 8

3^{-β…›(1-Ξ³)h} ≀ 3^{-3} < 1/8 for every natural h β‰₯ 2Q (p.fixed.geometry.one.grid.propagation). This is the drift half of the contraction.

theorem HCPolySupport.HighContrast.Multiscale.two_bigQ_mul_three_rpow_le_third (d : β„•) (hd : 2 ≀ d) (Ξ³ : ℝ) (hΞ³ : Ξ³ ∈ Set.Ico 0 1) :
2 * ↑(bigQ d Ξ³) * 3 ^ (-↑d * (↑(bigQ d Ξ³) - 1)) ≀ 1 / 3

2QΒ·3^{-d(Q-1)} ≀ 1/3 (p.fixed.geometry.one.grid.propagation).

theorem HCPolySupport.HighContrast.Multiscale.contraction_coefficient_lt (d : β„•) (hd : 2 ≀ d) (Ξ³ : ℝ) (hΞ³ : Ξ³ ∈ Set.Ico 0 1) (h : β„•) (hh : 2 * bigQ d Ξ³ ≀ h) :
3 ^ (-((1 - Ξ³) / 4) * ↑h) + 2 ^ (bigQ d Ξ³ - 1) * (3 ^ d * ↑(bigQ d Ξ³)) ^ bigQ d Ξ³ * 3 ^ (-(↑(bigQ d Ξ³) * ↑d / 2) * ↑h) < 1 / 8

The printed contraction coefficient 3^{-ΒΌ(1-Ξ³)h} + 2^{Q-1}(3^dQ)^Q 3^{-Qdh/2} ≀ 3^{-6} + Β½Β·3^{-12} < 1/8, for every natural h β‰₯ 2Q (p.fixed.geometry.one.grid.propagation).

theorem HCPolySupport.HighContrast.Multiscale.span_coefficient_le (d : β„•) (hd : 2 ≀ d) (Ξ³ : ℝ) (hΞ³ : Ξ³ ∈ Set.Ico 0 1) :
2 ^ (bigQ d Ξ³ - 1) * (3 ^ d * ↑(bigQ d Ξ³)) ^ bigQ d Ξ³ * 3 ^ (-↑d * ↑(bigQ d Ξ³) ^ 2) ≀ 1 / 8

The span-2Q coefficient 2^{Q-1}(3^dQ)^Q 3^{-dQΒ²} ≀ 1/8 used in Step 4, where the recurrence is applied with the fixed span 2Q so that the constants do not depend on the possibly larger h (p.fixed.geometry.one.grid.propagation).

theorem HCPolySupport.HighContrast.Multiscale.averaging_coefficient_le (d : β„•) (hd : 2 ≀ d) (Ξ³ : ℝ) (hΞ³ : Ξ³ ∈ Set.Ico 0 1) :
↑(bigQ d Ξ³) * 3 ^ (↑d / 2) * (1 + (2 * ↑d) ^ (↑(bigQ d Ξ³))⁻¹) ≀ 3 ^ d * ↑(bigQ d Ξ³)

QΒ·3^{d/2}(1 + (2d)^{1/Q}) ≀ 3^d Q for d β‰₯ 2, Q β‰₯ 2 (p.fixed.geometry.one.grid.propagation): the step that turns the recurrence's printed prefactor into the contraction's.

theorem HCPolySupport.HighContrast.Multiscale.exp_pred_mul_sub_one_le (N : β„•) (hN : 1 ≀ N) (t : ℝ) (ht : 0 ≀ t) :
Real.exp ((↑N - 1) * t) * (Real.exp t - 1) ≀ Real.exp (↑N * t) - 1

e^{(N-1)t}(e^t - 1) ≀ e^{Nt} - 1 for t β‰₯ 0 (p.fixed.geometry.one.grid.propagation).

e^x a ≀ a + e^x - 1 for 0 ≀ a ≀ 1 and x β‰₯ 0 (p.fixed.geometry.one.grid.propagation).

theorem HCPolySupport.HighContrast.Multiscale.profile_geometric_weight_sum_le (Ξ³ : ℝ) (hΞ³ : Ξ³ ∈ Set.Ico 0 1) (a b : β„€) :
βˆ‘ j ∈ Finset.Ico a b, 3 ^ (-((1 - Ξ³) / 4) * (↑b - 1 - ↑j)) ≀ (1 - 3 ^ (-((1 - Ξ³) / 4)))⁻¹

The geometric tail of the profile weights is bounded by a constant depending only on Ξ³ (p.fixed.geometry.one.grid.propagation, "the geometric sum is bounded by C").

Log-determinant loss algebra and the multiplicity bound #

Group B of the printed proof (p.fixed.geometry.one.grid.propagation): additivity of the log-determinant loss, its comparison with the synchronized loss, and conjunct 7, e.fixed.geometry.synchronized.multiplicity.

Part of the proof of HCPoly/Entry/Statements/OneGridPropagation.lean. Conventions of this group: q = 𝒬(π”ͺ) is Geometry.explicitRoundedGrid jStar metric at every loss, history, profile and drift; metric (not m) names the positive matrix, because m, n, mβ‚€ are generations; the source lower scale e.source.lower.scale is carried exactly by the HCPoly/Entry/OneGridPropagation.lean; this group's stronger internal algebraic and history lemmas omit an unused threshold, and this group's generic integration helpers take finiteness, integrability or measurability inputs that are proved at their actual use sites, not extra premises of the printed proposition.

Group B. Log-determinant loss algebra and the multiplicity bound #

(p.fixed.geometry.one.grid.propagation)

theorem HCPolySupport.HighContrast.Multiscale.logDetLoss_le_synchronizedLogDetLoss (d : β„•) (hd : 2 ≀ d) (Ξ³ : ℝ) (_hΞ³ : Ξ³ ∈ Set.Ico 0 1) (P : MeasureTheory.Measure (CoeffSpace d)) (E : BlockMat d) (Ξ¨ : ℝ β†’ ℝ) (K : ℝ) (S : CoeffSpace d β†’ ℝ) :
MeasureTheory.IsProbabilityMeasure P β†’ IsStationaryLaw P β†’ IsUnitRangeLaw P β†’ CoarseEllipticityDagger P Ξ³ E Ξ¨ K S β†’ βˆ€ (h : β„•), 1 ≀ h β†’ βˆ€ (jStar : β„•), 2 * d ≀ 3 ^ jStar β†’ βˆ€ (metric : Mat d), Matrix.PosDef metric β†’ βˆ€ (m j : β„€), ↑jStar + ↑h ≀ m β†’ m < j β†’ j ≀ m + ↑h β†’ detIncrement P (Geometry.explicitRoundedGrid jStar metric) (j - ↑h) (m + ↑h) ≀ synchCharge P (Geometry.explicitRoundedGrid jStar metric) (↑h) m

Ξ”^q_{j-h,m+h} ≀ Ξ”Μ‚^q_h(m) for m < j ≀ m+h (p.fixed.geometry.one.grid.propagation). Taking j = m+h gives the Ξ”^q_{m,m+h} ≀ Ξ”Μ‚^q_h(m) used in Step 3.

theorem HCPolySupport.HighContrast.Multiscale.synchronized_multiplicity (d : β„•) (hd : 2 ≀ d) (Ξ³ : ℝ) (hΞ³ : Ξ³ ∈ Set.Ico 0 1) (P : MeasureTheory.Measure (CoeffSpace d)) (E : BlockMat d) (Ξ¨ : ℝ β†’ ℝ) (K : ℝ) (S : CoeffSpace d β†’ ℝ) :
MeasureTheory.IsProbabilityMeasure P β†’ IsStationaryLaw P β†’ IsUnitRangeLaw P β†’ CoarseEllipticityDagger P Ξ³ E Ξ¨ K S β†’ βˆ€ (h : β„•), 2 * bigQ d Ξ³ ≀ h β†’ βˆ€ (jStar : β„•), 2 * d ≀ 3 ^ jStar β†’ βˆ€ (metric : Mat d), Matrix.PosDef metric β†’ βˆ€ (mβ‚€ : β„€), ↑jStar + ↑h ≀ mβ‚€ β†’ βˆ€ (Ksteps : β„•), 1 ≀ Ksteps β†’ βˆ‘ k ∈ Finset.range Ksteps, synchCharge P (Geometry.explicitRoundedGrid jStar metric) (↑h) (mβ‚€ + ↑k * ↑h) ≀ ↑h * detIncrement P (Geometry.explicitRoundedGrid jStar metric) (mβ‚€ + 1 - ↑h) (mβ‚€ + ↑Ksteps * ↑h)

Conjunct 7, e.fixed.geometry.synchronized.multiplicity (p.fixed.geometry.one.grid.propagation): expand the synchronized losses into adjacent increments, observe that the smallest index is at least j_* and that each increment occurs at most h times.