# script used to create shas.html from allbsd.* files:
from sage.all import isqrt
# 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]
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 = ''.join(["allbigsha/allbigsha.", str(n), "0000-", str(n), "9999"])
if verbose:
print("processing " + infilename)
infile = open(infilename)
for L in infile.readlines():
_, _, _, _, r, _, 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("\nTotals for all ranks: {} (total {})\n".format(total_tab, total))
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]))