source63 verify.py

from fractions import Fraction as Q
from pathlib import Path
import hashlib,json

ROOT=Path('/workspace/scratch/1b40da62dcbd')
OUT=ROOT/'c532_c631/prep'
sources={
 'File63':'project_sources/09-File_63.Scale_Neutral_480_483_490_Carrier-20260821-143254-.md',
 'File12':'project_sources/16-File_12.Calendrical_Physics-20260907-192303-.md',
 'File17':'project_sources/21-File_17.Prophetic_Time_Span_Anatomy-20260908-195734-.md',
 'C531_synthesis':'490d_How_Chronological_Families_Fit_Together_Draft_20260927.md',
}
inputs={
 'sources':{k:{'path':v,'sha256':hashlib.sha256((ROOT/v).read_bytes()).hexdigest()} for k,v in sources.items()},
 'operators':{'E':'25/23','P':'70/69','J':'300/299','phase_pair_ratio':'161/160'},
 'carrier':{'source':'File63 1.1–1.4','unit_pair':['40','161/4'],'multipliers':[3,10,12,30,300]},
 'dual_490':{'source':'File63 1.3 / File17 0.4','whole':'483','unexpanded_subset':'805/2','expanded_subset':'161/2','whole_P_output':'490','subset_E_output':'175/2','whole_E_output':'525'},
 'creation_cell':{'source':'File63 7.1–7.6','nodes':[14011,14008,14006,14004],'roles':['Day1','Day4','Day6','week_close'],'multipliers':[25,26,27],'carrier':483,'expected_rows':[[1936,1933,1931,1929],[1453,1450,1448,1446],[970,967,965,963]],'independent_roles':{'1936':'cumulative Jacob-death/Levi-birth upper','1933':'cumulative lower Levi branch','1929':'regular Jacob Bethel/Haran','1446':'Exodus','970':'precise Solomon accession','967':'ordinal Temple foundation'}},
 'flood_cells':{'source':'File63 8.1–8.3','envelope':[4836,4829],'core':[4834,4831],'carriers':[480,483,490],'multiplier':7,'expected_rows':[[1476,1469],[1455,1448],[1406,1399]],'core_carrier_multiple':6,'core_expected':[1936,1933]},
 'shem_jacob_square':{'source':'File63 8.5','heads':[5436,4836],'targets':[2006,1406],'horizontal':3430,'vertical':600},
 'birthlaw':{'source':'File63 4.1–4.4','nodes':[2166,2026,1366,1296,966],'roles':['regular Abraham birth','regular Isaac marriage','regular Joshua death','cumulative Joshua endpoint','cardinal Temple foundation'],'small_vector':[14,66,7,33],'scale':10,'exodus_partition':[720,480],'joshua_partition':[800,400]},
 'common_seeds':[{'source':'File63 9.7','seed':12558,'source_P':12740,'source_J':12600,'derived_E':13650},{'source':'File63 7.9','seed':143520,'source_P':145600,'source_J':144000,'derived_E':156000}],
 'actual_rounded_square':{'source':'File63 9.7; C531 3–4','actual':[14004,1446],'rounded':[14006,1406],'operator':'J','head_held':14004,'translation_BC':2},
 'phase_examples':{'source':'File63 7.6; File12 6','coordinate_rule':{'n':'Y','t':'Y-1/2'},'pairs':[[[1450,'t'],[1446,'n'],'7/2'],[[1449,'n'],[1446,'n'],'3'],[[970,'t'],[966,'n'],'7/2']]},
 'controls':['no primer work','no inverse-of-inverse','no new anchors','1486 excluded from current synthesis','generated nodes remain generated','phase and calendar units preserved']
}
E,P,J=Q(25,23),Q(70,69),Q(300,299)
checks=[]
def check(name,actual,expected):
 assert actual==expected,(name,actual,expected)
 checks.append({'name':name,'actual':str(actual),'expected':str(expected),'pass':True})
check('allocation_fraction',(P-1)/(E-1),Q(1,6))
check('dual_490_subset',Q(805,2)+E*Q(161,2),Q(490))
check('dual_490_whole',P*483,Q(490))
check('whole_E',E*483,Q(525))
check('day_volume_EP',525*336,490*360)
check('phase_then_P',Q(161,160)*P,Q(49,48))
check('phase_then_E',Q(161,160)*E,Q(35,32))
for m in inputs['carrier']['multipliers']:
 check('phase_excess_'+str(m),Q(161,4)*m-40*m,Q(m,4))
for k,row in zip(inputs['creation_cell']['multipliers'],inputs['creation_cell']['expected_rows']):
 actual=[x-k*483 for x in inputs['creation_cell']['nodes']]
 check('creation_cell_'+str(k),actual,row)
 check('creation_shape_'+str(k),[actual[i]-actual[i+1] for i in range(3)],[3,2,2])
for carrier,row in zip(inputs['flood_cells']['carriers'],inputs['flood_cells']['expected_rows']):
 check('flood_'+str(carrier),[x-7*carrier for x in inputs['flood_cells']['envelope']],row)
check('flood_core',[x-6*483 for x in inputs['flood_cells']['core']],[1936,1933])
check('shem_square_horizontal',[a-b for a,b in zip(inputs['shem_jacob_square']['heads'],inputs['shem_jacob_square']['targets'])],[3430,3430])
check('shem_square_vertical',[5436-4836,2006-1406],[600,600])
nodes=inputs['birthlaw']['nodes']
check('birthlaw_vector',[nodes[i]-nodes[i+1] for i in range(4)],[10*x for x in inputs['birthlaw']['small_vector']])
check('partition_transfer',[720+80,480-80],[800,400])
for row in inputs['common_seeds']:
 s=row['seed'];n=Q(s,897)
 check('common_seed_integral_'+str(s),n.denominator,1)
 for key,f in [('P',P),('J',J),('E',E)]:
  check('common_seed_'+str(s)+'_'+key,s*f,Q(row[('source_' if key!='E' else 'derived_')+key]))
 check('common_seed_volume_'+str(s),s*J*364,327600*n)
check('actual_rounded_span',J*(14004-1446),14006-1406)
check('actual_rounded_target',14004+J*(1446-14004)+2,Q(1406))
for a,b,span in inputs['phase_examples']['pairs']:
 qa=Q(a[0])-(Q(1,2) if a[1]=='t' else 0)
 qb=Q(b[0])-(Q(1,2) if b[1]=='t' else 0)
 check('phase_'+str(a)+'_'+str(b),qa-qb,Q(span))
(OUT/'source63_inputs.json').write_text(json.dumps(inputs,indent=2)+'\n')
(OUT/'source63_results.json').write_text(json.dumps({'status':'PASS','check_count':len(checks),'checks':checks},indent=2)+'\n')
print(json.dumps({'status':'PASS','check_count':len(checks),'files':['source63_inputs.json','source63_results.json']}))