validate.py

from pathlib import Path
from fractions import Fraction as F
import json,hashlib
import exact as X
HERE=Path(__file__).resolve().parent
d=json.loads((HERE/'diagnostics.json').read_text());p=json.loads((HERE/'source_packet.json').read_text());r=json.loads((HERE/'source_only_reproduction.json').read_text());checks={}
checks['sources_unchanged']=all(hashlib.sha256(Path(v['path']).read_bytes()).hexdigest()==v['sha256'] for v in p['sources'].values())
checks['questions_declared_and_bound']=d['question_file_sha256']==hashlib.sha256((HERE/'questions_declared.json').read_bytes()).hexdigest()
checks['old_pattern_values_reproduced_without_inputs']=r['derived_old_pattern_values']=={'Ishmael_Levi_birth_gap':161,'Levi_terminal_gap':161,'Covenant_Levi_span':299,'Covenant_Joseph_span':276}
checks['regular_all_expected']={k:r['regular_network'][k] for k in ['Ishmael_birth','Ishmael_death','Isaac_birth','Jacob_birth','Jacob_call','Jacob_death','Levi_birth','Levi_death','Joseph_birth','Joseph_death','terminal']}=={'Ishmael_birth':1865,'Ishmael_death':1728,'Isaac_birth':1851,'Jacob_birth':1791,'Jacob_call':1714,'Jacob_death':1644,'Levi_birth':1704,'Levi_death':1567,'Joseph_birth':1700,'Joseph_death':1590,'terminal':1406}
checks['root_original_literal'] = r['original_cumulative']==[2435,2260,2080,1933,1796,1663,1526]
checks['clutch_derived_from_sources']=r['derived_clutch']==14 and r['derived_half_clutch']==7 and not {'clutch','half','clutch14','half7'}&{x['name'] for x in p['source_local_parameters']}
expected_b=[2435,2260,2080,1933,1796,1663,1526]
checks['all_42_coordinate_values']=all(list(map(F,x['values'].values()))==[F(b)+x['phase_index']*F(7,2)-7*x['clutch_index'] for b in expected_b] for x in r['six_rails'])
M=[[F(x) for x in row] for row in d['claim_matrix']]
checks['all_nine_claims_hold']=all(x['source_value']==0 for x in d['claims'])
checks['five_null_vectors']=len(d['claim_nullspace'])==5 and all(all(sum(a*F(y[0]) for a,y in zip(row,v))==0 for row in M) for v in d['claim_nullspace'])
checks['single_claim_row_circuit']=len(d['claim_row_dependencies'])==1 and all(sum(M[i][j]*F(d['claim_row_dependencies'][0][i][0]) for i in range(9))==0 for j in range(13))
rows=d['claim_basis_indices'];cols=X.Matrix(M).rref()[1];minor=[[int(M[i][j]) for j in cols] for i in rows]
def det_mod(a,p=1000003):
 a=[[x%p for x in row] for row in a];det=1
 for j in range(len(a)):
  k=next((i for i in range(j,len(a)) if a[i][j]),None)
  if k is None:return 0
  if k!=j:a[k],a[j]=a[j],a[k];det=-det
  det=det*a[j][j]%p;iv=pow(a[j][j],p-2,p)
  for i in range(j+1,len(a)):
   f=a[i][j]*iv%p;a[i]=[(u-f*v)%p for u,v in zip(a[i],a[j])]
 return det%p
dm=det_mod(minor)
checks['claim_rank_eight_lower_bound_minor']=len(minor)==8 and dm!=0
checks['coordinate_incidence_unimodular']=abs(F(d['coordinate_map']['determinant']))==1
checks['clutch_observation_basis_unimodular']=abs(F(d['clutch_field']['observation_determinant']))==1
checks['clutch_33_dependencies_count']=len(d['clutch_field']['relation_matrix'])==33
FM=[]
for i in range(7):
 for j in range(2):
  for k in range(3):FM.append([int(z==i) for z in range(7)]+[j,-k])
checks['all_clutch_dependencies_exact']=all(all(sum(FM[i][j]*F(n[i][0]) for i in range(42))==0 for j in range(9)) for n in d['clutch_field']['relation_matrix'])
checks['phase_and_repeat_reductions']=d['clutch_field']['relation_half_equals_twice_phase']['rank']==8 and d['clutch_field']['repeated137_source_equality']['rank']==8 and d['clutch_field']['both_relations_rank']==7
join=d['endogenous_clutch_joins']
checks['forced_joins_zero']=join['Kohath_join_residual']==join['Levi_join_residual']==join['Jacob_join_equals_negative_first_147_condition']==0
checks['all_translated_nodes_Moses_independent']=join['translated_Moses_anchor_derivatives']==[0]*7 and join['Exodus_derivative_of_clutch']==0
checks['full_coupled_rank']=d['coupled_field']['rank']==d['coupled_field']['source_parameter_count']==19 and d['coupled_field']['new_scalar_clutch_parameter_count']==0
out={'status':'PASS' if all(checks.values()) else 'FAIL','checks':checks,'passed':sum(checks.values()),'total':len(checks),'rank8_minor_columns':list(cols),'rank8_minor_mod1000003':dm,'limits':'Validation of this preparation only; not numbered research or a statistical independence tally.'}
(HERE/'VALIDATION.json').write_text(json.dumps(out,indent=2)+'\n')
print(json.dumps({'status':out['status'],'passed':out['passed'],'total':out['total'],'failed':[k for k,v in checks.items() if not v]}))