forked from ciaran-schembri/Shimura
-
Notifications
You must be signed in to change notification settings - Fork 0
/
ShimDBaccess.m
89 lines (66 loc) · 2.6 KB
/
ShimDBaccess.m
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
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
//Merge Test
AttachSpec("spec");
SetLogFile("trying.log");
command := Sprintf("ls %o", "ShimDB");
ls := Pipe(command, "");
filenames := Split(ls, "\n");
for filename in filenames do
file:=Sprintf("ShimDB/%o",filename);
FP:=Read(file);
attr:=eval FP;
if attr`ShimGenus ge 2 then
filename;
X:=attr`ShimModel;
/*XG1:=GenusOneModel(X);
locally_sol:=IsLocallySolvable(XG1);
locally_sol;
if locally_sol eq true then
Points(X : Bound:=10000);
end if;
print "==================";*/
print X;
try
Y := SimplifiedModel(X);
RatAttempt(Y);
catch e "double check";
end try;
end if;
end for;
command := Sprintf("ls %o", "ShimDB");
ls := Pipe(command, "");
filenames := Split(ls, "\n");
atkinlehners:= [ eval Split(Split(file,"w")[2],">")[1] : file in filenames ];
ParallelSort(~atkinlehners, ~filenames);
for filename in filenames do
file:=Sprintf("ShimDB/%o",filename);
FP:=Read(file);
attr:=eval FP;
if attr`ShimGenus gt 1 and "s`ShimAtkinLehner" in FP then
if attr`ShimHasAdelicPoints eq true then
printf "& & $w_{%o}$ & %o & %o & %o & %o & %o %o & %o & $y^2 = %o$ \\\\ \n",
attr`ShimAtkinLehner, attr`ShimGenus, attr`ShimAnalyticRank, attr`ShimHasAdelicPoints,
attr`ShimRepresentsSurface, attr`ShimChabauty, attr`ShimPointSearch, attr`CMpoints, attr`ShimEquation;
end if;
end if;
end for;
Pipe("ls ShimDB", "");
filenames:= [ "../../Dropbox (Dartmouth College)/Shimura-Chabauty/code/%o"]
Sprintf("ls ../../Dropbox (Dartmouth College)/Shimura-Chabauty/code/%o", "ShimDatabase");
lst:=Pipe("ls ShimDB", "");
load "dir";
load "../../Dropbox (Dartmouth College)/Shimura-Chabauty/code/ShimDatabase/Shim-X(26,1)-<w2>-g1.m";
Open(filename);
if s`ShimGenus eq 0 then
printf "& & $w_{%o}$ & %o & %o & %o & %o & %o & %o & $%o$ \\\\ \n",
s`ShimAtkinLehner, s`ShimGenus, s`ShimInfinitelyManyPoints, s`ShimHasAdelicPoints,
s`ShimRepresentsSurface, "NA", s`ShimCMPoints, s`ShimEquation;
elif s`ShimGenus eq 1 then
printf "& & $w_{%o}$ & %o & $\\Z^%o + %o$ & %o & %o & %o & %o & $y^2 = %o$ \\\\ \n",
s`ShimAtkinLehner, s`ShimGenus, s`MordellWeil[1], s`MordellWeil[2], s`ShimHasAdelicPoints,
s`ShimRepresentsSurface, "NA", s`ShimCMPoints, s`ShimEquation;
else
printf "& & $w_{%o}$ & %o & %o & %o & %o & %o %o & %o & $y^2 = %o$ \\\\ \n",
s`ShimAtkinLehner, s`ShimGenus, s`ShimRankBounds, s`ShimHasAdelicPoints,
s`ShimRepresentsSurface, s`ShimChabauty, s`ShimPointSearch, s`CMpoints, s`ShimEquation;
end if;
Open("../../Dropbox (Dartmouth College)/Shimura-Chabauty/code/ShimDatabase/Shim-X(26,1)-<w2>-g1.m",p);