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.
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.
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.
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.
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)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. -/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
- Fearnhead–Nemeth–Oates–Sherlock, arXiv:2407.12751v1, §2.1.1 Component-wise updates and Gibbs moves; printed page 48 / PDF page 54 — Conditional-update motivation only; the generic finite-measure regular-conditional-law contract is not a numbered-source theorem assimilation.
- Chen–Štefankovič–Vigoda, arXiv:2307.13826v4, §1.3 Glauber dynamics/Gibbs sampler; printed/PDF page 5; support convention in §1.1 printed/PDF page 4 — Fixed-site update motivation, not full finite-binary-state algorithm or uniform-scan certification.
- Existing ASTIS declaration AutoSamplingTheory.TechnicalLemmas.Probability.HeatBath.heatBathSnd_apply — Exact public declaration and proof/construction read from the local source. Line numbers are current source anchors; no new source certification or Lean build performed.
- Focused tests (consumer evidence) — An example consumes exactly this pointwise equality for arbitrary x.
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.