Documentation

LeanPool.EllipticPDE.Regularity.LocalWeakForm

Localising the bilinear pairing to plain integrals #

The weak formulation a solution satisfies pairs H1amb Ω vectors through the bounded bilinear form Op.fullBilin. Every localised step of interior regularity instead works with plain integrals over a compact V ⊆ Ω against a test function supported in V, on Lp ℝ 2 (volume.restrict V) classes. This file is the bridge between the two.

The derivation is mechanical. A test function supported in V has its graph in H₀¹(Ω), so the weak formulation applies to it; EllipticCoeff.bilin_apply and FullEllipticOp.lowerBilin_apply expand the pairing into inner products; inner_actL_eq and inner_mulCoeffL_eq turn each of those into an integral over Ω; each integrand has a factor vanishing off tsupport v ⊆ V, so the integral shrinks to V; and the restriction of the whole-space extension of each coordinate agrees with that coordinate on V.

The localised identity is the hypothesis hLoc of differentiated_weakForm and differentiated_weakForm_div, which nothing discharged before. Discharging it makes the differentiated-equation identity of Evans, Partial Differential Equations (2nd ed.), §6.3.1, Theorem 2, reachable from the weak formulation itself.

Main declarations #

Localised plain-integral weak identity #

theorem EllipticPdes.Regularity.localWeakForm_of_fullBilin {d : ℕ} (Op : Sobolev.FullEllipticOp d) {Ω : Set (EuclideanSpace ℝ (Fin d))} (hΩm : MeasurableSet Ω) {V : Set (EuclideanSpace ℝ (Fin d))} (hVm : MeasurableSet V) (hVΩ : V ⊆ Ω) (u : ↥(Sobolev.H01 Ω)) (f : Sobolev.L2D Ω) (hu : ∀ (w : ↥(Sobolev.H01 Ω)), ((Op.fullBilin Ω) u) w = ∫ (x : EuclideanSpace ℝ (Fin d)) in Ω, ↑↑f x * ↑↑((↑w).ofLp 0) x) (v : EuclideanSpace ℝ (Fin d) → ℝ) (hvc : ContDiff ℝ (↑⊤) v) (hvcs : HasCompactSupport v) (hvV : tsupport v ⊆ V) :
((∑ i : Fin d, ∑ j : Fin d, ∫ (x : EuclideanSpace ℝ (Fin d)) in V, Op.a x i j * ↑↑(restrictL2 ((extendL2 hΩm) ((↑u).ofLp i.succ))) x * Sobolev.partialD j v x) + ∑ i : Fin d, ∫ (x : EuclideanSpace ℝ (Fin d)) in V, Op.b x i * ↑↑(restrictL2 ((extendL2 hΩm) ((↑u).ofLp i.succ))) x * v x) + ∫ (x : EuclideanSpace ℝ (Fin d)) in V, Op.c x * ↑↑(restrictL2 ((extendL2 hΩm) ((↑u).ofLp 0))) x * v x = ∫ (x : EuclideanSpace ℝ (Fin d)) in V, ↑↑(restrictL2 ((extendL2 hΩm) f)) x * v x

Localised weak formulation (Evans, Partial Differential Equations (2nd ed.), §6.3.1, Theorem 2). A weak solution u ∈ H₀¹(Ω) of L u = f, stated through the bilinear pairing Op.fullBilin over all of Ω, satisfies the plain-integral identity ∑_{i,j} ∫_V a_{ij}(∂ᵢu) ∂ⱼv + ∑_i ∫_V b_i (∂ᵢu) v + ∫_V c u v = ∫_V f v on any measurable V ⊆ Ω, for every test function v with tsupport v ⊆ V, with each coordinate read as the V-restriction of its whole-space extension by zero. This is the shape the differentiated-equation identity takes as its hypothesis hLoc.

Differentiated-equation identity from the weak formulation #

theorem EllipticPdes.Regularity.differentiated_weakForm_of_weakSolution {n : ℕ} (Op : Sobolev.FullEllipticOp (n + 1)) {Ω : Set (EuclideanSpace ℝ (Fin (n + 1)))} (hΩm : MeasurableSet Ω) (hΩo : IsOpen Ω) (hA : IsC2Coeff Op.toEllipticCoeff) (ℓ : Fin (n + 1)) (hb : ∀ (i : Fin (n + 1)), ContDiff ℝ 1 fun (x : EuclideanSpace ℝ (Fin (n + 1))) => Op.b x i) (hc : ContDiff ℝ 1 Op.c) (Mdb : Fin (n + 1) → ℝ) (hbdM : ∀ (i : Fin (n + 1)), ∀ᵐ (x : EuclideanSpace ℝ (Fin (n + 1))), |Sobolev.partialD ℓ (fun (y : EuclideanSpace ℝ (Fin (n + 1))) => Op.b y i) x| ≤ Mdb i) (Mdc : ℝ) (hcdM : ∀ᵐ (x : EuclideanSpace ℝ (Fin (n + 1))), |Sobolev.partialD ℓ Op.c x| ≤ Mdc) {V : Set (EuclideanSpace ℝ (Fin (n + 1)))} (hVc : IsCompact V) (hVΩ : V ⊆ Ω) :
∃ (C : ℝ), 0 ≤ C ∧ ∀ (u : ↥(Sobolev.H01 Ω)) (f : Sobolev.L2D Ω) (Df : ↥(MeasureTheory.Lp ℝ 2 (MeasureTheory.volume.restrict V))), HasWeakDerivOn V ℓ (restrictL2 ((extendL2 hΩm) f)) Df → (∀ (w : ↥(Sobolev.H01 Ω)), ((Op.fullBilin Ω) u) w = ∫ (x : EuclideanSpace ℝ (Fin (n + 1))) in Ω, ↑↑f x * ↑↑((↑w).ofLp 0) x) → ∃ (D2 : Fin (n + 1) → Fin (n + 1) → ↥(MeasureTheory.Lp ℝ 2 (MeasureTheory.volume.restrict V))), (∀ (k i : Fin (n + 1)), HasWeakDerivOn V k (restrictL2 ((extendL2 hΩm) ((↑u).ofLp i.succ))) (D2 k i)) ∧ (∀ (k i : Fin (n + 1)), ‖D2 k i‖ + ‖restrictL2 ((extendL2 hΩm) ((↑u).ofLp i.succ))‖ + ‖restrictL2 ((extendL2 hΩm) ((↑u).ofLp 0))‖ ≤ C * (‖f‖ + ‖(↑u).ofLp 0‖)) ∧ ∀ (φ : EuclideanSpace ℝ (Fin (n + 1)) → ℝ), ContDiff ℝ (↑⊤) φ → HasCompactSupport φ → tsupport φ ⊆ V → ∑ i : Fin (n + 1), ∑ j : Fin (n + 1), ∫ (x : EuclideanSpace ℝ (Fin (n + 1))) in V, Op.a x i j * ↑↑(D2 ℓ i) x * Sobolev.partialD j φ x = (((∫ (x : EuclideanSpace ℝ (Fin (n + 1))) in V, ↑↑Df x * φ x) - ∑ i : Fin (n + 1), ∫ (x : EuclideanSpace ℝ (Fin (n + 1))) in V, (Sobolev.partialD ℓ (fun (y : EuclideanSpace ℝ (Fin (n + 1))) => Op.b y i) x * ↑↑(restrictL2 ((extendL2 hΩm) ((↑u).ofLp i.succ))) x + Op.b x i * ↑↑(D2 ℓ i) x) * φ x) - ∫ (x : EuclideanSpace ℝ (Fin (n + 1))) in V, (Sobolev.partialD ℓ Op.c x * ↑↑(restrictL2 ((extendL2 hΩm) ((↑u).ofLp 0))) x + Op.c x * ↑↑(restrictL2 ((extendL2 hΩm) ((↑u).ofLp ℓ.succ))) x) * φ x) + ∑ i : Fin (n + 1), ∑ j : Fin (n + 1), ∫ (x : EuclideanSpace ℝ (Fin (n + 1))) in V, (Sobolev.partialD j (Sobolev.partialD ℓ fun (y : EuclideanSpace ℝ (Fin (n + 1))) => Op.a y i j) x * ↑↑(restrictL2 ((extendL2 hΩm) ((↑u).ofLp i.succ))) x + Sobolev.partialD ℓ (fun (y : EuclideanSpace ℝ (Fin (n + 1))) => Op.a y i j) x * ↑↑(D2 j i) x) * φ x

Differentiated-equation identity for a weak solution (Evans, Partial Differential Equations (2nd ed.), §6.3.1, Theorem 2). For a weak solution u ∈ H₀¹(Ω) of L u = f with C² principal coefficients and C¹ lower-order coefficients, and any compact V ⋐ Ω, the second weak derivatives of u on V exist and are bounded by the data, and for every direction ℓ in which the datum has a weak derivative Df, the pair (∂_ℓ∂ᵢu, Df) satisfies Evans' equation (34) against every test function supported in V.

Every hypothesis of differentiated_weakForm other than the weak ℓ-derivative of f is discharged here: the second derivatives come from interior_H2_estimate, the first derivatives from hasWeakDeriv_extendL2_of_mem_H01, and the localised weak identity from localWeakForm_of_fullBilin. The weak derivative of f cannot be dropped, since Evans (26) assumes f ∈ H^m(U) and his datum (36) contains D^α f.