Documentation

LeanPool.Superdiffusion.Algsuperdiff.Section4.Provider.ExcessDecay.HarmonicApproximationProvider

HarmonicApproximationProvider #

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

l.harmonic.approximation.good.scales, the provider #

What is proved #

One public theorem, harmonic_approximation_provider, and it carries no hypothesis beyond the anchor's own regime, funding line and geometry binders: for every dimension d there is a constant C > 0 such that, almost surely in the cutoff sample, every Dirichlet solution u on □_m with datum h and force g, and every weakly ã-harmonic competitor v = u − w on the child cube x + □_n (w ∈ H¹₀), satisfies two displays.

Here 𝓔_{n+3} = fluxCorrectedErrorRepresentative M L (n+3) ⟨s/8⟩ (τ_z ω) is the coarse-graining error factor read at the clause's own flux slot, and σ_n is the annealed comparator at the child scale.

The flux-corrected error factor on the flat leg #

The last leg is the manuscript's scale-weighted window average of ∇h, carried multiplied by 𝓔. Bounding 𝓔 · 1_𝒢 by a constant multiple of the event's own threshold is a monotone one-way step; the chain below declines it on that row, so the leg vanishes with the flux-corrected error instead of being priced against it. The row is re-proved rather than transferred, through AssemblyComposedErrorWeighted and ResidueRouteOneErrorWeighted (the two composition sites), ResiduePluginErrorWeighted and StitchDisplayLegsErrorWeighted (the residue lane's second consumer), AssemblyBoundaryClauseErrorWeighted and ProviderBoundaryErrorWeighted (the repackaging), and RebaseGeneralClauseErrorWeighted, RebaseJoinErrorWeighted, RebaseSealDistanceErrorWeighted (statement-only re-typings, the leg being pure nonnegative slack on the frontier-empty branch).

The last step, and what it costs #

exists_harmonicApproximation_provider_of_boundaryClauseReal_errorWeighted delivers the display with that leg still in normalized-L² form. The mean-zero fractional Poincaré inequality on the n+3 window (GradHSplit) turns it into the printed average leg plus a multiple of the display's own Gagliardo ∇h leg,

  ‖∇h‖_{L²(W')} ≤ ‖(∇h)_W‖ + C(d)·3^{ns}·[∇h]_{H^s(W')} ,

and 3^n · 3^{ns} = 3^{(1+s)n} is exactly that leg's weight, so no s-power is consumed. The remainder carries 𝓔, which is shed under the good-event indicator by the proved caps step ReindexSlot.ae_errorRepresentative_le_harmonicSlot_addThree (𝓔 ≤ C_𝓔/2 there). The whole price is the constant re-choice C := max ((1 + κ(d)·C_𝓔/2)·C₀) (max 2 C_𝓔), a function of d alone.

References #

theorem Algsuperdiff.Section4.Provider.ExcessDecay.harmonic_approximation_provider (d : ℕ) :
∃ (C : ℝ), 0 < C ∧ ∀ (M : Section3.ABKModel d), ∀ s ∈ Set.Icc (64 * M.gamma) 1, M.gamma ≤ C⁻¹ * Section3.Disorder.cstar M ^ 10 → M.gamma * |Real.log M.gamma| ^ 2 ≤ (s / 8).rpow (3 / 2) * Section3.Disorder.cstar M ^ 2 * (C⁻¹ * s ^ 4) → ∀ (hs : 0 < s) (L m n : ℤ), m ≤ L → n + 3 ≤ m → ∀ (x z : HCPolySupport.Vec d), x ∈ HCPolySupport.openCubeSet (HCPolySupport.originCube d m) → z ∈ HCPolySupport.openCubeSet (HCPolySupport.originCube d m) → (fun (y : HCPolySupport.Vec d) => x + y) '' HCPolySupport.openCubeSet (HCPolySupport.originCube d n) ⊆ (fun (y : HCPolySupport.Vec d) => z + y) '' HCPolySupport.openCubeSet (HCPolySupport.originCube d (n + 1)) ∩ HCPolySupport.openCubeSet (HCPolySupport.originCube d m) → ∀ᵐ (omega : Section3.Cutoff.CutoffSample d) ∂↑(Section3.Cutoff.cutoffSampleLaw M), ∀ (u h : HCPolySupport.H1Function (HCPolySupport.openCubeSet (HCPolySupport.originCube d m))) (g : HCPolySupport.Vec d → HCPolySupport.Vec d), Support.IsDirichletSolutionOn (Section3.Cutoff.coefficientCutoff M.nu L omega).toCoeffField (HCPolySupport.originCube d m) u h g → MeasureTheory.MemLp g 2 (Support.normalizedVolumeMeasureOn (HCPolySupport.openCubeSet (HCPolySupport.originCube d m))) → MeasureTheory.MemLp (HCPolySupport.Gagliardo.gagliardoKernel s 2 g) 2 (Support.normalizedGagliardoMeasureOn (HCPolySupport.openCubeSet (HCPolySupport.originCube d m))) → MeasureTheory.MemLp (HCPolySupport.Gagliardo.gagliardoKernel s 2 h.grad) 2 (Support.normalizedGagliardoMeasureOn (HCPolySupport.openCubeSet (HCPolySupport.originCube d m))) → ∀ (v : HCPolySupport.H1Function ((fun (y : HCPolySupport.Vec d) => x + y) '' HCPolySupport.openCubeSet (HCPolySupport.originCube d n))) (w : HCPolySupport.H10Function ((fun (y : HCPolySupport.Vec d) => x + y) '' HCPolySupport.openCubeSet (HCPolySupport.originCube d n))), Support.IsWeaklyHarmonicOn ((fun (y : HCPolySupport.Vec d) => x + y) '' HCPolySupport.openCubeSet (HCPolySupport.originCube d n)) v → (∀ (y : HCPolySupport.Vec d), v.toFun y = u.toFun y - w.toFun y) → (∀ (y : HCPolySupport.Vec d), v.grad y = u.grad y - w.grad y) → (Frozen.Section4.goodEventAt M (Support.cgEllipLowerConstant d) (n + 3) z ⟨s / 8, ⋯⟩ (C⁻¹ * s ^ 4)).indicator (fun (_omega' : Section3.Cutoff.CutoffSample d) => MeasureTheory.eLpNorm (fun (y : HCPolySupport.Vec d) => u.toFun y - v.toFun y) 2 (Support.normalizedVolumeMeasureOn ((fun (y : HCPolySupport.Vec d) => x + y) '' HCPolySupport.openCubeSet (HCPolySupport.originCube d n)))) omega ≤ ENNReal.ofReal (C * s.rpow (-4) * Support.fluxCorrectedErrorRepresentative M L (n + 3) ⟨s / 8, ⋯⟩ (Section3.Cutoff.translateCutoffSample z omega)) * (MeasureTheory.eLpNorm (fun (y : HCPolySupport.Vec d) => u.toFun y - HCPolySupport.volumeAverage ((fun (y' : HCPolySupport.Vec d) => z + y') '' HCPolySupport.openCubeSet (HCPolySupport.originCube d (n + 3)) ∩ HCPolySupport.openCubeSet (HCPolySupport.originCube d m)) u.toFun) 2 (Support.normalizedVolumeMeasureOn ((fun (y' : HCPolySupport.Vec d) => z + y') '' HCPolySupport.openCubeSet (HCPolySupport.originCube d (n + 3)) ∩ HCPolySupport.openCubeSet (HCPolySupport.originCube d m))) + ENNReal.ofReal (s.rpow (-(3 / 2)) * Real.rpow 3 ↑n * ‖HCPolySupport.volumeAverageVec ((fun (y' : HCPolySupport.Vec d) => z + y') '' HCPolySupport.openCubeSet (HCPolySupport.originCube d (n + 2)) ∩ HCPolySupport.openCubeSet (HCPolySupport.originCube d m)) h.grad‖)) + ENNReal.ofReal (C * s.rpow (-7) * (↑(Section3.Annealed.sigmaBar M n))⁻¹ * Real.rpow 3 ((1 + s) * ↑n)) * Support.normalizedGagliardoESeminormOn ((fun (y' : HCPolySupport.Vec d) => z + y') '' HCPolySupport.openCubeSet (HCPolySupport.originCube d (n + 3)) ∩ HCPolySupport.openCubeSet (HCPolySupport.originCube d m)) s g + ENNReal.ofReal (C * s.rpow (-6) * Real.rpow 3 ((1 + s) * ↑n)) * Support.normalizedGagliardoESeminormOn ((fun (y' : HCPolySupport.Vec d) => z + y') '' HCPolySupport.openCubeSet (HCPolySupport.originCube d (n + 3)) ∩ HCPolySupport.openCubeSet (HCPolySupport.originCube d m)) s h.grad + ENNReal.ofReal (C * s.rpow (-6) * Support.fluxCorrectedErrorRepresentative M L (n + 3) ⟨s / 8, ⋯⟩ (Section3.Cutoff.translateCutoffSample z omega) * Real.rpow 3 ↑n * ‖HCPolySupport.volumeAverageVec ((fun (y' : HCPolySupport.Vec d) => z + y') '' HCPolySupport.openCubeSet (HCPolySupport.originCube d (n + 2)) ∩ HCPolySupport.openCubeSet (HCPolySupport.originCube d m)) h.grad‖) ∧ ((fun (y' : HCPolySupport.Vec d) => z + y') '' HCPolySupport.openCubeSet (HCPolySupport.originCube d (n + 2)) ∩ frontier (HCPolySupport.openCubeSet (HCPolySupport.originCube d m)) = ∅ → (Frozen.Section4.goodEventAt M (Support.cgEllipLowerConstant d) (n + 2) z ⟨s / 8, ⋯⟩ (1 / 2)).indicator (fun (_omega' : Section3.Cutoff.CutoffSample d) => MeasureTheory.eLpNorm (fun (y : HCPolySupport.Vec d) => u.toFun y - v.toFun y) 2 (Support.normalizedVolumeMeasureOn ((fun (y : HCPolySupport.Vec d) => x + y) '' HCPolySupport.openCubeSet (HCPolySupport.originCube d n)))) omega ≤ ENNReal.ofReal (C * s.rpow (-4) * Support.fluxCorrectedErrorRepresentative M L (n + 2) ⟨s / 8, ⋯⟩ (Section3.Cutoff.translateCutoffSample z omega)) * MeasureTheory.eLpNorm (fun (y : HCPolySupport.Vec d) => u.toFun y - HCPolySupport.volumeAverage ((fun (y' : HCPolySupport.Vec d) => z + y') '' HCPolySupport.openCubeSet (HCPolySupport.originCube d (n + 2)) ∩ HCPolySupport.openCubeSet (HCPolySupport.originCube d m)) u.toFun) 2 (Support.normalizedVolumeMeasureOn ((fun (y' : HCPolySupport.Vec d) => z + y') '' HCPolySupport.openCubeSet (HCPolySupport.originCube d (n + 2)) ∩ HCPolySupport.openCubeSet (HCPolySupport.originCube d m))) + ENNReal.ofReal (C * s.rpow (-(19 / 2)) * (↑(Section3.Annealed.sigmaBar M n))⁻¹ * Real.rpow 3 ((1 + s) * ↑n)) * Support.normalizedGagliardoESeminormOn ((fun (y' : HCPolySupport.Vec d) => z + y') '' HCPolySupport.openCubeSet (HCPolySupport.originCube d (n + 2)) ∩ HCPolySupport.openCubeSet (HCPolySupport.originCube d m)) s g)

l.harmonic.approximation.good.scales, the provider — unconditional.

The harmonic-approximation display on good scales, with the general clause's flat leg carried as the manuscript's own scale-weighted window average of ∇h at the n+2 window, multiplied by the clause's own coarse-graining error factor 𝓔. Proved by running the boundary chain with the monotone error discard declined on that one row (ProviderBoundaryErrorWeighted, RebaseSealDistanceErrorWeighted) and then applying the fractional Poincaré split of GradHSplit.