Documentation

LeanPool.Superdiffusion.Algsuperdiff.Section5.Support.IterationStep

IterationStep #

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

One step of the perturbation iteration on y + □_n #

Section 5.1 compares the solution of the rough-field Dirichlet problem at scale m with the one at scale n ≤ m by treating k_m - k_n as a perturbation and iterating. One step of that iteration takes an H¹ function w on y + □_n with a continuous representative and produces

Both halves are assembled here. The Laplace half is the cube Schauder estimate Algsuperdiff.Frozen.External.cube_schauder at s = 1/2, unit coefficient and zero boundary datum, transported from the origin cube to y + □_n; the size of its right-hand side is the Hölder Leibniz rule of HolderCalculus; the antisymmetry of the perturbation rewrites the forcing field through SkewTransfer; and the next iterate exists by DirichletSolvability.

Main results #

References #

1. Reading a solution of the origin-cube problem on y + □_n #

The C^{0,α} bound transports from the origin cube to y + □_n.

A C^{0,α} bound depends on the field only through its values on the set.

2. The Laplace problem on y + □_n #

theorem Algsuperdiff.Section5.Support.exists_isDirichletSolutionAt_laplacian {d : ℕ} (hdim : 2 ≤ d) :
∃ (C : ℝ), 0 < C ∧ ∀ (y : HCPolySupport.Vec d) (n : ℤ) (G : HCPolySupport.Vec d → HCPolySupport.Vec d) (K : ℝ), Section4.Support.HolderSeminormBoundOn (cubeSetAt y n) (1 / 2) K G → ∃ (h : HCPolySupport.H1Function (cubeSetAt y n)), IsDirichletSolutionAt (fun (x : HCPolySupport.Vec d) => 1 • 1) y n h G ∧ (∀ x ∈ cubeSetAt y n, ‖h.grad x‖ ≤ C * Real.rpow 3 (↑n / 2) * K) ∧ Section4.Support.HolderSeminormBoundOn (cubeSetAt y n) (1 / 2) (C * K) h.grad

The Laplace problem in the cube y + □_n with a 1/2-Hölder forcing field.

For every G with [G]_{C^{0,1/2}(y+□_n)} ≤ K there is a zero-boundary weak solution h of -Δ h = ∇·G on y + □_n whose gradient obeys

  ‖∇h‖_{L^∞(y+□_n)} ≤ C 3^{n/2} K ,   [∇h]_{C^{0,1/2}(y+□_n)} ≤ C K ,

with C = C(d) independent of the scale, the centre and the forcing field.

This is the boundary Schauder estimate for the Poisson equation in a cube at s = 1/2, unit coefficient and zero boundary datum, read on y + □_n through the translation bridge of Translation.lean. The two displayed powers of the scale are exactly what the frozen estimate produces: its left-hand gauge is 3^{-n/2} ‖∇h‖_{L^∞} + [∇h]_{C^{0,1/2}} and its right-hand side at zero datum is C 3^{-n/2} · 3^{n/2} K = C K.

3. The size of the forcing field of the iteration #

theorem Algsuperdiff.Section5.Support.holderSeminormBoundOn_neg_smul {d : ℕ} {E : Type u_1} [NormedAddCommGroup E] [NormedSpace ℝ E] {U : Set (HCPolySupport.Vec d)} {alpha Kf Sf KG SG : ℝ} {f : HCPolySupport.Vec d → ℝ} {G : HCPolySupport.Vec d → E} (hKf : 0 ≤ Kf) (hSf : 0 ≤ Sf) (hf : Section4.Support.HolderSeminormBoundOn U alpha Kf f) (hfsup : ∀ x ∈ U, ‖f x‖ ≤ Sf) (hG : Section4.Support.HolderSeminormBoundOn U alpha KG G) (hGsup : ∀ x ∈ U, ‖G x‖ ≤ SG) :
Section4.Support.HolderSeminormBoundOn U alpha (Kf * SG + Sf * KG) fun (x : HCPolySupport.Vec d) => -f x • G x

The Hölder Leibniz rule in the form the iteration needs. The forcing field of the Laplace step is -w ∇·(k_m - k_n), a scalar field times a vector field.

4. One step of the iteration #

theorem Algsuperdiff.Section5.Support.exists_iterationStep {d : ℕ} (hdim : 2 ≤ d) :
∃ (C : ℝ), 0 < C ∧ ∀ (M : Section3.ABKModel d) (n : ℤ) (y : HCPolySupport.Vec d) (omega : Section3.Cutoff.CutoffSample d) (k : HCPolySupport.Vec d → HCPolySupport.Mat d), (∀ (p q : Fin d), ContDiff ℝ 1 fun (x : HCPolySupport.Vec d) => k x p q) → (∀ (x : HCPolySupport.Vec d), HCPolySupport.matTranspose (k x) = -k x) → ∀ (w : HCPolySupport.H1Function (cubeSetAt y n)) (wRep : HCPolySupport.Vec d → ℝ) (Kw Sw Kk Sk : ℝ), IsCubeRepresentative y n w wRep → Section4.Support.HolderSeminormBoundOn (cubeSetAt y n) (1 / 2) Kw wRep → (∀ x ∈ cubeSetAt y n, ‖wRep x‖ ≤ Sw) → Section4.Support.HolderSeminormBoundOn (cubeSetAt y n) (1 / 2) Kk (matFieldDiv k) → (∀ x ∈ cubeSetAt y n, ‖matFieldDiv k x‖ ≤ Sk) → ∃ (h : HCPolySupport.H1Function (cubeSetAt y n)) (w' : HCPolySupport.H1Function (cubeSetAt y n)), (IsDirichletSolutionAt (fun (x : HCPolySupport.Vec d) => 1 • 1) y n h fun (x : HCPolySupport.Vec d) => -HCPolySupport.matVecMul (k x) (w.grad x)) ∧ (∀ x ∈ cubeSetAt y n, ‖h.grad x‖ ≤ C * Real.rpow 3 (↑n / 2) * (Kw * Sk + Sw * Kk)) ∧ Section4.Support.HolderSeminormBoundOn (cubeSetAt y n) (1 / 2) (C * (Kw * Sk + Sw * Kk)) h.grad ∧ IsDirichletSolutionAt (Section3.Cutoff.coefficientCutoff M.nu n omega).toCoeffField y n w' h.grad

One step of the perturbation iteration.

Given a skew, entrywise C¹ matrix field k and an H¹ function w on y + □_n with a continuous representative wRep obeying [wRep] ≤ Kw and ‖wRep‖_{L^∞} ≤ Sw, and given the two sizes [∇·k] ≤ Kk, ‖∇·k‖_{L^∞} ≤ Sk of the divergence of the perturbation, the step produces

  • h, the zero-boundary solution of -Δ h = -∇·(w ∇·k) on y + □_n, stated in the equivalent antisymmetric form -∇·(k ∇w) that the antisymmetry transfer supplies, with

      ‖∇h‖_{L^∞(y+□_n)} ≤ C 3^{n/2} ([wRep] ‖∇·k‖_∞ + ‖wRep‖_∞ [∇·k]) ,
      [∇h]_{C^{0,1/2}(y+□_n)} ≤ C ([wRep] ‖∇·k‖_∞ + ‖wRep‖_∞ [∇·k]);
    
  • w', the next iterate: the solution of the rough-field problem at scale n with forcing field ∇h.

C = C(d) is the constant of the cube Schauder estimate; it does not depend on the scale, the centre, the sample, the perturbation or w.

5. The perturbation of the iteration is the shell increment #

Every entry of the shell increment k_q - k_p is C¹, the first hypothesis of the iteration step.

The shell increment k_q - k_p is pointwise antisymmetric, the second hypothesis of the iteration step.

6. The good cube event bounds the next iterate #

On the good cube event the localized regularity is at most 2 C_reg, and this bounds the next iterate.

  σ_n 3^{-n/2} [w']_{C^{0,1/2}(y+□_n)}
      ≤ 2 C_reg 3^{n/2} [∇h]_{C^{0,1/2}(y+□_n)} .

On the good cube event the next iterate has a continuous representative. The infimum defining the localized regularity is ⊤ when none exists.

theorem Algsuperdiff.Section5.Support.exists_iterationStep_of_mem_goodCubeEvent {d : ℕ} (hdim : 2 ≤ d) :
∃ (C : ℝ), 0 < C ∧ ∀ (M : Section3.ABKModel d) (Creg : ℝ) (n : ℤ) (y : HCPolySupport.Vec d) (ep : ℝ), ∀ omega ∈ goodCubeEvent M Creg n y ep, ∀ (k : HCPolySupport.Vec d → HCPolySupport.Mat d), (∀ (p q : Fin d), ContDiff ℝ 1 fun (x : HCPolySupport.Vec d) => k x p q) → (∀ (x : HCPolySupport.Vec d), HCPolySupport.matTranspose (k x) = -k x) → ∀ (w : HCPolySupport.H1Function (cubeSetAt y n)) (wRep : HCPolySupport.Vec d → ℝ) (Kw Sw Kk Sk : ℝ), IsCubeRepresentative y n w wRep → Section4.Support.HolderSeminormBoundOn (cubeSetAt y n) (1 / 2) Kw wRep → (∀ x ∈ cubeSetAt y n, ‖wRep x‖ ≤ Sw) → Section4.Support.HolderSeminormBoundOn (cubeSetAt y n) (1 / 2) Kk (matFieldDiv k) → (∀ x ∈ cubeSetAt y n, ‖matFieldDiv k x‖ ≤ Sk) → ∃ (h : HCPolySupport.H1Function (cubeSetAt y n)) (w' : HCPolySupport.H1Function (cubeSetAt y n)) (w'Rep : HCPolySupport.Vec d → ℝ), (IsDirichletSolutionAt (fun (x : HCPolySupport.Vec d) => 1 • 1) y n h fun (x : HCPolySupport.Vec d) => -HCPolySupport.matVecMul (k x) (w.grad x)) ∧ (∀ x ∈ cubeSetAt y n, ‖h.grad x‖ ≤ C * Real.rpow 3 (↑n / 2) * (Kw * Sk + Sw * Kk)) ∧ Section4.Support.HolderSeminormBoundOn (cubeSetAt y n) (1 / 2) (C * (Kw * Sk + Sw * Kk)) h.grad ∧ IsDirichletSolutionAt (Section3.Cutoff.coefficientCutoff M.nu n omega).toCoeffField y n w' h.grad ∧ IsCubeRepresentative y n w' w'Rep ∧ ENNReal.ofReal (↑(Section3.Annealed.sigmaBar M n) * Real.rpow 3 (-↑n / 2)) * holderSeminormOn (cubeSetAt y n) (1 / 2) w'Rep ≤ ENNReal.ofReal (2 * Creg) * ENNReal.ofReal (Real.rpow 3 (↑n / 2)) * holderSeminormOn (cubeSetAt y n) (1 / 2) h.grad

One step of the iteration on the good cube event. The step produces the Laplace solution h, the next iterate w', a continuous representative of w', and the estimate of the proof of Proposition 5.1 relating the two.