Skip to main content

Follow the source behind the result.

Search the paper sections, open their original context and follow the checks they actually cite. The knowledge atlas explains the result; this library locates its source.

v849_cfin_unique_cofinal_lean.py

verification/v849_cfin_unique_cofinal_lean.py

This is a registered verification module. Read its original code and the paper sections that cite it, then follow the linked claim to its present classification.

Read the publication code on this page

The registered Python source is shown as text. Displaying it does not execute the code.

"""v849 TFPT scope: finite compiler uniqueness, CFIN.UNIQUE.01.

Scope extraction from original module at ba2f0fc; original ID and
mathematical statements preserved. Original full source is retained
in the extraction provenance; no proof/status promotion.
"""
import contextlib

import io

import os

import re

import sys

import types

_HERE = os.path.dirname(os.path.abspath(__file__))

_SRC_CFIN = r'''
#!/usr/bin/env python3
# -*- coding: utf-8 -*-
r"""cfin_uniqueness_probe -- CFIN.UNIQUE.01: the missing canonicity
theorem for the four-bit compiler

      C_fin = (V, omega, q*, sigma, iota),

machine-decided by exhaustive finite computation (work package F1
side B, 2026-08-07 evening plan; the uniqueness companion of
NORMALFORM.CFIN.01 = finite_compiler_normal_form_probe.py,
NORMAL-FORM-ASSEMBLED 28/28, whose object and conventions are frozen
inputs here, not re-adjudicated).

FROZEN CLAIMS (frozen + SHA-hashed before first run):

THE ADMISSIBLE CLASS (the abstract compiler object, exactly the
certified structural properties of C_fin):
  a tuple (omega, q, sigma, iota) on the carrier V = F_2^4 with
  (A1) omega  : nondegenerate ALTERNATING bilinear form;
  (A2) q     : quadratic refinement of omega
               (q(v+w)+q(v)+q(w) = omega(v,w)) of ARF TYPE 1;
  (A3) sigma : order-3 operator (sigma^3 = id, sigma != id),
               symplectic for omega, preserving q, with the
               family-cycle fixed-class structure: EXACTLY 3 nonzero
               fixed classes (fixed space dim 2; the remaining 12
               classes then fall in 4 free 3-orbits, checked);
  (A4) iota  : linear parity lift V -> F_2^5 with image = the even
               five-slot code C_even(5), beta(v,w) =
               iota(v).iota(w) mod 2 == omega on all 256 cells,
               q_wt = wt(iota)/2 mod 2 == q, and sigma-equivariant:
               SOME slot permutation pi in S_5 satisfies
               iota(sigma v) = pi(iota(v)) (the S4.4 register
               intertwining of the normal form).
  ISOMORPHISM: (g, pi) in GL(4,2) x S_5 transporting all four data;
  after U0 (form transitivity) the classification group is
  G = Sp(4,2) x S_5, |G| = 720 x 120 = 86400, acting on (q, sigma,
  iota) at fixed omega = hbar (the J - I Gram convention of the
  normal form).

U0  FORM REDUCTION: census of all 64 alternating bilinear forms on
    F_2^4; the nondegenerate ones are ONE GL(4,2)-orbit containing
    hbar (so WLOG omega = hbar); |GL(4,2)| = 20160, |Sp(4,2)| = 720
    rebuilt by exhaustive enumeration.
U1  REFINEMENT REDUCTION: exactly 16 refinements of hbar (2^16
    brute force, the parent census); Arf census 6 + 10; Sp(4,2)
    transitive on the 6 Arf-1 refinements (so WLOG q = q*).
U2  SIGMA CENSUS: for each Arf-1 q, the admissible sigma set
    (order 3, symplectic, q-preserving, 3 nonzero fixed classes) is
    enumerated exactly; the free-orbit structure (4 x 3) checked;
    the canonical family 3-cycle is in the set for q = q*.
U3  IOTA CENSUS: for each admissible (q, sigma), the admissible
    iota set is enumerated exactly (DFS over basis images with dot
    constraints); the equivariant subcensus separately; the
    canonical iota is in the set.
U4  (a) UNIQUENESS: the FULL admissible tuple set at omega = hbar
    is enumerated; the G-orbit of the C_fin tuple is computed by
    applying all 86400 group elements; PASS = the orbit IS the
    admissible set (exactly one orbit) => every admissible compiler
    object on any carrier isomorphic to F_2^4 is isomorphic to
    C_fin (with U0/U1 as the reduction steps).  If more orbits:
    the typed orbit census.
U5  (b) THE AUTOMORPHISM GROUP: Aut(C_fin) = the joint stabilizer
    of (omega, q*, sigma, iota) in G, computed exactly: order
    (cross-checked against orbit-stabilizer |G| = |orbit| x |Aut|),
    element-order census, abelian/cyclic flags, identification
    against small groups (C6 vs S3 vs ...), the projection to
    Sp(4,2) (is pi determined by g?), and the joint stabilizer of
    (omega, q*, sigma, iota) IN Sp(4,2) (the assignment's reading).
U6  (c) UNIVERSAL PROPERTY CANDIDATE (measured, not narrated): the
    Gaussian-reduction diagram -- the E8 root census + deck action
    (rebuilt inline, parent S1 machinery) -- fixes (V, hbar, sigma)
    up to a MEASURED number of identifications (= the centralizer
    of sigma in Sp(4,2), cross-checked); the completions to a full
    C_fin object are measured: q-completions 4 -> 2 -> 1 under the
    frozen selector pins (sigma-invariance, q(A) = 1,
    q(F_Sigma) = 0), iota-completions counted and their orbit
    structure under the joint stabilizer computed.  TYPE WHAT
    HOLDS: initial/terminal (unique up to UNIQUE isomorphism) holds
    iff the completion is unique AND Aut(C_fin) = 1; otherwise the
    honest reading is "unique up to (non-unique) isomorphism, with
    |Aut(C_fin)| arrows".

C   CONTROLS (must fire; frozen fire rules):
    C1 DOT FORM: the non-alternating dot form admits ZERO quadratic
       refinements (2^16 brute force) -- admissible set EMPTY.
    C2 WRONG ARF: for EVERY Arf-0 refinement, the admissible iota
       set (even without any sigma condition) is EMPTY -- the even
       five-slot weight census (1, 10, 5) forces q_wt to be Arf-1.
    C3 NO SIGMA (trivial deck, parent C2): without the order-3
       datum the SELECTOR loses uniqueness: the pinned completion
       census on the Gaussian carrier becomes 16 -> 8 -> 4 != 1
       (fire iff != 1).  TYPED SIDE MEASUREMENT (predicted, frozen
       as a report line, not a fire rule): the BARE (q, iota) pair
       set without sigma stays ONE Sp x S_5 orbit -- sigma is
       load-bearing for CANONICITY (the selector), not for the bare
       orbit count; the measured value is printed either way.

KILLS (any one fires => typed):
  K1 carrier / group / census guard breaks       -> CFIN-UNDERDETERMINED
  K2 orbit-stabilizer arithmetic inconsistent    -> CFIN-UNDERDETERMINED
  K3 admissible set empty                        -> CFIN-UNDERDETERMINED
  K4 a control does not fire                     -> CFIN-UNDERDETERMINED

VERDICT (frozen enum, decision order):
  0. any kill                          -> CFIN-UNDERDETERMINED, exit 1
  1. exactly one G-orbit               -> CFIN-UNIQUE (+ |Aut| printed)
  2. N > 1 orbits                      -> CFIN-ORBITS-N (typed census)

FIREWALL: experiments/ probe; EXPLORATION ONLY -- writes nothing but
stdout; no verification/, paper, ledger, changelog or website
surface; no .md, no commits.  Exact integer arithmetic only; no
floats, no RNG, no fits.  NO physics claim, NO RH claim.  Runtime
cap 900 s.

Sources (read-only, machinery rebuilt inline):
finite_compiler_normal_form_probe.py (the object, conventions,
selector, parent controls C1/C2), verification/v752 + v774 (quotient
+ compiler), v775 (census conventions), v833 (Gaussian rebuild).

Run:
    experiments/tfpt-discovery/.venv/bin/python \
        experiments/tfpt-discovery/cfin_uniqueness_probe.py
"""
import hashlib
import itertools
import math
import time

T0 = time.time()
CHECKS = []
KILLS = []


def check(name, ok, detail="", kill=None):
    CHECKS.append((name, bool(ok)))
    if kill and not ok:
        KILLS.append(kill)
    print("  [%s] %s%s" % ("PASS" if ok else "FAIL", name,
                           ("  -- " + detail) if detail else ""), flush=True)
    return bool(ok)


def section(title):
    print("\n== %s ==  (t=%.1fs)" % (title, time.time() - T0), flush=True)


print("CFIN.UNIQUE.01 -- uniqueness + automorphisms of "
      "C_fin = (V, omega, q*, sigma, iota)")
print("FROZEN_SPEC SHA-256: %s"
      % hashlib.sha256(__doc__.encode("utf-8")).hexdigest())

# ------------------------------------------------------- bit machinery
W16 = [tuple(b) for b in itertools.product((0, 1), repeat=4)]
WIDX = {w: i for i, w in enumerate(W16)}
GJI = ((0, 1, 1, 1), (1, 0, 1, 1), (1, 1, 0, 1), (1, 1, 1, 0))
BASIS = [WIDX[tuple(1 if k == i else 0 for k in range(4))]
         for i in range(4)]
ADD = [[WIDX[tuple((a + b) % 2 for a, b in zip(W16[i], W16[j]))]
        for j in range(16)] for i in range(16)]


def form_table(mat):
    return [[sum(W16[i][r] * mat[r][c] * W16[j][c] for r in range(4)
                 for c in range(4)) % 2 for j in range(16)]
            for i in range(16)]


HB = form_table(GJI)
A_BIT, FSIG = (0, 0, 0, 1), (1, 1, 1, 0)
LA, LF = WIDX[A_BIT], WIDX[FSIG]
SIGP = tuple(WIDX[(w[2], w[0], w[1], w[3])] for w in W16)


def perm_of_images(g1, g2, g3, g4):
    """Label permutation of the linear map sending e_i -> g_i."""
    gs = (g1, g2, g3, g4)
    p = []
    for i in range(16):
        acc = 0
        for k in range(4):
            if W16[i][k]:
                acc = ADD[acc][gs[k]]
        p.append(acc)
    return tuple(p)


def inv_perm(p):
    q = [0] * 16
    for i, pi in enumerate(p):
        q[pi] = i
    return tuple(q)


# ======================================================================
section("U0: form reduction -- alternating forms are one GL orbit")
# ======================================================================
GL, SP = [], []
for imgs in itertools.product(range(1, 16), repeat=4):
    p = perm_of_images(*imgs)
    if len(set(p)) != 16:
        continue
    GL.append(p)
    if all(HB[p[BASIS[i]]][p[BASIS[j]]] == GJI[i][j]
           for i in range(4) for j in range(i + 1, 4)):
        SP.append(p)
check("U0.1 |GL(4,2)| = %d == 20160; |Sp(4,2)| = %d == 720 "
      "(exhaustive image enumeration)" % (len(GL), len(SP)),
      len(GL) == 20160 and len(SP) == 720, kill="K1")

alt_nondeg = set()
for bits in range(64):
    mat = [[0] * 4 for _ in range(4)]
    k = 0
    for i in range(4):
        for j in range(i + 1, 4):
            mat[i][j] = mat[j][i] = (bits >> k) & 1
            k += 1
    tab = form_table(mat)
    if all(any(tab[i][j] for j in range(16)) for i in range(1, 16)):
        alt_nondeg.add(tuple(tuple(mat[r]) for r in range(4)))
orbit_forms = set()
for p in GL:
    ip = inv_perm(p)
    mat = tuple(tuple(HB[ip[BASIS[i]]][ip[BASIS[j]]] for j in range(4))
                for i in range(4))
    orbit_forms.add(mat)
check("U0.2 nondegenerate alternating forms: %d of 64; GL-orbit of "
      "hbar has %d elements == the whole set (WLOG omega = hbar)"
      % (len(alt_nondeg), len(orbit_forms)),
      orbit_forms == alt_nondeg and GJI in alt_nondeg, kill="K1")

# ======================================================================
section("U1: refinement reduction -- Sp transitive on the 6 Arf-1")
# ======================================================================
def refinements_of_table(tab):
    out = []
    for mask in range(1 << 16):
        q = tuple((mask >> i) & 1 for i in range(16))
        ok = True
        for i in range(16):
            for j in range(16):
                if q[ADD[i][j]] ^ q[i] ^ q[j] != tab[i][j]:
                    ok = False
                    break
            if not ok:
                break
        if ok:
            out.append(q)
    return out


REFS = refinements_of_table(HB)
ARF1 = sorted(q for q in REFS if sum(1 for x in q if x == 0) == 6)
ARF0 = sorted(q for q in REFS if sum(1 for x in q if x == 0) == 10)
check("U1.1 EXACTLY 16 refinements of hbar; Arf census 6 + 10",
      len(REFS) == 16 and len(ARF1) == 6 and len(ARF0) == 10,
      kill="K1")

siginv = [q for q in REFS if all(q[SIGP[i]] == q[i] for i in range(16))]
cand_a = [q for q in siginv if q[LA] == 1]
cand = [q for q in cand_a if q[LF] == 0]
QSTAR = cand[0] if len(cand) == 1 else None
check("U1.2 selector 16 -> sigma-invariant %d -> q(A)=1: %d -> "
      "q(F_Sigma)=0: %d == 1 (the canonical q*, Arf 1)"
      % (len(siginv), len(cand_a), len(cand)),
      len(siginv) == 4 and len(cand_a) == 2 and len(cand) == 1
      and QSTAR in ARF1, kill="K1")

SP_INV = [inv_perm(p) for p in SP]
orb_q = {tuple(QSTAR[ip[i]] for i in range(16)) for ip in SP_INV}
check("U1.3 Sp(4,2) TRANSITIVE on Arf-1 refinements: orbit of q* has "
      "%d elements == the 6 Arf-1 forms (WLOG q = q*)" % len(orb_q),
      orb_q == set(ARF1), kill="K1")

# ======================================================================
section("U2: the admissible sigma census per Arf-1 refinement")
# ======================================================================
def sigmas_for(q):
    out = []
    for p in SP:
        if p == tuple(range(16)):
            continue
        p2 = tuple(p[p[i]] for i in range(16))
        if tuple(p[p2[i]] for i in range(16)) != tuple(range(16)):
            continue
        if sum(1 for i in range(1, 16) if p[i] == i) != 3:
            continue
        if any(q[p[i]] != q[i] for i in range(16)):
            continue
        orbs = {frozenset({i, p[i], p2[i]})
                for i in range(1, 16) if p[i] != i}
        assert len(orbs) == 4 and all(len(o) == 3 for o in orbs)
        out.append(p)
    return out


SIGS = {q: sigmas_for(q) for q in ARF1}
n_sig = sorted(len(v) for v in SIGS.values())
check("U2.1 admissible sigma per Arf-1 q: counts %s (order 3, "
      "symplectic, q-preserving, 3 nonzero fixed, 4 free 3-orbits "
      "checked); canonical family cycle IN the q* set"
      % n_sig, all(n > 0 for n in n_sig) and SIGP in SIGS[QSTAR],
      kill="K1")

STAB_Q = [k for k in range(720)
          if tuple(QSTAR[SP_INV[k][i]] for i in range(16)) == QSTAR]
conj_orbs = []
seen = set()
for s in SIGS[QSTAR]:
    if s in seen:
        continue
    o = set()
    for k in STAB_Q:
        p, ip = SP[k], SP_INV[k]
        o.add(tuple(p[s[ip[i]]] for i in range(16)))
    seen |= o
    conj_orbs.append(o)
check("U2.2 |Stab_Sp(q*)| = %d == 120; its conjugation orbits on the "
      "admissible sigma set: %s (sizes)"
      % (len(STAB_Q), sorted(len(o) for o in conj_orbs)),
      len(STAB_Q) == 120
      and sum(len(o) for o in conj_orbs) == len(SIGS[QSTAR]),
      kill="K1")

# ======================================================================
section("U3: the admissible iota census per (q, sigma)")
# ======================================================================
B5_EVEN = [b for b in itertools.product((0, 1), repeat=5)
           if sum(b) % 2 == 0]
S5 = list(itertools.permutations(range(5)))


def iotas_for(q, sigp=None, require_equiv=False):
    """All admissible parity lifts; DFS over the four basis images."""
    sols = []

    def rec(imgs):
        i = len(imgs)
        if i == 4:
            iot = []
            for li in range(16):
                acc = (0, 0, 0, 0, 0)
                for k in range(4):
                    if W16[li][k]:
                        acc = tuple((a + b) % 2
                                    for a, b in zip(acc, imgs[k]))
                iot.append(acc)
            if len(set(iot)) != 16:
                return
            if any((sum(x) // 2) % 2 != q[li]
                   for li, x in enumerate(iot)):
                return
            sols.append(tuple(iot))
            return
        for b in B5_EVEN:
            if any(sum(x * y for x, y in zip(b, imgs[j])) % 2
                   != GJI[i][j] for j in range(i)):
                continue
            rec(imgs + [b])

    rec([])
    if not require_equiv:
        return sols
    out = []
    for iot in sols:
        for pi in S5:
            if all(iot[sigp[li]] == tuple(iot[li][pi[s]]
                                          for s in range(5))
                   for li in range(16)):
                out.append(iot)
                break
    return out


IOTA_CANON = tuple(tuple(list(w) + [sum(w) % 2]) for w in W16)
iot_bare_qstar = iotas_for(QSTAR)
iot_eq_canon = iotas_for(QSTAR, SIGP, require_equiv=True)
check("U3.1 iota census for q = q*: %d bare admissible lifts (image "
      "= C_even(5), beta = hbar, q_wt = q*); %d sigma-equivariant "
      "for the canonical family cycle; the canonical iota is one of "
      "them" % (len(iot_bare_qstar), len(iot_eq_canon)),
      len(iot_bare_qstar) > 0 and len(iot_eq_canon) > 0
      and IOTA_CANON in iot_eq_canon
      and IOTA_CANON in iot_bare_qstar, kill="K1")

# ======================================================================
section("U4: (a) UNIQUENESS -- the full tuple census vs the G-orbit")
# ======================================================================
ADMISSIBLE = set()
for q in ARF1:
    for s in SIGS[q]:
        for iot in iotas_for(q, s, require_equiv=True):
            ADMISSIBLE.add((q, s, iot))
N_ADM = len(ADMISSIBLE)
CFIN = (QSTAR, SIGP, IOTA_CANON)
check("U4.1 admissible tuple census at omega = hbar: N = %d > 0; the "
      "C_fin tuple is admissible" % N_ADM,
      N_ADM > 0 and CFIN in ADMISSIBLE, kill="K3")

orbit = set()
stab_pairs = []
for k in range(720):
    p, ip = SP[k], SP_INV[k]
    q2 = tuple(QSTAR[ip[i]] for i in range(16))
    s2 = tuple(p[SIGP[ip[i]]] for i in range(16))
    base = [IOTA_CANON[ip[i]] for i in range(16)]
    for pi in S5:
        i2 = tuple(tuple(x[pi[s]] for s in range(5)) for x in base)
        t = (q2, s2, i2)
        orbit.add(t)
        if t == CFIN:
            stab_pairs.append((k, pi))
N_ORB = len(orbit)
N_AUT = len(stab_pairs)
check("U4.2 orbit-stabilizer: |orbit| x |Aut| = %d x %d = %d == "
      "|G| = 86400" % (N_ORB, N_AUT, N_ORB * N_AUT),
      N_ORB * N_AUT == 86400, kill="K2")
check("U4.3 the G-orbit is INSIDE the admissible set (transport "
      "preserves admissibility)", orbit <= ADMISSIBLE, kill="K2")

n_orbits = None
if orbit == ADMISSIBLE:
    n_orbits = 1
else:
    rest = set(ADMISSIBLE) - orbit
    n_orbits = 1
    while rest:
        rep = next(iter(rest))
        o2 = set()
        for k in range(720):
            p, ip = SP[k], SP_INV[k]
            q2 = tuple(rep[0][ip[i]] for i in range(16))
            s2 = tuple(p[rep[1][ip[i]]] for i in range(16))
            base = [rep[2][ip[i]] for i in range(16)]
            for pi in S5:
                i2 = tuple(tuple(x[pi[s]] for s in range(5))
                           for x in base)
                o2.add((q2, s2, i2))
        rest -= o2
        n_orbits += 1
check("U4.4 UNIQUENESS: admissible set = %d tuples, G-orbits = %d "
      "(PASS bar: exactly 1 -- every admissible compiler object is "
      "isomorphic to C_fin)" % (N_ADM, n_orbits), n_orbits == 1,
      kill=None)

# ======================================================================
section("U5: (b) Aut(C_fin) -- the joint stabilizer, exactly")
# ======================================================================
def perm_order(p):
    o, pp = 1, p
    ident = tuple(range(len(p)))
    while pp != ident:
        pp = tuple(p[x] for x in pp)
        o += 1
    return o


orders = sorted(math.lcm(perm_order(SP[k]), perm_order(pi))
                for k, pi in stab_pairs)
sp_parts = {k for k, _ in stab_pairs}
pi_per_g = {k: [pi for kk, pi in stab_pairs if kk == k]
            for k in sp_parts}
faithful_sp = (len(sp_parts) == N_AUT
               and all(len(v) == 1 for v in pi_per_g.values()))
abelian = True
for (k1, p1), (k2, p2) in itertools.combinations(stab_pairs, 2):
    g12 = tuple(SP[k1][SP[k2][i]] for i in range(16))
    g21 = tuple(SP[k2][SP[k1][i]] for i in range(16))
    p12 = tuple(p1[p2[s]] for s in range(5))
    p21 = tuple(p2[p1[s]] for s in range(5))
    if g12 != g21 or p12 != p21:
        abelian = False
        break
cyclic = max(orders) == N_AUT if N_AUT > 0 else False
small = {(1,): "1", (2,): "C2", (3,): "C3",
         (6, True, True): "C6", (6, False, False): "S3",
         (4, True, True): "C4", (4, True, False): "C2xC2",
         (12, True, True): "C12"}
if N_AUT in (1, 2, 3):
    iso = small.get((N_AUT,), "order %d" % N_AUT)
else:
    iso = small.get((N_AUT, abelian, cyclic),
                    "order %d (abelian=%s, cyclic=%s)"
                    % (N_AUT, abelian, cyclic))
check("U5.1 |Aut(C_fin)| = %d; element orders %s; abelian=%s, "
      "cyclic=%s -> Aut(C_fin) ~= %s"
      % (N_AUT, orders, abelian, cyclic, iso), N_AUT >= 1, kill="K2")
check("U5.2 the projection Aut -> Sp(4,2) is faithful and the slot "
      "permutation pi is DETERMINED by g (|Sp-part| = %d, one pi per "
      "g: %s) -- Aut(C_fin) IS the joint stabilizer of (omega, q*, "
      "sigma, iota) in Sp(4,2)"
      % (len(sp_parts), all(len(v) == 1 for v in pi_per_g.values())),
      faithful_sp, kill="K2")

# ======================================================================
section("U6: (c) the Gaussian-reduction diagram -- what holds")
# ======================================================================
G_NAIVE = ((1, 0, 0, 0, 0, 1, 1, 1), (0, 1, 0, 0, 1, 0, 1, 1),
           (0, 0, 1, 0, 1, 1, 0, 1), (0, 0, 0, 1, 1, 1, 1, 0))
C_NAIVE = frozenset(
    tuple(sum(m[k] * G_NAIVE[k][j] for k in range(4)) % 2
          for j in range(8))
    for m in itertools.product((0, 1), repeat=4))
PI_J = (1, 0, 3, 2, 5, 4, 7, 6)
PI_SIG = (2, 3, 4, 5, 0, 1, 6, 7)


def code_image(code, perm):
    return frozenset(tuple(c[perm[k]] for k in range(8)) for c in code)


placements = set()
for perm in itertools.permutations(range(8)):
    placements.add(code_image(C_NAIVE, perm))
both = [c for c in placements
        if code_image(c, PI_J) == c and code_image(c, PI_SIG) == c]
CSTAR = next(c for c in both if (1, 0, 1, 0, 1, 0, 1, 0) in c)

ROOTS = []
for k in range(8):
    for sgn in (2, -2):
        v = [0] * 8
        v[k] = sgn
        ROOTS.append(tuple(v))
for c in (c for c in CSTAR if sum(c) == 4):
    sup = [k for k in range(8) if c[k]]
    for signs in itertools.product((1, -1), repeat=4):
        v = [0] * 8
        for k, sgn in zip(sup, signs):
            v[k] = sgn
        ROOTS.append(tuple(v))


def J_vec(v):
    out = [0] * 8
    for k in range(4):
        out[2 * k] = -v[2 * k + 1]
        out[2 * k + 1] = v[2 * k]
    return tuple(out)


def sig_vec(v):
    return tuple(v[PI_SIG[k]] for k in range(8))


def in_L(v):
    return tuple(x % 2 for x in v) in CSTAR


def in_pi2L(v):
    w = [0] * 8
    for k in range(4):
        w[2 * k] = v[2 * k] + v[2 * k + 1]
        w[2 * k + 1] = v[2 * k + 1] - v[2 * k]
    if any(x % 2 for x in w):
        return False
    return in_L([x // 2 for x in w])


REPS = [(0,) * 8]
for r in ROOTS:
    if not any(in_pi2L(tuple(a - b for a, b in zip(r, rep)))
               for rep in REPS):
        REPS.append(r)


def label_of(v):
    for li, rep in enumerate(REPS):
        if in_pi2L(tuple(a - b for a, b in zip(v, rep))):
            return li
    raise AssertionError("vector not in L")


census = {}
for r in ROOTS:
    census[label_of(r)] = census.get(label_of(r), 0) + 1


def sig_label(li):
    return label_of(sig_vec(REPS[li]))


fixed_nz = [li for li in range(1, 16) if sig_label(li) == li]
check("U6.1 carrier guard (parent S1 rebuilt): 30 placements, 2 "
      "deck-stable, 16 classes, census 240 = 15 x 16, 3 nonzero "
      "sigma-fixed classes",
      len(placements) == 30 and len(both) == 2 and len(REPS) == 16
      and sorted(census.values()) == [16] * 15 and 0 not in census
      and len(fixed_nz) == 3, kill="K1")


def ip8(x, y):
    return sum(a * b for a, b in zip(x, y))


def hbar_vec(x, y):
    re2, im2 = ip8(x, y), ip8(x, J_vec(y))
    return ((re2 // 2) + (im2 // 2)) % 2


n_ident = 0
for o1 in range(1, 16):
    if o1 in fixed_nz:
        continue
    o2, o3 = sig_label(o1), sig_label(sig_label(o1))
    if len({o1, o2, o3}) != 3:
        continue
    r1, r2, r3 = REPS[o1], REPS[o2], REPS[o3]
    for anc in fixed_nz:
        ra = REPS[anc]
        bits_map = {}
        for bits in itertools.product((0, 1), repeat=4):
            v = tuple(bits[0] * a + bits[1] * b + bits[2] * c
                      + bits[3] * d
                      for a, b, c, d in zip(r1, r2, r3, ra))
            bits_map[label_of(v)] = bits
        if len(bits_map) != 16:
            continue
        ok = all(hbar_vec(REPS[lx], REPS[ly])
                 == HB[WIDX[bits_map[lx]]][WIDX[bits_map[ly]]]
                 for lx in range(16) for ly in range(16))
        ok = ok and all(
            WIDX[bits_map[sig_label(li)]]
            == SIGP[WIDX[bits_map[li]]] for li in range(16))
        if ok:
            n_ident += 1
centr = sum(1 for p in SP
            if tuple(p[SIGP[i]] for i in range(16))
            == tuple(SIGP[p[i]] for i in range(16)))
check("U6.2 identifications (V, hbar, sigma) -> (F_2^4, hb, sig): "
      "%d found == %d = |centralizer of sigma in Sp(4,2)| (the "
      "lattice fixes the carrier triple up to exactly the "
      "centralizer torsor)" % (n_ident, centr),
      n_ident == centr and n_ident > 0, kill="K1")

STAB_QS = [k for k in STAB_Q
           if tuple(SP[k][SIGP[SP_INV[k][i]]] for i in range(16))
           == SIGP]
compl_orbit = set()
point_stab = 0
for k in STAB_QS:
    p, ip = SP[k], SP_INV[k]
    base = [IOTA_CANON[ip[i]] for i in range(16)]
    for pi in S5:
        i2 = tuple(tuple(x[pi[s]] for s in range(5)) for x in base)
        compl_orbit.add(i2)
        if i2 == IOTA_CANON:
            point_stab += 1
n_compl = len(iot_eq_canon)
single_compl_orbit = compl_orbit == set(iot_eq_canon)
check("U6.3 completions over the fixed carrier triple: q-completions "
      "4 -> 2 -> 1 (selector, U1.2); iota-completions %d, forming "
      "%s single orbit under Stab(q*, sigma) x S_5 "
      "(|Stab_Sp(q*, sigma)| = %d)"
      % (n_compl, "a" if single_compl_orbit else "NOT a",
         len(STAB_QS)), single_compl_orbit, kill="K2")
strict_universal = (n_compl == 1 and N_AUT == 1)
check("U6.4 TYPED: strict initial/terminal (unique up to UNIQUE "
      "isomorphism) %s: completions = %d, |Aut(C_fin)| = %d -- the "
      "measured reading is 'unique up to isomorphism, with %d "
      "automorphism arrows' %s"
      % ("HOLDS" if strict_universal else "FAILS",
         n_compl, N_AUT, N_AUT,
         "(a one-object groupoid, not a terminal object)"
         if not strict_universal else ""),
      True, "measured, typed either way", kill=None)

# ======================================================================
section("C: controls (must fire)")
# ======================================================================
DOT = form_table(tuple(tuple(1 if i == j else 0 for j in range(4))
                       for i in range(4)))
refs_dot = refinements_of_table(DOT)
check("C1 FIRES: the non-alternating dot form admits %d == 0 "
      "refinements -- admissible set EMPTY (dropping alternation "
      "kills the object)" % len(refs_dot), len(refs_dot) == 0,
      kill="K4")

arf0_iotas = max(len(iotas_for(q)) for q in ARF0)
check("C2 FIRES: WRONG ARF: every Arf-0 refinement admits 0 "
      "admissible iota (max over the 10: %d) -- the C_even(5) "
      "weight census (1,10,5) forces Arf 1; admissible set EMPTY"
      % arf0_iotas, arf0_iotas == 0, kill="K4")

cand_nosig = [q for q in REFS if q[LA] == 1 and q[LF] == 0]
bare_pairs = set()
for q in ARF1:
    for iot in iotas_for(q):
        bare_pairs.add((q, iot))
pair_orbit = set()
for k in range(720):
    p, ip = SP[k], SP_INV[k]
    q2 = tuple(QSTAR[ip[i]] for i in range(16))
    base = [IOTA_CANON[ip[i]] for i in range(16)]
    for pi in S5:
        i2 = tuple(tuple(x[pi[s]] for s in range(5)) for x in base)
        pair_orbit.add((q2, i2))
pair_single = pair_orbit == bare_pairs
check("C3 FIRES: NO SIGMA (trivial deck): the pinned completion "
      "census becomes 16 -> %d -> %d != 1 (canonicity of q* DIES "
      "without the family 3-cycle)"
      % (len([q for q in REFS if q[LA] == 1]), len(cand_nosig)),
      len(cand_nosig) == 4 and len(cand_nosig) != 1, kill="K4")
check("C3' TYPED side measurement (report, not a fire rule): the "
      "bare (q, iota) pair set without sigma has %d elements and "
      "forms %s single Sp x S_5 orbit -- sigma is load-bearing for "
      "CANONICITY (the selector), %s for the bare orbit count"
      % (len(bare_pairs),
         "a" if pair_single else "NOT a",
         "not" if pair_single else "AND"),
      True, "measured, typed either way", kill=None)

# ======================================================================
section("VERDICT")
# ======================================================================
n_pass = sum(1 for _, ok in CHECKS if ok)
n_tot = len(CHECKS)
if KILLS or N_ADM == 0:
    VERDICT = "CFIN-UNDERDETERMINED"
elif n_orbits == 1:
    VERDICT = "CFIN-UNIQUE"
else:
    VERDICT = "CFIN-ORBITS-%d" % n_orbits
print("%d/%d checks passed" % (n_pass, n_tot))
print("VERDICT: %s%s" % (VERDICT,
                         "  (|Aut(C_fin)| = %d ~= %s)" % (N_AUT, iso)
                         if VERDICT == "CFIN-UNIQUE" else ""))
print("\nCENSUS SUMMARY (report only -- promotion separate):")
print("  * admissible tuples at omega = hbar: %d; G-orbits: %d"
      % (N_ADM, n_orbits))
print("  * |Aut(C_fin)| = %d ~= %s; element orders %s"
      % (N_AUT, iso, orders))
print("  * Gaussian diagram: carrier triple fixed up to %d "
      "identifications (= centralizer); q* unique via selector; "
      "%d iota-completions in %s orbit"
      % (n_ident, n_compl, "one" if single_compl_orbit else ">1"))
print("  * strict terminal-object reading: %s"
      % ("HOLDS" if strict_universal else
         "FAILS (unique up to NON-unique isomorphism)"))
print("Runtime: %.1f s" % (time.time() - T0))
print("ALL CHECKS PASSED" if n_pass == n_tot
      else "CHECKS FAILED: %d" % (n_tot - n_pass))
raise SystemExit(0 if (n_pass == n_tot
                       and VERDICT.startswith("CFIN-UNIQUE")) else 1)
'''

_PF_RE = re.compile(r"^\s*\[(PASS|FAIL)\]\s+(\S+)", re.M)

_VD_RE = re.compile(r"VERDICT[:\]]")

def _probe_file(name):
    cand = os.path.abspath(os.path.join(
        _HERE, os.pardir, "experiments", "tfpt-discovery", name + ".py"))
    return cand if os.path.isfile(cand) else None

def _census(out):
    marks = _PF_RE.findall(out)
    fails = sorted({tok for st, tok in marks if st == "FAIL"})
    verdict = ""
    for line in out.splitlines():
        if _VD_RE.search(line):
            verdict = line.strip()
    return len(marks), fails, verdict

def _exec_probe(name, src, run_entry=True):
    """Execute one embedded frozen probe source BYTE-EXACT in its own
    module namespace, registered in sys.modules under the probe's
    canonical import name; capture and re-emit stdout; return
    (stdout, exit_code, byte_equal_to_source_file_or_None)."""
    if src[:1] == "\n":
        src = src[1:]
    path = _probe_file(name)
    same = None
    if path is not None:
        with open(path, encoding="utf-8") as fh:
            same = (fh.read() == src)
    fname = path or os.path.abspath(__file__)
    mod = types.ModuleType(name)
    mod.__file__ = fname
    sys.modules[name] = mod
    buf = io.StringIO()
    code = 0
    with contextlib.redirect_stdout(buf):
        try:
            exec(compile(src, fname, "exec"), mod.__dict__)
            entry = mod.__dict__.get("main") or mod.__dict__.get("run")
            if run_entry and callable(entry):
                rc = entry()
                code = 0 if rc is None else int(rc)
        except SystemExit as exc:
            code = 0 if exc.code is None else int(exc.code)
        except Exception:                            # regression guard
            import traceback
            traceback.print_exc(file=sys.stdout)
            code = 99
    out = buf.getvalue()
    sys.stdout.write(out)
    sys.stdout.flush()
    return out, code, same

def _gate(name, out, code, same, exp_n, exp_fails, exp_verdict,
          exp_code, gates):
    n, fails, verdict = _census(out)
    ok = (n == exp_n and fails == list(exp_fails)
          and exp_verdict in verdict and code == exp_code
          and same is not False)
    gates.append(ok)
    prov = ("byte-exact vs experiments source" if same is True else
            "embedded copy (source file not present)" if same is None
            else "SOURCE MISMATCH")
    print("\n[%s] PATTERN GATE %s: %d checks (exp %d) | FAILs %s "
          "(exp %s) | exit %d (exp %d) | %s\n      verdict line: %s"
          % ("PASS" if ok else "FAIL", name, n, exp_n,
             ",".join(fails) if fails else "none",
             ",".join(exp_fails) if exp_fails else "none",
             code, exp_code, prov, verdict), flush=True)
    return ok

_PLAN = (
    ("cfin_uniqueness_probe", _SRC_CFIN, 22, (), "CFIN-UNIQUE", 0),
)

def run():
    gates = []
    for name, src, exp_n, exp_fails, exp_verdict, exp_code in _PLAN:
        out, code, same = _exec_probe(name, src)
        _gate(name, out, code, same, exp_n, exp_fails, exp_verdict, exp_code, gates)
    print("v849 compiler uniqueness: %d/%d pattern gates passed" % (sum(gates), len(gates)))
    return 0 if all(gates) else 1

if __name__ == "__main__":
    raise SystemExit(run())
Recorded source fingerprint

4cd1127dced960e9a72cafed0c77c57d866f714b4898a5c74a5f9ecbd40588a5

Paper sections that cite this module

What this source index covers

The public registry includes 116 document views, 6,711 section records and 986 physical LaTeX sources. Active publication documents, research material and archives retain their recorded roles. This is source navigation, not a new semantic review or proof run.