Documentation

LeanPool.Feige.SteinIdentity

Exponential Stein identity #

This file isolates the integration-by-parts step behind the local exponential transfer identity used in the proof of Theorem 2.1. The main theorem is a one-sided identity on the positive half-line. Its assumptions expose exactly the two weighted integrability conditions and the boundary condition at infinity needed for improper integration by parts.

theorem Feige.ExponentialStein.integrableOn_comp_mul_exp {φ : } ( : Continuous φ) (hbound : ∀ (x : ), φ x 1) (y a : ) :

A uniformly unit-bounded continuous function, composed with an affine path and multiplied by the exponential density, is integrable on the positive half-line.

theorem Feige.ExponentialStein.tendsto_comp_mul_exp_atTop {φ : } (hbound : ∀ (x : ), φ x 1) (y a : ) :
Filter.Tendsto (fun (e : ) => φ (y + a * e) * Real.exp (-e)) Filter.atTop (nhds 0)

The same exponential damping forces the affine composite to vanish at infinity.

theorem Feige.ExponentialStein.integral_Ioi_comp_exp_sub {φ φ' : } {y a : } ( : ∀ (x : ), HasDerivAt φ (φ' x) x) (hintφ : MeasureTheory.IntegrableOn (fun (e : ) => φ (y + a * e) * Real.exp (-e)) (Set.Ioi 0) MeasureTheory.volume) (hintφ' : MeasureTheory.IntegrableOn (fun (e : ) => a * φ' (y + a * e) * Real.exp (-e)) (Set.Ioi 0) MeasureTheory.volume) (hzero : Filter.Tendsto (fun (e : ) => φ (y + a * e) * Real.exp (-e)) (nhdsWithin 0 (Set.Ioi 0)) (nhds (φ y))) (htop : Filter.Tendsto (fun (e : ) => φ (y + a * e) * Real.exp (-e)) Filter.atTop (nhds 0)) :
( (e : ) in Set.Ioi 0, φ (y + a * e) * Real.exp (-e)) - φ y = (e : ) in Set.Ioi 0, a * φ' (y + a * e) * Real.exp (-e)

One-sided exponential integration by parts.

For a differentiable function φ, this says

∫₀∞ φ(y + a e)e⁻ᵉ de - φ(y) = a ∫₀∞ φ'(y + a e)e⁻ᵉ de.

The hypotheses are deliberately stated as the precise weighted integrability and boundary requirements. Bounded φ and bounded φ' satisfy them, as do many unbounded test functions.

theorem Feige.ExponentialStein.two_sided_integral_Ioi {φ φ' : } {y a b : } ( : ∀ (x : ), HasDerivAt φ (φ' x) x) (hintPlus : MeasureTheory.IntegrableOn (fun (e : ) => φ (y + a * e) * Real.exp (-e)) (Set.Ioi 0) MeasureTheory.volume) (hintDerivPlus : MeasureTheory.IntegrableOn (fun (e : ) => a * φ' (y + a * e) * Real.exp (-e)) (Set.Ioi 0) MeasureTheory.volume) (hzeroPlus : Filter.Tendsto (fun (e : ) => φ (y + a * e) * Real.exp (-e)) (nhdsWithin 0 (Set.Ioi 0)) (nhds (φ y))) (htopPlus : Filter.Tendsto (fun (e : ) => φ (y + a * e) * Real.exp (-e)) Filter.atTop (nhds 0)) (hintMinus : MeasureTheory.IntegrableOn (fun (e : ) => φ (y - b * e) * Real.exp (-e)) (Set.Ioi 0) MeasureTheory.volume) (hintDerivMinus : MeasureTheory.IntegrableOn (fun (e : ) => -b * φ' (y - b * e) * Real.exp (-e)) (Set.Ioi 0) MeasureTheory.volume) (hzeroMinus : Filter.Tendsto (fun (e : ) => φ (y - b * e) * Real.exp (-e)) (nhdsWithin 0 (Set.Ioi 0)) (nhds (φ y))) (htopMinus : Filter.Tendsto (fun (e : ) => φ (y - b * e) * Real.exp (-e)) Filter.atTop (nhds 0)) :
( (e : ) in Set.Ioi 0, φ (y + a * e) * Real.exp (-e)) - (e : ) in Set.Ioi 0, φ (y - b * e) * Real.exp (-e) = ( (e : ) in Set.Ioi 0, a * φ' (y + a * e) * Real.exp (-e)) + b * (e : ) in Set.Ioi 0, φ' (y - b * e) * Real.exp (-e)

Two-sided exponential Stein identity at a fixed value of y.

This is the two-sided identity before averaging over the law of Y. The first four analytic hypotheses are the one-sided requirements for y + aE; the next four are those for y - bE.

theorem Feige.ExponentialStein.integral_Ioi_mul_deriv_one_kink {u v u' v' : } {k a' b' : } (hk : 0 < k) (hu_cont : Continuous u) (hv_cont : Continuous v) (hu : ∀ (x : ), 0 < xx kHasDerivAt u (u' x) x) (hv : ∀ (x : ), 0 < xHasDerivAt v (v' x) x) (hu'_fin : IntervalIntegrable u' MeasureTheory.volume 0 k) (hv'_fin : IntervalIntegrable v' MeasureTheory.volume 0 k) (huv'_int : MeasureTheory.IntegrableOn (u * v') (Set.Ioi 0) MeasureTheory.volume) (hu'v_int : MeasureTheory.IntegrableOn (u' * v) (Set.Ioi 0) MeasureTheory.volume) (h_zero : u 0 * v 0 = a') (h_top : Filter.Tendsto (u * v) Filter.atTop (nhds b')) :
(x : ) in Set.Ioi 0, u x * v' x = b' - a' - (x : ) in Set.Ioi 0, u' x * v x

Improper integration by parts with one possible nondifferentiability point k in the positive half-line. The finite interval theorem does not require differentiability at its right endpoint, and the improper tail theorem does not require it at its left endpoint, so the two pieces join using continuity alone.