Documentation

LeanPool.HighContrastHomogenization.Entry.Response.Rows.AbstractCellPairingBound

High-contrast homogenization: Entry.Response.Rows.AbstractCellPairingBound #

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

The Abstract Cell-Pairing Row of the Cutoff-Mean Estimate #

For a stationary law above the threshold j_*, the annealed coarse block of a_- or a_+ is the same matrix at every aligned cell of a generation, so the flat average, over the generation's cells, of the square of each cell's own two-term annealed head collapses to the single head at the origin cell. Because the scale-s translations of the adapted grid are integral, the cell averages of the cutoff fluctuation over those cells share that common annealed value, and exact partition averaging turns their weighted average into the average of the fluctuation over the terminal cell, which vanishes since a cutoff has mean one. Combining the head collapse with the annealed Cauchy-Schwarz inequality that moves the expectation outside the flat average of the pathwise cell-pairing bound identifies the collapsed head with the depth-zero summand of L_s, assembling the abstract cell-pairing row bounded by √(L_s^∓) · √(2 tau^∓). The module serves the manuscript's p.response.transfer: it bounds the cell half of the cutoff-mean row by the source load L_s and the scale defect tau.

The flat average of the annealed heads collapses to the origin cell #

In the cell half of the cutoff-mean row of p.response.transfer the terminal cell is subdivided into the 3^{Hd} cells of the coarse generation s, and each contributes the square of the two-term head of its own annealed block. Above j_* every aligned cell of a given generation is an integer translate of the centred one and the law is translation invariant, so all those annealed blocks coincide: the flat average over the triadic index box collapses to the single head at the index 0, the depth-zero summand of the source load respSourceLoad. This is the step that makes the scale-s subdivision cost nothing.

The flat average of the annealed heads collapses, minus sign. For a stationary law P and the recentred coefficient a_- = a - g of p.response.transfer, at a scale j ≥ j_* the normalized sum over the 3^{nd} triadic subcells adaptedCellAtCenter (respGrid jStar F) j w of the squared two-term head of the corresponding annealed block is the head of the single cell at the index 0. Every subcell is an integer translate of the centred cell and the law is translation invariant, so the summand is independent of w.

The flat average of the annealed heads collapses, plus sign. For a stationary law P and the recentred coefficient a_+ = aᵀ + g of p.response.transfer, at a scale j ≥ j_* the normalized sum over the 3^{nd} triadic subcells adaptedCellAtCenter (respGrid jStar F) j w of the squared two-term head of the corresponding annealed block is the head of the single cell at the index 0. Every subcell is an integer translate of the centred cell and the law is translation invariant, so the summand is independent of w.

The cell half of the cutoff-mean row on the carriers #

The cell half of the two cutoff-mean rows of e.response.cutoff.estimate is the flat average, over the 3^{Hd} coarse cells of the terminal cell, of the cell mean of the cutoff times the crossed dual pairing of that cell's optimizer-mean defect. AK's direct full-dual pairing bounds each cell's pairing pathwise by the cell's own two-term head times the square root of twice the cell deficit; Cauchy--Schwarz with the expectation outside the flat average separates the two factors, the expectation meets the pathwise coarse blocks linearly through integral_avsum_pathwise_head_le, stationarity collapses the flat average of the annealed blocks to the single cell at the origin through avsum_sq_head_annealed_respCoeffMinus_eq, and sq_head_le_respSourceLoad identifies that collapsed head with the depth-zero summand of the source load L_s of p.response.transfer. The second factor is the annealed flat average of the cell deficits, which is twice the scale defect tau. This module assembles those landed pieces into the two carrier-level estimates, one for each sign of the recentred coefficient.

theorem HCPolySupport.HighContrast.Multiscale.abs_integral_avsum_weighted_pairing_le_respLsMinus {d : ℕ} (P : MeasureTheory.Measure (CoeffSpace d)) (hstat : IsStationaryLaw P) (jStar : ℕ) (F : BlockMat d) (H : ℕ) (s t : ℤ) (e : Vec d) (hjs : ↑jStar ≤ s) (c : (Fin d → ℤ) → ℝ) (hc : ∀ w ∈ triadicIndexBox d H, |c w| ≤ 1) (pairing D : (Fin d → ℤ) → CoeffSpace d → ℝ) (hsum : Summable (respSourceLoadSummand P jStar F s (respCoeffMinus F) (respYMinus P jStar F t e))) (hD : ∀ w ∈ triadicIndexBox d H, ∀ (a : CoeffSpace d), 0 ≤ D w a) (hbd : ∀ w ∈ triadicIndexBox d H, ∀ (a : CoeffSpace d), |pairing w a| ≤ (√(vecDot (respYMinus P jStar F t e).1 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffMinus F a)).upperLeft (respYMinus P jStar F t e).1)) + √(vecDot (respYMinus P jStar F t e).2 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffMinus F a)).lowerRight (respYMinus P jStar F t e).2))) * √(2 * D w a)) (hDval : ∫ (a : CoeffSpace d), (↑(triadicIndexBox d H).card)⁻¹ * ∑ w ∈ triadicIndexBox d H, 2 * D w a ∂P = 2 * respTauMinus P jStar F s t e) (hUL : ∀ w ∈ triadicIndexBox d H, ∀ (a : CoeffSpace d), 0 ≤ vecDot (respYMinus P jStar F t e).1 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffMinus F a)).upperLeft (respYMinus P jStar F t e).1)) (hLR : ∀ w ∈ triadicIndexBox d H, ∀ (a : CoeffSpace d), 0 ≤ vecDot (respYMinus P jStar F t e).2 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffMinus F a)).lowerRight (respYMinus P jStar F t e).2)) (hentUL : ∀ w ∈ triadicIndexBox d H, ∀ (i j : Fin d), MeasureTheory.Integrable (fun (a : CoeffSpace d) => (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffMinus F a)).upperLeft i j) P) (hentLR : ∀ w ∈ triadicIndexBox d H, ∀ (i j : Fin d), MeasureTheory.Integrable (fun (a : CoeffSpace d) => (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffMinus F a)).lowerRight i j) P) (hqUL : ∀ w ∈ triadicIndexBox d H, MeasureTheory.Integrable (fun (a : CoeffSpace d) => vecDot (respYMinus P jStar F t e).1 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffMinus F a)).upperLeft (respYMinus P jStar F t e).1)) P) (hqLR : ∀ w ∈ triadicIndexBox d H, MeasureTheory.Integrable (fun (a : CoeffSpace d) => vecDot (respYMinus P jStar F t e).2 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffMinus F a)).lowerRight (respYMinus P jStar F t e).2)) P) (hcross : ∀ w ∈ triadicIndexBox d H, MeasureTheory.Integrable (fun (a : CoeffSpace d) => √(vecDot (respYMinus P jStar F t e).1 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffMinus F a)).upperLeft (respYMinus P jStar F t e).1)) * √(vecDot (respYMinus P jStar F t e).2 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffMinus F a)).lowerRight (respYMinus P jStar F t e).2))) P) (hsq : ∀ w ∈ triadicIndexBox d H, MeasureTheory.Integrable (fun (a : CoeffSpace d) => (√(vecDot (respYMinus P jStar F t e).1 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffMinus F a)).upperLeft (respYMinus P jStar F t e).1)) + √(vecDot (respYMinus P jStar F t e).2 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffMinus F a)).lowerRight (respYMinus P jStar F t e).2))) ^ 2) P) (hmeas : ∀ (i k : Fin d), MeasureTheory.AEStronglyMeasurable (fun (a : CoeffSpace d) => (coarseBlockMatrix (adaptedCell (respGrid jStar F) s) (respCoeffMinus F a)).upperLeft i k) P) (hmeas' : ∀ (i k : Fin d), MeasureTheory.AEStronglyMeasurable (fun (a : CoeffSpace d) => (coarseBlockMatrix (adaptedCell (respGrid jStar F) s) (respCoeffMinus F a)).lowerRight i k) P) (hFint : MeasureTheory.Integrable (fun (a : CoeffSpace d) => (↑(triadicIndexBox d H).card)⁻¹ * ∑ w ∈ triadicIndexBox d H, (√(vecDot (respYMinus P jStar F t e).1 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffMinus F a)).upperLeft (respYMinus P jStar F t e).1)) + √(vecDot (respYMinus P jStar F t e).2 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffMinus F a)).lowerRight (respYMinus P jStar F t e).2))) ^ 2) P) (hDint : MeasureTheory.Integrable (fun (a : CoeffSpace d) => (↑(triadicIndexBox d H).card)⁻¹ * ∑ w ∈ triadicIndexBox d H, 2 * D w a) P) (hMidint : MeasureTheory.Integrable (fun (a : CoeffSpace d) => √((↑(triadicIndexBox d H).card)⁻¹ * ∑ w ∈ triadicIndexBox d H, (√(vecDot (respYMinus P jStar F t e).1 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffMinus F a)).upperLeft (respYMinus P jStar F t e).1)) + √(vecDot (respYMinus P jStar F t e).2 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffMinus F a)).lowerRight (respYMinus P jStar F t e).2))) ^ 2) * √((↑(triadicIndexBox d H).card)⁻¹ * ∑ w ∈ triadicIndexBox d H, 2 * D w a)) P) (hPint : MeasureTheory.Integrable (fun (a : CoeffSpace d) => (↑(triadicIndexBox d H).card)⁻¹ * ∑ w ∈ triadicIndexBox d H, c w * (√((√(vecDot (respYMinus P jStar F t e).1 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffMinus F a)).upperLeft (respYMinus P jStar F t e).1)) + √(vecDot (respYMinus P jStar F t e).2 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffMinus F a)).lowerRight (respYMinus P jStar F t e).2))) ^ 2) * √(2 * D w a))) P) (hPint' : MeasureTheory.Integrable (fun (a : CoeffSpace d) => (↑(triadicIndexBox d H).card)⁻¹ * ∑ w ∈ triadicIndexBox d H, c w * pairing w a) P) :
|∫ (a : CoeffSpace d), (↑(triadicIndexBox d H).card)⁻¹ * ∑ w ∈ triadicIndexBox d H, c w * pairing w a ∂P| ≤ √(respLsMinus P jStar F s t e) * √(2 * respTauMinus P jStar F s t e)

The cell half of the cutoff-mean row, minus sign. On the carriers of p.response.transfer, let pairing be the full-dual pairing of the optimizer-mean defect and let D be the scalar cell deficit. If each cell's pairing is controlled pathwise by the cell's own two-term head √(P · b_s^{1/2} P) + √(Q · (S_{*,s})^{-1/2} Q) of the pathwise coarse block of the recentred coefficient a_- = a - g times √(2 D), and the annealed flat average of 2 D is 2 tau^-, then, provided the pathwise coarse blocks are nonnegative, integrable and measurable as prescribed, the modulus of the sample expectation of the cutoff-weighted flat average of the pairing is at most √(L_s^-) * √(2 tau^-), where L_s^- is the source load of p.response.transfer at the recentred coefficient a_- and the dual variable Y^-. This is the cell half of the cutoff-mean row, with the expectation outside the flat average and no inversion of an annealed block and no Jensen step.

theorem HCPolySupport.HighContrast.Multiscale.abs_integral_avsum_weighted_pairing_le_respLsPlus {d : ℕ} (P : MeasureTheory.Measure (CoeffSpace d)) (hstat : IsStationaryLaw P) (jStar : ℕ) (F : BlockMat d) (H : ℕ) (s t : ℤ) (e : Vec d) (hjs : ↑jStar ≤ s) (c : (Fin d → ℤ) → ℝ) (hc : ∀ w ∈ triadicIndexBox d H, |c w| ≤ 1) (pairing D : (Fin d → ℤ) → CoeffSpace d → ℝ) (hsum : Summable (respSourceLoadSummand P jStar F s (respCoeffPlus F) (respYPlus P jStar F t e))) (hD : ∀ w ∈ triadicIndexBox d H, ∀ (a : CoeffSpace d), 0 ≤ D w a) (hbd : ∀ w ∈ triadicIndexBox d H, ∀ (a : CoeffSpace d), |pairing w a| ≤ (√(vecDot (respYPlus P jStar F t e).1 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffPlus F a)).upperLeft (respYPlus P jStar F t e).1)) + √(vecDot (respYPlus P jStar F t e).2 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffPlus F a)).lowerRight (respYPlus P jStar F t e).2))) * √(2 * D w a)) (hDval : ∫ (a : CoeffSpace d), (↑(triadicIndexBox d H).card)⁻¹ * ∑ w ∈ triadicIndexBox d H, 2 * D w a ∂P = 2 * respTauPlus P jStar F s t e) (hUL : ∀ w ∈ triadicIndexBox d H, ∀ (a : CoeffSpace d), 0 ≤ vecDot (respYPlus P jStar F t e).1 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffPlus F a)).upperLeft (respYPlus P jStar F t e).1)) (hLR : ∀ w ∈ triadicIndexBox d H, ∀ (a : CoeffSpace d), 0 ≤ vecDot (respYPlus P jStar F t e).2 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffPlus F a)).lowerRight (respYPlus P jStar F t e).2)) (hentUL : ∀ w ∈ triadicIndexBox d H, ∀ (i j : Fin d), MeasureTheory.Integrable (fun (a : CoeffSpace d) => (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffPlus F a)).upperLeft i j) P) (hentLR : ∀ w ∈ triadicIndexBox d H, ∀ (i j : Fin d), MeasureTheory.Integrable (fun (a : CoeffSpace d) => (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffPlus F a)).lowerRight i j) P) (hqUL : ∀ w ∈ triadicIndexBox d H, MeasureTheory.Integrable (fun (a : CoeffSpace d) => vecDot (respYPlus P jStar F t e).1 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffPlus F a)).upperLeft (respYPlus P jStar F t e).1)) P) (hqLR : ∀ w ∈ triadicIndexBox d H, MeasureTheory.Integrable (fun (a : CoeffSpace d) => vecDot (respYPlus P jStar F t e).2 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffPlus F a)).lowerRight (respYPlus P jStar F t e).2)) P) (hcross : ∀ w ∈ triadicIndexBox d H, MeasureTheory.Integrable (fun (a : CoeffSpace d) => √(vecDot (respYPlus P jStar F t e).1 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffPlus F a)).upperLeft (respYPlus P jStar F t e).1)) * √(vecDot (respYPlus P jStar F t e).2 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffPlus F a)).lowerRight (respYPlus P jStar F t e).2))) P) (hsq : ∀ w ∈ triadicIndexBox d H, MeasureTheory.Integrable (fun (a : CoeffSpace d) => (√(vecDot (respYPlus P jStar F t e).1 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffPlus F a)).upperLeft (respYPlus P jStar F t e).1)) + √(vecDot (respYPlus P jStar F t e).2 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffPlus F a)).lowerRight (respYPlus P jStar F t e).2))) ^ 2) P) (hmeas : ∀ (i k : Fin d), MeasureTheory.AEStronglyMeasurable (fun (a : CoeffSpace d) => (coarseBlockMatrix (adaptedCell (respGrid jStar F) s) (respCoeffPlus F a)).upperLeft i k) P) (hmeas' : ∀ (i k : Fin d), MeasureTheory.AEStronglyMeasurable (fun (a : CoeffSpace d) => (coarseBlockMatrix (adaptedCell (respGrid jStar F) s) (respCoeffPlus F a)).lowerRight i k) P) (hFint : MeasureTheory.Integrable (fun (a : CoeffSpace d) => (↑(triadicIndexBox d H).card)⁻¹ * ∑ w ∈ triadicIndexBox d H, (√(vecDot (respYPlus P jStar F t e).1 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffPlus F a)).upperLeft (respYPlus P jStar F t e).1)) + √(vecDot (respYPlus P jStar F t e).2 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffPlus F a)).lowerRight (respYPlus P jStar F t e).2))) ^ 2) P) (hDint : MeasureTheory.Integrable (fun (a : CoeffSpace d) => (↑(triadicIndexBox d H).card)⁻¹ * ∑ w ∈ triadicIndexBox d H, 2 * D w a) P) (hMidint : MeasureTheory.Integrable (fun (a : CoeffSpace d) => √((↑(triadicIndexBox d H).card)⁻¹ * ∑ w ∈ triadicIndexBox d H, (√(vecDot (respYPlus P jStar F t e).1 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffPlus F a)).upperLeft (respYPlus P jStar F t e).1)) + √(vecDot (respYPlus P jStar F t e).2 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffPlus F a)).lowerRight (respYPlus P jStar F t e).2))) ^ 2) * √((↑(triadicIndexBox d H).card)⁻¹ * ∑ w ∈ triadicIndexBox d H, 2 * D w a)) P) (hPint : MeasureTheory.Integrable (fun (a : CoeffSpace d) => (↑(triadicIndexBox d H).card)⁻¹ * ∑ w ∈ triadicIndexBox d H, c w * (√((√(vecDot (respYPlus P jStar F t e).1 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffPlus F a)).upperLeft (respYPlus P jStar F t e).1)) + √(vecDot (respYPlus P jStar F t e).2 (matVecMul (coarseBlockMatrix (adaptedCellAtCenter (respGrid jStar F) s w) (respCoeffPlus F a)).lowerRight (respYPlus P jStar F t e).2))) ^ 2) * √(2 * D w a))) P) (hPint' : MeasureTheory.Integrable (fun (a : CoeffSpace d) => (↑(triadicIndexBox d H).card)⁻¹ * ∑ w ∈ triadicIndexBox d H, c w * pairing w a) P) :
|∫ (a : CoeffSpace d), (↑(triadicIndexBox d H).card)⁻¹ * ∑ w ∈ triadicIndexBox d H, c w * pairing w a ∂P| ≤ √(respLsPlus P jStar F s t e) * √(2 * respTauPlus P jStar F s t e)

The cell half of the cutoff-mean row, plus sign. This is the Plus twin of abs_integral_avsum_weighted_pairing_le_respLsMinus: the same flat-average estimate with the recentred coefficient a_+ = aᵀ + g, the dual variable Y^+, the annealed head collapse avsum_sq_head_annealed_respCoeffPlus_eq, the scale defect tau^+ and the source load L_s^+ of p.response.transfer.

Cancellation of constant cell means after expectation #

In the centred cutoff decomposition of p.response.transfer, the constant cell means cancel after expectation. The scale-s translations of the grid are integral, so the cell averages of the cutoff fluctuation share a common annealed value. Exact partition averaging turns their average over the generation-(t - n) triadic subdivision of the terminal cell into the average of φ - 1 over that cell, and that average vanishes because the cutoff class has volume average one. The same partition then shows that any weighting of a family of cell functionals with a common annealed value and mean-zero normalized weights integrates to zero.

theorem HCPolySupport.HighContrast.Multiscale.avsum_volumeAverage_sub_one_eq_zero {d : ℕ} [NeZero d] {qq : Mat d} (hq : IsUnit qq) (t : ℤ) (n : ℕ) {φ : Vec d → ℝ} (hφ : IsResponseCutoff qq t φ) (hint : MeasureTheory.IntegrableOn (fun (x : Vec d) => φ x - 1) (adaptedCell qq t) MeasureTheory.volume) (hvol : (MeasureTheory.volume (adaptedCell qq t)).toReal ≠ 0) (hφint : MeasureTheory.IntegrableOn φ (adaptedCell qq t) MeasureTheory.volume) :
((↑(triadicIndexBox d n).card)⁻¹ * ∑ w ∈ triadicIndexBox d n, volumeAverage (adaptedCellAtCenter qq (t - ↑n) w) fun (x : Vec d) => φ x - 1) = 0

Cancellation of the cell means of the cutoff fluctuation. The normalized sum of the cell averages of φ - 1 over the 3 ^ (n * d) depth-n triadic subcells of the adapted cell HighContrast.adaptedCell qq t equals the average of φ - 1 over that cell, which is zero because φ has volume average one there. This is the first cancellation in p.response.transfer.

theorem HCPolySupport.HighContrast.Multiscale.integral_avsum_mul_eq_zero {d : ℕ} {α : Type u_1} [MeasurableSpace α] (P : MeasureTheory.Measure α) (n : ℕ) (c : (Fin d → ℤ) → ℝ) (G : (Fin d → ℤ) → α → ℝ) (m : ℝ) (hc : (↑(triadicIndexBox d n).card)⁻¹ * ∑ w ∈ triadicIndexBox d n, c w = 0) (hG : ∀ w ∈ triadicIndexBox d n, ∫ (a : α), G w a ∂P = m) (hint : ∀ w ∈ triadicIndexBox d n, MeasureTheory.Integrable (G w) P) :
∫ (a : α), (↑(triadicIndexBox d n).card)⁻¹ * ∑ w ∈ triadicIndexBox d n, c w * G w a ∂P = 0

Mean-zero weighting of a common annealed value. If each cell functional G w has the same annealed value m and the weights c w have normalized sum zero, then the annealed integral of the weighted sum vanishes. The common value is the one supplied by the integrality of the scale-s translations; this is the second cancellation in p.response.transfer.