# script used to create table.html from allbsd.* files:
#countrank.awk:
# cat curves.*0000-*9999 |
# gawk -v FIRST=$1 -v LAST=$2 'BEGIN{printf("Curve numbers by rank in the range %d...%d:\nrank:\t0\t1\t2\t3\t4\n",FIRST,LAST);}\
# ($1>=FIRST)&&($1<=LAST){r[$5]+=1;rt+=1;}\
# END {printf("number:\t%d\t%d\t%d\t%d\t%d\nTotal number: %d\n",
# r[0],r[1],r[2],r[3],r[4],rt);printf("\n| %d-%d | \n%d | \n",FIRST,LAST,rt);for(i=0;i<5;i++){printf("%d | \n",r[i]);}printf("
\n");}'
#countcurves.awk:
# cat allcurves.*0000-*9999 |
# gawk -v FIRST=$1 -v LAST=$2 'BEGIN{printf("Numbers of isogeny and isomorphism classes in the range %d...%d:\n",FIRST,LAST);ncu=0;ncl=0;}\
# ($1>=FIRST)&&($1<=LAST){;ncu+=1;if($3==1){ncl+=1;}}\
# END {printf("\n| %d-%d | \n%d | \n%d | \n
\n",FIRST,LAST,ncl,ncu);}'
# we do not call the output file "table.html" so we can compare the new
# version with the old
HTML_FILENAME = "newtable.html"
MAX_RANK = 4
def make_table(nmax=30, verbose=False):
total_tab = {}
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)]
range_total_all = [0 for n in range(nmax)]
total_all = 0
for n in range(nmax):
infilename = ''.join(["allcurves/allcurves.", str(n), "0000-", str(n), "9999"])
range_total_all[n] = len(open(infilename).readlines())
total_all += range_total_all[n]
for n in range(nmax):
infilename = ''.join(["curves/curves.", str(n), "0000-", str(n), "9999"])
if verbose:
print("processing {}".format(infilename))
infile = open(infilename)
for L in infile.readlines():
_, _, _, _, r, _, _ = L.split()
r = int(r)
total_tab[r] = total_tab.get(r, 0) + 1
range_tab[n][r] = range_tab[n].get(r, 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: {} (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("Elliptic Curves\n")
outfile.write("\n")
outfile.write("\n")
# Table 1 header
outfile.write("\n")
outfile.write("Number of isogeny classes of curves of conductor N < %s0000, sorted by rank\n"%nmax)
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 | r=0 | r=1 | r=2 | r=3 | r=4\n")
outfile.write(" |
\n")
outfile.write("\n")
#Table 1 footer
outfile.write("\n")
outfile.write("\n")
outfile.write("| 1-%s9999 | \n" % (nmax - 1))
outfile.write("%s | \n" % total)
for r in range(MAX_RANK+1):
outfile.write("%s | \n" % total_tab[r])
outfile.write("
\n")
outfile.write("\n")
#Table 1 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" % (n, n))
outfile.write("%s | \n" % range_total[n])
for r in range(MAX_RANK + 1):
outfile.write("%s | \n" % range_tab[n].get(r, 0))
outfile.write("
\n")
outfile.write("\n")
outfile.write("
\n")
outfile.write("\n")
outfile.write("
\n")
outfile.write("\n")
# Table 2 header
outfile.write("\n")
outfile.write("Total number of curves of conductor N < %s0000\n" % nmax)
outfile.write("
\n")
outfile.write("\n")
outfile.write("\n")
outfile.write("\n")
outfile.write("\n")
outfile.write("\n")
outfile.write("\n")
outfile.write("| N | \n")
outfile.write("# isogeny classes | \n")
outfile.write("# isomorphism classes | \n")
outfile.write("
\n")
outfile.write("\n")
# Table 2 footer
outfile.write("\n")
outfile.write("\n")
outfile.write("| 1-%s9999 | \n" % (nmax - 1))
outfile.write("%s | \n" % total)
outfile.write("%s | \n" % total_all)
outfile.write("
\n")
outfile.write("\n")
# Table 2 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" % (n, n))
outfile.write("%s | \n" % range_total[n])
outfile.write("%s | \n" % range_total_all[n])
outfile.write("
\n")
outfile.write("\n")
outfile.write("
\n")
# html footer
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.close()