Documentation

LeanPool.EllipticPDE.Existence.StrongMaximum

Hopf's lemma and the strong maximum principle #

Hopf's lemma: a C² subsolution on a ball, continuous on the closed ball, that is strictly below its value at a boundary point x₀ throughout the ball, has positive outward normal derivative at x₀. The barrier v = exp(-λ|x - y|²) - exp(-λ r²) is a subsolution on the annulus r/2 < |x - y| < r for λ large, vanishes on the outer sphere and is positive on the inner one, so u + ε v - u(x₀) is nonpositive on the boundary of the annulus for ε small and, by the weak maximum principle, on the annulus. Along the inward radius through x₀ the function u + ε v is therefore at most its value at x₀, and its one-sided derivative there, which is -∂_ν u(x₀) + ε ∂_ν(-v)(x₀), is nonpositive. The normal derivative of v is negative, which gives the strict inequality.

The strong maximum principle: a C² subsolution on a connected open set that attains its maximum at an interior point is constant. If not, the set where the function is below the maximum is open, nonempty, and has a frontier point inside the set; a small ball about a nearby point of it, of radius the distance to the level set of the maximum, lies in it and touches the level set at a point where Hopf's lemma gives a nonzero gradient, though the point is an interior maximum.

Main declarations #

References #

L. C. Evans, Partial Differential Equations (2nd ed.), §6.4.2 Lemma (Hopf's Lemma, p. 347) and Theorem 3 (p. 349); D. Gilbarg and N. S. Trudinger, Elliptic Partial Differential Equations of Second Order, §3.2 Lemma 3.4 (p. 34) and Theorem 3.5 (p. 35).

The barrier #

The squared distance to y as a sum of squares.

Equations
Instances For

    The squared distance is the squared norm of the difference.

    The derivative of the squared distance.

    noncomputable def EllipticPdes.Classical.barrierExp {d : ℕ} (lam : ℝ) (y x : EuclideanSpace ℝ (Fin d)) :

    The barrier's exponential part exp (-λ |x - y|²).

    Equations
    Instances For
      theorem EllipticPdes.Classical.barrierExp_pos {d : ℕ} (lam : ℝ) (y x : EuclideanSpace ℝ (Fin d)) :
      0 < barrierExp lam y x

      The exponential part is positive.

      theorem EllipticPdes.Classical.hasFDerivAt_barrierExp {d : ℕ} (lam : ℝ) (y x : EuclideanSpace ℝ (Fin d)) :
      HasFDerivAt (barrierExp lam y) (barrierExp lam y x • -lam • ∑ i : Fin d, (2 * (x.ofLp i - y.ofLp i)) • EuclideanSpace.proj i) x

      The derivative of the exponential part.

      theorem EllipticPdes.Classical.fderiv_barrierExp_apply {d : ℕ} (lam : ℝ) (y x ξ : EuclideanSpace ℝ (Fin d)) :
      (fderiv ℝ (barrierExp lam y) x) ξ = barrierExp lam y x * (-lam * ∑ i : Fin d, 2 * (x.ofLp i - y.ofLp i) * ξ.ofLp i)

      The value of the derivative of the exponential part on a vector.

      theorem EllipticPdes.Classical.partialD_barrierExp {d : ℕ} (lam : ℝ) (y : EuclideanSpace ℝ (Fin d)) (i : Fin d) (x : EuclideanSpace ℝ (Fin d)) :
      Sobolev.partialD i (barrierExp lam y) x = -2 * lam * (x.ofLp i - y.ofLp i) * barrierExp lam y x

      The first partials of the exponential part.

      theorem EllipticPdes.Classical.partialD_barrierExp_eq {d : ℕ} (lam : ℝ) (y : EuclideanSpace ℝ (Fin d)) (i : Fin d) :
      Sobolev.partialD i (barrierExp lam y) = fun (x : EuclideanSpace ℝ (Fin d)) => -2 * lam * ((x.ofLp i - y.ofLp i) * barrierExp lam y x)

      The first partials of the exponential part, as functions.

      theorem EllipticPdes.Classical.partialD_partialD_barrierExp {d : ℕ} (lam : ℝ) (y : EuclideanSpace ℝ (Fin d)) (i j : Fin d) (x : EuclideanSpace ℝ (Fin d)) :
      Sobolev.partialD i (Sobolev.partialD j (barrierExp lam y)) x = ((-2 * lam * if j = i then 1 else 0) + 4 * lam ^ 2 * (x.ofLp j - y.ofLp j) * (x.ofLp i - y.ofLp i)) * barrierExp lam y x

      The second partials of the exponential part.

      noncomputable def EllipticPdes.Classical.barrier {d : ℕ} (lam r : ℝ) (y x : EuclideanSpace ℝ (Fin d)) :

      The barrier exp (-λ |x - y|²) - exp (-λ r²).

      Equations
      Instances For

        The barrier is smooth.

        The partials of the barrier are those of its exponential part.

        theorem EllipticPdes.Classical.nondivOp_barrier {d : ℕ} (a : EuclideanSpace ℝ (Fin d) → Fin d → Fin d → ℝ) (b : EuclideanSpace ℝ (Fin d) → Fin d → ℝ) (c : EuclideanSpace ℝ (Fin d) → ℝ) (lam r : ℝ) (y x : EuclideanSpace ℝ (Fin d)) :
        nondivOp a b c (barrier lam r y) x = barrierExp lam y x * (2 * lam * ∑ i : Fin d, a x i i - 4 * lam ^ 2 * ∑ i : Fin d, ∑ j : Fin d, a x i j * (x.ofLp i - y.ofLp i) * (x.ofLp j - y.ofLp j) - 2 * lam * ∑ i : Fin d, b x i * (x.ofLp i - y.ofLp i)) + c x * barrier lam r y x

        Operator on the barrier.

        theorem EllipticPdes.Classical.nondivOp_barrier_nonpos {d : ℕ} (hd : 0 < d) {a : EuclideanSpace ℝ (Fin d) → Fin d → Fin d → ℝ} {b : EuclideanSpace ℝ (Fin d) → Fin d → ℝ} {c : EuclideanSpace ℝ (Fin d) → ℝ} {θ A B C : ℝ} (hθ : 0 < θ) {x y : EuclideanSpace ℝ (Fin d)} (hell : ∀ (ξ : Fin d → ℝ), θ * ∑ i : Fin d, ξ i ^ 2 ≤ ∑ i : Fin d, ∑ j : Fin d, a x i j * ξ i * ξ j) (ha : ∀ (i j : Fin d), |a x i j| ≤ A) (hb : ∀ (i : Fin d), |b x i| ≤ B) (hc0 : 0 ≤ c x) (hcC : c x ≤ C) {r : ℝ} (hr : 0 < r) (hq1 : r ^ 2 / 4 ≤ sqDist y x) (hq2 : sqDist y x ≤ r ^ 2) {lam : ℝ} (hlam : (2 * ↑d * A + B * (↑d + r ^ 2) + C) / (θ * r ^ 2) + 1 ≤ lam) :
        nondivOp a b c (barrier lam r y) x ≤ 0

        Barrier as a subsolution on the annulus for λ large: with the bounds on the coefficients and r²/4 ≤ |x - y|² ≤ r².

        Hopf's lemma #

        theorem EllipticPdes.Classical.deriv_nonpos_of_le_on_Ioo {φ : ℝ → ℝ} {φ' δ : ℝ} (hδ : 0 < δ) (hφ : HasDerivAt φ φ' 0) (hle : ∀ t ∈ Set.Ioo 0 δ, φ t ≤ φ 0) :
        φ' ≤ 0

        One-sided derivative at a right-sided maximum.

        theorem EllipticPdes.Classical.hopf_lemma_ball {d : ℕ} (hd : 0 < d) {y : EuclideanSpace ℝ (Fin d)} {r : ℝ} (hr : 0 < r) {a : EuclideanSpace ℝ (Fin d) → Fin d → Fin d → ℝ} {b : EuclideanSpace ℝ (Fin d) → Fin d → ℝ} {c : EuclideanSpace ℝ (Fin d) → ℝ} {θ A B C : ℝ} (hθ : 0 < θ) (hsymm : ∀ x ∈ Metric.ball y r, ∀ (i j : Fin d), a x i j = a x j i) (hell : ∀ x ∈ Metric.ball y r, ∀ (ξ : Fin d → ℝ), θ * ∑ i : Fin d, ξ i ^ 2 ≤ ∑ i : Fin d, ∑ j : Fin d, a x i j * ξ i * ξ j) (ha : ∀ x ∈ Metric.ball y r, ∀ (i j : Fin d), |a x i j| ≤ A) (hb : ∀ x ∈ Metric.ball y r, ∀ (i : Fin d), |b x i| ≤ B) (hc0 : ∀ x ∈ Metric.ball y r, 0 ≤ c x) (hcC : ∀ x ∈ Metric.ball y r, c x ≤ C) {u : EuclideanSpace ℝ (Fin d) → ℝ} (hu : ContDiffOn ℝ 2 u (Metric.ball y r)) (huc : ContinuousOn u (Metric.closedBall y r)) (hsub : ∀ x ∈ Metric.ball y r, nondivOp a b c u x ≤ 0) {x₀ : EuclideanSpace ℝ (Fin d)} (hx₀ : dist x₀ y = r) (hlt : ∀ x ∈ Metric.ball y r, u x < u x₀) (hcu : ∀ x ∈ Metric.ball y r, 0 ≤ c x * u x₀) (hdiff : DifferentiableAt ℝ u x₀) :
        0 < (fderiv ℝ u x₀) (x₀ - y)

        Hopf's lemma on a ball (Evans §6.4.2 Lemma, Gilbarg and Trudinger Lemma 3.4). A subsolution on a ball, continuous on the closed ball, strictly below its value at a point x₀ of the sphere throughout the ball, and differentiable at x₀, has positive derivative at x₀ in the outward radial direction x₀ - y. The zeroth-order coefficient is nonnegative and bounded, and c u(x₀) ≥ 0, which covers the clause c = 0 and the clause c ≥ 0 with u(x₀) ≥ 0.

        theorem EllipticPdes.Classical.hopf_lemma {d : ℕ} (hd : 0 < d) {U : Set (EuclideanSpace ℝ (Fin d))} {a : EuclideanSpace ℝ (Fin d) → Fin d → Fin d → ℝ} {b : EuclideanSpace ℝ (Fin d) → Fin d → ℝ} {θ A B : ℝ} (hθ : 0 < θ) (hsymm : ∀ x ∈ U, ∀ (i j : Fin d), a x i j = a x j i) (hell : ∀ x ∈ U, ∀ (ξ : Fin d → ℝ), θ * ∑ i : Fin d, ξ i ^ 2 ≤ ∑ i : Fin d, ∑ j : Fin d, a x i j * ξ i * ξ j) (ha : ∀ x ∈ U, ∀ (i j : Fin d), |a x i j| ≤ A) (hb : ∀ x ∈ U, ∀ (i : Fin d), |b x i| ≤ B) {c : EuclideanSpace ℝ (Fin d) → ℝ} {C : ℝ} (hc0 : ∀ x ∈ U, 0 ≤ c x) (hcC : ∀ x ∈ U, c x ≤ C) {u : EuclideanSpace ℝ (Fin d) → ℝ} (hu : ContDiffOn ℝ 2 u U) (huc : ContinuousOn u (closure U)) (hsub : ∀ x ∈ U, nondivOp a b c u x ≤ 0) {x₀ y : EuclideanSpace ℝ (Fin d)} {r : ℝ} (hr : 0 < r) (hball : Metric.ball y r ⊆ U) (hx₀ : dist x₀ y = r) (hlt : ∀ x ∈ U, u x < u x₀) (hcu : ∀ x ∈ U, 0 ≤ c x * u x₀) (hdiff : DifferentiableAt ℝ u x₀) :
        0 < (fderiv ℝ u x₀) (x₀ - y)

        Hopf's lemma (Evans §6.4.2 Lemma, Gilbarg and Trudinger Lemma 3.4). A subsolution on an open set, continuous on its closure, strictly below its value at a point x₀ throughout the set, and differentiable at x₀, has positive derivative at x₀ in the outward direction of any ball inside the set whose sphere passes through x₀. The zeroth-order coefficient is nonnegative and bounded with c u(x₀) ≥ 0, which covers both clauses of the sources.

        The strong maximum principle #

        theorem EllipticPdes.Classical.strong_maximum_principle {d : ℕ} (hd : 0 < d) {U : Set (EuclideanSpace ℝ (Fin d))} (hU : IsOpen U) (hUc : IsPreconnected U) {a : EuclideanSpace ℝ (Fin d) → Fin d → Fin d → ℝ} {b : EuclideanSpace ℝ (Fin d) → Fin d → ℝ} {c : EuclideanSpace ℝ (Fin d) → ℝ} {θ A B C : ℝ} (hθ : 0 < θ) (hsymm : ∀ x ∈ U, ∀ (i j : Fin d), a x i j = a x j i) (hell : ∀ x ∈ U, ∀ (ξ : Fin d → ℝ), θ * ∑ i : Fin d, ξ i ^ 2 ≤ ∑ i : Fin d, ∑ j : Fin d, a x i j * ξ i * ξ j) (ha : ∀ x ∈ U, ∀ (i j : Fin d), |a x i j| ≤ A) (hb : ∀ x ∈ U, ∀ (i : Fin d), |b x i| ≤ B) (hc0 : ∀ x ∈ U, 0 ≤ c x) (hcC : ∀ x ∈ U, c x ≤ C) {u : EuclideanSpace ℝ (Fin d) → ℝ} (hu : ContDiffOn ℝ 2 u U) (hsub : ∀ x ∈ U, nondivOp a b c u x ≤ 0) {x₀ : EuclideanSpace ℝ (Fin d)} (hx₀ : x₀ ∈ U) (hmax : ∀ x ∈ U, u x ≤ u x₀) (hcu : ∀ x ∈ U, 0 ≤ c x * u x₀) (x : EuclideanSpace ℝ (Fin d)) :
        x ∈ U → u x = u x₀

        Strong maximum principle (Evans §6.4.2 Theorem 3, Gilbarg and Trudinger Theorem 3.5). A subsolution, C² on a connected open set, that attains its maximum over the set at an interior point is constant on the set. The zeroth-order coefficient is nonnegative and bounded with c times the maximum nonnegative, which covers the clause c = 0 and the clause c ≥ 0 with a nonnegative maximum.

        theorem EllipticPdes.Classical.strong_maximum_principle_of_nonneg {d : ℕ} (hd : 0 < d) {U : Set (EuclideanSpace ℝ (Fin d))} (hU : IsOpen U) (hUc : IsPreconnected U) {a : EuclideanSpace ℝ (Fin d) → Fin d → Fin d → ℝ} {b : EuclideanSpace ℝ (Fin d) → Fin d → ℝ} {c : EuclideanSpace ℝ (Fin d) → ℝ} {θ A B C : ℝ} (hθ : 0 < θ) (hsymm : ∀ x ∈ U, ∀ (i j : Fin d), a x i j = a x j i) (hell : ∀ x ∈ U, ∀ (ξ : Fin d → ℝ), θ * ∑ i : Fin d, ξ i ^ 2 ≤ ∑ i : Fin d, ∑ j : Fin d, a x i j * ξ i * ξ j) (ha : ∀ x ∈ U, ∀ (i j : Fin d), |a x i j| ≤ A) (hb : ∀ x ∈ U, ∀ (i : Fin d), |b x i| ≤ B) (hc0 : ∀ x ∈ U, 0 ≤ c x) (hcC : ∀ x ∈ U, c x ≤ C) {u : EuclideanSpace ℝ (Fin d) → ℝ} (hu : ContDiffOn ℝ 2 u U) (hsub : ∀ x ∈ U, nondivOp a b c u x ≤ 0) {x₀ : EuclideanSpace ℝ (Fin d)} (hx₀ : x₀ ∈ U) (hmax : ∀ x ∈ U, u x ≤ u x₀) (hu0 : 0 ≤ u x₀) (x : EuclideanSpace ℝ (Fin d)) :
        x ∈ U → u x = u x₀

        Strong maximum principle with nonnegative zeroth-order coefficient (Evans §6.4.2 Theorem 3(ii)). With c ≥ 0, a subsolution that attains a nonnegative maximum at an interior point of a connected open set is constant on the set.