Documentation

LeanPool.Superdiffusion.Algsuperdiff.Section3.Provider.Diffusivity.ApproximateRecurrence.Closure.JunctionDischargeStepFiveDisplay

JunctionDischargeStepFiveDisplay #

Ported from scottnarmstrong/Superdiffusion at 7f20476385d4af418f110303085ad0740c9eba46. Module imports, support namespace and local option directives adapted for Lean Pool.

The Step-5 gauge junction: the averaged display and the endpoint #

ABK26, l.approximate.recurrence.formula, Step 5, together with e.what.nablaw.really.is and e.lower.bound.oscillations.

Closure.JunctionDischargeStepFive reduces Closure.Step5GaugeEndpoint to one binder, hdisplay:

  avsum_R E[ step5CellExcess ]
    <= - E|| grad w_D^{(K)} ||^2_{L2(cu_K)}
       + Fend h^2 shom_n^{-4} 3^{4 cgamma (n+h)} + cgamma^15 .

This module proves that display, and then the endpoint.

The three pieces of the display #

Carrier transports the caller must supply on hosc #

The proved Lean endpoint exists_gamma0_freshShellDirichlet_meshOscillation_le_gamma_pow_sixteen (LocalizationOscillationEndpointSixteen.lean) proves e.lower.bound.oscillations

Neither transport is performed here, and neither is assumed away: both live inside the single binder hosc, which is the manuscript's own display.

References #

e.what.nablaw.really.is #

The Cauchy--Schwarz step that feeds the fourth-moment endpoint #

theorem Algsuperdiff.Section3.Provider.Diffusivity.ApproximateRecurrence.Closure.integral_rpow_two_le_rpow_integral_rpow_four {Omega : Type u_1} [MeasurableSpace Omega] (mu : MeasureTheory.Measure Omega) [MeasureTheory.IsProbabilityMeasure mu] {X : Omega → ℝ} (hX : ∀ (omega : Omega), 0 ≤ X omega) (hXm : MeasureTheory.MemLp X 4 mu) :
∫ (omega : Omega), X omega ^ 2 ∂mu ≤ (∫ (omega : Omega), X omega ^ 4 ∂mu) ^ 2⁻¹

Lyapunov's inequality on a probability space, in the rpow spelling.

theorem Algsuperdiff.Section3.Provider.Diffusivity.ApproximateRecurrence.Closure.integral_mul_le_rpowFourRoot_mul_rpowFourRoot {Omega : Type u_1} [MeasurableSpace Omega] (mu : MeasureTheory.Measure Omega) [MeasureTheory.IsProbabilityMeasure mu] {X Y : Omega → ℝ} (hX : ∀ (omega : Omega), 0 ≤ X omega) (hY : ∀ (omega : Omega), 0 ≤ Y omega) (hXm : MeasureTheory.MemLp X 4 mu) (hYm : MeasureTheory.MemLp Y 4 mu) :
∫ (omega : Omega), X omega * Y omega ∂mu ≤ (∫ (omega : Omega), X omega ^ 4 ∂mu) ^ 4⁻¹ * (∫ (omega : Omega), Y omega ^ 4 ∂mu) ^ 4⁻¹

Cauchy--Schwarz in the sample, in the fourth-moment spelling.

The grid Cauchy--Schwarz for a per-cube second factor #

The one-sided step above is the shape the display needed while the forcing gauge of e.nablaw.oscillations was read at a single window z0. The gauge is sited at the cube's own base point, so the second factor of the cross term is itself a grid family, and the grid step must be Hoelder in both factors. The three lemmas below supply that.

Jensen for the concave square root, at the finite cube family.

Hoelder on the grid at a pair of fourth roots.

theorem Algsuperdiff.Section3.Provider.Diffusivity.ApproximateRecurrence.Closure.cubeFamilyAverage_integral_mul_le_gridFourthMomentRoot_mul_root {d : ℕ} {Omega : Type u_1} [MeasurableSpace Omega] (mu : MeasureTheory.Measure Omega) [MeasureTheory.IsProbabilityMeasure mu] (I : Finset (HCPolySupport.TriadicCube d)) (X Y : HCPolySupport.TriadicCube d → Omega → ℝ) (hX : ∀ (R : HCPolySupport.TriadicCube d) (omega : Omega), 0 ≤ X R omega) (hY : ∀ (R : HCPolySupport.TriadicCube d) (omega : Omega), 0 ≤ Y R omega) (hXm : ∀ R ∈ I, MeasureTheory.MemLp (X R) 4 mu) (hYm : ∀ R ∈ I, MeasureTheory.MemLp (Y R) 4 mu) :
(cubeFamilyAverage I fun (R : HCPolySupport.TriadicCube d) => ∫ (omega : Omega), X R omega * Y R omega ∂mu) ≤ gridFourthMomentRoot mu I X * gridFourthMomentRoot mu I Y

The two-sided grid/sample Cauchy--Schwarz that feeds the re-sited e.lower.bound.oscillations.

The per-cube cross defect and the grid identity #

The Step-5 per-cube cross defect: the cube average of the product h grad w, tested against e, minus the product of the two cube averages. It is the entire cost of replacing (h)_R (grad w)_R by (h grad w)_R.

Equations
  • One or more equations did not get rendered due to their size.
Instances For

    The grid form of the cross term: the descendant average of the products of the averages is the cube average of the product, minus the grid average of the defects. This is the exact tiling bookkeeping.

    The fourth-root gauge normalization #

    The two arithmetic folds of the display #

    The four legs of the averaged display #

    The averaged Step-5 display #

    theorem Algsuperdiff.Section3.Provider.Diffusivity.ApproximateRecurrence.Closure.step5CellDisplay_of_legs {d : ℕ} (M : ABKModel d) (n : ℤ) (h K j : ℕ) (e : HCPolySupport.Vec d) (nmesh : ℤ) (gauge : HCPolySupport.TriadicCube d → Cutoff.CutoffSample d → ℝ) {C1 C2 Cdef shellMoment gradMoment : ℝ} (wD : C(HCPolySupport.Vec d, HCPolySupport.Mat d) → HCPolySupport.H10Function (HCPolySupport.openCubeSet (HCPolySupport.originCube d ↑K))) (hgauge0 : ∀ (R : HCPolySupport.TriadicCube d) (omega : Cutoff.CutoffSample d), 0 ≤ gauge R omega) (hCdef : 0 ≤ Cdef) (hgate : Cdef * M.gamma ^ 15 ≤ 2) (hD : ∀ (omega : Cutoff.ShellSeq d), HCPolySupport.IsZeroTraceDirichletRhsWeakSolution (fun (x : HCPolySupport.Vec d) => 1) (HCPolySupport.openCubeSet (HCPolySupport.originCube d ↑K)) (alongIncrementPath n h wD omega) fun (x : HCPolySupport.Vec d) => -Corrector.streamForcing (↑(Annealed.sigmaBar M n))⁻¹ omega n (n + ↑h) e x) (hmemD : ∀ (omega : Cutoff.ShellSeq d), ∀ R ∈ HCPolySupport.descendantsAtDepth (HCPolySupport.originCube d ↑K) j, MeasureTheory.MemLp (fun (x : HCPolySupport.Vec d) => (alongIncrementPath n h wD omega).grad x) 2 (HCPolySupport.normalizedCubeMeasure R)) (hintD : ∀ (omega : Cutoff.ShellSeq d), MeasureTheory.IntegrableOn (fun (x : HCPolySupport.Vec d) => HCPolySupport.vecNormSq ((alongIncrementPath n h wD omega).grad x)) (HCPolySupport.cubeSet (HCPolySupport.originCube d ↑K)) MeasureTheory.volume) (hintP : ∀ (omega : Cutoff.ShellSeq d), MeasureTheory.IntegrableOn (fun (x : HCPolySupport.Vec d) => HCPolySupport.vecDot e (HCPolySupport.matVecMul ((Cutoff.finiteShellIncrement omega n (n + ↑h)).toFun x) ((alongIncrementPath n h wD omega).grad x))) (HCPolySupport.cubeSet (HCPolySupport.originCube d ↑K)) MeasureTheory.volume) (hsampN : ∀ R ∈ HCPolySupport.descendantsAtDepth (HCPolySupport.originCube d ↑K) j, MeasureTheory.Integrable (fun (omega : Cutoff.CutoffSample d) => HCPolySupport.vecNormSq (HCPolySupport.cubeAverageVec R fun (x : HCPolySupport.Vec d) => (alongIncrementPath n h wD ↑omega).grad x)) ↑(Cutoff.cutoffSampleLaw M)) (hsampX : ∀ R ∈ HCPolySupport.descendantsAtDepth (HCPolySupport.originCube d ↑K) j, MeasureTheory.Integrable (fun (omega : Cutoff.CutoffSample d) => HCPolySupport.vecDot e (HCPolySupport.matVecMul (freshShellCubeAverage R (↑omega) n (n + ↑h)) (HCPolySupport.cubeAverageVec R fun (x : HCPolySupport.Vec d) => (alongIncrementPath n h wD ↑omega).grad x))) ↑(Cutoff.cutoffSampleLaw M)) (hsampS : ∀ R ∈ HCPolySupport.descendantsAtDepth (HCPolySupport.originCube d ↑K) j, MeasureTheory.Integrable (fun (omega : Cutoff.CutoffSample d) => HCPolySupport.vecNormSq (HCPolySupport.matVecMul (freshShellCubeAverage R (↑omega) n (n + ↑h)) (HCPolySupport.cubeAverageVec R fun (x : HCPolySupport.Vec d) => (alongIncrementPath n h wD ↑omega).grad x))) ↑(Cutoff.cutoffSampleLaw M)) (hsampDef : ∀ R ∈ HCPolySupport.descendantsAtDepth (HCPolySupport.originCube d ↑K) j, MeasureTheory.Integrable (fun (omega : Cutoff.CutoffSample d) => step5CrossDefect R (↑omega) n (n + ↑h) e fun (x : HCPolySupport.Vec d) => (alongIncrementPath n h wD ↑omega).grad x) ↑(Cutoff.cutoffSampleLaw M)) (hsampG : MeasureTheory.Integrable (fun (omega : Cutoff.CutoffSample d) => HCPolySupport.cubeAverage (HCPolySupport.originCube d ↑K) fun (x : HCPolySupport.Vec d) => HCPolySupport.vecNormSq ((alongIncrementPath n h wD ↑omega).grad x)) ↑(Cutoff.cutoffSampleLaw M)) (hdefect : ∀ R ∈ HCPolySupport.descendantsAtDepth (HCPolySupport.originCube d ↑K) j, ∀ (omega : Cutoff.CutoffSample d), (step5CrossDefect R (↑omega) n (n + ↑h) e fun (x : HCPolySupport.Vec d) => (alongIncrementPath n h wD ↑omega).grad x) ≤ Cdef * (meshOscillationCell nmesh (fun (x : HCPolySupport.Vec d) => (alongIncrementPath n h wD ↑omega).grad x) R * gauge R omega)) (hmemOsc : ∀ R ∈ HCPolySupport.descendantsAtDepth (HCPolySupport.originCube d ↑K) j, MeasureTheory.MemLp (fun (omega : Cutoff.CutoffSample d) => meshOscillationCell nmesh (fun (x : HCPolySupport.Vec d) => (alongIncrementPath n h wD ↑omega).grad x) R) 4 ↑(Cutoff.cutoffSampleLaw M)) (hmemGauge : ∀ R ∈ HCPolySupport.descendantsAtDepth (HCPolySupport.originCube d ↑K) j, MeasureTheory.MemLp (fun (omega : Cutoff.CutoffSample d) => gauge R omega) 4 ↑(Cutoff.cutoffSampleLaw M)) (hosc : (gridFourthMomentRoot (↑(Cutoff.cutoffSampleLaw M)) (HCPolySupport.descendantsAtDepth (HCPolySupport.originCube d ↑K) j) fun (R : HCPolySupport.TriadicCube d) (omega : Cutoff.CutoffSample d) => meshOscillationCell nmesh (fun (x : HCPolySupport.Vec d) => (alongIncrementPath n h wD ↑omega).grad x) R) + (↑(Annealed.sigmaBar M n))⁻¹ * gridFourthMomentRoot (↑(Cutoff.cutoffSampleLaw M)) (HCPolySupport.descendantsAtDepth (HCPolySupport.originCube d ↑K) j) gauge ≤ M.gamma ^ 15) (hC1 : 0 ≤ C1) (hgrad0 : 0 ≤ gradMoment) (hprod : (HCPolySupport.descendantsAverage (HCPolySupport.originCube d ↑K) j fun (R : HCPolySupport.TriadicCube d) => ∫ (omega : Cutoff.CutoffSample d), HCPolySupport.vecNormSq (HCPolySupport.matVecMul (freshShellCubeAverage R (↑omega) n (n + ↑h)) (HCPolySupport.cubeAverageVec R fun (x : HCPolySupport.Vec d) => (alongIncrementPath n h wD ↑omega).grad x)) ∂↑(Cutoff.cutoffSampleLaw M)) ≤ shellMoment * gradMoment) (hshell : shellMoment ≤ C1 * ↑h * 3 ^ (2 * M.gamma * ↑(n + ↑h))) (hgradM : gradMoment ≤ C2 * (↑(Annealed.sigmaBar M n) ^ 2)⁻¹ * ↑h * 3 ^ (2 * M.gamma * ↑(n + ↑h))) :
    (HCPolySupport.descendantsAverage (HCPolySupport.originCube d ↑K) j fun (R : HCPolySupport.TriadicCube d) => ∫ (omega : Cutoff.CutoffSample d), step5CellExcess (Annealed.sigmaBar M n) R n (n + ↑h) e (alongIncrementPath n h wD) ↑omega ∂↑(Cutoff.cutoffSampleLaw M)) ≤ -dirichletCubeEnergy M n h K wD + C1 * C2 * ↑h ^ 2 * (↑(Annealed.sigmaBar M n) ^ 4)⁻¹ * 3 ^ (4 * M.gamma * ↑(n + ↑h)) + M.gamma ^ 15

    The averaged Step-5 display, at the closure's grid carrier and the finite localization cube. This is exactly the binder hdisplay of step5GaugeEndpoint_of_cellDisplay.

    The endpoint #

    theorem Algsuperdiff.Section3.Provider.Diffusivity.ApproximateRecurrence.Closure.step5GaugeEndpoint_of_correctors_gauge {d : ℕ} (M : ABKModel d) (n : ℤ) (h K j : ℕ) {e : HCPolySupport.Vec d} (he : HCPolySupport.vecNormSq e = 1) (nmesh : ℤ) (gauge : HCPolySupport.TriadicCube d → Cutoff.CutoffSample d → ℝ) {C1 C2 Cdef shellMoment gradMoment : ℝ} (wD : C(HCPolySupport.Vec d, HCPolySupport.Mat d) → HCPolySupport.H10Function (HCPolySupport.openCubeSet (HCPolySupport.originCube d ↑K))) (wN : C(HCPolySupport.Vec d, HCPolySupport.Mat d) → HCPolySupport.H1MeanZeroFunction (HCPolySupport.openCubeSet (HCPolySupport.originCube d ↑K))) (hgauge0 : ∀ (R : HCPolySupport.TriadicCube d) (omega : Cutoff.CutoffSample d), 0 ≤ gauge R omega) (hCdef : 0 ≤ Cdef) (hgate : Cdef * M.gamma ^ 15 ≤ 2) (hD : ∀ (omega : Cutoff.ShellSeq d), HCPolySupport.IsZeroTraceDirichletRhsWeakSolution (fun (x : HCPolySupport.Vec d) => 1) (HCPolySupport.openCubeSet (HCPolySupport.originCube d ↑K)) (alongIncrementPath n h wD omega) fun (x : HCPolySupport.Vec d) => -Corrector.streamForcing (↑(Annealed.sigmaBar M n))⁻¹ omega n (n + ↑h) e x) (hN : ∀ (omega : Cutoff.ShellSeq d), HCPolySupport.IsMeanZeroNeumannRhsWeakSolution (fun (x : HCPolySupport.Vec d) => 1) (HCPolySupport.openCubeSet (HCPolySupport.originCube d ↑K)) (alongIncrementPath n h wN omega) fun (x : HCPolySupport.Vec d) => -Corrector.streamForcing (↑(Annealed.sigmaBar M n))⁻¹ omega n (n + ↑h) 0 x) (hmemD : ∀ (omega : Cutoff.ShellSeq d), ∀ R ∈ HCPolySupport.descendantsAtDepth (HCPolySupport.originCube d ↑K) j, MeasureTheory.MemLp (fun (x : HCPolySupport.Vec d) => (alongIncrementPath n h wD omega).grad x) 2 (HCPolySupport.normalizedCubeMeasure R)) (hintD : ∀ (omega : Cutoff.ShellSeq d), MeasureTheory.IntegrableOn (fun (x : HCPolySupport.Vec d) => HCPolySupport.vecNormSq ((alongIncrementPath n h wD omega).grad x)) (HCPolySupport.cubeSet (HCPolySupport.originCube d ↑K)) MeasureTheory.volume) (hintP : ∀ (omega : Cutoff.ShellSeq d), MeasureTheory.IntegrableOn (fun (x : HCPolySupport.Vec d) => HCPolySupport.vecDot e (HCPolySupport.matVecMul ((Cutoff.finiteShellIncrement omega n (n + ↑h)).toFun x) ((alongIncrementPath n h wD omega).grad x))) (HCPolySupport.cubeSet (HCPolySupport.originCube d ↑K)) MeasureTheory.volume) (hsampN : ∀ R ∈ HCPolySupport.descendantsAtDepth (HCPolySupport.originCube d ↑K) j, MeasureTheory.Integrable (fun (omega : Cutoff.CutoffSample d) => HCPolySupport.vecNormSq (HCPolySupport.cubeAverageVec R fun (x : HCPolySupport.Vec d) => (alongIncrementPath n h wD ↑omega).grad x)) ↑(Cutoff.cutoffSampleLaw M)) (hsampX : ∀ R ∈ HCPolySupport.descendantsAtDepth (HCPolySupport.originCube d ↑K) j, MeasureTheory.Integrable (fun (omega : Cutoff.CutoffSample d) => HCPolySupport.vecDot e (HCPolySupport.matVecMul (freshShellCubeAverage R (↑omega) n (n + ↑h)) (HCPolySupport.cubeAverageVec R fun (x : HCPolySupport.Vec d) => (alongIncrementPath n h wD ↑omega).grad x))) ↑(Cutoff.cutoffSampleLaw M)) (hsampS : ∀ R ∈ HCPolySupport.descendantsAtDepth (HCPolySupport.originCube d ↑K) j, MeasureTheory.Integrable (fun (omega : Cutoff.CutoffSample d) => HCPolySupport.vecNormSq (HCPolySupport.matVecMul (freshShellCubeAverage R (↑omega) n (n + ↑h)) (HCPolySupport.cubeAverageVec R fun (x : HCPolySupport.Vec d) => (alongIncrementPath n h wD ↑omega).grad x))) ↑(Cutoff.cutoffSampleLaw M)) (hsampDef : ∀ R ∈ HCPolySupport.descendantsAtDepth (HCPolySupport.originCube d ↑K) j, MeasureTheory.Integrable (fun (omega : Cutoff.CutoffSample d) => step5CrossDefect R (↑omega) n (n + ↑h) e fun (x : HCPolySupport.Vec d) => (alongIncrementPath n h wD ↑omega).grad x) ↑(Cutoff.cutoffSampleLaw M)) (hsampG : MeasureTheory.Integrable (fun (omega : Cutoff.CutoffSample d) => HCPolySupport.cubeAverage (HCPolySupport.originCube d ↑K) fun (x : HCPolySupport.Vec d) => HCPolySupport.vecNormSq ((alongIncrementPath n h wD ↑omega).grad x)) ↑(Cutoff.cutoffSampleLaw M)) (hdefect : ∀ R ∈ HCPolySupport.descendantsAtDepth (HCPolySupport.originCube d ↑K) j, ∀ (omega : Cutoff.CutoffSample d), (step5CrossDefect R (↑omega) n (n + ↑h) e fun (x : HCPolySupport.Vec d) => (alongIncrementPath n h wD ↑omega).grad x) ≤ Cdef * (meshOscillationCell nmesh (fun (x : HCPolySupport.Vec d) => (alongIncrementPath n h wD ↑omega).grad x) R * gauge R omega)) (hmemOsc : ∀ R ∈ HCPolySupport.descendantsAtDepth (HCPolySupport.originCube d ↑K) j, MeasureTheory.MemLp (fun (omega : Cutoff.CutoffSample d) => meshOscillationCell nmesh (fun (x : HCPolySupport.Vec d) => (alongIncrementPath n h wD ↑omega).grad x) R) 4 ↑(Cutoff.cutoffSampleLaw M)) (hmemGauge : ∀ R ∈ HCPolySupport.descendantsAtDepth (HCPolySupport.originCube d ↑K) j, MeasureTheory.MemLp (fun (omega : Cutoff.CutoffSample d) => gauge R omega) 4 ↑(Cutoff.cutoffSampleLaw M)) (hosc : (gridFourthMomentRoot (↑(Cutoff.cutoffSampleLaw M)) (HCPolySupport.descendantsAtDepth (HCPolySupport.originCube d ↑K) j) fun (R : HCPolySupport.TriadicCube d) (omega : Cutoff.CutoffSample d) => meshOscillationCell nmesh (fun (x : HCPolySupport.Vec d) => (alongIncrementPath n h wD ↑omega).grad x) R) + (↑(Annealed.sigmaBar M n))⁻¹ * gridFourthMomentRoot (↑(Cutoff.cutoffSampleLaw M)) (HCPolySupport.descendantsAtDepth (HCPolySupport.originCube d ↑K) j) gauge ≤ M.gamma ^ 15) (hC1 : 0 ≤ C1) (hgrad0 : 0 ≤ gradMoment) (hprod : (HCPolySupport.descendantsAverage (HCPolySupport.originCube d ↑K) j fun (R : HCPolySupport.TriadicCube d) => ∫ (omega : Cutoff.CutoffSample d), HCPolySupport.vecNormSq (HCPolySupport.matVecMul (freshShellCubeAverage R (↑omega) n (n + ↑h)) (HCPolySupport.cubeAverageVec R fun (x : HCPolySupport.Vec d) => (alongIncrementPath n h wD ↑omega).grad x)) ∂↑(Cutoff.cutoffSampleLaw M)) ≤ shellMoment * gradMoment) (hshell : shellMoment ≤ C1 * ↑h * 3 ^ (2 * M.gamma * ↑(n + ↑h))) (hgradM : gradMoment ≤ C2 * (↑(Annealed.sigmaBar M n) ^ 2)⁻¹ * ↑h * 3 ^ (2 * M.gamma * ↑(n + ↑h))) :
    Step5GaugeEndpoint M n h K j e (C1 * C2) wD wN

    Closure.Step5GaugeEndpoint, from the two correctors and the named legs, at a per-cube forcing gauge. This is the Step-5 mirror of Closure.JunctionDischargeStepFour.step4GaugeEndpoint_of_correctors.

    The gauge family gauge R omega is the right factor of e.nablaw.oscillations read at the cube R's own base point; the producer instantiates it at `3^nmesh

    • shellDerivNormSum n (n+h) (triadicCubeShift R) omega`.