Documentation

LeanPool.HighContrastHomogenization.Provider.Response.ProfileWeakFullLp

High-contrast homogenization: Provider.Response.ProfileWeakFullLp #

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

Full annealed-centered weak estimates #

The random optimizer average is replaced by its annealed center before the pointwise diagonal estimate is integrated. A sample map is retained so the result applies directly to constant-skew recentering.

theorem HCPolySupport.HighContrast.Response.eLpNorm_profilePrimalWeakRoot_sample_le {d : ℕ} [NeZero d] {P : MeasureTheory.Measure (CoeffSpace d)} {q : Mat d} (hq : Matrix.PosDef q) (t : ℤ) (H : ℕ) {m : Mat d} (hm : Matrix.PosDef m) {E : BlockMat d} (hE : IsSymmetricBlockMat E) (hEpd : Book.Ch02.BlockPosDef E) {rho alpha : ℝ} (hrho : 0 < rho) (hrho1 : rho < 1) (halpha : alpha = (1 - rho) / 2) (sample : CoeffSpace d → CoeffSpace d) (p r : Vec d) (hU : AEMeasurable (fun (a : CoeffSpace d) => diagonalWeakCellSum q t H (1 / 2) E (sample a)) P) (hV : AEMeasurable (fun (a : CoeffSpace d) => diagonalWeakAverageSum q t H (1 / 2) rho E (sample a)) P) (hM : AEMeasurable (fun (a : CoeffSpace d) => diagonalWeakMaximum rho q t E (sample a)) P) (hEn : AEMeasurable (fun (a : CoeffSpace d) => diagonalWeakEnergy hq t (sample a) p r) P) (havg : AEMeasurable (fun (a : CoeffSpace d) => blockCellAverage (adaptedCell q t) (diagonalWeakState hq t (sample a) p r)) P) (hfinite : ∀ᵐ (a : CoeffSpace d) ∂P, diagonalWeakMaximum rho q t E (sample a) ≠ ⊤) (hroot : MeasureTheory.AEStronglyMeasurable (profilePrimalWeakRoot m hq t sample p r (profilePrimalCenter P hq t sample p r)) P) :
MeasureTheory.eLpNorm (profilePrimalWeakRoot m hq t sample p r (profilePrimalCenter P hq t sample p r)) 2 P ≤ ENNReal.ofReal (16 * diagonalWeakMetricFactor m E * diagonalWeakLoadMinus E p r) * (MeasureTheory.eLpNorm (fun (a : CoeffSpace d) => diagonalWeakCellSum q t H (1 / 2) E (sample a)) 2 P + MeasureTheory.eLpNorm (fun (a : CoeffSpace d) => diagonalWeakAverageSum q t H (1 / 2) rho E (sample a)) 2 P) + ENNReal.ofReal (16 * diagonalWeakMetricFactor m E / (2 * alpha)) * ((profileBadEnergy P (fun (a : CoeffSpace d) => diagonalWeakMaximum rho q t E (sample a)) fun (a : CoeffSpace d) => diagonalWeakEnergy hq t (sample a) p r) + profileGoodEnergy P alpha H (fun (a : CoeffSpace d) => diagonalWeakMaximum rho q t E (sample a)) fun (a : CoeffSpace d) => diagonalWeakEnergy hq t (sample a) p r) + ENNReal.ofReal constantSeminormCoefficient * profilePrimalCenterVariance P m hq t sample p r

The full primal weak root for an arbitrary coefficient sample map has the recent, optimizer-event, and exact centering-variance bound.

theorem HCPolySupport.HighContrast.Response.eLpNorm_profileAdjointWeakRoot_sample_le {d : ℕ} [NeZero d] {P : MeasureTheory.Measure (CoeffSpace d)} {q : Mat d} (hq : Matrix.PosDef q) (t : ℤ) (H : ℕ) {m : Mat d} (hm : Matrix.PosDef m) {E : BlockMat d} (hE : IsSymmetricBlockMat E) (hEpd : Book.Ch02.BlockPosDef E) {rho alpha : ℝ} (hrho : 0 < rho) (hrho1 : rho < 1) (halpha : alpha = (1 - rho) / 2) (sample : CoeffSpace d → CoeffSpace d) (p r : Vec d) (hU : AEMeasurable (fun (a : CoeffSpace d) => diagonalWeakCellSum q t H (1 / 2) E (sample a)) P) (hV : AEMeasurable (fun (a : CoeffSpace d) => diagonalWeakAverageSum q t H (1 / 2) rho E (sample a)) P) (hM : AEMeasurable (fun (a : CoeffSpace d) => diagonalWeakMaximum rho q t E (sample a)) P) (hEn : AEMeasurable (fun (a : CoeffSpace d) => diagonalWeakAdjointEnergy hq t (sample a) p r) P) (havg : AEMeasurable (fun (a : CoeffSpace d) => blockCellAverage (adaptedCell q t) (diagonalWeakAdjointState hq t (sample a) p r)) P) (hfinite : ∀ᵐ (a : CoeffSpace d) ∂P, diagonalWeakMaximum rho q t E (sample a) ≠ ⊤) (hroot : MeasureTheory.AEStronglyMeasurable (profileAdjointWeakRoot m hq t sample p r (profileAdjointCenter P hq t sample p r)) P) :
MeasureTheory.eLpNorm (profileAdjointWeakRoot m hq t sample p r (profileAdjointCenter P hq t sample p r)) 2 P ≤ ENNReal.ofReal (16 * diagonalWeakMetricFactor m E * diagonalWeakLoadPlus E p r) * (MeasureTheory.eLpNorm (fun (a : CoeffSpace d) => diagonalWeakCellSum q t H (1 / 2) E (sample a)) 2 P + MeasureTheory.eLpNorm (fun (a : CoeffSpace d) => diagonalWeakAverageSum q t H (1 / 2) rho E (sample a)) 2 P) + ENNReal.ofReal (16 * diagonalWeakMetricFactor m E / (2 * alpha)) * ((profileBadEnergy P (fun (a : CoeffSpace d) => diagonalWeakMaximum rho q t E (sample a)) fun (a : CoeffSpace d) => diagonalWeakAdjointEnergy hq t (sample a) p r) + profileGoodEnergy P alpha H (fun (a : CoeffSpace d) => diagonalWeakMaximum rho q t E (sample a)) fun (a : CoeffSpace d) => diagonalWeakAdjointEnergy hq t (sample a) p r) + ENNReal.ofReal constantSeminormCoefficient * profileAdjointCenterVariance P m hq t sample p r

The full adjoint weak root retains its independent plus load under the same arbitrary sample map.

The released split level #

theorem HCPolySupport.HighContrast.Response.eLpNorm_profilePrimalWeakRoot_sample_at_level_le {d : ℕ} [NeZero d] {P : MeasureTheory.Measure (CoeffSpace d)} {q : Mat d} (hq : Matrix.PosDef q) (t : ℤ) (H : ℕ) {m : Mat d} (hm : Matrix.PosDef m) {E : BlockMat d} (hE : IsSymmetricBlockMat E) (hEpd : Book.Ch02.BlockPosDef E) {rho alpha lev : ℝ} (hrho : 0 < rho) (hrho1 : rho < 1) (hlev : 0 < lev) (halpha : alpha = (1 - rho) / 2) (sample : CoeffSpace d → CoeffSpace d) (p r : Vec d) (hU : AEMeasurable (fun (a : CoeffSpace d) => diagonalWeakCellSum q t H (1 / 2) E (sample a)) P) (hV : AEMeasurable (fun (a : CoeffSpace d) => diagonalWeakAverageSum q t H (1 / 2) rho E (sample a)) P) (hM : AEMeasurable (fun (a : CoeffSpace d) => diagonalWeakMaximum rho q t E (sample a)) P) (hEn : AEMeasurable (fun (a : CoeffSpace d) => diagonalWeakEnergy hq t (sample a) p r) P) (havg : AEMeasurable (fun (a : CoeffSpace d) => blockCellAverage (adaptedCell q t) (diagonalWeakState hq t (sample a) p r)) P) (hfinite : ∀ᵐ (a : CoeffSpace d) ∂P, diagonalWeakMaximum rho q t E (sample a) ≠ ⊤) (hroot : MeasureTheory.AEStronglyMeasurable (profilePrimalWeakRoot m hq t sample p r (profilePrimalCenter P hq t sample p r)) P) :
MeasureTheory.eLpNorm (profilePrimalWeakRoot m hq t sample p r (profilePrimalCenter P hq t sample p r)) 2 P ≤ ENNReal.ofReal (recentConstantAtLevel lev * diagonalWeakMetricFactor m E * diagonalWeakLoadMinus E p r) * (MeasureTheory.eLpNorm (fun (a : CoeffSpace d) => diagonalWeakCellSum q t H (1 / 2) E (sample a)) 2 P + MeasureTheory.eLpNorm (fun (a : CoeffSpace d) => diagonalWeakAverageSum q t H (1 / 2) rho E (sample a)) 2 P) + ENNReal.ofReal (maxGroupConstantAtLevel lev * diagonalWeakMetricFactor m E / (2 * alpha)) * ((profileBadEnergyAt P (ENNReal.ofReal lev) (fun (a : CoeffSpace d) => diagonalWeakMaximum rho q t E (sample a)) fun (a : CoeffSpace d) => diagonalWeakEnergy hq t (sample a) p r) + profileGoodEnergyAt P (ENNReal.ofReal lev) alpha H (fun (a : CoeffSpace d) => diagonalWeakMaximum rho q t E (sample a)) fun (a : CoeffSpace d) => diagonalWeakEnergy hq t (sample a) p r) + ENNReal.ofReal constantSeminormCoefficient * profilePrimalCenterVariance P m hq t sample p r

The full primal weak root at a released split level.

theorem HCPolySupport.HighContrast.Response.eLpNorm_profileAdjointWeakRoot_sample_at_level_le {d : ℕ} [NeZero d] {P : MeasureTheory.Measure (CoeffSpace d)} {q : Mat d} (hq : Matrix.PosDef q) (t : ℤ) (H : ℕ) {m : Mat d} (hm : Matrix.PosDef m) {E : BlockMat d} (hE : IsSymmetricBlockMat E) (hEpd : Book.Ch02.BlockPosDef E) {rho alpha lev : ℝ} (hrho : 0 < rho) (hrho1 : rho < 1) (hlev : 0 < lev) (halpha : alpha = (1 - rho) / 2) (sample : CoeffSpace d → CoeffSpace d) (p r : Vec d) (hU : AEMeasurable (fun (a : CoeffSpace d) => diagonalWeakCellSum q t H (1 / 2) E (sample a)) P) (hV : AEMeasurable (fun (a : CoeffSpace d) => diagonalWeakAverageSum q t H (1 / 2) rho E (sample a)) P) (hM : AEMeasurable (fun (a : CoeffSpace d) => diagonalWeakMaximum rho q t E (sample a)) P) (hEn : AEMeasurable (fun (a : CoeffSpace d) => diagonalWeakAdjointEnergy hq t (sample a) p r) P) (havg : AEMeasurable (fun (a : CoeffSpace d) => blockCellAverage (adaptedCell q t) (diagonalWeakAdjointState hq t (sample a) p r)) P) (hfinite : ∀ᵐ (a : CoeffSpace d) ∂P, diagonalWeakMaximum rho q t E (sample a) ≠ ⊤) (hroot : MeasureTheory.AEStronglyMeasurable (profileAdjointWeakRoot m hq t sample p r (profileAdjointCenter P hq t sample p r)) P) :
MeasureTheory.eLpNorm (profileAdjointWeakRoot m hq t sample p r (profileAdjointCenter P hq t sample p r)) 2 P ≤ ENNReal.ofReal (recentConstantAtLevel lev * diagonalWeakMetricFactor m E * diagonalWeakLoadPlus E p r) * (MeasureTheory.eLpNorm (fun (a : CoeffSpace d) => diagonalWeakCellSum q t H (1 / 2) E (sample a)) 2 P + MeasureTheory.eLpNorm (fun (a : CoeffSpace d) => diagonalWeakAverageSum q t H (1 / 2) rho E (sample a)) 2 P) + ENNReal.ofReal (maxGroupConstantAtLevel lev * diagonalWeakMetricFactor m E / (2 * alpha)) * ((profileBadEnergyAt P (ENNReal.ofReal lev) (fun (a : CoeffSpace d) => diagonalWeakMaximum rho q t E (sample a)) fun (a : CoeffSpace d) => diagonalWeakAdjointEnergy hq t (sample a) p r) + profileGoodEnergyAt P (ENNReal.ofReal lev) alpha H (fun (a : CoeffSpace d) => diagonalWeakMaximum rho q t E (sample a)) fun (a : CoeffSpace d) => diagonalWeakAdjointEnergy hq t (sample a) p r) + ENNReal.ofReal constantSeminormCoefficient * profileAdjointCenterVariance P m hq t sample p r

The full adjoint weak root at a released split level.