Skip to content

W7 · CP-4 follow-up: the NaN-under-both draws and all-transect non-finite gradients — singular systems or a formula-independent overflow? #164

Description

@Jammy2211

From phase_08_regularization/RESULTS.md (human call 2026-08-24: adopt slogdet, chase the rest).

Questions

  1. The 32 (A100) / 20 (CPU) delaunay_adapt_split draws that are NaN under both cholesky and slogdet: are the curvature+reg matrices genuinely singular (λ⁴ fragility, feat(jax_compile): make warm compile machine-identifiable via cache_state #104 — eigen-spectrum at those draws), or does the NaN arise upstream (mapping matrix, data-vector overflow, the sqrt(dual_area) grad hazard) so no log-det formula could help? Replay from slogdet_ab/draws/*.npz with per-term instrumentation.
  2. The λ-transect gradients are non-finite under both arms on all 128 draws on both tiers (while prior/truth-bar draws mostly are finite): identify the term (is it the transect's log-coefficient axis itself?) — this may be a driver artefact rather than a likelihood hazard.
  3. The marginal-band disagreement (up to 9,619 nats on A100 vs 2.27 on CPU): tier-dependent — characterise the coefficient range and whether it is fp64 LU pivoting on GPU.

Deliverable

RESULTS.md addendum with attribution per draw class; a recommendation on whether the PyAutoArray default should flip (feeds W9).

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions