Documentation

LeanPool.Superdiffusion.Algsuperdiff.Section24.Sensitivity.Provider.ResponseUnconditional.Provider

Provider #

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

The unconditional response estimate #

Source: ABK26, the J half of l.J.sensitivity.no.conditions.

Provider.ResponseUnconditional.Closure reduces the frozen target Algsuperdiff.Frozen.Section24.responseJ_sensitivity_unconditional to a single analytic input, the per-mesoscopic-cube conditional J-sensitivity at the normalized loading p = (μσ₀)^{-1/2} e, q = (μσ₀)^{1/2} e. Provider.ResponseUnconditional.CubeSensitivity proves exactly that input at a general loading. This module instantiates the loading and assembles the provider.

At the normalized loading

|p|² = (μσ₀)^{-1} ,      |p·q| = |e|² = 1 ,

which is precisely the pairing of the two K₁ groups of the frozen statement, and at a cube of depth H below the unit cube |R|² = 3^{-2H}, the damping factor of the frozen statement.

The value budget is Provider.ResponseUnconditional.valueBudget, whose depth-H descendant average over the unit cube is exactly ‖h‖²_{L²(□₀)}.

The normalized loading #

The per-mesoscopic-cube input #

theorem Algsuperdiff.Section24.Sensitivity.Provider.ResponseUnconditional.mesoscale_cube_sensitivity {d : ℕ} (dimension : 2 ≤ d) :
∃ (C₀ : ℝ) (K₁ : ℝ), 1 ≤ C₀ ∧ 0 ≤ K₁ ∧ ∀ (a : HCPolySupport.Book.Ch02.CoeffOn (HCPolySupport.Book.Ch02.cubeDomain (HCPolySupport.originCube d 0))) (h : Frozen.Section24.UnitCubeSkewW2Infinity d) (F G : HCPolySupport.Book.Ch02.TriadicCoeffFamily d), (F.coeffOn (HCPolySupport.originCube d 0)).AEEq a → (G.coeffOn (HCPolySupport.originCube d 0)).AEEq (Frozen.Section24.perturbCoeffOn (HCPolySupport.Book.Ch02.cubeDomain (HCPolySupport.originCube d 0)) a h.toLInfSkewMatrixFieldOn 1) → ∀ (s μ σ0 : ℝ) (e : HCPolySupport.Vec d), 0 < s → s ≤ 1 / 4 → 0 < μ → 0 < σ0 → HCPolySupport.Book.Ch02.vecNorm e = 1 → ∃ (V : HCPolySupport.TriadicCube d → ℝ), (∀ (R : HCPolySupport.TriadicCube d), 0 ≤ V R) ∧ HCPolySupport.descendantsAverage (HCPolySupport.originCube d 0) (Mesoscale.mesoscaleDepth s (1 + C₀ * (h.gradientW1Infinity * (Frozen.Section24.unitCubeLambda s (HCPolySupport.Book.Ch02.MultiscaleExponent.finite 2) a)⁻¹))) V ≤ h.valueL2 ^ 2 ∧ ∀ R ∈ HCPolySupport.descendantsAtDepth (HCPolySupport.originCube d 0) (Mesoscale.mesoscaleDepth s (1 + C₀ * (h.gradientW1Infinity * (Frozen.Section24.unitCubeLambda s (HCPolySupport.Book.Ch02.MultiscaleExponent.finite 2) a)⁻¹))), HCPolySupport.Book.Ch02.responseJ (HCPolySupport.Book.Ch02.cubeDomain R) (G.coeffOn R) ((√(μ * σ0))⁻¹ • e) (√(μ * σ0) • e) ≤ 3 * HCPolySupport.Book.Ch02.responseJ (HCPolySupport.Book.Ch02.cubeDomain R) (F.coeffOn R) ((√(μ * σ0))⁻¹ • e) (√(μ * σ0) • e) + K₁ * ((μ * σ0)⁻¹ * (HCPolySupport.Book.Ch02.lambdaSq R (3 / 8) (HCPolySupport.Book.Ch02.MultiscaleExponent.finite 2) F)⁻¹ * (V R + Real.rpow 3 (-(2 * ↑(Mesoscale.mesoscaleDepth s (1 + C₀ * (h.gradientW1Infinity * (Frozen.Section24.unitCubeLambda s (HCPolySupport.Book.Ch02.MultiscaleExponent.finite 2) a)⁻¹))))) * h.gradientW1Infinity ^ 2)) + K₁ * (Real.rpow 3 (-(2 * ↑(Mesoscale.mesoscaleDepth s (1 + C₀ * (h.gradientW1Infinity * (Frozen.Section24.unitCubeLambda s (HCPolySupport.Book.Ch02.MultiscaleExponent.finite 2) a)⁻¹))))) * (HCPolySupport.Book.Ch02.lambdaSq R (3 / 8) (HCPolySupport.Book.Ch02.MultiscaleExponent.finite 2) F)⁻¹ ^ 2 * h.gradientW1Infinity ^ 2)

The per-mesoscopic-cube conditional J-sensitivity at the normalized loading, in exactly the shape consumed by responseJ_sensitivity_unconditional_of_mesoscale_cube_sensitivity.

The provider #

theorem Algsuperdiff.Section24.Sensitivity.Provider.ResponseUnconditional.responseJ_sensitivity_unconditional {d : ℕ} (dimension : 2 ≤ d) :
∃ (C : ℝ), 0 < C ∧ ∀ (a : HCPolySupport.Book.Ch02.CoeffOn (HCPolySupport.Book.Ch02.cubeDomain (HCPolySupport.originCube d 0))) (h : Frozen.Section24.UnitCubeSkewW2Infinity d) (s t μ σ0 : ℝ) (e : HCPolySupport.Vec d), 0 < s → s ≤ 1 / 4 → 0 < t → t ≤ 1 / 4 → 0 < μ → 0 < σ0 → HCPolySupport.Book.Ch02.vecNorm e = 1 → HCPolySupport.Book.Ch02.responseJ (HCPolySupport.Book.Ch02.cubeDomain (HCPolySupport.originCube d 0)) (Frozen.Section24.perturbCoeffOn (HCPolySupport.Book.Ch02.cubeDomain (HCPolySupport.originCube d 0)) a h.toLInfSkewMatrixFieldOn 1) ((√(μ * σ0))⁻¹ • e) (√(μ * σ0) • e) ≤ C * μ⁻¹ * (1 + h.gradientW1Infinity * (Frozen.Section24.unitCubeLambda s (HCPolySupport.Book.Ch02.MultiscaleExponent.finite 2) a)⁻¹).rpow (2 * t / (1 - 2 * s)) * Frozen.Section24.unitCubeHomogenizationError t (HCPolySupport.Book.Ch02.MultiscaleExponent.finite 2) (HCPolySupport.Book.Ch02.MultiscaleExponent.finite 2) a (HCPolySupport.scalarMatrix σ0) ^ 2 + C * μ⁻¹ * σ0 * (Frozen.Section24.unitCubeLambda s (HCPolySupport.Book.Ch02.MultiscaleExponent.finite 2) a)⁻¹ * (1 + h.gradientW1Infinity * (Frozen.Section24.unitCubeLambda s (HCPolySupport.Book.Ch02.MultiscaleExponent.finite 2) a)⁻¹).rpow (2 * s / (1 - 2 * s)) * (σ0⁻¹ ^ 2 * h.valueL2 ^ 2 + (μ - 1) ^ 2) + C * μ⁻¹ ^ 2 * σ0⁻¹ ^ 2 * h.gradientW1Infinity ^ 2 + C * min 1 (h.gradientW1Infinity * (Frozen.Section24.unitCubeLambda s (HCPolySupport.Book.Ch02.MultiscaleExponent.finite 2) a)⁻¹) ^ 2

Unconditional four-term sensitivity estimate for the response functional.

This is the provider form of the frozen target Algsuperdiff.Frozen.Section24.responseJ_sensitivity_unconditional, with exactly the frozen binders.