Samplinglib
Lean gate not recorded for this source state main · 0e31a3cda412
ASTIS mathematical exposition

The exact output law is a retained Dirac factor times a conditional law

AutoSamplingTheory.TechnicalLemmas.Probability.HeatBath.heatBathSnd_apply · theorem · Teaching coverage

Statement

For every input pair x under the finite-measure conditional-existence hypotheses, the measure returned by heatBathSnd μ at x is exactly the ordinary product of the Dirac measure at x's first coordinate and the selected conditional law of the second coordinate given that value. This is a pointwise identity for the chosen kernel, not pointwise uniqueness of conditional distributions.

\[H_\mu(x)=\delta_{x_1}\otimes q(x_1)\qquad(x\in\alpha\times\beta).\]

All objects and hypotheses

  • {α β : Type*} — Two possibly universe-polymorphic coordinate types α and β.
  • [MeasurableSpace α] [MeasurableSpace β] — Both carry measurable-space structures; α × β has the product sigma-algebra.
  • [StandardBorelSpace β] — The resampled space β is Standard Borel, supplying the regular-conditional-distribution existence contract.
  • [Nonempty β] — β is nonempty, so the selected conditional kernel can be totalized as a probability law even on marginal-null fibers.
  • (μ : Measure (α × β)) [IsFiniteMeasure μ] — μ is any finite joint measure. It need not be normalized and may be zero. α need not be Standard Borel or nonempty.
  • (x : α × β) — An arbitrary input x=(a,b), with no positivity or support condition.

Notation and interpretation

m

The first marginal, represented by μ.map Prod.fst.

\[m=(\operatorname{fst})_\#\mu\]
q

The particular regular conditional probability kernel selected by Mathlib, not a pointwise unique conditional density.

\[q(a)=\operatorname{condDistrib}(\operatorname{snd}\mid\operatorname{fst};\mu)(a)\]
μ ⊗ₘ q

Composition-product of a measure and a kernel, for measurable A ⊆ α × β. It is a joint measure, unlike ordinary kernel composition which forgets the intermediate state.

\[(m\otimes_m q)(A)=\int_\alpha q(a,\{b:(a,b)\in A\})\,m(da)\]
K ∘ₘ μ

The measure after one kernel application. Invariant means exactly K ∘ₘ μ = μ.

\[(K\circ_m\mu)(A)=\int K(x,A)\,\mu(dx)\]
  • All sets used to test equality of measures are measurable unless explicitly stated otherwise.
  • All measure values, reciprocals and scalar multiplications in the Lean conditional formula are in ℝ≥0∞. Integral displays of kernel probabilities denote nonnegative extended integrals, not Bochner integrals with unstated integrability hypotheses.
  • The notation r_i, F_i, e_i, q_i and K_i in the teaching text is explanatory mathematical shorthand, not new public Lean definitions.
  • Markovness quantifies over every input; a selected conditional version's law characterization is only marginal-a.e. unless a positive-atom/fiber theorem applies.
  • Finite target mass is not probability normalization. Fin (n+1) counts sites, not states. No law is assumed to have a density.

Mathematical proof

1. Replace deterministic precomposition by input substitution

Composing a kernel with a deterministic map evaluates that kernel at the mapped input. The identity is between kernels and does not require a singleton-measurable retained space.

\[(L\circ_k D_f)(x)=L(f(x)),\qquad f=\operatorname{fst}.\]
Corresponding Lean step

Kernel.comp_deterministic_eq_comap turns the composition into a comap; Kernel.comap_apply evaluates the comap.

2. Evaluate the product kernel

The product kernel evaluated at x1 is the product of its two measure values at x1. Markovness supplies the s-finiteness needed by the kernel-product evaluation machinery.

\[(\operatorname{Id}_\alpha\times_k q)(x_1)=\operatorname{Id}_\alpha(x_1)\otimes q(x_1).\]
Corresponding Lean step

Kernel.prod_apply.

3. Evaluate the identity kernel

The identity kernel returns the Dirac measure at its input. Substituting this gives the desired measure identity. The x2 component never enters the expression.

\[\operatorname{Id}_\alpha(x_1)=\delta_{x_1},\qquad H_\mu(x)=\delta_{x_1}\otimes q(x_1).\]
Corresponding Lean step

Kernel.id_apply completes the rw chain.

Lean statement · heatBathSnd_apply

Once a starting pair is fixed, the update's whole law can be written down compactly: a point mass on the retained coordinate, multiplied by the conditional law on the updated coordinate. The equality is valid even at a null conditioning value because it refers to the version already selected.

Braces mark parameters Lean can infer; square brackets request structures such as a measurable space or probability measure. Named hypotheses are mathematical premises, not facts established by this declaration. Section parameters are described in the mathematical hypotheses above; the module link retains their exact source context.

@[simp] theorem heatBathSnd_apply (μ : Measure (α × β)) [IsFiniteMeasure μ]
    (x : α × β) :
    heatBathSnd μ x = (Measure.dirac x.1).prod (condDistrib Prod.snd Prod.fst μ x.1)

Exact module and namespace context

Lean proof · heatBathSnd_apply

`@[simp]` marks the proven equality as a standard simplification rule for later uses; it does not weaken the statement. `rw [heatBathSnd, ...]` performs a sequence of equality rewrites, each corresponding to one of the evaluation steps. `x.1` means the first projection of x. `Measure.prod` is a product of two measures, whereas the construction used the product `×ₖ` of two kernels. Equality of this retained marginal law with a Dirac measure is not automatically literal almost-everywhere equality of retained values on an arbitrary coarse sigma-algebra; the coordinate-retention theorem makes its singleton assumption explicit.

Braces mark parameters Lean can infer; square brackets request structures such as a measurable space or probability measure. Named hypotheses are mathematical premises, not facts established by this declaration. Section parameters are described in the mathematical hypotheses above; the module link retains their exact source context.

@[simp] theorem heatBathSnd_apply (μ : Measure (α × β)) [IsFiniteMeasure μ]
    (x : α × β) :
    heatBathSnd μ x = (Measure.dirac x.1).prod (condDistrib Prod.snd Prod.fst μ x.1) := by
  rw [heatBathSnd, Kernel.comp_deterministic_eq_comap, Kernel.comap_apply,
    Kernel.prod_apply, Kernel.id_apply]

/-- A second-coordinate heat-bath update leaves its finite joint target invariant.

The zero measure is allowed. This theorem does not imply irreducibility,
reversibility, convergence from another initial law, or any mixing rate. -/

Exact module and namespace context

Scope and omitted-condition boundaries

  • These are existing source-neutral measure/kernel interfaces, not newly proved claims in this research packet. Reading source and historical tests is not a fresh Lean build or independent verification.
  • Finite target measure means μ(univ) < ∞, not μ(univ) = 1; zero and finite unnormalized targets are allowed unless a positive-fiber hypothesis rules them out.
  • A selected regular conditional distribution is a probability kernel at every conditioning value. Its conditional-law characterization is determined only almost everywhere under the conditioning marginal; null fibers have no asserted feasible/Gibbs-support property.
  • No reversibility, irreducibility, ergodicity, convergence from a different initial law, mixing rate, random scan, scan-order correctness, executable conditional sampling or cost bound follows from these declarations alone.

Source and reuse

ASTIS parents called

Mathlib API called (external library)

  • ProbabilityTheory.Kernel.comp_deterministic_eq_comap
  • ProbabilityTheory.Kernel.comap_apply
  • ProbabilityTheory.Kernel.prod_apply
  • ProbabilityTheory.Kernel.id_apply

Mathematical sources

Older frontier cells and the coordinate card say the printed copy-index repair is pending. The later positive-fiber card and source-support audit record separate acceptance of the exact one-index correction as ASTIS-REPAIR-20260908-GlauberCopyIndex. The old PositiveFiberUpdate audit remains possible-source-error for its limited page-5 evidence window. Do not overwrite that history, call the original PDF corrected, call a null-fiber start a counterexample to the positive-support source convention, or infer full source fidelity. This packet reads local evidence only; it performs no new PDF/source certification.

ASTIS prose is not a quotation or a source-equivalence certificate. Definitions and aliases are explained as constructions, not counted as new mathematical proofs.