s1003.py

from research import *
s=begin(1003,'Exact rational field normal form','How can the full File46 head field be retained when some images are nonintegral?',{'anchor':14006,'source_mesh':10,'indices':[92,89,46,43,0]},['File46','C931 technical','C1002'])
p=s['inputs'];heads=[p['anchor']+p['source_mesh']*n for n in p['indices']];mesh=E*p['source_mesh'];images=[p['anchor']+mesh*n for n in p['indices']]
a=artifact('model/key_full_field_normal_form.json',json.dumps(clean({'anchor':p['anchor'],'indices':p['indices'],'source_mesh':10,'image_mesh':mesh,'source_heads':heads,'images':images,'status':'exact comparison mesh, not a new historical calendar'}),indent=2)+'\n')
finish(s,{'field':a,'image_mesh':mesh,'integral_indices':[n for n,y in zip(p['indices'],images) if y.denominator==1]},'The entire labelled index pattern survives on the exact mesh250/23. Only indices divisible23 give integer images; keeping fractions preserves the comparison without repairing dates.','Test the complete componentwise calendar assignment field.',{'mesh':mesh==F(250,23),'integral_subset':[n for n,y in zip(p['indices'],images) if y.denominator==1]==[92,46,0],'source_heads':heads==[14926,14896,14466,14436,14006]})