Flow chart of with uniform coefficient bounds
ProvedOAI.SmoothLocal.Flow.exists_height_flow_uniform_coefficientsNotation. is the model square; a metric is a smooth positive matrix field ; is a scalar height; is its covariant Hessian; is the Gaussian curvature; ; is the Hessian quotient; is the model curvature; the patch box and the central box are the squares and ; the sheared coordinates are , with sheared height ; and are the class width and speed of (4.2); is the test metric (4.12) built from ; a bounded admissible height of bound is an admissible height with (4.1).
Let be smooth on with bound on , , , and on , with , as before. There is such that the chart is smooth on its domain, , , stays in and on , and with an explicit two-sided bound by for the bound of . Moreover the chart coefficients of the differentiated equation have explicit uniform bounds, an explicit positive lower bound, and for every ,
which is the structure (3.10).
Formalization Note. This is the upstream Lean declaration OAI.SmoothLocal.Flow.exists_height_flow_uniform_coefficients in lean/OAI/Geometry/IsometricImmersion/Flows/HeightFlowUniformBounds.lean, stated verbatim; the Lean statement is authoritative for the precise hypotheses. It is an intermediate result of OpenAI's formalization.
import Mathlib import Definitions.Def_IsometricImmersion import Definitions.Def_SmoothLocalIsometricSupport namespace OAI noncomputable section open Set Filter Function open scoped ContDiff Topology namespace SmoothLocal.Flow open SmoothLocal.Geometry SmoothLocal.ODE SmoothLocal.Weighted
theorem exists_height_flow_uniform_coefficients
{g : MetricField} {z : Coord → ℝ} {U : Set Coord} {G Z d c e0 : ℝ}
(hg : SmoothPositiveOn g U) (hU : IsOpen U) (hSU : modelSquare ⊆ U)
(hz : ContDiffOn ℝ ∞ z U) (hG : 0 ≤ G) (hZ : 0 ≤ Z) (hd : 0 < d) (hc : 0 < c) (he0 : 0 < e0)
(hgB : ∀ i j : Fin 2, CoordinateBound (fun p => g p i j) modelSquare 4 G)
(hzB : CoordinateBound z modelSquare 5 Z)
(hdet : ∀ p ∈ modelSquare, d ≤ |(g p).det|)
(hyy : ∀ p ∈ modelSquare, c ≤ |covHessian g z p 1 1|)
(hEfloor : ∀ p ∈ modelSquare, e0 ≤ heightEnergy g z p)
(hsmall : ∀ p ∈ modelSquare, |hessianQuotient g z p| ≤ (1 : ℝ) / 100)
(hD : ∀ p ∈ modelSquare, (covHessian g z p).det = gaussianCurvature g p * heightEnergy g z p) :
∃ Y : ℝ → ℝ → ℝ,
ContDiffOn ℝ ∞ (capChart Y) capChartDomain ∧
ContinuousOn (uncurry Y) (Icc (-2 : ℝ) 2 ×ˢ Icc (-2 : ℝ) 2) ∧
ContDiffOn ℝ ∞ (fun p : ℝ × ℝ => Y p.2 p.1) (pairRectangle 2 (-2) 2) ∧
(∀ s ∈ Ioo (-2 : ℝ) 2, ∀ t ∈ Ioo (-2 : ℝ) 2, 0 < deriv (fun r => Y r t) s) ∧
(∀ s ∈ Icc (-2 : ℝ) 2, Y s 0 = s) ∧
(∀ s ∈ Icc (-2 : ℝ) 2, ∀ t ∈ Icc (-2 : ℝ) 2,
HasDerivWithinAt (Y s) (-hessianQuotient g z (coordinatePoint t (Y s t))) (Icc (-2 : ℝ) 2) t) ∧
(∀ s ∈ Icc (-2 : ℝ) 2, ∀ t ∈ Icc (-2 : ℝ) 2, Y s t ∈ Icc (-3 : ℝ) 3) ∧
(∀ s ∈ Icc (-2 : ℝ) 2, ∀ t ∈ Icc (-2 : ℝ) 2, |Y s t - s| ≤ (1 : ℝ) / 50) ∧
(∀ p ∈ capChartDomain, |capFlowHeight Y p - p 1| ≤ (1 : ℝ) / 50) ∧
(∀ p ∈ capChartDomain,
0 < capChartRho Y p ∧
1 / Real.exp (2 * heightQuotientJetBound G Z d c) ≤ capChartRho Y p ∧
capChartRho Y p ≤ 1 / Real.exp (-2 * heightQuotientJetBound G Z d c)) ∧
ContDiffOn ℝ ∞ (heightChartG1 g z Y) capChartDomain ∧
CoordinateBound (heightChartG1 g z Y) capChartDomain 1 (heightG1JetBound G Z d c) ∧
(∀ p ∈ capChartDomain, heightG1Lower G Z d c e0 ≤ heightChartG1 g z Y p) ∧
(∀ ell : ℕ,
ContDiffOn ℝ ∞ (heightChartB g z Y ell) capChartDomain ∧
CoordinateBound (heightChartB g z Y ell) capChartDomain 1 (heightBJetBound G Z d c ell) ∧
ContDiffOn ℝ ∞ (heightChartRemainder g z Y ell) capChartDomain ∧
∀ p ∈ capChartDomain,
|heightChartRemainder g z Y ell p| ≤ heightRemainderBound G Z d c e0 ell ∧
heightChartC g z Y ell p = ((ell : ℝ) + 1) * coordPartial 1 (heightChartA g z Y) p +
heightChartA g z Y p * heightChartRemainder g z Y ell p) := by
sorry
end SmoothLocal.Flow
end
end OAI