| |
| """Measure recovered original-asset outputs against conventional geometry formulas.""" |
| import hashlib,json,math,struct |
| from pathlib import Path |
| from analyze_qvm_acceleration import word,number |
| from analyze_qvm import write_json |
|
|
| def f(v):return struct.unpack('<f',struct.pack('<f',v))[0] |
| def conventional_basis(t): |
| a=list(map(number,t['sample_words64_108']));left=[a[3:6]+[0],a[6:9]+[0],a[9:12]+[0],[0,0,0,1]];right=[[number(v) for v in t['initial_record'][3+4*r:7+4*r]] for r in range(4)];out=[] |
| for row in left: |
| for col in range(4): |
| value=f(row[0]*right[0][col]) |
| for k in range(1,4):value=f(value+f(row[k]*right[k][col])) |
| out.append(word(value)) |
| return out |
|
|
| def main(): |
| points=[];matrices=[];sources={} |
| for stem in ['md3-original-leaf','md3-original-hierarchy','md3-original-model-hit']: |
| path=Path('docs/research')/(stem+'.json');data=path.read_bytes();sources[stem]=hashlib.sha256(data).hexdigest();rows=json.loads(data)['probes'] |
| for row_index,p in enumerate(rows): |
| leafs=[p] if stem=='md3-original-leaf' else [r['effect'] for r in (p['leaf_reports'] if stem=='md3-original-hierarchy' else p['hierarchies'][0]['leaf_reports'])] |
| for leaf_index,r in enumerate(leafs): |
| origin=list(map(number,r['input']['origin']));direction=list(map(number,r['input']['direction'])) |
| for record in range(r['before'][16],r['after'][16]): |
| values=r['after'][17+5*record:22+5*record];t=number(values[4]);expected=[word(f(x+f(t*y))) for x,y in zip(origin,direction)];distance=math.sqrt(sum((number(x)-number(y))**2 for x,y in zip(values[1:4],expected)));points.append({'source':stem,'row':row_index,'leaf':leaf_index,'record':record,'category':values[0],'parameter':values[4],'source_point':values[1:4],'conventional_ray_point':expected,'distance':distance}) |
| if stem=='md3-original-hierarchy': |
| for index,t in enumerate(p['transforms']): |
| actual=t['final_record'][3:19];expected=conventional_basis(t);matrices.append({'source':stem,'row':row_index,'transform':index,'source_basis':actual,'conventional_basis':expected,'different_indices':[j for j,(a,b) in enumerate(zip(actual,expected)) if a!=b]}) |
| write_json(Path('docs/research/md3-geometry-audit.json'),{'source_report_hashes':sources,'points':points,'matrices':matrices,'summary':{'point_comparisons':len(points),'noncoincident_points':sum(p['distance']>1e-5 for p in points),'largest_point_distance':max(p['distance'] for p in points),'fourth_row_differences':sum(bool(p['different_indices']) for p in matrices)},'limits':'Numeric comparison of recovered QVM outputs on original assets against conventional formulas. Neither formula is executed in the original game here. Differences are evidence requiring native/reference validation, not authorization to change source behavior.'}) |
| if __name__=='__main__':main() |
|
|
