# script used to create shas.html from allbsd.* files:
# we do not call the output file "shas.html" so we can compare the new
# version with the old
HTML_FILENAME = "newshas.html"
MAX_RANK = 4
SHA_LIST = range(2,35)+[37,41,43,47,50,75]
from sage.all import isqrt
def make_rankshatable(nmax=30, verbose=False):
total_tab = {}
rank_tab = [{} for r in range(MAX_RANK+1)]
range_tab = [{} for n in range(nmax)]
total = 0
rank_total = [0 for r in range(MAX_RANK+1)]
range_total = [0 for n in range(nmax)]
for n in range(nmax):
infilename = "allbigsha/allbigsha."+str(n)+"0000-"+str(n)+"9999"
if verbose:
print("processing "+infilename)
infile = open(infilename)
for L in infile.readlines():
N, cl, num, ainvs, r, t, S = L.split()
r = int(r)
S = int(S)
s = int(isqrt(S))
total_tab[s] = total_tab.get(s,0)+1
rank_tab[r][s] = rank_tab[r].get(s,0)+1
range_tab[n][s] = range_tab[n].get(s,0)+1
total +=1
rank_total[r] +=1
range_total[n] +=1
infile.close()
if verbose:
print("Totals for range {}0000-{}9999: {} (total {})".format(n,n,range_tab[n],range_total))
if verbose:
print
print("Totals for all ranks: {} (total {})".format(total_tab,total))
print
outfilename = HTML_FILENAME
outfile = open(outfilename, mode='w')
# header info for html file
outfile.write('\n')
outfile.write("\n")
outfile.write("\n")
outfile.write('\n')
outfile.write("Analytic Shas\n")
outfile.write("\n")
outfile.write("\n")
#
# First table: by ranges, not subdivided by rank
#
outfile.write("\n")
outfile.write("Number of curves with nontrivial (analytic) Sha, by conductor range\n")
outfile.write("
\n")
outfile.write("\n")
outfile.write("Details in files allbigsha.00000-09999 etc.\n")
outfile.write("
\n")
# Table head
outfile.write("\n")
outfile.write("\n")
outfile.write("\n")
outfile.write("\n")
outfile.write("\n")
outfile.write("\n")
outfile.write("\n")
outfile.write("\n")
outfile.write("\n")
outfile.write("\n")
outfile.write('\n')
outfile.write("| N | all>1")
for s in SHA_LIST:
outfile.write(" | %s2"%s)
outfile.write("\n |
\n")
outfile.write("\n")
# Table foot
outfile.write("\n")
outfile.write("\n")
outfile.write("| 1-%s9999 | \n"%str(nmax-1))
outfile.write("%s | \n"%total)
for s in SHA_LIST:
outfile.write("%s | \n"%total_tab.get(s,0))
outfile.write("\n
\n")
outfile.write('\n')
outfile.write("| N | all>1")
for s in SHA_LIST:
outfile.write(" | %s2"%s)
outfile.write("\n |
\n")
outfile.write("\n")
# Table body
outfile.write("\n")
for n in range(nmax):
outfile.write("\n")
if n==0:
outfile.write("| 1-9999 | \n")
else:
outfile.write("%s0000-%s9999 | \n"%(str(n),str(n)))
outfile.write("%s | \n"%range_total[n])
for s in SHA_LIST:
outfile.write("%s | \n"%range_tab[n].get(s,' '))
outfile.write("
\n")
outfile.write("\n")
outfile.write("
\n")
outfile.write("
\n")
#
# Second table: by ranks, not subdivided by range
#
outfile.write("\n")
outfile.write("Number of curves with nontrivial (analytic) Sha, by rank\n")
outfile.write("
\n")
outfile.write("\n")
outfile.write("\n")
outfile.write("\n")
outfile.write("\n")
outfile.write("\n")
outfile.write("\n")
outfile.write("\n")
outfile.write("\n")
outfile.write("\n")
outfile.write("\n")
outfile.write('\n')
outfile.write("| | all>1")
for s in SHA_LIST:
outfile.write(" | %s2"%s)
outfile.write("\n |
\n")
outfile.write("\n")
# Table foot
outfile.write("\n")
outfile.write("\n")
outfile.write("| all ranks | \n")
outfile.write("%s | \n"%total)
for s in SHA_LIST:
outfile.write("%s | \n"%total_tab.get(s,0))
outfile.write("\n
\n")
outfile.write('\n')
outfile.write("| | all>1")
for s in SHA_LIST:
outfile.write(" | %s2"%s)
outfile.write("\n |
\n")
outfile.write("\n")
# Table body
outfile.write("\n")
for r in range(MAX_RANK+1):
if rank_total[r]>0:
outfile.write("\n")
outfile.write("| r=%s | \n"%str(r))
outfile.write("%s | \n"%rank_total[r])
for s in SHA_LIST:
outfile.write("%s | \n"%rank_tab[r].get(s,' '))
outfile.write("
\n")
outfile.write("\n")
outfile.write("
\n")
outfile.write("
\n")
# footer info for html file
outfile.write("\n")
outfile.write(' 
\n')
outfile.write("
\n")
outfile.write("\n")
outfile.write("\n")
if verbose:
for r in range(MAX_RANK+1):
if rank_total[r]>0:
print("Totals for rank {}: {} (total {})".format(r,rank_tab[r],rank_total[r]))