Infinite-dimensional operator ridgelet transform

8.7. Manuscript Section 7🔗

8.7.1. Example 7.1 — Closed-form transform and its filtered network (ex:closed-form)🔗

Blueprint node: Proposition 6.2.1. Status: verified (all 11 Lean theorems verified).

Formalization note. Standing hypotheses: W self-adjoint, ⟨Wx,x⟩ ≥ 0, injective; Q^{1/2} is data S with IsPositiveSqrt S Q; M = S * W * S with HasSummableTrace; f_W = gaussianTarget W, D = fredholmDet (S*W*S), κ_W = gaussianKappa S W = ⟨S (I+M)⁻¹ S ξ, ξ⟩ with Ring.inverse. i_a: 𝒢_Q f_W = D^{-1/2} e^{-κ_W/2}; i_b: R_ρ f_W(a,c) = D^{-1/2}(ρ*φ_{κ_W(a)})(c) with gaussianSmooth (convolution with Mathlib's gaussianReal 0 v, δ_0 for v = 0), stated for every (a,c); i_c: f_W ∈ 𝒟_α (MemSpectralCore); i_d: non-cylindricity for W of infinite rank (IsCylindrical allows any finite-rank linear L, continuous or not, as in the manuscript's definition). ii_a: G = 𝒢_Q f_W regular along rays for every frequency window of ρ; ii_b: T_α f_W is represented by g_G (frame-operator identity of thm:C(iii), for the element of 𝒟_α equal a.e. to f_W); ii_c: eq:filtered-gaussian-target with P^{1/2} as data R, S_W = gaussianTargetResolvent S W, Σ_s = mixtureLayerCovariance R S_W s, det(I + 2s P^{1/2} S_W P^{1/2}) = fredholmDet ((2s) • (R S_W R)), the s-integral over Ioi 0. iii: β is a tempered distribution that is a Lipschitz non-polynomial function b (IsTemperedFunction, LipschitzWith L b, ¬IsPolynomialFun b); iii_a: R_ρ f_W = γ_G pointwise; iii_b: ∫ (1+‖a‖²+|c|²)|R_ρ f_W| dλ_α < ∞ (finite variation and second moment); iii_c: S_β[R_ρ f_W λ_α] = C_{β,ρ} g_G (integralNetworkDensity, temperedAdmissibilityConst); iii_d: the rate eq:spectral-barron in the conventions of Section 6, exactly as thm_E_iv with γ = R_ρ f_W: the Bochner expectation over sampleLaw n (densityLaw λ_α γ) of ‖densitySampledNetwork β λ_α γ θ − C g_G‖_{C(K)} (compactSupNorm) is at most 8V/√n (|β(0)| + Lip(β) R_K M₂) with V = densityWeight λ_α γ, R_K = compactRadius K, M₂² = secondMoment (densityLaw λ_α γ), n ≥ 1 (the earlier local Examples.polarSample/normalizedLaw/supNormOn/compactRadius/secondMoment and the lower-integral form were replaced by the Sampling/Defs objects).

OperatorRidgelet.Paper.ex_closed_form_i_a, theorem in Challenge/Examples.lean:

/-- **Example [ex:closed-form]** Closed-form transform and its filtered network.  For `W` bounded,
positive, injective, self-adjoint with `M = Q^{1/2} W Q^{1/2}` trace class,
`𝒢_Q f_W(ξ) = D^{-1/2} e^{-κ_W(ξ)/2}` with `D = det(I+M)`. -/
theorem ex_closed_form_i_a (hH : ¬ FiniteDimensional ℝ H) {P Q : H →L[ℝ] H}
    (hP : IsTraceClassCovariance P) (hQ : IsTraceClassCovariance Q) {N : ℝ → Measure H}
    (hN : IsCenteredGaussianLayers P N) {α : ℝ} (hα : 0 < α) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (W S : H →L[ℝ] H)
    (hW : IsSelfAdjoint W) (hW0 : ∀ x, 0 ≤ ⟪W x, x⟫) (hWi : Function.Injective W)
    (hS : IsPositiveSqrt S Q) (hM : HasSummableTrace (S * W * S)) :
    ∀ ξ : H, gaussFourier μ (gaussianTarget W) ξ =
      (((Real.sqrt (fredholmDet (S * W * S)))⁻¹ : ℝ) : ℂ) *
        Complex.exp (-((gaussianKappa S W ξ / 2 : ℝ) : ℂ)) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_closed_form_i_b, theorem in Challenge/Examples.lean:

/-- **Example [ex:closed-form]** Closed-form transform and its filtered network.  For every
band-pass `ρ`, `R_ρ f_W(a,c) = D^{-1/2} (ρ * φ_{κ_W(a)})(c)`. -/
theorem ex_closed_form_i_b (hH : ¬ FiniteDimensional ℝ H) {P Q : H →L[ℝ] H}
    (hP : IsTraceClassCovariance P) (hQ : IsTraceClassCovariance Q) {N : ℝ → Measure H}
    (hN : IsCenteredGaussianLayers P N) {α : ℝ} (hα : 0 < α) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (ρ : SchwartzMap ℝ ℝ)
    (hρ : IsBandPass ρ) (W S : H →L[ℝ] H) (hW : IsSelfAdjoint W) (hW0 : ∀ x, 0 ≤ ⟪W x, x⟫)
    (hWi : Function.Injective W) (hS : IsPositiveSqrt S Q) (hM : HasSummableTrace (S * W * S)) :
    ∀ p : H × ℝ, ridgelet μ ρ (gaussianTarget W) p =
      (((Real.sqrt (fredholmDet (S * W * S)))⁻¹ * gaussianSmooth ρ (gaussianKappa S W p.1) p.2 :
        ℝ) : ℂ) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_closed_form_i_c, theorem in Challenge/Examples.lean:

/-- **Example [ex:closed-form]** Closed-form transform and its filtered network.  In particular
`f_W ∈ 𝒟_α` for every `α > 0`. -/
theorem ex_closed_form_i_c (hH : ¬ FiniteDimensional ℝ H) {P Q : H →L[ℝ] H}
    (hP : IsTraceClassCovariance P) (hQ : IsTraceClassCovariance Q) {N : ℝ → Measure H}
    (hN : IsCenteredGaussianLayers P N) {α : ℝ} (hα : 0 < α) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (W S : H →L[ℝ] H)
    (hW : IsSelfAdjoint W) (hW0 : ∀ x, 0 ≤ ⟪W x, x⟫) (hWi : Function.Injective W)
    (hS : IsPositiveSqrt S Q) (hM : HasSummableTrace (S * W * S)) :
    MemSpectralCore μ (gaussianMixture N α) (gaussianTarget W) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_closed_form_i_d, theorem in Challenge/Examples.lean:

omit [SecondCountableTopology H] [MeasurableSpace H] [BorelSpace H] in
/-- **Example [ex:closed-form]** Closed-form transform and its filtered network.  `f_W` is not
cylindrical when `W` has infinite rank. -/
theorem ex_closed_form_i_d (W : H →L[ℝ] H) (hW : IsSelfAdjoint W) (hW0 : ∀ x, 0 ≤ ⟪W x, x⟫)
    (hWi : Function.Injective W) (hrank : HasInfiniteRank (W : H →ₗ[ℝ] H)) :
    ¬ IsCylindrical (gaussianTarget W) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_closed_form_ii_a, theorem in Challenge/Examples.lean:

/-- **Example [ex:closed-form]** Closed-form transform and its filtered network.  For every
band-pass `ρ`, the density `G = 𝒢_Q f_W` is regular along rays. -/
theorem ex_closed_form_ii_a (hH : ¬ FiniteDimensional ℝ H) {P Q : H →L[ℝ] H}
    (hP : IsTraceClassCovariance P) (hQ : IsTraceClassCovariance Q) {N : ℝ → Measure H}
    (hN : IsCenteredGaussianLayers P N) {α : ℝ} (hα : 0 < α) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (ρ : SchwartzMap ℝ ℝ)
    (hρ : IsBandPass ρ) (W S : H →L[ℝ] H) (hW : IsSelfAdjoint W) (hW0 : ∀ x, 0 ≤ ⟪W x, x⟫)
    (hWi : Function.Injective W) (hS : IsPositiveSqrt S Q) (hM : HasSummableTrace (S * W * S)) :
    ∀ I : Set ℝ, IsFrequencyWindow ρ I →
      IsRegularAlongRays (gaussianMixture N α) I (gaussFourier μ (gaussianTarget W)) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_closed_form_ii_b, theorem in Challenge/Examples.lean:

/-- **Example [ex:closed-form]** Closed-form transform and its filtered network.  The image
`T_α f_W` is represented by the bounded continuous function `g_G`, `G = 𝒢_Q f_W`:
`T_α f_W [g] = ∫ g_G(x) conj(g(x)) μ_Q(dx)` for `g ∈ 𝒟_α`. -/
theorem ex_closed_form_ii_b (hH : ¬ FiniteDimensional ℝ H) {P Q : H →L[ℝ] H}
    (hP : IsTraceClassCovariance P) (hQ : IsTraceClassCovariance Q) {N : ℝ → Measure H}
    (hN : IsCenteredGaussianLayers P N) {α : ℝ} (hα : 0 < α) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (W S : H →L[ℝ] H)
    (hW : IsSelfAdjoint W) (hW0 : ∀ x, 0 ≤ ⟪W x, x⟫) (hWi : Function.Injective W)
    (hS : IsPositiveSqrt S Q) (hM : HasSummableTrace (S * W * S)) :
    ∀ fW : spectralCore μ (gaussianMixture N α), (fW : H → ℂ) =ᵐ[μ] gaussianTarget W →
      ∀ g : spectralCore μ (gaussianMixture N α),
        frameOperator μ (gaussianMixture N α) (spectralEmbed μ (gaussianMixture N α) fW)
            (spectralEmbed μ (gaussianMixture N α) g) =
          ∫ x, spectralTarget (gaussianMixture N α) (gaussFourier μ (gaussianTarget W)) x *
            (starRingEnd ℂ) ((g : Lp ℂ 2 μ) x) ∂μ := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_closed_form_ii_c, theorem in Challenge/Examples.lean:

/-- **Example [ex:closed-form]** Closed-form transform and its filtered network.  With
`S_W = Q^{1/2}(I+M)⁻¹Q^{1/2}`, `R = P^{1/2}`, and
`Σ_s = 2s P^{1/2}(I + 2s P^{1/2} S_W P^{1/2})⁻¹ P^{1/2}`, the representing function is
`g_G(x) = D^{-1/2} ∫₀^∞ det(I + 2s P^{1/2} S_W P^{1/2})^{-1/2} exp(-½⟨Σ_s x,x⟩) s^{α/2-1} ds`
(`eq:filtered-gaussian-target`). -/
theorem ex_closed_form_ii_c (hH : ¬ FiniteDimensional ℝ H) {P Q : H →L[ℝ] H}
    (hP : IsTraceClassCovariance P) (hQ : IsTraceClassCovariance Q) {N : ℝ → Measure H}
    (hN : IsCenteredGaussianLayers P N) {α : ℝ} (hα : 0 < α) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (W S R : H →L[ℝ] H)
    (hW : IsSelfAdjoint W) (hW0 : ∀ x, 0 ≤ ⟪W x, x⟫) (hWi : Function.Injective W)
    (hS : IsPositiveSqrt S Q) (hM : HasSummableTrace (S * W * S)) (hR : IsPositiveSqrt R P) :
    ∀ x : H, spectralTarget (gaussianMixture N α) (gaussFourier μ (gaussianTarget W)) x =
      (((Real.sqrt (fredholmDet (S * W * S)))⁻¹ : ℝ) : ℂ) *
        ∫ s in Set.Ioi (0 : ℝ),
          (((Real.sqrt (fredholmDet ((2 * s) • (R * gaussianTargetResolvent S W * R))))⁻¹ *
            Real.exp (-⟪mixtureLayerCovariance R (gaussianTargetResolvent S W) s x, x⟫ / 2) *
            s ^ (α / 2 - 1) : ℝ) : ℂ) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_closed_form_iii_a, theorem in Challenge/Examples.lean:

/-- **Example [ex:closed-form]** Closed-form transform and its filtered network.  The ridgelet
coefficient of `f_W` is the coefficient `γ_G` of its density `G = 𝒢_Q f_W`:
`R_ρ f_W = γ_G`. -/
theorem ex_closed_form_iii_a (hH : ¬ FiniteDimensional ℝ H) {P Q : H →L[ℝ] H}
    (hP : IsTraceClassCovariance P) (hQ : IsTraceClassCovariance Q) {N : ℝ → Measure H}
    (hN : IsCenteredGaussianLayers P N) {α : ℝ} (hα : 0 < α) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (ρ : SchwartzMap ℝ ℝ)
    (hρ : IsBandPass ρ) (W S : H →L[ℝ] H) (hW : IsSelfAdjoint W) (hW0 : ∀ x, 0 ≤ ⟪W x, x⟫)
    (hWi : Function.Injective W) (hS : IsPositiveSqrt S Q) (hM : HasSummableTrace (S * W * S)) :
    ridgelet μ ρ (gaussianTarget W) = coefficientFormula ρ (gaussFourier μ (gaussianTarget W)) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_closed_form_iii_b, theorem in Challenge/Examples.lean:

/-- **Example [ex:closed-form]** Closed-form transform and its filtered network.  The ridgelet
coefficient `R_ρ f_W` has finite variation and second moment:
`∫ (1 + ‖a‖² + |c|²) |R_ρ f_W(a,c)| λ_α(da,dc) < ∞`. -/
theorem ex_closed_form_iii_b (hH : ¬ FiniteDimensional ℝ H) {P Q : H →L[ℝ] H}
    (hP : IsTraceClassCovariance P) (hQ : IsTraceClassCovariance Q) {N : ℝ → Measure H}
    (hN : IsCenteredGaussianLayers P N) {α : ℝ} (hα : 0 < α) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (ρ : SchwartzMap ℝ ℝ)
    (hρ : IsBandPass ρ) (W S : H →L[ℝ] H) (hW : IsSelfAdjoint W) (hW0 : ∀ x, 0 ≤ ⟪W x, x⟫)
    (hWi : Function.Injective W) (hS : IsPositiveSqrt S Q) (hM : HasSummableTrace (S * W * S)) :
    Integrable (fun p : H × ℝ => (1 + ‖p.1‖ ^ 2 + |p.2| ^ 2) * ‖ridgelet μ ρ (gaussianTarget W) p‖)
      (parameterMeasure (gaussianMixture N α)) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_closed_form_iii_c, theorem in Challenge/Examples.lean:

/-- **Example [ex:closed-form]** Closed-form transform and its filtered network.  For every
real, globally Lipschitz, non-polynomial `β` (including ReLU), the integral network
`S_β[R_ρ f_W λ_α]` equals `C^{(α)}_{β,ρ} g_G`. -/
theorem ex_closed_form_iii_c (hH : ¬ FiniteDimensional ℝ H) {P Q : H →L[ℝ] H}
    (hP : IsTraceClassCovariance P) (hQ : IsTraceClassCovariance Q) {N : ℝ → Measure H}
    (hN : IsCenteredGaussianLayers P N) {α : ℝ} (hα : 0 < α) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (ρ : SchwartzMap ℝ ℝ)
    (hρ : IsBandPass ρ) (W S : H →L[ℝ] H) (hW : IsSelfAdjoint W) (hW0 : ∀ x, 0 ≤ ⟪W x, x⟫)
    (hWi : Function.Injective W) (hS : IsPositiveSqrt S Q) (hM : HasSummableTrace (S * W * S))
    (β : TemperedDistribution ℝ ℂ) (b : ℝ → ℝ) (hβ : IsTemperedFunction β b) {L : ℝ≥0}
    (hb : LipschitzWith L b) (hbp : ¬ IsPolynomialFun b) :
    integralNetworkDensity (fun t => (b t : ℂ)) (parameterMeasure (gaussianMixture N α))
        (ridgelet μ ρ (gaussianTarget W)) =
      fun x => temperedAdmissibilityConst α β ρ *
        spectralTarget (gaussianMixture N α) (gaussFourier μ (gaussianTarget W)) x := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_closed_form_iii_d, theorem in Challenge/Examples.lean:

/-- **Example [ex:closed-form]** Closed-form transform and its filtered network.  For every
real, globally Lipschitz, non-polynomial `β`, the sampled network `eq:polar-network` of
`R_ρ f_W λ_α` (with `V = ‖R_ρ f_W‖_{L¹(λ_α)}` and samples from `p = |R_ρ f_W| λ_α / V`, as in
Theorem `thm:E`(iv)) converges to `C^{(α)}_{β,ρ} g_G` at the rate `n^{-1/2}` in `C(K)`, as in
`eq:spectral-barron`: `E‖f_n - C g_G‖_{C(K)} ≤ 8V n^{-1/2} (|β(0)| + Lip(β) R_K M₂)`. -/
theorem ex_closed_form_iii_d (hH : ¬ FiniteDimensional ℝ H) {P Q : H →L[ℝ] H}
    (hP : IsTraceClassCovariance P) (hQ : IsTraceClassCovariance Q) {N : ℝ → Measure H}
    (hN : IsCenteredGaussianLayers P N) {α : ℝ} (hα : 0 < α) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (ρ : SchwartzMap ℝ ℝ)
    (hρ : IsBandPass ρ) (W S : H →L[ℝ] H) (hW : IsSelfAdjoint W) (hW0 : ∀ x, 0 ≤ ⟪W x, x⟫)
    (hWi : Function.Injective W) (hS : IsPositiveSqrt S Q) (hM : HasSummableTrace (S * W * S))
    (β : TemperedDistribution ℝ ℂ) (b : ℝ → ℝ) (hβ : IsTemperedFunction β b) {L : ℝ≥0}
    (hb : LipschitzWith L b) (hbp : ¬ IsPolynomialFun b) (K : Set H) (hK : IsCompact K)
    (n : ℕ) (hn : 0 < n) :
    ∫ θ, compactSupNorm K (fun x =>
          densitySampledNetwork (fun t => (b t : ℂ)) (parameterMeasure (gaussianMixture N α))
              (ridgelet μ ρ (gaussianTarget W)) θ x -
            temperedAdmissibilityConst α β ρ *
              spectralTarget (gaussianMixture N α) (gaussFourier μ (gaussianTarget W)) x)
        ∂sampleLaw n (densityLaw (parameterMeasure (gaussianMixture N α))
          (ridgelet μ ρ (gaussianTarget W))) ≤
      8 * densityWeight (parameterMeasure (gaussianMixture N α))
            (ridgelet μ ρ (gaussianTarget W)) / Real.sqrt n *
        (|b 0| + (L : ℝ) * compactRadius K *
          Real.sqrt (secondMoment (densityLaw (parameterMeasure (gaussianMixture N α))
            (ridgelet μ ρ (gaussianTarget W))))) := by

Status: verified by comparator.

8.7.2. Example 7.2 — ReLU and Gaussian networks with Gaussian parameters (ex:gaussian-parameter)🔗

Blueprint node: Proposition 6.3.2. Status: verified (all 8 Lean theorems verified).

Formalization note. H = ℓ²(ℕ) with Q e_j = q_j e_j is stated on an abstract H with a Hilbert basis e : HilbertBasis ℕ ℝ H, q_j > 0 summable, and Q (e j) = q j • e j (the manuscript's setting up to the unitary identification with ℓ²); the parameter law is μ with IsCenteredGaussian Q μ; F_Q = gaussianParameterReLU μ, Φ_Q = gaussianParameterGauss μ with Φ = gaussianFun = e^{-u²/2} of Tempered/Defs (the manuscript's Gaussian activation, not the distribution function gaussianCdf; the duplicate Examples definition gaussianAct was removed). i, ii: the closed forms; iii: the ReLU integral over μ ⊗ db with φ'' = gaussianActDeriv2; iv, v: neither target is cylindrical; vi: 'a Gaussian-activation network with finite coefficient measure is also a ReLU network with finite coefficient measure, with all parameter moments finite' is stated for a Y-valued coefficient density γ with respect to a σ-finite λ, integrable with all moments finite: vi_a the integral-network identity with the hinge coefficient measure hingeCoefficientMeasure λ γ (a Mathlib VectorMeasure), vi_b its finiteness, vi_c its moments. The hinge pushforward identity vi_a retains BorelSpace H, the manuscript's standing compatibility of the measurable and Hilbert topologies; continuity alone would not make the parameter map measurable for an arbitrary measurable structure.

OperatorRidgelet.Paper.ex_gaussian_parameter_i, theorem in Challenge/Examples.lean:

/-- **Example [ex:gaussian-parameter]** ReLU and Gaussian networks with Gaussian parameters.
`F_Q(x) = ∫ ReLU(⟨a,x⟩) 𝒩(0,Q)(da) = √(⟨Qx,x⟩/2π)`. -/
theorem ex_gaussian_parameter_i {Q : H →L[ℝ] H} (e : HilbertBasis ℕ ℝ H) (q : ℕ → ℝ)
    (hq : ∀ j, 0 < q j) (hqs : Summable q) (hQe : ∀ j, Q (e j) = q j • e j) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) :
    ∀ x : H, gaussianParameterReLU μ x = Real.sqrt (⟪Q x, x⟫ / (2 * Real.pi)) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_gaussian_parameter_ii, theorem in Challenge/Examples.lean:

/-- **Example [ex:gaussian-parameter]** ReLU and Gaussian networks with Gaussian parameters.
`Φ_Q(x) = ∫ Φ(⟨a,x⟩) 𝒩(0,Q)(da) = (1 + ⟨Qx,x⟩)^{-1/2}`. -/
theorem ex_gaussian_parameter_ii {Q : H →L[ℝ] H} (e : HilbertBasis ℕ ℝ H) (q : ℕ → ℝ)
    (hq : ∀ j, 0 < q j) (hqs : Summable q) (hQe : ∀ j, Q (e j) = q j • e j) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) :
    ∀ x : H, gaussianParameterGauss μ x = (Real.sqrt (1 + ⟪Q x, x⟫))⁻¹ := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_gaussian_parameter_iii, theorem in Challenge/Examples.lean:

/-- **Example [ex:gaussian-parameter]** ReLU and Gaussian networks with Gaussian parameters.
`Φ_Q(x) = ∫∫ ReLU(⟨a,x⟩ - b) φ''(b) 𝒩(0,Q)(da) db` with `φ''(b) = (b² - 1) e^{-b²/2}`. -/
theorem ex_gaussian_parameter_iii {Q : H →L[ℝ] H} (e : HilbertBasis ℕ ℝ H) (q : ℕ → ℝ)
    (hq : ∀ j, 0 < q j) (hqs : Summable q) (hQe : ∀ j, Q (e j) = q j • e j) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) :
    ∀ x : H, gaussianParameterGauss μ x =
      ∫ p : H × ℝ, relu (⟪p.1, x⟫ - p.2) * gaussianActDeriv2 p.2 ∂(μ.prod volume) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_gaussian_parameter_iv, theorem in Challenge/Examples.lean:

/-- **Example [ex:gaussian-parameter]** ReLU and Gaussian networks with Gaussian parameters.
`F_Q` is not cylindrical. -/
theorem ex_gaussian_parameter_iv {Q : H →L[ℝ] H} (e : HilbertBasis ℕ ℝ H) (q : ℕ → ℝ)
    (hq : ∀ j, 0 < q j) (hqs : Summable q) (hQe : ∀ j, Q (e j) = q j • e j) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) :
    ¬ IsCylindrical (gaussianParameterReLU μ) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_gaussian_parameter_v, theorem in Challenge/Examples.lean:

/-- **Example [ex:gaussian-parameter]** ReLU and Gaussian networks with Gaussian parameters.
`Φ_Q` is not cylindrical. -/
theorem ex_gaussian_parameter_v {Q : H →L[ℝ] H} (e : HilbertBasis ℕ ℝ H) (q : ℕ → ℝ)
    (hq : ∀ j, 0 < q j) (hqs : Summable q) (hQe : ∀ j, Q (e j) = q j • e j) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) :
    ¬ IsCylindrical (gaussianParameterGauss μ) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_gaussian_parameter_vi_a, theorem in Challenge/Examples.lean:

omit [CompleteSpace H] [SecondCountableTopology H] in
/-- **Example [ex:gaussian-parameter]** ReLU and Gaussian networks with Gaussian parameters.  A
Gaussian-activation network with finite coefficient measure `γ λ` (with all parameter moments
finite) is the ReLU network with the coefficient measure `Γ' = (a,c,b) ↦ (a, c-b)`-pushforward
of `φ''(b) γ(a,c) λ(da,dc) db`. -/
theorem ex_gaussian_parameter_vi_a {Y : Type*} [NormedAddCommGroup Y] [InnerProductSpace ℂ Y]
    [CompleteSpace Y] (lam : Measure (H × ℝ)) [SigmaFinite lam] (γ : H × ℝ → Y)
    (hγ : Integrable γ lam)
    (hmom : ∀ k : ℕ, Integrable (fun θ : H × ℝ => (1 + ‖θ.1‖ + |θ.2|) ^ k * ‖γ θ‖) lam) :
    integralNetworkDensity (fun t => (gaussianFun t : ℂ)) lam γ =
      integralNetwork (fun t => (relu t : ℂ)) (hingeCoefficientMeasure lam γ) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_gaussian_parameter_vi_b, theorem in Challenge/Examples.lean:

omit [CompleteSpace H] [SecondCountableTopology H] [BorelSpace H] in
/-- **Example [ex:gaussian-parameter]** ReLU and Gaussian networks with Gaussian parameters.  The
ReLU coefficient measure `Γ'` of a Gaussian-activation network with finite coefficient measure
is finite. -/
theorem ex_gaussian_parameter_vi_b {Y : Type*} [NormedAddCommGroup Y] [InnerProductSpace ℂ Y]
    [CompleteSpace Y] (lam : Measure (H × ℝ)) [SigmaFinite lam] (γ : H × ℝ → Y)
    (hγ : Integrable γ lam)
    (hmom : ∀ k : ℕ, Integrable (fun θ : H × ℝ => (1 + ‖θ.1‖ + |θ.2|) ^ k * ‖γ θ‖) lam) :
    IsFiniteMeasure (hingeCoefficientMeasure lam γ).variation := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_gaussian_parameter_vi_c, theorem in Challenge/Examples.lean:

omit [CompleteSpace H] [SecondCountableTopology H] [BorelSpace H] in
/-- **Example [ex:gaussian-parameter]** ReLU and Gaussian networks with Gaussian parameters.  The
ReLU coefficient measure `Γ'` of a Gaussian-activation network with finite coefficient measure
has all parameter moments finite. -/
theorem ex_gaussian_parameter_vi_c {Y : Type*} [NormedAddCommGroup Y] [InnerProductSpace ℂ Y]
    [CompleteSpace Y] (lam : Measure (H × ℝ)) [SigmaFinite lam] (γ : H × ℝ → Y)
    (hγ : Integrable γ lam)
    (hmom : ∀ k : ℕ, Integrable (fun θ : H × ℝ => (1 + ‖θ.1‖ + |θ.2|) ^ k * ‖γ θ‖) lam) :
    ∀ k : ℕ, ∫⁻ θ : H × ℝ, ENNReal.ofReal ((1 + ‖θ.1‖ + |θ.2|) ^ k)
      ∂(hingeCoefficientMeasure lam γ).variation < ⊤ := by

Status: verified by comparator.

8.7.3. Corollary 7.3 — Discretization of the Gaussian-parameter ReLU network (cor:relu-discretization)🔗

Blueprint node: Corollary 6.3.3. Status: verified (its Lean theorem is verified).

Formalization note. Stated in the conventions of Section 6 (Sampling/Defs), as the case V = 1, c = 0, h = 1 of thm:lipschitz-barron used in the manuscript's proof: the samples are θ_j = (a_j, 0) with law ι_# 𝒩(0,Q) = μ.map (a ↦ (a, 0)), F_{Q,n} = sampledNetwork ReLU 1 1 θ = n⁻¹ Σ ReLU(⟨a_j,·⟩), ‖·‖_{C(K)} = compactSupNorm K, R_K = compactRadius K, tr Q = traceOf Q (Transform/Defs), the expectation over n ≥ 1 samples is the Bochner integral over sampleLaw n, and F_Q = gaussianParameterReLU μ (real valued, coerced to ℂ). The earlier local Examples.gaussianReLUSample/supNormOn/compactRadius and the lower-integral form were replaced.

OperatorRidgelet.Paper.cor_relu_discretization, theorem in Challenge/Examples.lean:

/-- **Corollary [cor:relu-discretization]** Discretization of the Gaussian-parameter ReLU
network.  For `a_1, …, a_n` independent with law `𝒩(0,Q)` (the samples `θ_j = (a_j, 0)` of the
law `ι_# 𝒩(0,Q)`, `ι(a) = (a, 0)`) and `F_{Q,n}(x) = n⁻¹ ∑_j ReLU(⟨a_j,x⟩)` (the sampled network
`eq:polar-network` with `V = 1` and phase `1`), every compact `K ⊆ H` satisfies
`E‖F_{Q,n} - F_Q‖_{C(K)} ≤ 8 R_K √(tr Q) / √n`. -/
theorem cor_relu_discretization {Q : H →L[ℝ] H} (hQ : IsTraceClassCovariance Q) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (K : Set H) (hK : IsCompact K)
    (n : ℕ) (hn : 0 < n) :
    ∫ θ, compactSupNorm K (fun x =>
          sampledNetwork (fun t => (relu t : ℂ)) 1 (fun _ => (1 : ℂ)) θ x -
            (gaussianParameterReLU μ x : ℂ))
        ∂sampleLaw n (μ.map fun a => (a, (0 : ℝ))) ≤
      8 * compactRadius K * Real.sqrt (traceOf Q) / Real.sqrt n := by

Status: verified by comparator.

8.7.4. Example 7.4 — Neural-operator layer as an integral network (ex:operator-layer)🔗

Blueprint node: Proposition 6.4.2. Status: verified (all 27 Lean theorems verified).

Formalization note. Setting: (Ω,m) a finite measure space, a : Ω → H, b : Ω → Y with IsLayerData m a b (Borel, ‖a_y‖ bounded, ∫‖b_y‖ < ∞), β continuous of polynomial growth (HasPolynomialGrowth), ℱ = operatorLayer m a b β, F_φ = layerObservable = ⟨ℱ(·),φ⟩_Y = inner ℂ φ (ℱ x), w_φ = layerWeight, ‖A‖_∞ = layerSupNorm a, A = layerA m a : H →ₗ (Ω →ₘ[m] ℝ), Γ = layerMeasure m a b (VectorMeasure map of withDensityᵥ), S_y = layerCovariance Q a y. i_a: ℱ = integralNetwork β Γ; i_b: ‖Γ‖_TV ≤ ∫⁻‖b_y‖; i_c: ∫(‖a‖²+c²)d|Γ| ≤ ‖A‖²_∞ ‖Γ‖_TV (second parameter moment ≤ ‖A‖²_∞); i_d: the L²(ζ;Y) rate of cor:vector-rates for the polar sampled network polarSampledNetwork β Γ of Γ = layerMeasure m a b with samples from polarLaw Γ (Sampling/Defs, exactly the objects of cor:vector-rates), as a Bochner expectation over sampleLaw n, with the manuscript's constant 2(∫‖b_y‖ dm)²/n(|β(0)|²+Lip(β)²(1+∫‖x‖²dζ)‖A‖²_∞) (which dominates the constant of cor:vector-rates by i_b and i_c); i_e: the C(K) bound of thm:lipschitz-barron for the polar sampled network of the scalar measure Γ_φ = layerMeasure m a w_φ = ι_#(w_φ m) with samples from polarLaw Γ_φ, with the constant 8‖w_φ‖_{L¹(m)}/√n(|β(0)|+Lip(β)R_K‖A‖_∞), n ≥ 1 (the earlier local Examples.layerSampleVec/layerSampleScalar, which sampled y ∈ Ω from normalizedLaw, and the lower-integral form were replaced by the Sampling/Defs objects the manuscript's proof invokes). ii (β = Φ = gaussianFun): ii_a F_φ ∈ 𝒟_α; ii_b, ii_c the two formulas of eq:operator-layer-transform (with gaussianSmooth); ii_d S_y ≥ (1+‖Q‖‖A‖²_∞)⁻¹Q; ii_e regularity along rays; ii_f the frame-operator representation of thm:C(iii); ii_g R_ρ F_φ = γ_G; ii_h finite variation and moments; ii_i the synthesis identity with any Lipschitz non-polynomial β'; ii_j the rate of eq:spectral-barron in the conventions of Section 6, as thm_E_iv with γ = R_ρ F_φ (densitySampledNetwork, densityLaw, densityWeight, secondMoment, compactSupNorm, compactRadius); ii_k–ii_p the same for ℱ as a Y-valued target (membership in 𝒟_α(Y), 𝒢_Q ℱ, regularity, R_ρ ℱ = γ, moments, synthesis). iii_a: the ReLU double integral over m ⊗ db; iii_b: ℱ = integralNetwork ReLU of layerHingeMeasure; iii_c, iii_d: finiteness and moments of that measure. iv: A of infinite rank (HasInfiniteRank (layerA m a)), w_φ > 0 a.e. (real positive: 0 < re, im = 0) ⇒ F_φ not cylindrical. The pushforward identities i_a/iii_b and sampling rates i_d/i_e retain BorelSpace H. The integrated sampling rate i_d also retains SecondCountableTopology H and SecondCountableTopology Y for joint measurability, as provided by the standing separability assumption. i_f: the uniform bound E||F_N - F||_{C(K;Y)} <= B_1(4|beta(0)| + 8 Lip(beta) R_K ||A||_infty)/sqrt N of eq:operator-layer-uniform for the whole output function, new in the 2026-09-13 revision, for the polar sampled network of Gamma = layerMeasure m a b itself with B_1 = int ||b_y|| dm; the manuscript renames the output kernel b_y -> v_y, and i_e remains the scalar observable with ||w_phi||_{L^1(m)} in place of B_1.

OperatorRidgelet.Paper.ex_operator_layer_i_a, theorem in Challenge/Examples.lean:

omit [CompleteSpace H] [SecondCountableTopology H] [SecondCountableTopology Y] in
/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  Network
structure: `ℱ = S_β[Γ]` with the `Y`-valued measure `Γ = ι_#(b_y m(dy))`, `ι(y) = (a_y, 0)`. -/
theorem ex_operator_layer_i_a (m : Measure Ω) [IsFiniteMeasure m] (a : Ω → H) (b : Ω → Y)
    (hL : IsLayerData m a b) (β : ℝ → ℝ) (hβc : Continuous β) (hβp : HasPolynomialGrowth β) :
    operatorLayer m a b β = integralNetwork (fun t => (β t : ℂ)) (layerMeasure m a b) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_i_b, theorem in Challenge/Examples.lean:

omit [CompleteSpace H] [SecondCountableTopology H] [BorelSpace H] [SecondCountableTopology Y] in
/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  The total
variation of `Γ` is at most `∫ ‖b_y‖ m(dy)`. -/
theorem ex_operator_layer_i_b (m : Measure Ω) [IsFiniteMeasure m] (a : Ω → H) (b : Ω → Y)
    (hL : IsLayerData m a b) :
    totalVariation (layerMeasure m a b) ≤ ∫⁻ y, ‖b y‖ₑ ∂m := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_i_c, theorem in Challenge/Examples.lean:

omit [CompleteSpace H] [SecondCountableTopology H] [BorelSpace H] [SecondCountableTopology Y] in
/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  The second
parameter moment of `Γ` is at most `‖A‖_∞²`: `∫ (‖a‖² + c²) d|Γ| ≤ ‖A‖_∞² ‖Γ‖_TV`. -/
theorem ex_operator_layer_i_c (m : Measure Ω) [IsFiniteMeasure m] (a : Ω → H) (b : Ω → Y)
    (hL : IsLayerData m a b) :
    ∫⁻ θ : H × ℝ, ENNReal.ofReal (‖θ.1‖ ^ 2 + |θ.2| ^ 2) ∂(layerMeasure m a b).variation ≤
      ENNReal.ofReal (layerSupNorm a ^ 2) * totalVariation (layerMeasure m a b) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_i_d, theorem in Challenge/Examples.lean:

omit [CompleteSpace H] in
/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  Hence
Corollary `cor:vector-rates` gives width-`n` networks approximating `ℱ` at the rate `n^{-1/2}` in
`L²(ζ;Y)`: for globally Lipschitz `β`, the polar sampled network of `Γ` (samples from
`p = |Γ|/V`) satisfies, by (b) and (c),
`E‖f_n - ℱ‖²_{L²(ζ;Y)} ≤ 2 (∫‖b_y‖ m(dy))² n⁻¹ (|β(0)|² + Lip(β)² (1 + ∫‖x‖² dζ) ‖A‖_∞²)`. -/
theorem ex_operator_layer_i_d (m : Measure Ω) [IsFiniteMeasure m] (a : Ω → H) (b : Ω → Y)
    (hL : IsLayerData m a b) (β : ℝ → ℝ) {L : ℝ≥0} (hβ : LipschitzWith L β) (ζ : Measure H)
    [IsProbabilityMeasure ζ] (hζ : Integrable (fun x : H => ‖x‖ ^ 2) ζ) (n : ℕ) (hn : 0 < n) :
    ∫ θ, (∫ x, ‖polarSampledNetwork (fun t => (β t : ℂ)) (layerMeasure m a b) θ x -
            operatorLayer m a b β x‖ ^ 2 ∂ζ)
        ∂sampleLaw n (polarLaw (layerMeasure m a b)) ≤
      2 * (∫ y, ‖b y‖ ∂m) ^ 2 / n *
        (|β 0| ^ 2 + (L : ℝ) ^ 2 * (1 + ∫ x, ‖x‖ ^ 2 ∂ζ) * layerSupNorm a ^ 2) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_i_e, theorem in Challenge/Examples.lean:

omit [CompleteSpace H] [SecondCountableTopology H] [SecondCountableTopology Y] in
/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  For each
`φ ∈ Y` and globally Lipschitz `β`, Theorem `thm:lipschitz-barron` applied to
`F_φ = S_β[Γ_φ]`, `Γ_φ = ι_#(w_φ m)`, gives for the polar sampled network `F_{φ,n}` of `Γ_φ`
`E‖F_{φ,n} - F_φ‖_{C(K)} ≤ 8 ‖w_φ‖_{L¹(m)} n^{-1/2} (|β(0)| + Lip(β) R_K ‖A‖_∞)`. -/
theorem ex_operator_layer_i_e (m : Measure Ω) [IsFiniteMeasure m] (a : Ω → H) (b : Ω → Y)
    (hL : IsLayerData m a b) (β : ℝ → ℝ) {L : ℝ≥0} (hβ : LipschitzWith L β) (φ : Y) (K : Set H)
    (hK : IsCompact K) (n : ℕ) (hn : 0 < n) :
    ∫ θ, compactSupNorm K (fun x =>
          polarSampledNetwork (fun t => (β t : ℂ)) (layerMeasure m a (layerWeight b φ)) θ x -
            layerObservable m a b β φ x)
        ∂sampleLaw n (polarLaw (layerMeasure m a (layerWeight b φ))) ≤
      8 * (∫ y, ‖layerWeight b φ y‖ ∂m) / Real.sqrt n *
        (|β 0| + (L : ℝ) * compactRadius K * layerSupNorm a) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_i_f, theorem in Challenge/Examples.lean:

/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  For globally
Lipschitz `β`, Theorem `thm:lipschitz-barron` applied to `F = S_β[Γ]`, `Γ = ι_#(v_y m)`, bounds
the error of the polar sampled network uniformly on the whole output function:
`E‖F_N - F‖_{C(K;Y)} ≤ B₁ N^{-1/2} (4|β(0)| + 8 Lip(β) R_K ‖A‖_∞)` with
`B₁ = ∫ ‖v_y‖ m(dy)` (`eq:operator-layer-uniform`). -/
theorem ex_operator_layer_i_f (m : Measure Ω) [IsFiniteMeasure m] (a : Ω → H) (b : Ω → Y)
    (hL : IsLayerData m a b) (β : ℝ → ℝ) {L : ℝ≥0} (hβ : LipschitzWith L β) (K : Set H)
    (hK : IsCompact K) (n : ℕ) (hn : 0 < n) :
    ∫ θ, compactSupNorm K (fun x =>
          polarSampledNetwork (fun t => (β t : ℂ)) (layerMeasure m a b) θ x -
            operatorLayer m a b β x)
        ∂sampleLaw n (polarLaw (layerMeasure m a b)) ≤
      (∫ y, ‖b y‖ ∂m) / Real.sqrt n *
        (4 * |β 0| + 8 * (L : ℝ) * compactRadius K * layerSupNorm a) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_ii_a, theorem in Challenge/Examples.lean:

/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  Gaussian
activation `β = Φ`: `F_φ ∈ 𝒟_α` for every `α > 0`. -/
theorem ex_operator_layer_ii_a (hH : ¬ FiniteDimensional ℝ H) {P Q : H →L[ℝ] H}
    (hP : IsTraceClassCovariance P) (hQ : IsTraceClassCovariance Q) {N : ℝ → Measure H}
    (hN : IsCenteredGaussianLayers P N) {α : ℝ} (hα : 0 < α) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (m : Measure Ω) [IsFiniteMeasure m]
    (a : Ω → H) (b : Ω → Y) (hL : IsLayerData m a b) (φ : Y) :
    MemSpectralCore μ (gaussianMixture N α) (layerObservable m a b gaussianFun φ) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_ii_b, theorem in Challenge/Examples.lean:

omit [SecondCountableTopology H] [BorelSpace H] [SecondCountableTopology Y] in
/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  Gaussian
activation: with `σ_y² = ⟨Qa_y,a_y⟩` and `S_y = Q - (1+σ_y²)⁻¹ (Qa_y) ⊗ (Qa_y)`,
`𝒢_Q F_φ(ξ) = ∫ w_φ(y) (1+σ_y²)^{-1/2} e^{-⟨S_yξ,ξ⟩/2} m(dy)` (`eq:operator-layer-transform`). -/
theorem ex_operator_layer_ii_b {Q : H →L[ℝ] H} (hQ : IsTraceClassCovariance Q) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (m : Measure Ω) [IsFiniteMeasure m]
    (a : Ω → H) (b : Ω → Y) (hL : IsLayerData m a b) (φ : Y) :
    ∀ ξ : H, gaussFourier μ (layerObservable m a b gaussianFun φ) ξ =
      ∫ y, layerWeight b φ y * (((Real.sqrt (1 + ⟪Q (a y), a y⟫))⁻¹ *
        Real.exp (-⟪layerCovariance Q a y ξ, ξ⟫ / 2) : ℝ) : ℂ) ∂m := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_ii_c, theorem in Challenge/Examples.lean:

omit [SecondCountableTopology H] [BorelSpace H] [SecondCountableTopology Y] in
/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  Gaussian
activation: for every band-pass `ρ`,
`R_ρ F_φ(a,c) = ∫ w_φ(y) (1+σ_y²)^{-1/2} (ρ * φ_{⟨S_ya,a⟩})(c) m(dy)`
(`eq:operator-layer-transform`). -/
theorem ex_operator_layer_ii_c {Q : H →L[ℝ] H} (hQ : IsTraceClassCovariance Q) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (ρ : SchwartzMap ℝ ℝ)
    (hρ : IsBandPass ρ) (m : Measure Ω) [IsFiniteMeasure m] (a : Ω → H) (b : Ω → Y)
    (hL : IsLayerData m a b) (φ : Y) :
    ∀ p : H × ℝ, ridgelet μ ρ (layerObservable m a b gaussianFun φ) p =
      ∫ y, layerWeight b φ y * (((Real.sqrt (1 + ⟪Q (a y), a y⟫))⁻¹ *
        gaussianSmooth ρ ⟪layerCovariance Q a y p.1, p.1⟫ p.2 : ℝ) : ℂ) ∂m := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_ii_d, theorem in Challenge/Examples.lean:

omit [SecondCountableTopology H] [MeasurableSpace H] [BorelSpace H]
  [SecondCountableTopology Y] in
/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  Gaussian
activation: `S_y ≥ (1 + ‖Q‖ ‖A‖_∞²)⁻¹ Q`. -/
theorem ex_operator_layer_ii_d {Q : H →L[ℝ] H} (hQ : IsTraceClassCovariance Q) (m : Measure Ω)
    [IsFiniteMeasure m] (a : Ω → H) (b : Ω → Y) (hL : IsLayerData m a b) :
    ∀ (y : Ω) (ξ : H),
      (1 + ‖Q‖ * layerSupNorm a ^ 2)⁻¹ * ⟪Q ξ, ξ⟫ ≤ ⟪layerCovariance Q a y ξ, ξ⟫ := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_ii_e, theorem in Challenge/Examples.lean:

/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  Gaussian
activation: consequently `𝒢_Q F_φ` is regular along rays, for every band-pass `ρ`. -/
theorem ex_operator_layer_ii_e (hH : ¬ FiniteDimensional ℝ H) {P Q : H →L[ℝ] H}
    (hP : IsTraceClassCovariance P) (hQ : IsTraceClassCovariance Q) {N : ℝ → Measure H}
    (hN : IsCenteredGaussianLayers P N) {α : ℝ} (hα : 0 < α) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (ρ : SchwartzMap ℝ ℝ)
    (hρ : IsBandPass ρ) (m : Measure Ω) [IsFiniteMeasure m] (a : Ω → H) (b : Ω → Y)
    (hL : IsLayerData m a b) (φ : Y) :
    ∀ I : Set ℝ, IsFrequencyWindow ρ I →
      IsRegularAlongRays (gaussianMixture N α) I
        (gaussFourier μ (layerObservable m a b gaussianFun φ)) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_ii_f, theorem in Challenge/Examples.lean:

/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  Gaussian
activation: the reconstruction formulas of Theorem `thm:C` hold for `F_φ`; in particular
`T_α F_φ` is represented by `g_G`, `G = 𝒢_Q F_φ`: `T_α F_φ [g] = ∫ g_G(x) conj(g(x)) μ_Q(dx)`
for `g ∈ 𝒟_α`. -/
theorem ex_operator_layer_ii_f (hH : ¬ FiniteDimensional ℝ H) {P Q : H →L[ℝ] H}
    (hP : IsTraceClassCovariance P) (hQ : IsTraceClassCovariance Q) {N : ℝ → Measure H}
    (hN : IsCenteredGaussianLayers P N) {α : ℝ} (hα : 0 < α) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (m : Measure Ω) [IsFiniteMeasure m]
    (a : Ω → H) (b : Ω → Y) (hL : IsLayerData m a b) (φ : Y) :
    ∀ F : spectralCore μ (gaussianMixture N α),
      (F : H → ℂ) =ᵐ[μ] layerObservable m a b gaussianFun φ →
      ∀ g : spectralCore μ (gaussianMixture N α),
        frameOperator μ (gaussianMixture N α) (spectralEmbed μ (gaussianMixture N α) F)
            (spectralEmbed μ (gaussianMixture N α) g) =
          ∫ x, spectralTarget (gaussianMixture N α)
              (gaussFourier μ (layerObservable m a b gaussianFun φ)) x *
            (starRingEnd ℂ) ((g : Lp ℂ 2 μ) x) ∂μ := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_ii_g, theorem in Challenge/Examples.lean:

omit [SecondCountableTopology H] [BorelSpace H] [SecondCountableTopology Y] in
/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  Gaussian
activation: the ridgelet coefficient of `F_φ` is the coefficient `γ_G` of `G = 𝒢_Q F_φ`. -/
theorem ex_operator_layer_ii_g {Q : H →L[ℝ] H} (hQ : IsTraceClassCovariance Q) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (ρ : SchwartzMap ℝ ℝ)
    (hρ : IsBandPass ρ) (m : Measure Ω) [IsFiniteMeasure m] (a : Ω → H) (b : Ω → Y)
    (hL : IsLayerData m a b) (φ : Y) :
    ridgelet μ ρ (layerObservable m a b gaussianFun φ) =
      coefficientFormula ρ (gaussFourier μ (layerObservable m a b gaussianFun φ)) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_ii_h, theorem in Challenge/Examples.lean:

/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  Gaussian
activation: the ridgelet coefficient `R_ρ F_φ` has finite variation and moments,
`∫ (1 + ‖a‖² + |c|²) |R_ρ F_φ(a,c)| λ_α(da,dc) < ∞`. -/
theorem ex_operator_layer_ii_h (hH : ¬ FiniteDimensional ℝ H) {P Q : H →L[ℝ] H}
    (hP : IsTraceClassCovariance P) (hQ : IsTraceClassCovariance Q) {N : ℝ → Measure H}
    (hN : IsCenteredGaussianLayers P N) {α : ℝ} (hα : 0 < α) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (ρ : SchwartzMap ℝ ℝ)
    (hρ : IsBandPass ρ) (m : Measure Ω) [IsFiniteMeasure m] (a : Ω → H) (b : Ω → Y)
    (hL : IsLayerData m a b) (φ : Y) :
    Integrable (fun p : H × ℝ =>
        (1 + ‖p.1‖ ^ 2 + |p.2| ^ 2) * ‖ridgelet μ ρ (layerObservable m a b gaussianFun φ) p‖)
      (parameterMeasure (gaussianMixture N α)) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_ii_i, theorem in Challenge/Examples.lean:

/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  Gaussian
activation: the ridgelet coefficient `R_ρ F_φ` synthesizes, with any real Lipschitz
non-polynomial `β'`, the target `C^{(α)}_{β',ρ} T_α F_φ` (represented by `g_G`). -/
theorem ex_operator_layer_ii_i (hH : ¬ FiniteDimensional ℝ H) {P Q : H →L[ℝ] H}
    (hP : IsTraceClassCovariance P) (hQ : IsTraceClassCovariance Q) {N : ℝ → Measure H}
    (hN : IsCenteredGaussianLayers P N) {α : ℝ} (hα : 0 < α) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (ρ : SchwartzMap ℝ ℝ)
    (hρ : IsBandPass ρ) (m : Measure Ω) [IsFiniteMeasure m] (a : Ω → H) (b : Ω → Y)
    (hL : IsLayerData m a b) (φ : Y) (β' : TemperedDistribution ℝ ℂ) (b' : ℝ → ℝ)
    (hβ' : IsTemperedFunction β' b') {L : ℝ≥0} (hb' : LipschitzWith L b')
    (hb'p : ¬ IsPolynomialFun b') :
    integralNetworkDensity (fun t => (b' t : ℂ)) (parameterMeasure (gaussianMixture N α))
        (ridgelet μ ρ (layerObservable m a b gaussianFun φ)) =
      fun x => temperedAdmissibilityConst α β' ρ *
        spectralTarget (gaussianMixture N α)
          (gaussFourier μ (layerObservable m a b gaussianFun φ)) x := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_ii_j, theorem in Challenge/Examples.lean:

/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  Gaussian
activation: the sampled network `eq:polar-network` of `R_ρ F_φ λ_α` with a real Lipschitz
non-polynomial `β'` (with `V = ‖R_ρ F_φ‖_{L¹(λ_α)}` and samples from `p = |R_ρ F_φ| λ_α / V`,
as in Theorem `thm:E`(iv)) converges to `C^{(α)}_{β',ρ} g_G` at the finite-width rate of
`eq:spectral-barron`. -/
theorem ex_operator_layer_ii_j (hH : ¬ FiniteDimensional ℝ H) {P Q : H →L[ℝ] H}
    (hP : IsTraceClassCovariance P) (hQ : IsTraceClassCovariance Q) {N : ℝ → Measure H}
    (hN : IsCenteredGaussianLayers P N) {α : ℝ} (hα : 0 < α) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (ρ : SchwartzMap ℝ ℝ)
    (hρ : IsBandPass ρ) (m : Measure Ω) [IsFiniteMeasure m] (a : Ω → H) (b : Ω → Y)
    (hL : IsLayerData m a b) (φ : Y) (β' : TemperedDistribution ℝ ℂ) (b' : ℝ → ℝ)
    (hβ' : IsTemperedFunction β' b') {L : ℝ≥0} (hb' : LipschitzWith L b')
    (hb'p : ¬ IsPolynomialFun b') (K : Set H) (hK : IsCompact K) (n : ℕ) (hn : 0 < n) :
    ∫ θ, compactSupNorm K (fun x =>
          densitySampledNetwork (fun t => (b' t : ℂ)) (parameterMeasure (gaussianMixture N α))
              (ridgelet μ ρ (layerObservable m a b gaussianFun φ)) θ x -
            temperedAdmissibilityConst α β' ρ *
              spectralTarget (gaussianMixture N α)
                (gaussFourier μ (layerObservable m a b gaussianFun φ)) x)
        ∂sampleLaw n (densityLaw (parameterMeasure (gaussianMixture N α))
          (ridgelet μ ρ (layerObservable m a b gaussianFun φ))) ≤
      8 * densityWeight (parameterMeasure (gaussianMixture N α))
            (ridgelet μ ρ (layerObservable m a b gaussianFun φ)) / Real.sqrt n *
        (|b' 0| + (L : ℝ) * compactRadius K *
          Real.sqrt (secondMoment (densityLaw (parameterMeasure (gaussianMixture N α))
            (ridgelet μ ρ (layerObservable m a b gaussianFun φ))))) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_ii_k, theorem in Challenge/Examples.lean:

/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  The same
holds for `ℱ` itself as a `Y`-valued target: `ℱ ∈ 𝒟_α(Y)` for every `α > 0`. -/
theorem ex_operator_layer_ii_k (hH : ¬ FiniteDimensional ℝ H) {P Q : H →L[ℝ] H}
    (hP : IsTraceClassCovariance P) (hQ : IsTraceClassCovariance Q) {N : ℝ → Measure H}
    (hN : IsCenteredGaussianLayers P N) {α : ℝ} (hα : 0 < α) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (m : Measure Ω) [IsFiniteMeasure m]
    (a : Ω → H) (b : Ω → Y) (hL : IsLayerData m a b) :
    MemSpectralCoreVec μ (gaussianMixture N α) (operatorLayer m a b gaussianFun) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_ii_l, theorem in Challenge/Examples.lean:

omit [SecondCountableTopology H] [BorelSpace H] [SecondCountableTopology Y] in
/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  The same
holds for `ℱ` itself as a `Y`-valued target:
`𝒢_Q ℱ(ξ) = ∫ (1+σ_y²)^{-1/2} e^{-⟨S_yξ,ξ⟩/2} b_y m(dy)`. -/
theorem ex_operator_layer_ii_l {Q : H →L[ℝ] H} (hQ : IsTraceClassCovariance Q) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (m : Measure Ω) [IsFiniteMeasure m]
    (a : Ω → H) (b : Ω → Y) (hL : IsLayerData m a b) :
    ∀ ξ : H, gaussFourierVec μ (operatorLayer m a b gaussianFun) ξ =
      ∫ y, (((Real.sqrt (1 + ⟪Q (a y), a y⟫))⁻¹ *
        Real.exp (-⟪layerCovariance Q a y ξ, ξ⟫ / 2) : ℝ) : ℂ) • b y ∂m := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_ii_m, theorem in Challenge/Examples.lean:

/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  The same
holds for `ℱ` itself as a `Y`-valued target: `𝒢_Q ℱ` is regular along rays for every band-pass
`ρ`. -/
theorem ex_operator_layer_ii_m (hH : ¬ FiniteDimensional ℝ H) {P Q : H →L[ℝ] H}
    (hP : IsTraceClassCovariance P) (hQ : IsTraceClassCovariance Q) {N : ℝ → Measure H}
    (hN : IsCenteredGaussianLayers P N) {α : ℝ} (hα : 0 < α) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (ρ : SchwartzMap ℝ ℝ)
    (hρ : IsBandPass ρ) (m : Measure Ω) [IsFiniteMeasure m] (a : Ω → H) (b : Ω → Y)
    (hL : IsLayerData m a b) :
    ∀ I : Set ℝ, IsFrequencyWindow ρ I →
      IsRegularAlongRays (gaussianMixture N α) I
        (gaussFourierVec μ (operatorLayer m a b gaussianFun)) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_ii_n, theorem in Challenge/Examples.lean:

omit [SecondCountableTopology H] [BorelSpace H] [SecondCountableTopology Y] in
/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  The same
holds for `ℱ` itself as a `Y`-valued target: `R_ρ ℱ = γ_{𝒢_Q ℱ}` for every band-pass `ρ`. -/
theorem ex_operator_layer_ii_n {Q : H →L[ℝ] H} (hQ : IsTraceClassCovariance Q) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (ρ : SchwartzMap ℝ ℝ)
    (hρ : IsBandPass ρ) (m : Measure Ω) [IsFiniteMeasure m] (a : Ω → H) (b : Ω → Y)
    (hL : IsLayerData m a b) :
    ridgeletVec μ ρ (operatorLayer m a b gaussianFun) =
      coefficientFormulaVec ρ (gaussFourierVec μ (operatorLayer m a b gaussianFun)) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_ii_o, theorem in Challenge/Examples.lean:

/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  The same
holds for `ℱ` itself as a `Y`-valued target: `R_ρ ℱ` has finite variation and moments. -/
theorem ex_operator_layer_ii_o (hH : ¬ FiniteDimensional ℝ H) {P Q : H →L[ℝ] H}
    (hP : IsTraceClassCovariance P) (hQ : IsTraceClassCovariance Q) {N : ℝ → Measure H}
    (hN : IsCenteredGaussianLayers P N) {α : ℝ} (hα : 0 < α) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (ρ : SchwartzMap ℝ ℝ)
    (hρ : IsBandPass ρ) (m : Measure Ω) [IsFiniteMeasure m] (a : Ω → H) (b : Ω → Y)
    (hL : IsLayerData m a b) :
    Integrable (fun p : H × ℝ =>
        (1 + ‖p.1‖ ^ 2 + |p.2| ^ 2) * ‖ridgeletVec μ ρ (operatorLayer m a b gaussianFun) p‖)
      (parameterMeasure (gaussianMixture N α)) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_ii_p, theorem in Challenge/Examples.lean:

/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  The same
holds for `ℱ` itself as a `Y`-valued target: `R_ρ ℱ` synthesizes, with any real Lipschitz
non-polynomial `β'`, the `Y`-valued target `C^{(α)}_{β',ρ} g_{𝒢_Q ℱ}`. -/
theorem ex_operator_layer_ii_p (hH : ¬ FiniteDimensional ℝ H) {P Q : H →L[ℝ] H}
    (hP : IsTraceClassCovariance P) (hQ : IsTraceClassCovariance Q) {N : ℝ → Measure H}
    (hN : IsCenteredGaussianLayers P N) {α : ℝ} (hα : 0 < α) (μ : Measure H)
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian Q μ) (ρ : SchwartzMap ℝ ℝ)
    (hρ : IsBandPass ρ) (m : Measure Ω) [IsFiniteMeasure m] (a : Ω → H) (b : Ω → Y)
    (hL : IsLayerData m a b) (β' : TemperedDistribution ℝ ℂ) (b' : ℝ → ℝ)
    (hβ' : IsTemperedFunction β' b') {L : ℝ≥0} (hb' : LipschitzWith L b')
    (hb'p : ¬ IsPolynomialFun b') :
    integralNetworkDensity (fun t => (b' t : ℂ)) (parameterMeasure (gaussianMixture N α))
        (ridgeletVec μ ρ (operatorLayer m a b gaussianFun)) =
      fun x => temperedAdmissibilityConst α β' ρ •
        spectralTarget (gaussianMixture N α)
          (gaussFourierVec μ (operatorLayer m a b gaussianFun)) x := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_iii_a, theorem in Challenge/Examples.lean:

omit [CompleteSpace H] [SecondCountableTopology H] [BorelSpace H] [SecondCountableTopology Y] in
/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  ReLU form:
by `eq:gaussian-parameter-closed-forms`, the Gaussian-activation layer is
`ℱ(x) = ∫∫ b_y φ''(b) ReLU(⟨a_y,x⟩ - b) m(dy) db`. -/
theorem ex_operator_layer_iii_a (m : Measure Ω) [IsFiniteMeasure m] (a : Ω → H) (b : Ω → Y)
    (hL : IsLayerData m a b) :
    ∀ x : H, operatorLayer m a b gaussianFun x =
      ∫ p : Ω × ℝ, ((gaussianActDeriv2 p.2 * relu (⟪a p.1, x⟫ - p.2) : ℝ) : ℂ) • b p.1
        ∂(m.prod volume) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_iii_b, theorem in Challenge/Examples.lean:

omit [CompleteSpace H] [SecondCountableTopology H] [SecondCountableTopology Y] in
/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  ReLU form:
the Gaussian-activation layer is the ReLU network with the coefficient measure
`(y,b) ↦ (a_y, -b)`-pushforward of `φ''(b) b_y m(dy) db`. -/
theorem ex_operator_layer_iii_b (m : Measure Ω) [IsFiniteMeasure m] (a : Ω → H) (b : Ω → Y)
    (hL : IsLayerData m a b) :
    operatorLayer m a b gaussianFun =
      integralNetwork (fun t => (relu t : ℂ)) (layerHingeMeasure m a b) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_iii_c, theorem in Challenge/Examples.lean:

omit [CompleteSpace H] [SecondCountableTopology H] [BorelSpace H] [SecondCountableTopology Y] in
/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  ReLU form:
the ReLU coefficient measure is finite. -/
theorem ex_operator_layer_iii_c (m : Measure Ω) [IsFiniteMeasure m] (a : Ω → H) (b : Ω → Y)
    (hL : IsLayerData m a b) :
    IsFiniteMeasure (layerHingeMeasure m a b).variation := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_iii_d, theorem in Challenge/Examples.lean:

omit [CompleteSpace H] [SecondCountableTopology H] [BorelSpace H] [SecondCountableTopology Y] in
/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.  ReLU form:
the ReLU coefficient measure has all moments finite. -/
theorem ex_operator_layer_iii_d (m : Measure Ω) [IsFiniteMeasure m] (a : Ω → H) (b : Ω → Y)
    (hL : IsLayerData m a b) :
    ∀ k : ℕ, ∫⁻ θ : H × ℝ, ENNReal.ofReal ((1 + ‖θ.1‖ + |θ.2|) ^ k)
      ∂(layerHingeMeasure m a b).variation < ⊤ := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_operator_layer_iv, theorem in Challenge/Examples.lean:

omit [CompleteSpace H] [SecondCountableTopology H] [MeasurableSpace H] [BorelSpace H]
  [SecondCountableTopology Y] in
/-- **Example [ex:operator-layer]** Neural-operator layer as an integral network.
Non-cylindricity: if `A` has infinite rank, `β = Φ`, and `w_φ > 0` `m`-almost everywhere, then
`F_φ` is not cylindrical. -/
theorem ex_operator_layer_iv (m : Measure Ω) [IsFiniteMeasure m] (a : Ω → H) (b : Ω → Y)
    (hL : IsLayerData m a b) (φ : Y) (hA : HasInfiniteRank (layerA m a))
    (hw : ∀ᵐ y ∂m, 0 < (layerWeight b φ y).re ∧ (layerWeight b φ y).im = 0) :
    ¬ IsCylindrical (layerObservable m a b gaussianFun φ) := by

Status: verified by comparator.

8.7.5. Example 7.5 — Periodic convolution layer (ex:convolution)🔗

Blueprint node: Proposition 6.4.4. Status: verified (all 13 Lean theorems verified).

Formalization note. The 2026-09-13 revision added the hypothesis sigma = sigma_Gauss to the non-cylindricity clause, which is the hypothesis ex_convolution_x already carries (gaussianFun): the manuscript was corrected to the formalized statement, and the item is in sync.

OperatorRidgelet.Paper.ex_convolution_i, theorem in Challenge/Examples.lean:

/-- **Example [ex:convolution]** Periodic convolution layer.  With `a_y = k(y - ·)`,
`⟨a_y, x⟩ = (k * x)(y)`. -/
theorem ex_convolution_i (d : ℕ) (k : TorusL2 d) :
    ∀ (x : TorusL2 d) (y : Torus d),
      ⟪convDirection k y, x⟫ = ∫ t, k (y - t) * x t ∂torusHaar d := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_convolution_ii, theorem in Challenge/Examples.lean:

/-- **Example [ex:convolution]** Periodic convolution layer.  With `a_y = k(y - ·)` and
`b_y = ψ(· - y)`, the layer is `ℱ(x) = ψ * β(k * x)`. -/
theorem ex_convolution_ii (d : ℕ) (k ψ : TorusL2 d) (β : ℝ → ℝ) (hβc : Continuous β)
    (hβp : HasPolynomialGrowth β) :
    ∀ x : TorusL2 d,
      ⇑(operatorLayer (torusHaar d) (convDirection k) (convOutput ψ) β x)
        =ᵐ[torusHaar d] fun t =>
          ∫ y, ((ψ (t - y) * β ⟪convDirection k y, x⟫ : ℝ) : ℂ) ∂torusHaar d := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_convolution_iii, theorem in Challenge/Examples.lean:

/-- **Example [ex:convolution]** Periodic convolution layer.  `‖A‖_∞ = ‖k‖₂`. -/
theorem ex_convolution_iii (d : ℕ) (k : TorusL2 d) :
    layerSupNorm (convDirection k) = ‖k‖ := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_convolution_iv, theorem in Challenge/Examples.lean:

/-- **Example [ex:convolution]** Periodic convolution layer.  `∫ ‖b_y‖ dy = ‖ψ‖₂`. -/
theorem ex_convolution_iv (d : ℕ) (ψ : TorusL2 d) :
    ∫ y, ‖convOutput ψ y‖ ∂torusHaar d = ‖ψ‖ := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_convolution_v, theorem in Challenge/Examples.lean:

/-- **Example [ex:convolution]** Periodic convolution layer.  The convolution layer satisfies the
standing hypotheses of the neural-operator layer (`y ↦ a_y`, `y ↦ b_y` are continuous and
bounded into `L²`), so Example `ex:operator-layer` applies. -/
theorem ex_convolution_v (d : ℕ) (k ψ : TorusL2 d) :
    IsLayerData (torusHaar d) (convDirection k) (convOutput ψ) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_convolution_vi, theorem in Challenge/Examples.lean:

/-- **Example [ex:convolution]** Periodic convolution layer.  `ℱ` commutes with all translations
of `𝕋^d`: `ℱ(τ_z x) = τ_z ℱ(x)`. -/
theorem ex_convolution_vi (d : ℕ) (k ψ : TorusL2 d) (β : ℝ → ℝ) (hβc : Continuous β)
    (hβp : HasPolynomialGrowth β) :
    ∀ (z : Torus d) (x : TorusL2 d),
      operatorLayer (torusHaar d) (convDirection k) (convOutput ψ) β (torusTranslate d z x) =
        torusTranslateC d z
          (operatorLayer (torusHaar d) (convDirection k) (convOutput ψ) β x) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_convolution_vii, theorem in Challenge/Examples.lean:

/-- **Example [ex:convolution]** Periodic convolution layer.  `ℱ` commutes with every isometry
`σ` of `𝕋^d` (an isometric automorphism of the group, measure preserving) that fixes `k` and
`ψ`: `ℱ(x ∘ σ) = ℱ(x) ∘ σ`. -/
theorem ex_convolution_vii (d : ℕ) (k ψ : TorusL2 d) (β : ℝ → ℝ) (hβc : Continuous β)
    (hβp : HasPolynomialGrowth β) :
    ∀ (σ : Torus d ≃+ Torus d), Isometry σ →
      ∀ hσ : MeasurePreserving σ (torusHaar d) (torusHaar d),
        (fun t => k (σ t)) =ᵐ[torusHaar d] k → (fun t => ψ (σ t)) =ᵐ[torusHaar d] ψ →
        ∀ x : TorusL2 d,
          operatorLayer (torusHaar d) (convDirection k) (convOutput ψ) β
              (Lp.compMeasurePreserving σ hσ x) =
            Lp.compMeasurePreserving σ hσ
              (operatorLayer (torusHaar d) (convDirection k) (convOutput ψ) β x) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_convolution_viii, theorem in Challenge/Examples.lean:

/-- **Example [ex:convolution]** Periodic convolution layer.  If `k̂(n) ≠ 0` for infinitely many
`n ∈ ℤ^d`, then `A` has infinite rank. -/
theorem ex_convolution_viii (d : ℕ) (k : TorusL2 d)
    (hk : Set.Infinite {n : Fin d → ℤ | torusFourierCoeff (fun t => (k t : ℂ)) n ≠ 0}) :
    HasInfiniteRank (layerA (torusHaar d) (convDirection k)) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_convolution_ix, theorem in Challenge/Examples.lean:

/-- **Example [ex:convolution]** Periodic convolution layer.  With `φ ≡ 1` the observable is
`F_1(x) = ψ̂(0) ∫ β((k * x)(y)) dy`. -/
theorem ex_convolution_ix (d : ℕ) (k ψ : TorusL2 d) (β : ℝ → ℝ) (hβc : Continuous β)
    (hβp : HasPolynomialGrowth β) :
    ∀ x : TorusL2 d,
      layerObservable (torusHaar d) (convDirection k) (convOutput ψ) β (torusOne d) x =
        torusFourierCoeff (fun t => (ψ t : ℂ)) 0 *
          ∫ y, ((β ⟪convDirection k y, x⟫ : ℝ) : ℂ) ∂torusHaar d := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_convolution_x, theorem in Challenge/Examples.lean:

/-- **Example [ex:convolution]** Periodic convolution layer.  If `k̂(n) ≠ 0` for infinitely many
`n` and `ψ̂(0) ≠ 0`, then the observable `F_1` of the Gaussian-activation layer is not
cylindrical. -/
theorem ex_convolution_x (d : ℕ) (k ψ : TorusL2 d)
    (hk : Set.Infinite {n : Fin d → ℤ | torusFourierCoeff (fun t => (k t : ℂ)) n ≠ 0})
    (hψ : torusFourierCoeff (fun t => (ψ t : ℂ)) 0 ≠ 0) :
    ¬ IsCylindrical
      (layerObservable (torusHaar d) (convDirection k) (convOutput ψ) gaussianFun
        (torusOne d)) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_convolution_xi, theorem in Challenge/Examples.lean:

/-- **Example [ex:convolution]** Periodic convolution layer.  For `s > d/2`, the operator
`(I - Δ)^{-s}` is injective, positive, self-adjoint, and trace class. -/
theorem ex_convolution_xi (d : ℕ) (s : ℝ) (hs : (d : ℝ) / 2 < s) :
    IsTraceClassCovariance (besselOperator d s) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_convolution_xii, theorem in Challenge/Examples.lean:

/-- **Example [ex:convolution]** Periodic convolution layer.  `(I - Δ)^{-s}` is translation
invariant: `(I - Δ)^{-s} τ_z = τ_z (I - Δ)^{-s}`. -/
theorem ex_convolution_xii (d : ℕ) (s : ℝ) (hs : (d : ℝ) / 2 < s) :
    ∀ (z : Torus d) (x : TorusL2 d),
      besselOperator d s (torusTranslate d z x) = torusTranslate d z (besselOperator d s x) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_convolution_xiii, theorem in Challenge/Examples.lean:

/-- **Example [ex:convolution]** Periodic convolution layer.  With `Q = P = (I - Δ)^{-s}`,
`s > d/2`, the transform is equivariant: `R_ρ[f ∘ τ_z](a,c) = R_ρ f(τ_z a, c)` for every
translation `τ_z` and every `f ∈ L¹(μ_Q)`. -/
theorem ex_convolution_xiii (d : ℕ) (s : ℝ) (hs : (d : ℝ) / 2 < s)
    [MeasurableSpace (TorusL2 d)] [BorelSpace (TorusL2 d)] (μ : Measure (TorusL2 d))
    [IsProbabilityMeasure μ] (hμ : IsCenteredGaussian (besselOperator d s) μ)
    (ρ : SchwartzMap ℝ ℝ) (f : TorusL2 d → ℂ) (hf : Integrable f μ) :
    ∀ (z : Torus d) (p : TorusL2 d × ℝ),
      ridgelet μ ρ (fun x => f (torusTranslate d z x)) p =
        ridgelet μ ρ f (torusTranslate d z p.1, p.2) := by

Status: verified by comparator.

8.7.6. Example 7.6 — Dirichlet solution operator with a pointwise nonlinearity (ex:dirichlet)🔗

Blueprint node: Proposition 6.4.6. Status: verified (all 11 Lean theorems verified).

Formalization note. Ω = (0,1) is the subtype UnitOpenInterval with Lebesgue measure volume (a probability measure), H = UnitL2 = Lp ℝ 2 volume, Y = UnitL2C (complex), g = dirichletKernel, 𝖦 = dirichletOperator (the integral operator with kernel g, chosen through its defining property), λ_n = dirichletEigenvalue, e_n = dirichletEigenfunction (e_0 = 0), a_y = dirichletDirection y = g(y,·), b_y = dirichletOutput y, 𝖦_N = dirichletReLUTruncation N (spectralReLUNetwork over n ∈ [1,N]). i: (𝖦x)(y) = ∫ g(y,t)x(t)dt a.e.; ii: for continuous x, u = 𝖦x has u(0) = u(1) = 0, is differentiable on (0,1), and u'' = u - x there; iii: eigenpairs for n ≥ 1; iv: IsTraceClassCovariance 𝖦 (injective, positive, self-adjoint, trace class); v: infinite rank; vi: sup_y ‖g(y,·)‖₂ < ∞ (= ‖A‖_∞ by definition); vii: IsLayerData (ex:operator-layer applies); viii: ℱ(x) = 𝖦β(𝖦x); ix: the exact ReLU network as a HasSum; x: ‖𝖦x - 𝖦_N x‖ ≤ λ_{N+1}‖x‖; xi: λ_{N+1} ≤ π⁻²(N+1)⁻² (the O(N^{-2}) rate). The remark that with Q = P = 𝖦 the input measure is the law of the solution with white-noise source is not formalized. Prose corrections in the 2026-09-13 revision: the layer is described as the nonlinear correction in the Picard update u_1 = L_D^{-1}x + F(x) for -u'' + u = sigma(u) + x; with Q = P = L_D^{-1} the input measure is the Gaussian field with covariance L_D^{-1} (applying L_D^{-1} to white noise would give L_D^{-2}); and the earlier remark that no rate is available for the nonlinear layer is replaced by a reference to the uniform sampling estimate of ex:operator-layer. The mathematical content of the formalized statements is unchanged.

OperatorRidgelet.Paper.ex_dirichlet_i, theorem in Challenge/Examples.lean:

/-- **Example [ex:dirichlet]** Dirichlet solution operator with a pointwise nonlinearity.  `𝖦` is
the integral operator with kernel `g`: `(𝖦x)(y) = ∫₀¹ g(y,t) x(t) dt`. -/
theorem ex_dirichlet_i :
    ∀ x : UnitL2, (dirichletOperator x : UnitOpenInterval → ℝ) =ᵐ[volume]
      fun y : UnitOpenInterval => ∫ t : UnitOpenInterval, dirichletKernel y t * x t := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_dirichlet_ii, theorem in Challenge/Examples.lean:

/-- **Example [ex:dirichlet]** Dirichlet solution operator with a pointwise nonlinearity.  For a
continuous source `x`, `u = 𝖦x` solves `-u'' + u = x` on `(0,1)` with `u(0) = u(1) = 0`. -/
theorem ex_dirichlet_ii (x : ℝ → ℝ) (hx : Continuous x) :
    dirichletSolution x 0 = 0 ∧ dirichletSolution x 1 = 0 ∧
      (∀ y ∈ Set.Ioo (0 : ℝ) 1, DifferentiableAt ℝ (dirichletSolution x) y) ∧
      ∀ y ∈ Set.Ioo (0 : ℝ) 1,
        HasDerivAt (deriv (dirichletSolution x)) (dirichletSolution x y - x y) y := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_dirichlet_iii, theorem in Challenge/Examples.lean:

/-- **Example [ex:dirichlet]** Dirichlet solution operator with a pointwise nonlinearity.  `𝖦` has
the eigenpairs `λ_n = (1 + π²n²)⁻¹`, `e_n(t) = √2 sin(nπt)`, `n ≥ 1`. -/
theorem ex_dirichlet_iii :
    ∀ n : ℕ, 1 ≤ n →
      dirichletOperator (dirichletEigenfunction n) =
        dirichletEigenvalue n • dirichletEigenfunction n := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_dirichlet_iv, theorem in Challenge/Examples.lean:

/-- **Example [ex:dirichlet]** Dirichlet solution operator with a pointwise nonlinearity.  `𝖦` is
injective, positive, self-adjoint, and trace class. -/
theorem ex_dirichlet_iv : IsTraceClassCovariance dirichletOperator := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_dirichlet_v, theorem in Challenge/Examples.lean:

/-- **Example [ex:dirichlet]** Dirichlet solution operator with a pointwise nonlinearity.  `𝖦` has
infinite rank. -/
theorem ex_dirichlet_v : HasInfiniteRank (dirichletOperator : UnitL2 →ₗ[ℝ] UnitL2) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_dirichlet_vi, theorem in Challenge/Examples.lean:

/-- **Example [ex:dirichlet]** Dirichlet solution operator with a pointwise nonlinearity.
`‖A‖_∞ ≤ sup_y ‖g(y,·)‖₂ < ∞`: the directions `a_y = g(y,·)` are bounded in `L²(0,1)`. -/
theorem ex_dirichlet_vi :
    BddAbove (Set.range fun y : UnitOpenInterval => ‖dirichletDirection y‖) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_dirichlet_vii, theorem in Challenge/Examples.lean:

/-- **Example [ex:dirichlet]** Dirichlet solution operator with a pointwise nonlinearity.  With
`a_y = b_y = g(y,·)` the standing hypotheses of the neural-operator layer hold, so Example
`ex:operator-layer` applies. -/
theorem ex_dirichlet_vii : IsLayerData volume dirichletDirection dirichletOutput := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_dirichlet_viii, theorem in Challenge/Examples.lean:

/-- **Example [ex:dirichlet]** Dirichlet solution operator with a pointwise nonlinearity.  With
`a_y = b_y = g(y,·)` the layer is `ℱ(x) = 𝖦 β(𝖦x)`. -/
theorem ex_dirichlet_viii (β : ℝ → ℝ) (hβc : Continuous β) (hβp : HasPolynomialGrowth β) :
    ∀ x : UnitL2, operatorLayer volume dirichletDirection dirichletOutput β x =
      Complex.ofRealCLM.compLp
        (dirichletOperator (toLpOrZero 2 volume fun t => β (dirichletOperator x t))) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_dirichlet_ix, theorem in Challenge/Examples.lean:

/-- **Example [ex:dirichlet]** Dirichlet solution operator with a pointwise nonlinearity.  `𝖦` is
the exact ReLU network `𝖦x = ∑_n λ_n e_n [ReLU(⟨e_n,x⟩) - ReLU(-⟨e_n,x⟩)]`. -/
theorem ex_dirichlet_ix :
    ∀ x : UnitL2, HasSum
      (fun n : ℕ => (dirichletEigenvalue n *
        (relu ⟪dirichletEigenfunction n, x⟫ - relu (-⟪dirichletEigenfunction n, x⟫))) •
          dirichletEigenfunction n)
      (dirichletOperator x) := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_dirichlet_x, theorem in Challenge/Examples.lean:

/-- **Example [ex:dirichlet]** Dirichlet solution operator with a pointwise nonlinearity.  The
`2N`-neuron truncation has error `‖𝖦x - 𝖦_N x‖ ≤ λ_{N+1} ‖x‖`, uniformly on bounded sets. -/
theorem ex_dirichlet_x :
    ∀ (n : ℕ) (x : UnitL2),
      ‖dirichletOperator x - dirichletReLUTruncation n x‖ ≤ dirichletEigenvalue (n + 1) * ‖x‖ := by

Status: verified by comparator.

OperatorRidgelet.Paper.ex_dirichlet_xi, theorem in Challenge/Examples.lean:

/-- **Example [ex:dirichlet]** Dirichlet solution operator with a pointwise nonlinearity.  The
truncation error is `O(N^{-2})`: `λ_{N+1} ≤ π⁻² (N+1)⁻²`. -/
theorem ex_dirichlet_xi :
    ∀ n : ℕ, dirichletEigenvalue (n + 1) ≤ (Real.pi ^ 2 * ((n : ℝ) + 1) ^ 2)⁻¹ := by

Status: verified by comparator.