Documentation

TauCeti.Analysis.Convex.MongeAmpere.Basic

The Aleksandrov Monge–Ampère measure of a convex function #

Let E be a finite-dimensional real inner product space with an additive Haar measure μ, and let f : E → EReal be convex (convex real epigraph, never ⊥), with D the interior of its effective domain {x | f x ≠ ⊤}. The subgradient image of a set s is ∂f(s) = ⋃ x ∈ s, ∂f(x) (TauCeti.subgradientImage), the union of the subdifferentials of f for the inner product. The Aleksandrov Monge–Ampère measure of f is

MA_f(s) = μ(∂f(s ∩ D))

for Borel s. It gives the weak (Aleksandrov) sense in which a merely convex function solves the Monge–Ampère equation det D²f = ν. The same formula defines a measure for every μ absolutely continuous with respect to a Haar measure, such as the weighted measure s ↦ ∫_{∂f(s ∩ D)} g for μ = volume.withDensity g.

That s ↦ μ(∂f(s ∩ D)) is a measure rests on two facts.

A finite convex function u on an open convex set Ω is covered by extending it by ⊤ off Ω; then D = Ω and the subdifferential is the set of y with u x + ⟪x' - x, y⟫ ≤ u x' for all x' ∈ Ω (TauCeti.mongeAmpereMeasure_ite_apply).

For a convex function that is twice differentiable on D, the Aleksandrov measure is the classical one, MA_f = det (D²f) dx on D (TauCeti.mongeAmpereMeasure_eq_withDensity_det). On D the subdifferential is the gradient, so MA_f(s) = μ(∇f(s ∩ D)). This is computed by the change of variables formula along the gradient (TauCeti.withDensity_image_gradient_eq_setLIntegral): on a set where the subdifferential is the single gradient ∇f x and has a positive derivative H x, the gradient is injective where H x is nondegenerate (TauCeti.apply_sub_eq_zero_of_forall_mem_subdifferential_norm_sub_le) and maps the set where H x is degenerate to a null set. The Hessian of a convex function is positive (TauCeti.isPositive_of_hasFDerivAt_gradient).

The Aleksandrov maximum principle bounds a convex function by its Monge–Ampère mass: if u is convex on a bounded open convex Ω, continuous on the closure and nonnegative on the frontier, then max (-u x₀) 0 ^ n ≤ C diam(Ω) ^ (n - 1) dist(x₀, ∂Ω) MA_u(Ω) for x₀ ∈ Ω, with C = 2 ^ (n + 1) / μ (ball 0 1) (TauCeti.ofReal_neg_pow_mul_addHaar_ball_le_mul_mongeAmpereMeasure). Every slope p of an affine function through (x₀, u x₀) lying below u on the frontier is a subgradient of u at some point of Ω (TauCeti.exists_forall_add_inner_le_of_forall_frontier). When u ≥ 0 on the frontier and u x₀ < 0, these slopes include a convex set containing a ball of radius -u x₀ / diam Ω and a point of norm -u x₀ / dist(x₀, ∂Ω), and the volume of such a set is bounded below by Convex.ofReal_dist_mul_pow_mul_addHaar_ball_le.

Main definitions #

Main statements #

References #

Almost no point is a subgradient at two distinct points. For any f : E → EReal on a finite-dimensional real inner product space, the set of y lying in the subdifferentials of f at two distinct points is null for every additive Haar measure.

Subgradient images of disjoint sets are almost disjoint. For any f : E → EReal on a finite-dimensional real inner product space, the subgradient images of disjoint sets meet in a null set for every additive Haar measure.

Subgradient images of Borel subsets of the domain are null measurable. For a convex f : E → EReal on a finite-dimensional real inner product space, the subgradient image of the part of a measurable set inside the interior of the effective domain is null measurable for every additive Haar measure. If f takes the value ⊥, then no point has a subgradient and the image is empty.

The Aleksandrov Monge–Ampère measure of f : E → EReal with respect to a measure μ that is absolutely continuous with respect to an additive Haar measure. When f is convex (convex real epigraph, never ⊥), it is the measure with MA_f(s) = μ(∂f(s ∩ D)) for measurable s, where D is the interior of the effective domain and ∂f the subgradient image for the inner product (mongeAmpereMeasure_apply). For an additive Haar μ this is the Aleksandrov measure; for μ = ν.withDensity g it is the weighted measure ∫_{∂f(s ∩ D)} g dν. If f is not convex, or μ is not absolutely continuous with respect to a Haar measure, it is 0.

Instances For
    theorem TauCeti.mongeAmpereMeasure_apply {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] {f : E → EReal} [FiniteDimensional ℝ E] [MeasurableSpace E] [BorelSpace E] (μ : MeasureTheory.Measure E) {ν : MeasureTheory.Measure E} [ν.IsAddHaarMeasure] (hf : Convex ℝ {p : E × ℝ | f p.1 ≤ ↑p.2}) (hbot : ∀ (x : E), f x ≠ ⊥) (hμ : μ.AbsolutelyContinuous ν) {s : Set E} (hs : MeasurableSet s) :
    (mongeAmpereMeasure μ f) s = μ (subgradientImage (innerₗ E) f (s ∩ interior {x : E | f x ≠ ⊤}))

    The Aleksandrov Monge–Ampère measure of a convex function evaluated on a measurable set s is the measure of the subgradient image of the part of s in the interior of the effective domain.

    The Aleksandrov Monge–Ampère measure of a function whose real epigraph is not convex is 0.

    The Aleksandrov Monge–Ampère measure of a function taking the value ⊥ is 0.

    The Aleksandrov Monge–Ampère measure with respect to a measure that is not absolutely continuous with respect to an additive Haar measure is 0.

    The Aleksandrov Monge–Ampère measure is additive in the reference measure, among measures absolutely continuous with respect to an additive Haar measure.

    @[simp]

    The Aleksandrov Monge–Ampère measure is concentrated on the interior of the effective domain.

    theorem TauCeti.mongeAmpereMeasure_lt_top {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] {f : E → EReal} [FiniteDimensional ℝ E] [MeasurableSpace E] [BorelSpace E] (μ : MeasureTheory.Measure E) {ν : MeasureTheory.Measure E} [ν.IsAddHaarMeasure] [MeasureTheory.IsFiniteMeasureOnCompacts μ] (hf : Convex ℝ {p : E × ℝ | f p.1 ≤ ↑p.2}) (hbot : ∀ (x : E), f x ≠ ⊥) (hμ : μ.AbsolutelyContinuous ν) {K : Set E} (hK : IsCompact K) (hKD : K ⊆ interior {x : E | f x ≠ ⊤}) :

    The Aleksandrov Monge–Ampère measure of a convex function, with respect to a measure that is finite on compact sets, is finite on every compact subset of the interior of the effective domain.

    The Aleksandrov Monge–Ampère measure of a convex function, with respect to a measure that is finite on compact sets, is a locally finite measure on the interior of the effective domain.

    The Aleksandrov Monge–Ampère measure with respect to a measure that is finite on compact sets is σ-finite. For a convex function it is carried by the interior of the effective domain, which is a countable union of compact sets of finite Monge–Ampère measure; otherwise it is 0.

    theorem TauCeti.mongeAmpereMeasure_ite_apply {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] [FiniteDimensional ℝ E] [MeasurableSpace E] [BorelSpace E] (μ : MeasureTheory.Measure E) {ν : MeasureTheory.Measure E} [ν.IsAddHaarMeasure] {Ω : Set E} [DecidablePred fun (x : E) => x ∈ Ω] {u : E → ℝ} (hΩ : IsOpen Ω) (hu : ConvexOn ℝ Ω u) (hμ : μ.AbsolutelyContinuous ν) {s : Set E} (hs : MeasurableSet s) :
    (mongeAmpereMeasure μ fun (x : E) => if x ∈ Ω then ↑(u x) else ⊤) s = μ (⋃ x ∈ s ∩ Ω, {y : E | ∀ x' ∈ Ω, u x + inner ℝ (x' - x) y ≤ u x'})

    The Aleksandrov Monge–Ampère measure of a finite convex function on an open set. Let u be convex on an open set Ω, extended by ⊤ off Ω. On a measurable set s, its Monge–Ampère measure is the measure of the set of y that are subgradients of u relative to Ω at some point of s ∩ Ω, i.e. satisfy u x + ⟪x' - x, y⟫ ≤ u x' for every x' ∈ Ω.

    theorem TauCeti.mongeAmpereMeasure_ite_eq_restrict {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] [FiniteDimensional ℝ E] [MeasurableSpace E] [BorelSpace E] (μ : MeasureTheory.Measure E) {ν : MeasureTheory.Measure E} [ν.IsAddHaarMeasure] {Ω : Set E} [DecidablePred fun (x : E) => x ∈ Ω] {u : E → ℝ} {U : Set E} [DecidablePred fun (x : E) => x ∈ U] (hΩ : IsOpen Ω) (hU : IsOpen U) (hUc : Convex ℝ U) (hUΩ : U ⊆ Ω) (hu : ConvexOn ℝ Ω u) (hμ : μ.AbsolutelyContinuous ν) :
    (mongeAmpereMeasure μ fun (x : E) => if x ∈ U then ↑(u x) else ⊤) = (mongeAmpereMeasure μ fun (x : E) => if x ∈ Ω then ↑(u x) else ⊤).restrict U

    The Aleksandrov Monge–Ampère measure is local. Let u be convex on an open set Ω, and let U ⊆ Ω be open and convex. The Monge–Ampère measure of u restricted to U (extended by ⊤ off U) is the restriction to U of the Monge–Ampère measure of u on Ω: a subgradient of u relative to the neighbourhood U of a point is a subgradient relative to Ω.

    theorem TauCeti.withDensity_restrict_le_mongeAmpereMeasure {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] [FiniteDimensional ℝ E] [MeasurableSpace E] [BorelSpace E] (μ : MeasureTheory.Measure E) {ν : MeasureTheory.Measure E} [ν.IsAddHaarMeasure] {Ω : Set E} [DecidablePred fun (x : E) => x ∈ Ω] {u : E → ℝ} {U : Set E} [DecidablePred fun (x : E) => x ∈ U] {f g : E → ℝ} (hΩ : IsOpen Ω) (hU : IsOpen U) (hUc : Convex ℝ U) (hUΩ : U ⊆ Ω) (hu : ConvexOn ℝ Ω u) (hμ : μ.AbsolutelyContinuous ν) (hg : ∀ x ∈ U, g x ≤ f x) (hMA : ((μ.restrict Ω).withDensity fun (x : E) => ENNReal.ofReal (f x)) ≤ mongeAmpereMeasure μ fun (x : E) => if x ∈ Ω then ↑(u x) else ⊤) :
    ((μ.restrict U).withDensity fun (x : E) => ENNReal.ofReal (g x)) ≤ mongeAmpereMeasure μ fun (x : E) => if x ∈ U then ↑(u x) else ⊤

    Density below an Aleksandrov subsolution on an open convex set. Let u be convex on an open set Ω and U ⊆ Ω open and convex. If f dμ ≤ MA_u^Ω on Ω, where MA_u^Ω is the Monge–Ampère measure of u extended by ⊤ off Ω, then every density g ≤ f on U satisfies g dμ ≤ MA_u^U on U, where MA_u^U is the Monge–Ampère measure of u extended by ⊤ off U.

    theorem TauCeti.mongeAmpereMeasure_le_withDensity_restrict {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] [FiniteDimensional ℝ E] [MeasurableSpace E] [BorelSpace E] (μ : MeasureTheory.Measure E) {ν : MeasureTheory.Measure E} [ν.IsAddHaarMeasure] {Ω : Set E} [DecidablePred fun (x : E) => x ∈ Ω] {u : E → ℝ} {U : Set E} [DecidablePred fun (x : E) => x ∈ U] {f g : E → ℝ} (hΩ : IsOpen Ω) (hU : IsOpen U) (hUc : Convex ℝ U) (hUΩ : U ⊆ Ω) (hu : ConvexOn ℝ Ω u) (hμ : μ.AbsolutelyContinuous ν) (hg : ∀ x ∈ U, f x ≤ g x) (hMA : (mongeAmpereMeasure μ fun (x : E) => if x ∈ Ω then ↑(u x) else ⊤) ≤ (μ.restrict Ω).withDensity fun (x : E) => ENNReal.ofReal (f x)) :
    (mongeAmpereMeasure μ fun (x : E) => if x ∈ U then ↑(u x) else ⊤) ≤ (μ.restrict U).withDensity fun (x : E) => ENNReal.ofReal (g x)

    Density above an Aleksandrov supersolution on an open convex set. Let u be convex on an open set Ω and U ⊆ Ω open and convex. If MA_u^Ω ≤ f dμ on Ω, where MA_u^Ω is the Monge–Ampère measure of u extended by ⊤ off Ω, then every density g ≥ f on U satisfies MA_u^U ≤ g dμ on U, where MA_u^U is the Monge–Ampère measure of u extended by ⊤ off U.

    Change of variables along the gradient #

    theorem TauCeti.injOn_gradient_inter_setOf_det_ne_zero {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] {f : E → EReal} [FiniteDimensional ℝ E] (hf : Convex ℝ {p : E × ℝ | f p.1 ≤ ↑p.2}) {H : E → E →L[ℝ] E} {t : Set E} (hP : ∀ x ∈ t, subdifferential (innerₗ E) f x = {gradient (fun (x' : E) => (f x').toReal) x} ∧ ∀ ε > 0, ∀ᶠ (x' : E) in nhds x, ∀ y' ∈ subdifferential (innerₗ E) f x', ‖y' - gradient (fun (x' : E) => (f x').toReal) x - (H x) (x' - x)‖ ≤ ε * ‖x' - x‖) :
    Set.InjOn (gradient fun (x' : E) => (f x').toReal) (t ∩ {x : E | (H x).det ≠ 0})

    The gradient of a convex f : E → EReal is injective on the points of a set t at which the subdifferential of f is the single gradient ∇F x, with the expansion y' = ∇F x + H x (x' - x) + o(‖x' - x‖), and det (H x) ≠ 0. If two such points x₁ ≠ x₂ had the same gradient, it would be a subgradient along the segment between them, and the expansion at x₁ would vanish in the direction x₂ - x₁ (apply_sub_eq_zero_of_forall_mem_subdifferential_norm_sub_le).

    theorem TauCeti.withDensity_image_gradient_eq_setLIntegral {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] {f : E → EReal} [FiniteDimensional ℝ E] [MeasurableSpace E] [BorelSpace E] (μ : MeasureTheory.Measure E) [μ.IsAddHaarMeasure] (hf : Convex ℝ {p : E × ℝ | f p.1 ≤ ↑p.2}) {H : E → E →L[ℝ] E} {t : Set E} (htZ : MeasurableSet (t ∩ {x : E | (H x).det ≠ 0})) (hP : ∀ x ∈ t, subdifferential (innerₗ E) f x = {gradient (fun (x' : E) => (f x').toReal) x} ∧ (H x).IsPositive ∧ ∀ ε > 0, ∀ᶠ (x' : E) in nhds x, ∀ y' ∈ subdifferential (innerₗ E) f x', ‖y' - gradient (fun (x' : E) => (f x').toReal) x - (H x) (x' - x)‖ ≤ ε * ‖x' - x‖) (g : E → ENNReal) :
    (μ.withDensity g) ((gradient fun (x' : E) => (f x').toReal) '' t) = ∫⁻ (x : E) in t, ENNReal.ofReal (H x).det * g (gradient (fun (x' : E) => (f x').toReal) x) ∂μ

    Change of variables along the gradient of a convex function. Let f : E → EReal be convex, with real representative F x = (f x).toReal, and let t be a set at each point of which the subdifferential of f is the single gradient ∇F x, with the expansion y' = ∇F x + H x (x' - x) + o(‖x' - x‖) for a positive H x, and suppose the part of t where det (H x) ≠ 0 is measurable. Then ∫_{∇F(t)} g dμ = ∫_t det (H x) g (∇F x) dμ for every additive Haar measure μ. The gradient is injective on the nondegenerate part of t, by apply_sub_eq_zero_of_forall_mem_subdifferential_norm_sub_le, and maps the degenerate part to a null set.

    theorem TauCeti.mongeAmpereMeasure_eq_withDensity_det {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] {f : E → EReal} [FiniteDimensional ℝ E] [MeasurableSpace E] [BorelSpace E] (μ : MeasureTheory.Measure E) [μ.IsAddHaarMeasure] (hf : Convex ℝ {p : E × ℝ | f p.1 ≤ ↑p.2}) (hbot : ∀ (x : E), f x ≠ ⊥) (hd : ∀ x ∈ interior {x : E | f x ≠ ⊤}, DifferentiableAt ℝ (fun (x' : E) => (f x').toReal) x) (hd₂ : ∀ x ∈ interior {x : E | f x ≠ ⊤}, DifferentiableAt ℝ (gradient fun (x' : E) => (f x').toReal) x) :
    mongeAmpereMeasure μ f = (μ.restrict (interior {x : E | f x ≠ ⊤})).withDensity fun (x : E) => ENNReal.ofReal (fderiv ℝ (gradient fun (x' : E) => (f x').toReal) x).det

    The Monge–Ampère measure of a twice differentiable convex function is det (D²f) dx. Let f : E → EReal be convex and never ⊥, with real representative F x = (f x).toReal, and suppose that on the interior D of the effective domain both F and its gradient ∇ F are differentiable. Then for every additive Haar measure μ, the Aleksandrov Monge–Ampère measure of f is μ restricted to D with density det (D(∇ F)), the determinant of the Hessian.

    theorem TauCeti.mongeAmpereMeasure_ite_eq_withDensity_det {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] [FiniteDimensional ℝ E] [MeasurableSpace E] [BorelSpace E] (μ : MeasureTheory.Measure E) [μ.IsAddHaarMeasure] {Ω : Set E} [DecidablePred fun (x : E) => x ∈ Ω] {u : E → ℝ} (hΩ : IsOpen Ω) (hu : ConvexOn ℝ Ω u) (hd : DifferentiableOn ℝ u Ω) (hd₂ : DifferentiableOn ℝ (gradient u) Ω) :
    (mongeAmpereMeasure μ fun (x : E) => if x ∈ Ω then ↑(u x) else ⊤) = (μ.restrict Ω).withDensity fun (x : E) => ENNReal.ofReal (fderiv ℝ (gradient u) x).det

    The Monge–Ampère measure of a twice differentiable convex function on an open set. If u is convex on an open set Ω, and u and its gradient are differentiable on Ω, then the Aleksandrov Monge–Ampère measure of u (extended by ⊤ off Ω) is det (D²u) dx on Ω: it is μ restricted to Ω with density det (D(∇ u)), for every additive Haar measure μ. This applies in particular to every convex u of class C² on Ω.

    The Aleksandrov maximum principle #

    theorem TauCeti.ofReal_neg_pow_mul_addHaar_ball_le_mul_mongeAmpereMeasure {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] [FiniteDimensional ℝ E] [MeasurableSpace E] [BorelSpace E] {Ω : Set E} {u : E → ℝ} [Nontrivial E] (μ : MeasureTheory.Measure E) [μ.IsAddHaarMeasure] [DecidablePred fun (x : E) => x ∈ Ω] (hΩo : IsOpen Ω) (hΩ : Bornology.IsBounded Ω) (hu : ConvexOn ℝ Ω u) (hc : ContinuousOn u (closure Ω)) (hfr : ∀ x ∈ frontier Ω, 0 ≤ u x) {x₀ : E} (hx₀ : x₀ ∈ Ω) :
    ENNReal.ofReal (-u x₀) ^ Module.finrank ℝ E * μ (Metric.ball 0 1) ≤ 2 ^ (Module.finrank ℝ E + 1) * ENNReal.ofReal (Metric.diam Ω ^ (Module.finrank ℝ E - 1) * Metric.infDist x₀ Ωᶜ) * (mongeAmpereMeasure μ fun (x : E) => if x ∈ Ω then ↑(u x) else ⊤) Ω

    The Aleksandrov maximum principle. Let Ω be a bounded open convex subset of a real inner product space E of finite dimension n ≥ 1, and let u be convex on Ω, continuous on its closure, and nonnegative on its frontier. Then at every x₀ ∈ Ω,

    max (-u x₀) 0 ^ n * μ (ball 0 1) ≤ 2 ^ (n + 1) * diam Ω ^ (n - 1) * dist(x₀, ∂Ω) * MA_u(Ω),

    where dist(x₀, ∂Ω) is the distance from x₀ to the complement of Ω and MA_u is the Aleksandrov Monge–Ampère measure of u (extended by ⊤ off Ω) with respect to the additive Haar measure μ. When MA_u(Ω) is finite, this bounds the negative part max (-u x₀) 0 above on Ω and shows that max (-u x₀) 0 = O(dist(x₀, ∂Ω) ^ (1 / n)) as x₀ approaches the boundary.