writeCMtoDB.m

AttachSpec("spec");
SetDebugOnError(true);


for filename in ShimDBFilenames() do
  s:=ShimDBRead(filename : version:=2);
  W:=s`ShimAtkinLehner;
  if #W eq 2 then
    if s`ShimGenus ge 1 and Type(s`ShimRationalPoints) eq SetEnum
      and s`ShimRationalPoints ne {} then
        D:=s`ShimDiscriminant; N:=s`ShimLevel;
      //if D ne 6 or [D,N,W[2]] notin [[21,2,42],[26,3,78], [26,3,39],
        //[6,17,51],[6,19,38],[6,19,6],[6,29,87],[6,31,62]] then
        if N eq 1 and [D,W[2]] notin [[35,5],[51,17]] then //and [D,W[2]] in [[134,67],[159,53],[206,103],[38,19],[51,17],[62,31],[69,23],[86,43],[87,29],[93,93],[94,47]] then
        D; N; W;
        s`ShimRationalPoints;
        cmpr,cmun:=RationalCMPoints(D,N,W[2]);
        print "CM proven:"; cmpr;
        print "CM unproven:"; cmun;
        print "=====================";
      end if;
    end if;
  end if;
end for;



for i in [1..#GYList()] do
  item:=GYList()[i];
  D:=item[1];
  N:=item[2];
  for W in AllAtkinLehners(D,N) do
    if #W eq 2 then
      s:=ShimDBRecord(D,N : version:=2);
      if s`ShimGenus ge 1 and Type(s`ShimRationalPoints) eq SetEnum
        and s`ShimRationalPoints ne {} then
        s`ShimRationalPoints;
      end if;
    end if;
  end for;
end for;

        //if W ne [1,6,13,78] then
          printf "%o %o %o\n", disc,level,W;
          MakeShimDBObject(disc,level,W);
          attr:=ShimDBRecord(disc,level,W : version:=1);
        //end if;
        //C:=ShimuraCurveQuotient(disc,level,W);
        //printf "Genus = %o\n", Genus(C);
        //DefiningEquations(C);
      end for;
    end if;
  //end if;
end for;

ProfilePrintByTotalTime(:Max:=20);


command := Sprintf("ls %o", "ShimDB-v1");
ls := Pipe(command, "");
filenames := Split(ls, "\n");

for filename in filenames do
  file:=Sprintf("ShimDB-v1/%o",filename);
  file;
  FP:=Read(file);
  attr:=eval FP;
  if attr`ShimAtkinLehner ne [1] then
  //new_attributes:= [ <"ShimRationalPoints","\"{}\"">, <"ShimTest", "\"NA\""> ];
    X:=attr`ShimModel;
    ShimRationalPoints, ShimPointsProvedCorrect, ShimPointsNotes := RationalPointsAnyGenus(X);
    ShimPointsEverywhereLocally:= HasAdelicPointsAnyGenus(X);
    if Type(ShimRationalPoints) eq SetEnum then
      points:=Set([ Eltseq(P) : P in ShimRationalPoints]);
    elif Type(ShimRationalPoints) eq MonStgElt then
      points:=Sprintf("\"%o\"",ShimRationalPoints);
    else
      assert Type(ShimRationalPoints) eq BoolElt;
      assert Genus(X) eq 0;
      points:=ShimRationalPoints;
    end if;

    //if Genus(X) ge 2 then
      //ShimRationalPoints:= Sprintf("s`ShimModel!%o", ShimRationalPoints);
    new_attributes:=
    [
    <"ShimRationalPoints", Sprint(points)>,
    <"ShimPointsProvedCorrect", Sprint(ShimPointsProvedCorrect)>,
    <"ShimPointsNotes", Sprintf("\"%o\"",ShimPointsNotes)>,
    <"ShimPointsEverywhereLocally", Sprint(ShimPointsEverywhereLocally)>
    ];
    ShimAddAttributes(filename, new_attributes);
  end if;
end for;



command := Sprintf("ls %o", "ShimDB-v2");
ls := Pipe(command, "");
filenames := Split(ls, "\n");

points_unproven:=[];
Hasse_violations:=[];
for filename in filenames do
  file:=Sprintf("ShimDB-v2/%o",filename);
  filename;
  FP:=Read(file);
  attr:=eval FP;
  if attr`ShimGenus ne 0 and Type(attr`ShimRationalPoints) eq SetEnum
    and attr`ShimRationalPoints eq {} and attr`ShimPointsEverywhereLocally then
      Append(~Hasse_violations, filename);
  end if;

  if attr`ShimPointsProvedCorrect eq false then
    Append(~points_unproven, filename);
  end if;

end for;

points_unproven;
Hasse_violations;

for filename in filenames do
  file:=Sprintf("ShimDB-v2/%o",filename);
  filename;
  FP:=Read(file);
  s:=eval FP;
  proj:=s`ShimProjectionEquations;
  quotient_points:=s`ShimRationalPoints;
  W:=s`ShimAtkinLehner;
  discriminant:=s`ShimDiscriminant;
  PointsRepresentatingPQMSurface(proj,quotient_points,W,discriminant);
end for;




/* Hasse_violations:=

[ Shim-X(119,1)-g4-[1,7].m, Shim-X(14,5)-g1-[1,5].m, Shim-X(39,2)-g3-[1,78].m, Shim-X(6,29)-g2-[1,6].m, Shim-X(6,37)-g2-[1,3].m, Shim-X(87,1)-g2-[1,3].m, Shim-X(93,1)-g3-[1,3].m ]

points_unproven:=
[ Shim-X(10,19)-g2-[1,190].m, Shim-X(10,23)-g3-[1,2,23,46].m, Shim-X(119,1)-g5-[1,17].m, Shim-X(134,1)-g3-[1,2].m, Shim-X(134,1)-g3-[1,67].m, Shim-X(159,1)-g5-[1,53].m, Shim-X(206,1)-g5-[1,103].m, Shim-X(87,1)-g3-[1,29].m, Shim-X(93,1)-g2-[1,93].m ]*/