Documentation

LeanPool.HighContrastHomogenization.Entry.Response.Direct.OptimizerEnergyDefectRow

High-contrast homogenization: Entry.Response.Direct.OptimizerEnergyDefectRow #

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

The terminal-optimizer replacement row, assembled on the response's own data #

The (φ-1)-weighted optimizer energy on the terminal cell splits, over the subdivision into the aligned cells of the coarse scale, into an oscillation part — costing the cutoff class's Lipschitz gain times the terminal response — and a subcell-mean part, the mean-zero weighting of the subcell half-energies; combining the two gives the first error row of p.response.transfer. This file instantiates that row on the same carriers the cutoff estimate uses: the difference between the cutoff half-energy of the terminal optimizer and the terminal response, integrated over the law of coefficients. The pathwise identity behind that difference, cutoffHalfEnergyAux - respJ equal to half the (φ-1)-weighted optimizer energy, holds only for a coefficient that is pointwise elliptic on the cell, so it is transported here across an almost-everywhere equality to the two response coefficients respCoeffMinus and respCoeffPlus, each elliptic only almost everywhere.

The terminal-optimizer replacement row of p.response.transfer #

The (φ - 1)-weighted optimizer energy on the terminal cell splits, over the subdivision into the aligned cells of the coarse scale, into its oscillation part and its subcell-mean part. The oscillation part costs the Lipschitz gain of the cutoff class times the terminal response. The subcell-mean part is the mean-zero weighting of the subcell half-energies, and the abstract row bounds its expectation by τ + 2 √(τ cJ), where τ is the scale defect and cJ the common annealed subcell response. Since cJ = EJ + τ and √(τ (EJ + τ)) ≤ √(τ EJ) + τ, the two halves combine into the printed shape 6 τ + 4 √(τ EJ) + 2 K EJ.

Paper: p.response.transfer.

theorem HCPolySupport.HighContrast.Multiscale.abs_integral_le_row1_of_parts {α : Type u_1} [MeasurableSpace α] (P : MeasureTheory.Measure α) {d : ℕ} (n : ℕ) (c : (Fin d → ℤ) → ℝ) (Ecell J D : (Fin d → ℤ) → α → ℝ) (Wtot Jt : α → ℝ) (K cJ τ EJ : ℝ) (hK : 0 ≤ K) (hosc : ∀ (a : α), |Wtot a - (↑(triadicIndexBox d n).card)⁻¹ * ∑ w ∈ triadicIndexBox d n, c w * Ecell w a| ≤ K * (2 * Jt a)) (hcmp : ∀ w ∈ triadicIndexBox d n, ∀ (a : α), |1 / 2 * Ecell w a - J w a| ≤ D w a + 2 * √(J w a * D w a)) (hc0 : (↑(triadicIndexBox d n).card)⁻¹ * ∑ w ∈ triadicIndexBox d n, c w = 0) (hc1 : ∀ w ∈ triadicIndexBox d n, |c w| ≤ 1) (hJ0 : ∀ w ∈ triadicIndexBox d n, ∀ (a : α), 0 ≤ J w a) (hD0 : ∀ w ∈ triadicIndexBox d n, ∀ (a : α), 0 ≤ D w a) (hJt0 : ∀ (a : α), 0 ≤ Jt a) (hJint : ∀ w ∈ triadicIndexBox d n, MeasureTheory.Integrable (J w) P) (hDint : ∀ w ∈ triadicIndexBox d n, MeasureTheory.Integrable (D w) P) (hEint : ∀ w ∈ triadicIndexBox d n, MeasureTheory.Integrable (Ecell w) P) (hWint : MeasureTheory.Integrable Wtot P) (hJtint : MeasureTheory.Integrable Jt P) (hJval : ∀ w ∈ triadicIndexBox d n, ∫ (a : α), J w a ∂P = cJ) (hDval : ∫ (a : α), (↑(triadicIndexBox d n).card)⁻¹ * ∑ w ∈ triadicIndexBox d n, D w a ∂P = τ) (hEJ : ∫ (a : α), Jt a ∂P = EJ) (hcJ : cJ = EJ + τ) :
|∫ (a : α), Wtot a ∂P| ≤ 6 * τ + 4 * √(τ * EJ) + 2 * K * EJ

The terminal-optimizer replacement row. Suppose the terminal response Wtot differs from the normalized mean-zero weighting N ∑ c w Ecell w of the subcell half-energies, where N = |triadicIndexBox d n|⁻¹, by at most K times twice the terminal response Jt. Suppose further that on each subcell the half-energy (1/2) Ecell and the subcell response J satisfy the one-parameter comparison |(1/2) Ecell - J| ≤ D + 2 √(J D) with nonnegative J and deficit D, that the weights c are bounded by 1 and average to 0, that ∫ J w = cJ is the same on every subcell, and that the annealed flat average of the deficits D is the scale defect τ. When cJ = EJ + τ with EJ = ∫ Jt, the expectation of the terminal response is at most 6 τ + 4 √(τ EJ) + 2 K EJ. This is the first error row of p.response.transfer assembled from its oscillation and subcell-mean halves.

The terminal-optimizer replacement row on the carriers of the cutoff estimate #

This assembles the first error row of p.response.transfer on the same carriers the cutoff estimate uses. The left-hand side is the difference between the cutoff half-energy of the terminal optimizer and the terminal response, integrated over the law of coefficients. The three geometric inputs — the pathwise identification of that difference with the (φ - 1)-weighted optimizer energy, the oscillation split of that energy over the subdivision into the aligned coarse cells, and the subcell comparison of the half-energy with the subcell response — enter as hypotheses in the shape in which the modules that prove them produce them. The conclusion is the printed bound

C (τ^- + (τ^- E[J_t^-])^{1/2} + 3^{-H} E[J_t^-])

with the explicit dimensional constant C = max 3 (32 d^2 * responseCutoffProfileConst).

Paper: p.response.transfer.

theorem HCPolySupport.HighContrast.Multiscale.abs_integral_cutoffHalfEnergy_sub_respJ_le_rowMinus {d : ℕ} (P : MeasureTheory.Measure (CoeffSpace d)) (jStar : ℕ) (F : BlockMat d) (H : ℕ) (s t : ℤ) (e : Vec d) (φ : Vec d → ℝ) (uM : (a : CoeffSpace d) → AHarmonicFunction (respCoeffMinus F a) (respCell jStar F t)) (hid : ∀ (a : CoeffSpace d), cutoffHalfEnergyAux (respCell jStar F t) φ (respCoeffMinus F a) (uM a) - respJ (respGrid jStar F) t (respP (respMean P jStar F t) e) (respqMinus P jStar F t e) (respCoeffMinus F a) = 1 / 2 * volumeAverage (respCell jStar F t) fun (x : Vec d) => (φ x - 1) * scalarVariationEnergyIntegrand (respCoeffMinus F a) (uM a) x) (hosc : ∀ (a : CoeffSpace d), |(volumeAverage (respCell jStar F t) fun (x : Vec d) => (φ x - 1) * scalarVariationEnergyIntegrand (respCoeffMinus F a) (uM a) x) - (↑(triadicIndexBox d H).card)⁻¹ * ∑ w ∈ triadicIndexBox d H, (volumeAverage (adaptedCellAtCenter (respGrid jStar F) s w) fun (x : Vec d) => φ x - 1) * volumeAverage (adaptedCellAtCenter (respGrid jStar F) s w) (scalarVariationEnergyIntegrand (respCoeffMinus F a) (uM a))| ≤ 32 * ↑d ^ 2 * responseCutoffProfileConst * 3 ^ (-↑H) * (2 * respJ (respGrid jStar F) t (respP (respMean P jStar F t) e) (respqMinus P jStar F t e) (respCoeffMinus F a))) (hcmp : ∀ w ∈ triadicIndexBox d H, ∀ (a : CoeffSpace d), |1 / 2 * volumeAverage (adaptedCellAtCenter (respGrid jStar F) s w) (scalarVariationEnergyIntegrand (respCoeffMinus F a) (uM a)) - ResponseJ (adaptedCellAtCenter (respGrid jStar F) s w) (respP (respMean P jStar F t) e) (respqMinus P jStar F t e) (respCoeffMinus F a)| ≤ ResponseJ (adaptedCellAtCenter (respGrid jStar F) s w) (respP (respMean P jStar F t) e) (respqMinus P jStar F t e) (respCoeffMinus F a) - volumeAverage (adaptedCellAtCenter (respGrid jStar F) s w) (scalarResponseIntegrand (respCell jStar F t) (respCoeffMinus F a) (respP (respMean P jStar F t) e) (respqMinus P jStar F t e) (uM a)) + 2 * √(ResponseJ (adaptedCellAtCenter (respGrid jStar F) s w) (respP (respMean P jStar F t) e) (respqMinus P jStar F t e) (respCoeffMinus F a) * (ResponseJ (adaptedCellAtCenter (respGrid jStar F) s w) (respP (respMean P jStar F t) e) (respqMinus P jStar F t e) (respCoeffMinus F a) - volumeAverage (adaptedCellAtCenter (respGrid jStar F) s w) (scalarResponseIntegrand (respCell jStar F t) (respCoeffMinus F a) (respP (respMean P jStar F t) e) (respqMinus P jStar F t e) (uM a))))) (hc0 : ((↑(triadicIndexBox d H).card)⁻¹ * ∑ w ∈ triadicIndexBox d H, volumeAverage (adaptedCellAtCenter (respGrid jStar F) s w) fun (x : Vec d) => φ x - 1) = 0) (hc1 : ∀ w ∈ triadicIndexBox d H, |volumeAverage (adaptedCellAtCenter (respGrid jStar F) s w) fun (x : Vec d) => φ x - 1| ≤ 1) (hJ0 : ∀ w ∈ triadicIndexBox d H, ∀ (a : CoeffSpace d), 0 ≤ ResponseJ (adaptedCellAtCenter (respGrid jStar F) s w) (respP (respMean P jStar F t) e) (respqMinus P jStar F t e) (respCoeffMinus F a)) (hD0 : ∀ w ∈ triadicIndexBox d H, ∀ (a : CoeffSpace d), 0 ≤ ResponseJ (adaptedCellAtCenter (respGrid jStar F) s w) (respP (respMean P jStar F t) e) (respqMinus P jStar F t e) (respCoeffMinus F a) - volumeAverage (adaptedCellAtCenter (respGrid jStar F) s w) (scalarResponseIntegrand (respCell jStar F t) (respCoeffMinus F a) (respP (respMean P jStar F t) e) (respqMinus P jStar F t e) (uM a))) (hJt0 : ∀ (a : CoeffSpace d), 0 ≤ respJ (respGrid jStar F) t (respP (respMean P jStar F t) e) (respqMinus P jStar F t e) (respCoeffMinus F a)) (hJint : ∀ w ∈ triadicIndexBox d H, MeasureTheory.Integrable (fun (a : CoeffSpace d) => ResponseJ (adaptedCellAtCenter (respGrid jStar F) s w) (respP (respMean P jStar F t) e) (respqMinus P jStar F t e) (respCoeffMinus F a)) P) (hDint : ∀ w ∈ triadicIndexBox d H, MeasureTheory.Integrable (fun (a : CoeffSpace d) => ResponseJ (adaptedCellAtCenter (respGrid jStar F) s w) (respP (respMean P jStar F t) e) (respqMinus P jStar F t e) (respCoeffMinus F a) - volumeAverage (adaptedCellAtCenter (respGrid jStar F) s w) (scalarResponseIntegrand (respCell jStar F t) (respCoeffMinus F a) (respP (respMean P jStar F t) e) (respqMinus P jStar F t e) (uM a))) P) (hEint : ∀ w ∈ triadicIndexBox d H, MeasureTheory.Integrable (fun (a : CoeffSpace d) => volumeAverage (adaptedCellAtCenter (respGrid jStar F) s w) (scalarVariationEnergyIntegrand (respCoeffMinus F a) (uM a))) P) (hWint : MeasureTheory.Integrable (fun (a : CoeffSpace d) => volumeAverage (respCell jStar F t) fun (x : Vec d) => (φ x - 1) * scalarVariationEnergyIntegrand (respCoeffMinus F a) (uM a) x) P) (hJtint : MeasureTheory.Integrable (fun (a : CoeffSpace d) => respJ (respGrid jStar F) t (respP (respMean P jStar F t) e) (respqMinus P jStar F t e) (respCoeffMinus F a)) P) (hJval : ∀ w ∈ triadicIndexBox d H, ∫ (a : CoeffSpace d), ResponseJ (adaptedCellAtCenter (respGrid jStar F) s w) (respP (respMean P jStar F t) e) (respqMinus P jStar F t e) (respCoeffMinus F a) ∂P = ∫ (a : CoeffSpace d), respJ (respGrid jStar F) s (respP (respMean P jStar F t) e) (respqMinus P jStar F t e) (respCoeffMinus F a) ∂P) (hDval : ∫ (a : CoeffSpace d), (↑(triadicIndexBox d H).card)⁻¹ * ∑ w ∈ triadicIndexBox d H, (ResponseJ (adaptedCellAtCenter (respGrid jStar F) s w) (respP (respMean P jStar F t) e) (respqMinus P jStar F t e) (respCoeffMinus F a) - volumeAverage (adaptedCellAtCenter (respGrid jStar F) s w) (scalarResponseIntegrand (respCell jStar F t) (respCoeffMinus F a) (respP (respMean P jStar F t) e) (respqMinus P jStar F t e) (uM a))) ∂P = respTauMinus P jStar F s t e) :
|∫ (a : CoeffSpace d), cutoffHalfEnergyAux (respCell jStar F t) φ (respCoeffMinus F a) (uM a) - respJ (respGrid jStar F) t (respP (respMean P jStar F t) e) (respqMinus P jStar F t e) (respCoeffMinus F a) ∂P| ≤ max 3 (32 * ↑d ^ 2 * responseCutoffProfileConst) * (respTauMinus P jStar F s t e + √(respTauMinus P jStar F s t e * respEJMinus P jStar F t e) + 3 ^ (-↑H) * respEJMinus P jStar F t e)

The terminal-optimizer replacement row on the cutoff carriers. Assume the pathwise identification hid of the cutoff half-energy minus the terminal response with the (φ - 1)-weighted optimizer energy; the oscillation split hosc of that weighted energy over the aligned cells of the coarse scale, costing 32 d^2 responseCutoffProfileConst 3^{-H} times twice the terminal response; and the subcell comparison hcmp of half the subcell energy with the subcell response through its nonnegative deficit. Assume the cutoff fluctuation has mean-zero normalized subcell means hc0 and is bounded by 1 on each subcell hc1, the subcell responses and deficits are nonnegative and integrable, the subcell responses have a common annealed value hJval, and the annealed flat mean of the deficits is the scale defect hDval. Then the law-integral of the cutoff half-energy minus the terminal response is at most max 3 (32 d^2 responseCutoffProfileConst) times τ^- + (τ^- E[J_t^-])^{1/2} + 3^{-H} E[J_t^-]. This is the first error row of p.response.transfer on the carriers of the cutoff estimate.

The cutoff energy defect at an almost-everywhere elliptic coefficient #

The response identity cutoffHalfEnergy_sub_respJ_eq identifies the cutoff energy defect cutoffHalfEnergyAux - respJ with half the (φ - 1)-weighted optimizer energy, but only for a coefficient that is pointwise elliptic on the cell. The two coefficient fields of the response problem, respCoeffMinus F a and respCoeffPlus F a, are elliptic only almost everywhere: a point of CoeffSpace d is an a.e. class, and no representative is elliptic pointwise.

The identity is nevertheless invariant under an a.e. replacement of the coefficient, because carrying a harmonic function along the replacement leaves its gradient — and hence both the cutoff half-energy and the pathwise response — unchanged. Running the pointwise identity at the elliptic representative supplied by exists_elliptic_representative_respCell_respCoeffMinus and transporting back gives the defect identity verbatim for the carrier coefficient and the given maximizer:

theorem HCPolySupport.HighContrast.Multiscale.cutoffHalfEnergyAux_sub_respJ_eq_respCoeffMinus {d jStar : ℕ} (hjStar : 2 * d ≤ 3 ^ jStar) {F : BlockMat d} (hm : Matrix.PosDef (explicitCanonicalMetric F)) (t : ℤ) (φ : Vec d → ℝ) (hφ : IsResponseCutoff (respGrid jStar F) t φ) (p r : Vec d) (a : CoeffSpace d) (u : AHarmonicFunction (respCoeffMinus F a) (respCell jStar F t)) (hmax : IsResponseMaximizer (respCell jStar F t) p r (respCoeffMinus F a) u) :
cutoffHalfEnergyAux (respCell jStar F t) φ (respCoeffMinus F a) u - respJ (respGrid jStar F) t p r (respCoeffMinus F a) = 1 / 2 * volumeAverage (respCell jStar F t) fun (x : Vec d) => (φ x - 1) * scalarVariationEnergyIntegrand (respCoeffMinus F a) u x

The cutoff energy defect of the minus family, with no pointwise ellipticity hypothesis.

theorem HCPolySupport.HighContrast.Multiscale.cutoffHalfEnergyAux_sub_respJ_eq_respCoeffPlus {d jStar : ℕ} (hjStar : 2 * d ≤ 3 ^ jStar) {F : BlockMat d} (hm : Matrix.PosDef (explicitCanonicalMetric F)) (t : ℤ) (φ : Vec d → ℝ) (hφ : IsResponseCutoff (respGrid jStar F) t φ) (p r : Vec d) (a : CoeffSpace d) (u : AHarmonicFunction (respCoeffPlus F a) (respCell jStar F t)) (hmax : IsResponseMaximizer (respCell jStar F t) p r (respCoeffPlus F a) u) :
cutoffHalfEnergyAux (respCell jStar F t) φ (respCoeffPlus F a) u - respJ (respGrid jStar F) t p r (respCoeffPlus F a) = 1 / 2 * volumeAverage (respCell jStar F t) fun (x : Vec d) => (φ x - 1) * scalarVariationEnergyIntegrand (respCoeffPlus F a) u x

The adjoint twin.