1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
|
# Script to fix Lr1 and Sha values when we had forgotten to divide by r! (361725<=N<=370000)
squares = [n*n for n in srange(1,100)]
def fix_bsd(infilename):
"""
If rank is >1 check to see if Lr1 and Sha need to be divided by r!
Input line fields:
conductor iso number ainvs rank torsion tamagawa_product real_period special_value regulator sha_an
Sample input line:
11 a 1 [0,-1,1,-10,-20] 0 5 5 1.2692093042795534217 0.25384186085591068434 1 1.00000000000000000000
"""
infile = file(infilename,mode='r')
outfile = file(infilename+".fixed", mode='w')
for line in file(infilename).readlines():
data = line.split()
r=int(data[4])
s = round(float(data[10]))
if not s in squares:
print("fixing %s%s%s: rank = %s, Sha = %s" % (data[0],data[1],data[2],data[4],data[10]))
rfact = factorial(r)
data[10]=str(RR(data[10])/rfact)
data[8]=str(RR(data[8])/rfact)
print("fixed, new Sha = %s" % data[10])
line = ' '.join(data)+'\n'
outfile.write(line)
infile.close()
outfile.close()
|
