File indexing completed on 2024-11-25 02:30:05
0001 import sys
0002 if len(sys.argv) !=3 :
0003 print("usage:" + sys.argv[0] + " bla bla")
0004 exit
0005
0006 mm=0.
0007 n=0
0008 ave=0
0009 tot=0
0010 ln=0
0011 ml=0
0012 file1 = open(sys.argv[1])
0013 file2 = open(sys.argv[2])
0014 for l1 in file1.readlines() :
0015 ln+=1
0016 l2 = file2.readline()
0017 a1 =l1.split('(')
0018 a2 = l2.split('(')
0019 if len(a1)<2 :
0020 continue
0021 p1 = a1[1].split(',')
0022 p1[2] = p1[2].rstrip(') \n')
0023 p2 = a2[1].split(',')
0024 p2[2] = p2[2].rstrip(') \n')
0025 f1 = [float(i) for i in p1]
0026 f2 = [float(i) for i in p2]
0027 diff = []
0028 tot+=1
0029 for i in range(len(f1)):
0030 diff.append(f1[i]-f2[i])
0031 for x in diff:
0032 if (abs(x)>mm):
0033 mm=x
0034 ml=ln
0035 if (x!=0) :
0036 n+=1
0037 ave+=abs(x)
0038 print(tot, mm, n, ave/n)
0039
0040 ml = ml-5
0041 file1.seek(0)
0042 file2.seek(0)
0043 l1 = file1.readlines()
0044 l2 = file2.readlines()
0045 for i in range(10) :
0046 print(l1[ml+i].rstrip('\n'))
0047 print(l2[ml+i].rstrip('\n'))
