Documentation

LeanPool.ViscositySolutionTheory.ViscositySolns.Analysis.SemiconvexJensen.MatrixConclusion

Ported for Lean Pool: imports relocated and code adapted to Lean/Mathlib v4.34.0.

Aleksandrov--Jensen theorem for semiconvex functions (MatrixConclusion) #

Part of the Aleksandrov--Jensen development assembling the localized Jensen contact-set theorem with Aleksandrov differentiability. Split from SemiconvexJensen.lean; see the umbrella module docstring.

theorem ViscositySolns.exists_strictifiedQuadraticObjective_twoSidedJet_sequence_norm_bounds_hessian {n : ℕ} (hJensen : JensenContactSetPositiveMeasureOnClosedBallTheorem n) (hAleksandrov : AleksandrovSecondDifferentiabilityByJetsOnClosedBallTheorem n) {g : Point n → ℝ} {B : Hessian n} (hB : Matrix.IsHermitian B) (hconv : ConvexOn ℝ Set.univ g) (hmax : IsMaxOn (semiconvexQuadraticObjective g B) Set.univ 0) :
∃ (xSeq : ℕ → Point n) (aSeq : ℕ → Point n) (XSeq : ℕ → Hessian n) (JSeq : ℕ → Jet n), (∀ (k : ℕ), ‖xSeq k‖ ≤ 1 / (↑k + 1)) ∧ (∀ (k : ℕ), ‖aSeq k‖ < 1 / (↑k + 1)) ∧ Filter.Tendsto xSeq Filter.atTop (nhds 0) ∧ Filter.Tendsto aSeq Filter.atTop (nhds 0) ∧ (∀ (k : ℕ), Matrix.IsHermitian (XSeq k)) ∧ (∀ (k : ℕ), Matrix.IsHermitian (JSeq k).hessian) ∧ (∀ (k : ℕ), HasSecondOrderJet (strictifiedQuadraticObjective g B) (xSeq k) (aSeq k) (XSeq k)) ∧ (∀ (k : ℕ), JSeq k = { gradient := aSeq k, hessian := XSeq k } + coordinateQuarticJetAt (xSeq k) + quadraticModelJetAt 0 0 B (xSeq k)) ∧ (∀ (k : ℕ), JSeq k ∈ ClosedSuperjet Set.univ g (xSeq k)) ∧ (∀ (k : ℕ), (JSeq k).hessian ≤ (coordinateQuarticJetAt (xSeq k)).hessian + B) ∧ ∀ (k : ℕ), HasSecondOrderJet g (xSeq k) (JSeq k).gradient (JSeq k).hessian

Selected data along the sequence 1 / (k + 1), retaining ordinary two-sided jets of the original function.

In quantified mathematical form, under the hypotheses of the localized Aleksandrov--Jensen selected-data theorem, there exist sequences x_k, a_k, X_k, and J_k such that, for every natural number k,

  • ‖x_k‖ ≤ 1 / (k + 1);
  • ‖a_k‖ < 1 / (k + 1);
  • X_k and (J_k).hessian are Hermitian;
  • (a_k, X_k) is an ordinary two-sided second-order jet of x ↦ g x - (1 / 2) * ⟪B x, x⟫ - (∑ i, x_i^2)^2 at x_k;
  • J_k is obtained from (a_k, X_k) by adding the coordinate-quartic jet at x_k and the recentered quadratic jet at x_k;
  • J_k belongs to the closed superjet of g at x_k;
  • (J_k).hessian ≤ B + (coordinate-quartic Hessian at x_k);
  • (J_k).gradient and (J_k).hessian form an ordinary two-sided second-order jet of g at x_k.

Moreover x_k -> 0 and a_k -> 0 as k -> ∞.

theorem ViscositySolns.exists_strictifiedQuadraticObjective_closedSuperjet_sequence_hessian_subsequence {n : ℕ} (hJensen : JensenContactSetPositiveMeasureOnClosedBallTheorem n) (hAleksandrov : AleksandrovSecondDifferentiabilityByJetsOnClosedBallTheorem n) {g : Point n → ℝ} {B : Hessian n} (hB : Matrix.IsHermitian B) (hconv : ConvexOn ℝ Set.univ g) (hmax : IsMaxOn (semiconvexQuadraticObjective g B) Set.univ 0) :
∃ (xSeq : ℕ → Point n) (aSeq : ℕ → Point n) (XSeq : ℕ → Hessian n) (JSeq : ℕ → Jet n) (Z : Hessian n) (φ : ℕ → ℕ), (∀ (k : ℕ), ‖xSeq k‖ ≤ 1 / (↑k + 1)) ∧ (∀ (k : ℕ), ‖aSeq k‖ < 1 / (↑k + 1)) ∧ Filter.Tendsto xSeq Filter.atTop (nhds 0) ∧ Filter.Tendsto aSeq Filter.atTop (nhds 0) ∧ (∀ (k : ℕ), Matrix.IsHermitian (XSeq k)) ∧ (∀ (k : ℕ), Matrix.IsHermitian (JSeq k).hessian) ∧ (∀ (k : ℕ), HasSecondOrderJet (strictifiedQuadraticObjective g B) (xSeq k) (aSeq k) (XSeq k)) ∧ (∀ (k : ℕ), JSeq k = { gradient := aSeq k, hessian := XSeq k } + coordinateQuarticJetAt (xSeq k) + quadraticModelJetAt 0 0 B (xSeq k)) ∧ (∀ (k : ℕ), JSeq k ∈ ClosedSuperjet Set.univ g (xSeq k)) ∧ (∀ (k : ℕ), (JSeq k).hessian ≤ (coordinateQuarticJetAt (xSeq k)).hessian + B) ∧ StrictMono φ ∧ Filter.Tendsto (fun (k : ℕ) => (JSeq (φ k)).hessian) Filter.atTop (nhds Z)

Selected closed-superjet data together with a convergent Hessian subsequence.

In standard mathematical terms, under the hypotheses of the localized Aleksandrov--Jensen selected-data theorem, one may choose selected points x_k, first-order vectors a_k, auxiliary Hessians X_k, and closed superjets J_k ∈ \overline J^{2,+}_{R^n} g(x_k) as in exists_strictifiedQuadraticObjective_closedSuperjet_sequence_norm_bounds_hessian. Moreover, the Hessian matrices (J_k).hessian have a convergent subsequence.

theorem ViscositySolns.exists_strictifiedQuadraticObjective_closedSuperjet_sequence_hessian_subsequence_order {n : ℕ} (hJensen : JensenContactSetPositiveMeasureOnClosedBallTheorem n) (hAleksandrov : AleksandrovSecondDifferentiabilityByJetsOnClosedBallTheorem n) {g : Point n → ℝ} {B : Hessian n} (hB : Matrix.IsHermitian B) (hconv : ConvexOn ℝ Set.univ g) (hmax : IsMaxOn (semiconvexQuadraticObjective g B) Set.univ 0) :
∃ (xSeq : ℕ → Point n) (aSeq : ℕ → Point n) (XSeq : ℕ → Hessian n) (JSeq : ℕ → Jet n) (Z : Hessian n) (φ : ℕ → ℕ), (∀ (k : ℕ), ‖xSeq k‖ ≤ 1 / (↑k + 1)) ∧ (∀ (k : ℕ), ‖aSeq k‖ < 1 / (↑k + 1)) ∧ Filter.Tendsto xSeq Filter.atTop (nhds 0) ∧ Filter.Tendsto aSeq Filter.atTop (nhds 0) ∧ (∀ (k : ℕ), Matrix.IsHermitian (XSeq k)) ∧ (∀ (k : ℕ), Matrix.IsHermitian (JSeq k).hessian) ∧ (∀ (k : ℕ), HasSecondOrderJet (strictifiedQuadraticObjective g B) (xSeq k) (aSeq k) (XSeq k)) ∧ (∀ (k : ℕ), JSeq k = { gradient := aSeq k, hessian := XSeq k } + coordinateQuarticJetAt (xSeq k) + quadraticModelJetAt 0 0 B (xSeq k)) ∧ (∀ (k : ℕ), JSeq k ∈ ClosedSuperjet Set.univ g (xSeq k)) ∧ (∀ (k : ℕ), (JSeq k).hessian ≤ (coordinateQuarticJetAt (xSeq k)).hessian + B) ∧ StrictMono φ ∧ Filter.Tendsto (fun (k : ℕ) => (JSeq (φ k)).hessian) Filter.atTop (nhds Z) ∧ Matrix.IsHermitian Z ∧ 0 ≤ Z ∧ Z ≤ B

Selected closed-superjet data with a Hessian subsequential limit satisfying the matrix-order bounds required in the Aleksandrov--Jensen theorem.

In standard mathematical terms, the selected Hessians admit a subsequential limit Z, and this limit is Hermitian, nonnegative, and bounded above by the quadratic matrix B.

theorem ViscositySolns.tendsto_strictifiedQuadraticObjective_closedSuperjet_sequence_gradient_zero {n : ℕ} {xSeq aSeq : ℕ → Point n} {XSeq : ℕ → Hessian n} {JSeq : ℕ → Jet n} {B : Hessian n} {φ : ℕ → ℕ} (hx : Filter.Tendsto xSeq Filter.atTop (nhds 0)) (ha : Filter.Tendsto aSeq Filter.atTop (nhds 0)) (hφ : StrictMono φ) (hJ : ∀ (k : ℕ), JSeq k = { gradient := aSeq k, hessian := XSeq k } + coordinateQuarticJetAt (xSeq k) + quadraticModelJetAt 0 0 B (xSeq k)) :
Filter.Tendsto (fun (k : ℕ) => (JSeq (φ k)).gradient) Filter.atTop (nhds 0)

The gradient components of the selected closed-superjet sequence converge to zero along every strictly increasing subsequence.

In standard mathematical terms, if J_k = (a_k, X_k) + coordinateQuarticJetAt x_k + quadraticModelJetAt 0 0 B x_k, x_k -> 0, and a_k -> 0, then (J_{φ(k)}).gradient -> 0 for every strictly increasing φ : Nat -> Nat.

theorem ViscositySolns.tendsto_jet_of_tendsto_gradient_hessian {n : ℕ} {ι : Type u_1} {l : Filter ι} {JSeq : ι → Jet n} {J : Jet n} (hgrad : Filter.Tendsto (fun (i : ι) => (JSeq i).gradient) l (nhds J.gradient)) (hhess : Filter.Tendsto (fun (i : ι) => (JSeq i).hessian) l (nhds J.hessian)) :
Filter.Tendsto JSeq l (nhds J)

Convergence of jets follows from convergence of their gradient and Hessian components.

theorem ViscositySolns.exists_strictifiedQuadraticObjective_closedSuperjet_sequence_subsequence_graph_tendsto {n : ℕ} (hJensen : JensenContactSetPositiveMeasureOnClosedBallTheorem n) (hAleksandrov : AleksandrovSecondDifferentiabilityByJetsOnClosedBallTheorem n) {g : Point n → ℝ} (hg : Continuous g) {B : Hessian n} (hB : Matrix.IsHermitian B) (hconv : ConvexOn ℝ Set.univ g) (hmax : IsMaxOn (semiconvexQuadraticObjective g B) Set.univ 0) :
∃ (xSeq : ℕ → Point n) (aSeq : ℕ → Point n) (XSeq : ℕ → Hessian n) (JSeq : ℕ → Jet n) (Z : Hessian n) (φ : ℕ → ℕ), (∀ (k : ℕ), ‖xSeq k‖ ≤ 1 / (↑k + 1)) ∧ (∀ (k : ℕ), ‖aSeq k‖ < 1 / (↑k + 1)) ∧ Filter.Tendsto xSeq Filter.atTop (nhds 0) ∧ Filter.Tendsto aSeq Filter.atTop (nhds 0) ∧ (∀ (k : ℕ), Matrix.IsHermitian (XSeq k)) ∧ (∀ (k : ℕ), Matrix.IsHermitian (JSeq k).hessian) ∧ (∀ (k : ℕ), HasSecondOrderJet (strictifiedQuadraticObjective g B) (xSeq k) (aSeq k) (XSeq k)) ∧ (∀ (k : ℕ), JSeq k = { gradient := aSeq k, hessian := XSeq k } + coordinateQuarticJetAt (xSeq k) + quadraticModelJetAt 0 0 B (xSeq k)) ∧ (∀ (k : ℕ), JSeq k ∈ ClosedSuperjet Set.univ g (xSeq k)) ∧ (∀ (k : ℕ), (JSeq k).hessian ≤ (coordinateQuarticJetAt (xSeq k)).hessian + B) ∧ StrictMono φ ∧ Filter.Tendsto (fun (k : ℕ) => (JSeq (φ k)).hessian) Filter.atTop (nhds Z) ∧ Matrix.IsHermitian Z ∧ 0 ≤ Z ∧ Z ≤ B ∧ Filter.Tendsto (fun (k : ℕ) => xSeq (φ k)) Filter.atTop (nhds 0) ∧ Filter.Tendsto (fun (k : ℕ) => g (xSeq (φ k))) Filter.atTop (nhds (g 0)) ∧ Filter.Tendsto (fun (k : ℕ) => JSeq (φ k)) Filter.atTop (nhds { gradient := 0, hessian := Z })

Selected closed-superjet data with convergence of the graph coordinates needed for the closed-superjet graph argument.

In standard mathematical terms, the selected subsequence may be chosen so that x_{φ(k)} -> 0, g(x_{φ(k)}) -> g(0), and J_{φ(k)} -> (0, Z), where Z is Hermitian and satisfies 0 ≤ Z ≤ B.

theorem ViscositySolns.exists_strictifiedQuadraticObjective_twoSidedJet_sequence_subsequence_graph_tendsto {n : ℕ} (hJensen : JensenContactSetPositiveMeasureOnClosedBallTheorem n) (hAleksandrov : AleksandrovSecondDifferentiabilityByJetsOnClosedBallTheorem n) {g : Point n → ℝ} (hg : Continuous g) {B : Hessian n} (hB : Matrix.IsHermitian B) (hconv : ConvexOn ℝ Set.univ g) (hmax : IsMaxOn (semiconvexQuadraticObjective g B) Set.univ 0) :
∃ (xSeq : ℕ → Point n) (aSeq : ℕ → Point n) (XSeq : ℕ → Hessian n) (JSeq : ℕ → Jet n) (Z : Hessian n) (φ : ℕ → ℕ), (∀ (k : ℕ), ‖xSeq k‖ ≤ 1 / (↑k + 1)) ∧ (∀ (k : ℕ), ‖aSeq k‖ < 1 / (↑k + 1)) ∧ Filter.Tendsto xSeq Filter.atTop (nhds 0) ∧ Filter.Tendsto aSeq Filter.atTop (nhds 0) ∧ (∀ (k : ℕ), Matrix.IsHermitian (XSeq k)) ∧ (∀ (k : ℕ), Matrix.IsHermitian (JSeq k).hessian) ∧ (∀ (k : ℕ), HasSecondOrderJet (strictifiedQuadraticObjective g B) (xSeq k) (aSeq k) (XSeq k)) ∧ (∀ (k : ℕ), JSeq k = { gradient := aSeq k, hessian := XSeq k } + coordinateQuarticJetAt (xSeq k) + quadraticModelJetAt 0 0 B (xSeq k)) ∧ (∀ (k : ℕ), JSeq k ∈ ClosedSuperjet Set.univ g (xSeq k)) ∧ (∀ (k : ℕ), (JSeq k).hessian ≤ (coordinateQuarticJetAt (xSeq k)).hessian + B) ∧ StrictMono φ ∧ Filter.Tendsto (fun (k : ℕ) => (JSeq (φ k)).hessian) Filter.atTop (nhds Z) ∧ Matrix.IsHermitian Z ∧ 0 ≤ Z ∧ Z ≤ B ∧ Filter.Tendsto (fun (k : ℕ) => xSeq (φ k)) Filter.atTop (nhds 0) ∧ Filter.Tendsto (fun (k : ℕ) => g (xSeq (φ k))) Filter.atTop (nhds (g 0)) ∧ Filter.Tendsto (fun (k : ℕ) => JSeq (φ k)) Filter.atTop (nhds { gradient := 0, hessian := Z }) ∧ ∀ (k : ℕ), HasSecondOrderJet g (xSeq k) (JSeq k).gradient (JSeq k).hessian

Selected closed-superjet data with graph convergence and ordinary two-sided jets of the original function.

In quantified mathematical form, under the hypotheses of exists_strictifiedQuadraticObjective_closedSuperjet_sequence_subsequence_graph_tendsto, the selected sequences may also be chosen so that, for every natural number k, (J_k).gradient and (J_k).hessian form an ordinary two-sided second-order jet of g at x_k.

theorem ViscositySolns.exists_semiconvex_twoSidedJet_sequence_subsequence_graph_tendsto {n : ℕ} (hJensen : JensenContactSetPositiveMeasureOnClosedBallTheorem n) (hAleksandrov : AleksandrovSecondDifferentiabilityByJetsOnClosedBallTheorem n) {lambda : ℝ} (_hlambda : 0 ≤ lambda) {f : Point n → ℝ} (hf : Continuous f) {B : Hessian n} (hB : Matrix.IsHermitian B) (hsemi : CoordinateSemiconvexOn lambda Set.univ f) (hmax : IsMaxOn (semiconvexQuadraticObjective f B) Set.univ 0) :
∃ (xSeq : ℕ → Point n) (JSeq : ℕ → Jet n) (Z : Hessian n) (φ : ℕ → ℕ), (∀ (k : ℕ), Matrix.IsHermitian (JSeq k).hessian) ∧ (∀ (k : ℕ), HasSecondOrderJet f (xSeq k) (JSeq k).gradient (JSeq k).hessian) ∧ StrictMono φ ∧ Filter.Tendsto (fun (k : ℕ) => (JSeq (φ k)).hessian) Filter.atTop (nhds Z) ∧ Matrix.IsHermitian Z ∧ -lambda • 1 ≤ Z ∧ Z ≤ B ∧ { gradient := 0, hessian := Z } ∈ ClosedSuperjet Set.univ f 0

Selected ordinary two-sided jets for a semiconvex function.

In quantified mathematical form, let f : R^n -> R be continuous and suppose that f satisfies the coordinate semiconvexity inequality on all of R^n with constant lambda ≥ 0. Let B be a symmetric matrix, and suppose that

x ↦ f x - (1 / 2) * ⟪B x, x⟫

has a maximum at 0 on all of R^n. Then there are points x_k, jets J_k, a matrix Z, and a strictly increasing map φ : Nat -> Nat such that:

  • for every k, J_k is an ordinary two-sided second-order jet of f at x_k;
  • the Hessians (J_{φ(k)}).hessian converge to Z;
  • (0, Z) belongs to the closed superjet of f at 0;
  • -lambda I ≤ Z ≤ B.

The convex-function form of the Aleksandrov--Jensen matrix theorem.

In quantified mathematical form, this says that for every continuous convex function g : R^n -> R and every Hermitian matrix B, if x ↦ g x - (1 / 2) * ⟪B x, x⟫ has a maximum at 0 on all of R^n, then there exists a Hermitian matrix Z such that (0, Z) ∈ \overline J^{2,+}_{R^n} g(0) and 0 ≤ Z ≤ B.

Equations
  • One or more equations did not get rendered due to their size.
Instances For

    The localized Jensen contact-set theorem and localized Aleksandrov second-differentiability theorem imply the convex Aleksandrov--Jensen matrix theorem.

    The localized Jensen contact-set theorem has already been proved in this file, so the convex Aleksandrov--Jensen matrix theorem follows from the localized Aleksandrov second-differentiability theorem alone.

    The Aleksandrov--Jensen semiconvex matrix theorem in the convexified form.

    In quantified mathematical form, this is the statement: for every continuous f : R^n -> R and every Hermitian matrix B, if 0 ≤ lambda, if x ↦ f x + (lambda / 2) * ∑ i, x_i^2 is convex on all of R^n, and if x ↦ f x - (1 / 2) * ⟪B x, x⟫ has a maximum at 0 on all of R^n, then there exists a Hermitian matrix X such that (0, X) ∈ \overline J^{2,+}_{R^n} f(0) and -lambda I ≤ X ≤ B.

    Equations
    Instances For

      The convex-function Aleksandrov--Jensen theorem implies the semiconvex Aleksandrov--Jensen theorem for each constant lambda.

      In quantified mathematical form, apply the convex theorem to g x = f x + (lambda / 2) * ∑ i, x_i^2 and the matrix B + lambda I. The closed-superjet conclusion for g is then transformed into the corresponding closed-superjet conclusion for f by subtracting the quadratic jet.

      The convex-function Aleksandrov--Jensen theorem implies the family of semiconvex theorem statements for all constants.

      The public Aleksandrov--Jensen theorem target implies the coordinate semiconvex formulation used in the comparison files.

      The coordinate semiconvex formulation and the public convexified Aleksandrov--Jensen formulation are equivalent.

      theorem ViscositySolns.AleksandrovJensenSemiconvexMatrixTheoremOn.compactSupConvolution_of_isCompact {n : ℕ} {lambda : ℝ} (hAJ : AleksandrovJensenSemiconvexMatrixTheoremOn lambda) (hlambda : 0 ≤ lambda) {K : Set (Point n)} {v : Point n → ℝ} (hKne : K.Nonempty) (hKcompact : IsCompact K) (hv : UpperSemicontinuousOn v K) {B : Hessian n} (hB : Matrix.IsHermitian B) (hmax : IsMaxOn (semiconvexQuadraticObjective (fun (ξ : Point n) => compactSupConvolution lambda K v ξ) B) Set.univ 0) :
      SemiconvexMatrixConclusion lambda (fun (ξ : Point n) => compactSupConvolution lambda K v ξ) B

      Apply the public Aleksandrov--Jensen theorem target to a compact sup-convolution.