Commit ec01aa15 authored by Peter Bruin's avatar Peter Bruin

fix counting of basis elements found

parent 3d3d2221
......@@ -198,19 +198,21 @@ find_bases(GEN J, unsigned long l, GEN proj, int tries) {
Q = jacobian_Frob_polynomial(J, P, centerlift(gel(proj, i)), 1);
if (jacobian_is_zero(J, Q))
continue;
found++;
if (j == 0) {
R = jacobian_Frob(J, Q, 1);
if (is_multiple(J, Q, R, l))
if (is_multiple(J, Q, R, l)) {
gel(bases, i) = mkvec(Q);
else {
gel(bases, i) = mkvec2(Q, R);
found++;
} else {
gel(bases, i) = mkvec2(Q, R);
found += 2;
}
} else {
R = gmael(bases, i, 1);
if (!is_multiple(J, R, Q, l))
if (!is_multiple(J, R, Q, l)) {
gel(bases, i) = mkvec2(R, Q);
found++;
}
}
}
if (gc_needed(av, 1))
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment