Gradient bound for the comparison difference, with constants independent of
ProvedOAI.SmoothLocal.Perturbation.exists_actual_uniform_comparison_gradientNotation. 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 and positive on an open with bound and on , and fix . There is with and small such that for every there are such that for and all large : for every patch perturbation , flow data and bounded admissible height of bound with the usual bounds and within of the test metric, the Taylor approximation of the sheared height has residual of size with weighted moment at most , the source satisfies , and for , ,
Formalization Note. This is the upstream Lean declaration OAI.SmoothLocal.Perturbation.exists_actual_uniform_comparison_gradient in lean/OAI/Geometry/IsometricImmersion/Comparison/ActualUniformComparisonGradient.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 MeasureTheory open scoped ContDiff Topology Matrix Matrix.Norms.Elementwise namespace SmoothLocal.Perturbation open SmoothLocal.Geometry SmoothLocal.Pulse SmoothLocal.HighEquation SmoothLocal.Flow open SmoothLocal.ODE SmoothLocal.Weighted SmoothLocal.Hyperbolic SmoothLocal.Taylor
theorem exists_actual_uniform_comparison_gradient
{gStar : MetricField} {V : Set Coord}
(hgStar : SmoothPositiveOn gStar V) (hV : IsOpen V) (hSV : modelSquare ⊆ V)
{G d kappa q0 : ℝ} (M : ℕ) (hG : 0 ≤ G) (hd : 0 < d)
(hgStarB : ∀ i j k, k ≤ 2 → ∀ p ∈ modelSquare,
‖iteratedFDeriv ℝ k (fun q => gStar q i j) p‖ ≤ G)
(hdStar : ∀ p ∈ modelSquare, d ≤ (gStar p).det)
(hkappa : 0 < kappa) (hM : 0 < M) (hq0 : |q0| ≤ 1/20) :
∃ r : ℝ, 0 < r ∧ r < 1/2 ∧ boundedClassWidth kappa M*r ≤ 1/20 ∧
heightQuotientJetBound G (M : ℝ) d (1/(M : ℝ))*
(r+107*(boundedClassWidth kappa M*r)/100) ≤ 9/(100*boundedClassWidth kappa M) ∧
∀ N : ℕ, 2 < N → ∃ Cres Csource Cgradient : ℝ,
0 ≤ Cres ∧ 0 < Csource ∧ 0 ≤ Cgradient ∧
∀ delta : ℝ, 0 < delta → delta ≤ 1/2 →
∀ᶠ tau : ℕ in atTop,
∀ (g0 : MetricField) (eta : metricPatchSet g0 kappa) (U W : Set Coord)
(z : Coord → ℝ) (Y : ℝ → ℝ → ℝ),
SmoothPositiveOn (perturbedMetric g0 eta.val) U → IsOpen U →
CapInductionHeight (perturbedMetric g0 eta.val) U (M : ℝ) (1/(M : ℝ)) (1/(M : ℝ)) z →
CapInductionFlow (perturbedMetric g0 eta.val) U G (M : ℝ) d (1/(M : ℝ)) (1/(M : ℝ)) kappa z Y W →
BoundedAdmissibleHeight (perturbedMetric g0 eta.val) M z →
(∀ i j k, k ≤ 4 → ∀ p ∈ modelSquare,
‖iteratedFDeriv ℝ k (fun q => perturbedMetric g0 eta.val q i j) p‖ ≤ G) →
(∀ p ∈ modelSquare, d ≤ |(perturbedMetric g0 eta.val p).det|) →
|hessianQuotient (perturbedMetric g0 eta.val) z 0-q0| ≤ 1/(100*boundedClassWidth kappa M) →
(∀ i j : Fin 2, ∀ k ≤ tau, ∀ p ∈ modelSquare,
‖iteratedFDeriv ℝ k (fun q => perturbedMetric g0 eta.val q i j-
testMetric gStar q0 (boundedClassWidth kappa M*r/16) N delta (tau : ℝ) q i j) p‖ ≤
metricApproximationAccuracy tau) →
let zs := heightInShearCoordinates z q0
let gs := metricInShearCoordinates gStar q0
let P := taylorApproximation gs (-delta/(tau : ℝ))
(heightCauchyValue zs (-delta/(tau : ℝ)))
(heightCauchyVelocity zs (-delta/(tau : ℝ))) N
ContDiffOn ℝ ∞ P (spatialStrip (Ioo (-(boundedClassWidth kappa M*r)) (boundedClassWidth kappa M*r))) ∧
(∀ x ∈ Ioo (-(boundedClassWidth kappa M*r)) (boundedClassWidth kappa M*r),
P ![x,-delta/(tau : ℝ)] = zs ![x,-delta/(tau : ℝ)] ∧
coordPartial 1 P ![x,-delta/(tau : ℝ)] = coordPartial 1 zs ![x,-delta/(tau : ℝ)]) ∧
(∀ x, ∃ Q : Polynomial ℝ, Q.natDegree ≤ N+1 ∧ ∀ t, Q.eval t = P ![x,t]) ∧
(∀ x : ℝ, |x| ≤ boundedClassWidth kappa M*r/2 →
∀ t ∈ Icc (-delta/(tau : ℝ)) (delta/(tau : ℝ)),
(boundedClassSpeed kappa M)^2/(8*(M : ℝ)) ≤ |covHessian gs P ![x,t] 0 0| ∧
|qResidual gs P ![x,t]| ≤ Cres/(tau : ℝ)^N) ∧
ContinuousOn (qResidual gs P)
(pulseStrip (boundedClassWidth kappa M*r/16) delta (tau : ℝ)) ∧
|pulseWeightedMoment (boundedClassWidth kappa M*r/16) delta (tau : ℝ) (qResidual gs P)| ≤
(2*Cres*(∫ x : ℝ, axisBump (boundedClassWidth kappa M*r/16) x))*delta/(tau : ℝ)^(N+1) ∧
(∀ p ∈ pulseStrip (boundedClassWidth kappa M*r/2) delta (tau : ℝ),
|actualComparisonSource gStar (perturbedMetric g0 eta.val) z q0 P p| ≤
Csource*(tau : ℝ)^2/(tau : ℝ)^N) ∧
(∀ theta ∈ Icc (-(delta/(tau : ℝ))) (delta/(tau : ℝ)), ∀ i : Fin 2,
Real.sqrt (∫ x in Icc (-(boundedClassWidth kappa M*r/16)) (boundedClassWidth kappa M*r/16),
(coordPartial i (comparisonDifference P zs) (boxPoint x theta))^2) ≤
Cgradient*delta*(tau : ℝ)/(tau : ℝ)^N) := by
sorry
end SmoothLocal.Perturbation
end
end OAI