Documentation

LeanPool.HighContrastHomogenization.Entry.Response.Cutoff.CutoffDualityBoundInputs

High-contrast homogenization: Entry.Response.Cutoff.CutoffDualityBoundInputs #

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

The negative-Besov duality bound, assembled on the reference cube #

This file serves the cutoff estimate e.response.cutoff.estimate. On the adapted cell U_t = q ⋄_t, it identifies the cutoff-weighted pairing of a harmonic optimizer's gradient and flux defects, after integration by parts and the change of variables x = q y, with minus the centred cube average of the potential against the flux defect contracted with the cutoff's gradient (AK.HC Lemma A.1, (A.4)). It reads the generic negative-Besov duality bound of the cutoff-product estimate on the reference cube, records the Besov scale weight of that cube at every scale, and shows that coordinatewise square-integrability of a doubled field on the adapted cell supplies the measurability and the two slotwise pullback bounds the duality estimate needs. Instantiating these at the centred doubled optimizer field, it assembles the resulting coefficient and seminorm bounds into the shape the final inequality consumes: a dimension-only constant times the scale weight 3^{-t} times the square of the scale-average seminorm.

The cutoff pairing as a centred cube average #

On the adapted cell U_t = q ⋄_t the cutoff-weighted pairing of the gradient defect of a harmonic optimizer against its flux defect is, after integration by parts, the negative of the normalized average of the centred potential against the flux defect contracted with the gradient of the cutoff. Composing that identity with the change of variables x = q y and with the adjointness of the pulled-back pairing rewrites it on the reference cube originCube d t, where the first slot of the pairing is the pulled-back flux defect q⁻¹ (b ∇v - q₀) and the second is the gradient of the pulled-back cutoff.

Because the flux defect is weakly divergence free, the pairing against the cutoff gradient has zero average, so centring the potential changes nothing. The resulting identity is the exact shape in which the negative-Besov duality bound e.response.cutoff.estimate (AK.HC Lemma A.1, (A.4)) is stated, and it is the identity half of the pathwise cutoff bound: no inequality is involved.

theorem HCPolySupport.HighContrast.Multiscale.volumeAverage_cutoff_pairing_eq_neg_cubeAverage_centered {d jStar : ℕ} (hjStar : 2 * d ≤ 3 ^ jStar) {F : BlockMat d} (hm : Matrix.PosDef (explicitCanonicalMetric F)) (t : ℤ) {b : CoeffField d} (u : AHarmonicFunction b (adaptedCell (respGrid jStar F) t)) (Y : BlockVec d) {φ : Vec d → ℝ} (hφ : ContDiff ℝ (↑⊤) φ) (hφ_compact : HasCompactSupport φ) (hφ_sub : tsupport φ ⊆ adaptedCell (respGrid jStar F) t) (hflux : MemVectorL2 (adaptedCell (respGrid jStar F) t) fun (x : Vec d) => matVecMul (b x) (u.toH1.grad x)) (hint1 : MeasureTheory.IntegrableOn (fun (x : Vec d) => φ x * vecDot (matVecMul (b x) (u.toH1.grad x) - Y.2) (u.toH1.grad x - Y.1)) (adaptedCell (respGrid jStar F) t) MeasureTheory.volume) (hint2 : MeasureTheory.IntegrableOn (fun (x : Vec d) => (u.toH1.toFun x - vecDot Y.1 x) * vecDot (matVecMul (b x) (u.toH1.grad x) - Y.2) fun (j : Fin d) => (fderiv ℝ φ x) (basisVec j)) (adaptedCell (respGrid jStar F) t) MeasureTheory.volume) (hcube : MeasureTheory.IntegrableOn (fun (y : Vec d) => vecDot (matVecMul (respGrid jStar F)⁻¹ (optimizerField b u (matVecMul (respGrid jStar F) y) - Y).2) (scalarCutoffGradientField (fun (z : Vec d) => φ (matVecMul (respGrid jStar F) z)) y)) (cubeSet (originCube d t)) MeasureTheory.volume) (hcubeProd : MeasureTheory.IntegrableOn (fun (y : Vec d) => (respCenteredPullbackH1 hjStar hm t b u Y).toFun y * vecDot (matVecMul (respGrid jStar F)⁻¹ (optimizerField b u (matVecMul (respGrid jStar F) y) - Y).2) (scalarCutoffGradientField (fun (z : Vec d) => φ (matVecMul (respGrid jStar F) z)) y)) (cubeSet (originCube d t)) MeasureTheory.volume) :
(volumeAverage (respCell jStar F t) fun (x : Vec d) => φ x * vecDot ((optimizerField b u x).1 - Y.1) ((optimizerField b u x).2 - Y.2)) = -cubeAverage (originCube d t) fun (y : Vec d) => vecDot (matVecMul (respGrid jStar F)⁻¹ (optimizerField b u (matVecMul (respGrid jStar F) y) - Y).2) (((respCenteredPullbackH1 hjStar hm t b u Y).toFun y - cubeAverage (originCube d t) fun (z : Vec d) => (respCenteredPullbackH1 hjStar hm t b u Y).toFun z) • scalarCutoffGradientField (fun (z : Vec d) => φ (matVecMul (respGrid jStar F) z)) y)

The cutoff-weighted pairing of the gradient defect of a harmonic optimizer against its flux defect, normalized over the adapted cell U_t = q ⋄_t, equals the negative of the normalized average over the reference cube originCube d t of the Euclidean pairing of the pulled-back flux defect q⁻¹ (b ∇v - q₀) against the centred pulled-back potential multiplied into the gradient of the pulled-back cutoff. With q the selected grid and the potential centred at its cube average, this is the identity form of the cutoff estimate e.response.cutoff.estimate (AK.HC Lemma A.1, (A.4)) on which the negative-Besov duality bound is applied.

The negative-Besov duality bound at the pulled-back cutoff data #

The generic cutoff-product duality bound abs_cubeAverage_vecDot_centered_scalar_cutoff_le_scaledWeakNormProduct in namespace Book.Ch05.Section53.JUpperBoundWeakNorms is the direct full-dual pairing form of AK.HC Lemma A.1, (A.4). This file reads it on the reference cube after the linear change of variables of the cutoff argument: with q = respGrid jStar F the centred potential is the pulled-back response respCenteredPullbackH1, the flux is the pulled-back flux defect q⁻¹ (b ∇v - q0), and the weight field is the gradient of the pulled-back cutoff scalarCutoffGradientField (φ ∘ q). Every analytic hypothesis of the generic bound — the membrane and smoothness data of the weight field, the derivative bound, the square integrability of the flux, and the two negative-Besov seminorm bounds — is carried through unchanged at the pulled-back data.

theorem HCPolySupport.HighContrast.Multiscale.abs_cubeAverage_pullback_pairing_le {d jStar : ℕ} (hjStar : 2 * d ≤ 3 ^ jStar) {F : BlockMat d} (hm : Matrix.PosDef (explicitCanonicalMetric F)) (t : ℤ) (b : CoeffField d) (u : AHarmonicFunction b (adaptedCell (respGrid jStar F) t)) (Y : BlockVec d) {φ : Vec d → ℝ} {B gradWeak fluxWeak : ℝ} (hB : 0 ≤ B) (hξLp : MeasureTheory.MemLp (scalarCutoffGradientField fun (y : Vec d) => φ (matVecMul (respGrid jStar F) y)) ⊤ (normalizedCubeMeasure (originCube d t))) (hξ : ∀ (i : Fin d), ContDiff ℝ ↑⊤ fun (x : Vec d) => scalarCutoffGradientField (fun (y : Vec d) => φ (matVecMul (respGrid jStar F) y)) x i) (hderiv : ∀ (i : Fin d), ∀ z ∈ cubeSet (originCube d t), ‖fderiv ℝ (fun (x : Vec d) => scalarCutoffGradientField (fun (y : Vec d) => φ (matVecMul (respGrid jStar F) y)) x i) z‖ ≤ B) (hflux : MeasureTheory.MemLp (fun (y : Vec d) => matVecMul (respGrid jStar F)⁻¹ (optimizerField b u (matVecMul (respGrid jStar F) y) - Y).2) 2 (normalizedCubeMeasure (originCube d t))) (hgradWeak : ∀ (N : ℕ), cubeBesovNegativeVectorPartialSeminorm (originCube d t) (1 / 2) N (respCenteredPullbackH1 hjStar hm t b u Y).grad ≤ gradWeak) (hfluxWeak : ∀ (N : ℕ), (cubeBesovNegativeVectorPartialSeminorm (originCube d t) (1 / 2) N fun (y : Vec d) => matVecMul (respGrid jStar F)⁻¹ (optimizerField b u (matVecMul (respGrid jStar F) y) - Y).2) ≤ fluxWeak) :
|cubeAverage (originCube d t) fun (x : Vec d) => vecDot (matVecMul (respGrid jStar F)⁻¹ (optimizerField b u (matVecMul (respGrid jStar F) x) - Y).2) (((respCenteredPullbackH1 hjStar hm t b u Y).toFun x - cubeAverage (originCube d t) fun (y : Vec d) => (respCenteredPullbackH1 hjStar hm t b u Y).toFun y) • scalarCutoffGradientField (fun (y : Vec d) => φ (matVecMul (respGrid jStar F) y)) x)| ≤ (2 * cubeScaleFactor (originCube d t) * B + 3 * cubeLpNorm (originCube d t) ⊤ (scalarCutoffGradientField fun (y : Vec d) => φ (matVecMul (respGrid jStar F) y))) * (Book.Ch01.Legacy.fullVectorPoincareConstant (originCube d t) * 3 ^ (↑d + 1) * ↑(Fintype.card (Fin d))) * (↑(Fintype.card (Fin d)) * (3 ^ (↑d + (1 - 1 / 2)) * cubeBesovScaleWeight (-(1 - 1 / 2 - 1 / 2)) (originCube d t))) * (cubeBesovScaleWeight (-(1 / 2)) (originCube d t) * gradWeak * (cubeBesovScaleWeight (-(1 / 2)) (originCube d t) * fluxWeak))

The negative-Besov duality bound (AK.HC Lemma A.1, (A.4)) evaluated on the reference cube openCubeSet (originCube d t) at the pulled-back cutoff data. With q = respGrid jStar F, the potential is the centred pulled-back response respCenteredPullbackH1, the flux is the pulled-back flux defect q⁻¹ (b ∇v - q0), and the weight field is the gradient of the pulled-back cutoff φ ∘ q; the two exponents take the symmetric value s = t = 1/2.

The conclusion is the generic bound with its four coefficient lets expanded: the scaled gradient and flux seminorms, the cutoff-product coefficient built from the derivative bound B and the cutoff gradient, and the flux coefficient carrying the scale weights.

Constant arithmetic of the pathwise cutoff bound #

The negative-Besov duality bound of the cutoff argument yields, on a reference cube, the product of two coefficient bounds and two pulled-back seminorm bounds. This file assembles those factors into the single shape consumed by the annealed estimate: a dimensional constant times the scale weight 3 ^ (-t) times the square of the scale-average seminorm.

theorem HCPolySupport.HighContrast.Multiscale.mul_le_const_mul_rpow_mul_sq {d : ℕ} (t : ℤ) {gradCoeff fluxCoeff scaledGrad scaledFlux G f A K L : ℝ} (hA : 0 ≤ A) (hG : 0 ≤ G) (hf : 0 ≤ f) (hK : 0 ≤ K) (hL : 0 ≤ L) (hgc : 0 ≤ gradCoeff) (hfc : 0 ≤ fluxCoeff) (hsg : 0 ≤ scaledGrad) (hsf : 0 ≤ scaledFlux) (hgradCoeff : gradCoeff ≤ K * 3 ^ (-↑t)) (hfluxCoeff : fluxCoeff ≤ L) (hscaledGrad : scaledGrad ≤ G * A) (hscaledFlux : scaledFlux ≤ f * A) (hGf : G * f ≤ 3 * ↑d ^ 2) :
gradCoeff * fluxCoeff * (scaledGrad * scaledFlux) ≤ K * L * (3 * ↑d ^ 2) * (3 ^ (-↑t) * A ^ 2)

The product of the two coefficient bounds and the two pulled-back seminorm bounds of the pathwise cutoff estimate is bounded by the dimensional constant K * L * 3 * d ^ 2 times the scale weight 3 ^ (-t) and the squared scale-average seminorm A ^ 2. The coefficient bounds are gradCoeff ≤ K * 3 ^ (-t) and fluxCoeff ≤ L, the seminorm bounds are scaledGrad ≤ G * A and scaledFlux ≤ f * A, and the two Frobenius factors satisfy G * f ≤ 3 * d ^ 2; all quantities are nonnegative. This is the form in which the pathwise cutoff bound e.response.cutoff.estimate feeds the annealed one.

The Besov scale weight of the reference cube #

The Besov scale weight cubeBesovScaleWeight s Q = (cubeScaleFactor Q) ^ (-s) assigns to a triadic cube Q the geometric weight 3^{-s · scale Q}, and the scale factor is cubeScaleFactor Q = 3 ^ (scale Q). On the reference cube originCube d t, whose scale is t, both collapse to elementary powers of 3: the scale factor is 3 ^ t and the weight is 3^{-s t}. In particular the exponent -1/2 weight is 3^{t/2} and the exponent 0 weight is 1. These are the weight evaluations used when the cutoff support of the response is expanded in the cube Besov scale.

The scale factor of the reference cube originCube d t is 3 ^ t, the integer power being read as a real power: cubeScaleFactor (originCube d t) = 3 ^ (t : ℝ).

The Besov scale weight of exponent s on the reference cube originCube d t equals 3^{-s t}: cubeBesovScaleWeight s (originCube d t) = 3 ^ (-s * (t : ℝ)).

The exponent -1/2 Besov scale weight of the reference cube originCube d t is the positive power 3^{t/2}: cubeBesovScaleWeight (-(1/2)) (originCube d t) = 3 ^ ((t : ℝ) / 2).

The exponent 0 Besov scale weight of the reference cube originCube d t is 1: cubeBesovScaleWeight 0 (originCube d t) = 1.

The pulled-back cube seminorms of a square-integrable doubled field #

The slotwise pullback bounds of the cutoff argument, partialSeminorm_pullback_fst_le and partialSeminorm_pullback_snd_le, take as side conditions the integrability of the squared Euclidean length of the metric-scaled doubled field together with its a.e. strong measurability. This file supplies both from the single hypothesis that all 2d coordinates of the field are square integrable on the adapted cell.

Paper: the cutoff estimate e.response.cutoff.estimate and its pullback to the reference cube.

theorem HCPolySupport.HighContrast.Multiscale.memLp_two_coords_metricScaled {d : ℕ} {V : Set (Vec d)} (m : Mat d) {X : Vec d → BlockVec d} (h1 : ∀ (i : Fin d), MeasureTheory.MemLp (fun (x : Vec d) => (X x).1 i) 2 (MeasureTheory.volume.restrict V)) (h2 : ∀ (i : Fin d), MeasureTheory.MemLp (fun (x : Vec d) => (X x).2 i) 2 (MeasureTheory.volume.restrict V)) :
(∀ (i : Fin d), MeasureTheory.MemLp (fun (x : Vec d) => matVecMul (matSqrt m) (X x).1 i) 2 (MeasureTheory.volume.restrict V)) ∧ ∀ (i : Fin d), MeasureTheory.MemLp (fun (x : Vec d) => matVecMul (matSqrt m)⁻¹ (X x).2 i) 2 (MeasureTheory.volume.restrict V)

The metric scaling of a doubled field, by m^{1/2} in the first slot and m^{-1/2} in the second, carries square integrable coordinates to square integrable coordinates.

theorem HCPolySupport.HighContrast.Multiscale.partialSeminorm_pullback_fst_le_of_coords {d : ℕ} [NeZero d] (q : Mat d) (hq : IsUnit q) (m : Mat d) (hm : Matrix.PosDef m) (t : ℤ) (N : ℕ) (X : Vec d → BlockVec d) (h1 : ∀ (i : Fin d), MeasureTheory.MemLp (fun (x : Vec d) => (X x).1 i) 2 (MeasureTheory.volume.restrict (adaptedCell q t))) (h2 : ∀ (i : Fin d), MeasureTheory.MemLp (fun (x : Vec d) => (X x).2 i) 2 (MeasureTheory.volume.restrict (adaptedCell q t))) :

The slotwise pullback bound for the gradient slot, with the square-integrability side conditions supplied coordinatewise.

theorem HCPolySupport.HighContrast.Multiscale.partialSeminorm_pullback_snd_le_of_coords {d : ℕ} [NeZero d] (q : Mat d) (hq : IsUnit q) (m : Mat d) (hm : Matrix.PosDef m) (t : ℤ) (N : ℕ) (X : Vec d → BlockVec d) (h1 : ∀ (i : Fin d), MeasureTheory.MemLp (fun (x : Vec d) => (X x).1 i) 2 (MeasureTheory.volume.restrict (adaptedCell q t))) (h2 : ∀ (i : Fin d), MeasureTheory.MemLp (fun (x : Vec d) => (X x).2 i) 2 (MeasureTheory.volume.restrict (adaptedCell q t))) :

The slotwise pullback bound for the flux slot, with the square-integrability side conditions supplied coordinatewise.

The two negative-Besov inputs of the duality bound #

The slotwise pullback bounds partialSeminorm_pullback_fst_le_of_coords and partialSeminorm_pullback_snd_le_of_coords hold for an arbitrary doubled field X, subject to coordinatewise square integrability on the adapted cell. This file instantiates them at the centred doubled optimizer field X = (∇v, b∇v) − Y. The first slot is then the gradient of the centred pulled-back potential, and the second the pulled-back flux defect, while the square-integrability hypotheses are supplied by the elliptic representative of the coefficient.

Paper: the cutoff estimate e.response.cutoff.estimate and the scale-average seminorm bound of e.response.weak.estimate.

theorem HCPolySupport.HighContrast.Multiscale.cubeBesov_respCenteredPullbackH1_grad_le {d : ℕ} [NeZero d] {jStar : ℕ} (hjStar : 2 * d ≤ 3 ^ jStar) {F : BlockMat d} (hm : Matrix.PosDef (explicitCanonicalMetric F)) (t : ℤ) {lam Lam : ℝ} {b f : CoeffField d} (hEll : IsEllipticFieldOn lam Lam (respCell jStar F t) f) (hbf : b =ᵐ[volumeMeasureOn (respCell jStar F t)] f) (u : AHarmonicFunction b (respCell jStar F t)) (Y : BlockVec d) (N : ℕ) :

The negative-Besov seminorm on the reference cube of the gradient of the centred pulled-back potential is bounded by the scale-average seminorm of the metric-scaled doubled optimizer state, with the metric factor qᵀ m^{-1/2}. This is the gradient slot of the duality bound, obtained by applying partialSeminorm_pullback_fst_le_of_coords to the centred doubled optimizer field and rewriting its first slot through the gradient of the pulled-back centred potential.

theorem HCPolySupport.HighContrast.Multiscale.cubeBesov_pullback_fluxDefect_le {d : ℕ} [NeZero d] {jStar : ℕ} (hjStar : 2 * d ≤ 3 ^ jStar) {F : BlockMat d} (hm : Matrix.PosDef (explicitCanonicalMetric F)) (t : ℤ) {lam Lam : ℝ} {b f : CoeffField d} (hEll : IsEllipticFieldOn lam Lam (respCell jStar F t) f) (hbf : b =ᵐ[volumeMeasureOn (respCell jStar F t)] f) (u : AHarmonicFunction b (respCell jStar F t)) (Y : BlockVec d) (N : ℕ) :

The negative-Besov seminorm on the reference cube of the pulled-back flux defect is bounded by the scale-average seminorm of the metric-scaled doubled optimizer state, with the metric factor q⁻¹ m^{1/2}. This is the flux slot of the duality bound, obtained by applying partialSeminorm_pullback_snd_le_of_coords to the centred doubled optimizer field.