High-contrast homogenization: #
Provider.PolynomialHomogenization.Root.Corrector.C1CanonicalProjectionAbsoluteDecay
Imported from the Apache-2.0 HighContrastHomogenization development at commit
7a13dbcd8d6609264a713373f5c69ceeac870472.
Absolute decay of the canonical finite projection #
The same-projection decay is combined with the scale-uniform energy bounds for the joint corrector and the finite affine solution. The resulting bound depends only on the slope and retains the real decay exponent.
theorem
HCPolySupport.HighContrast.Root.exists_canonicalFiniteCorrectorAbsoluteDecayConstants
(d : ℕ)
[NeZero d]
(s eta : ℝ)
(hs : 0 < s)
(hs_lt : s < 1 / 2)
(hetaHalf : 1 / 2 ≤ eta)
(heta1 : eta < 1)
:
∃ (C : ℝ) (c : ℝ),
0 < C ∧ c ∈ Set.Ioo 0 1 ∧ ∀ (a : Book.Ch02.TriadicCoeffFamily d) (delta : ℝ) (n0 : ℤ),
delta ∈ Set.Ioc 0 c →
ScalarIdentityGoodTail a s delta n0 →
∀ (hCauchy : FiniteAffineCorrectionLocalCauchy a) (Phi : Vec d → NormalizedLocalH1Carrier d),
IsFiniteAffineCorrectionJointLocalEquation a Phi →
∀ (n m : ℕ),
n0 ≤ ↑n →
∀ (hnm2 : n + 2 ≤ m) (hT : Function.Injective ⇑(finiteTrialHarmonicGradientLinearMap a ⋯)),
have P := finiteCorrectorWeightedProjection a hCauchy ⋯ hT;
∀ (e : Vec d),
∀ q ∈ Finset.Icc n m,
(weightedGradNorm (a.coeffOn (originCube d ↑q)).toCoeffField (openCubeSet (originCube d ↑q))
fun (x : Vec d) =>
e + (Phi e).globalGradientRepresentative x - (finiteAffineSolution a (↑m) (P e)).toH1.grad x) ≤ ENNReal.ofReal (C * Real.rpow 3 (-eta * ↑(m - q)) * euclideanNorm e)
The canonical base projection approximates the joint affine-corrector gradient at every intermediate cube with a slope-proportional bound.