Skip to content
Snippets Groups Projects
Commit 7acb3d40 authored by Gonzalo Tornaría's avatar Gonzalo Tornaría
Browse files

simplify test for 'hard' field for C353

parent 3eea1b7d
No related branches found
No related tags found
No related merge requests found
t0 := Cputime();
AttachSpec("mod.spec"); AttachSpec("mod.spec");
debug:=false; debug:=false;
...@@ -73,8 +75,8 @@ procedure solve(f) ...@@ -73,8 +75,8 @@ procedure solve(f)
// finding the embeding R1 in R2 takes a long time otherwise // finding the embeding R1 in R2 takes a long time otherwise
SetClassGroupBounds("GRH"); SetClassGroupBounds("GRH");
time R2 := AbelianExtension(NumberField(Factorization(ChangeRing(f,R1))[1,1])); R2 := AbelianExtension(NumberField(Factorization(ChangeRing(f,R1))[1,1]));
time R2 := NumberField(R2); R2 := NumberField(R2);
//// ////
...@@ -83,26 +85,17 @@ procedure solve(f) ...@@ -83,26 +85,17 @@ procedure solve(f)
"\nQuadratic extension given by"; "\nQuadratic extension given by";
g181a; g181a;
time proj, phis := ModPhis(m, R2); proj, phis := ModPhis(m, R2);
"\nUsing rho(Frob)";
[ModPrimes(m, R2, proj, phi : debug:=debug, check_frob) : phi in phis];
"\nUsing rho(Frob)";
[ModPrimes(m, R2, proj, phi : debug:=debug,
check_frob, check_residual := residual_traces) : phi in phis];
"\nUsing the cycle decomposition of rho(Frob)"; "\nUsing the cycle decomposition of rho(Frob)";
[ModPrimes(m, R2, proj, phi : debug:=debug, [ModPrimes(m, R2, proj, phi : debug:=debug,
check_residual:=residual_traces) : phi in phis]; check_residual:=residual_traces) : phi in phis];
"\nUsing the order of rho(Frob)"; "\nUsing the actual rho(Frob)";
[ModPrimes(m, R2, proj, phi : debug:=debug, only_order, [ModPrimes(m, R2, proj, phi : debug:=debug, check_frob,
check_residual:=residual_traces) : phi in phis]; check_residual:=residual_traces) : phi in phis];
"\nNOT Using the order of rho(Frob)"; "\n";
[ModPrimes(m, R2, proj, phi : debug:=debug, no_order,
check_residual := residual_traces) : phi in phis];
end procedure; end procedure;
...@@ -111,3 +104,5 @@ solve(g137b); ...@@ -111,3 +104,5 @@ solve(g137b);
solve(g181a); solve(g181a);
solve(g181b); solve(g181b);
printf "\nTotal time: %o s\n", Cputime(t0);
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment