X-Git-Url: https://ginac.de/CLN/cln.git//cln.git?a=blobdiff_plain;ds=sidebyside;f=src%2Fnumtheory%2Fcl_nt_sqrtmodp.cc;h=5f804edc929da1a90ec5db70c9e23c4e02a1dee6;hb=3af2cde18b3aabed4c808b0113daa81c2263b0bd;hp=3727ff9851cf1b9afaa7b3b5662996ab34d31efe;hpb=d6d3cb2bdee716a84015b4e3428b3775e929dc58;p=cln.git diff --git a/src/numtheory/cl_nt_sqrtmodp.cc b/src/numtheory/cl_nt_sqrtmodp.cc index 3727ff9..5f804ed 100644 --- a/src/numtheory/cl_nt_sqrtmodp.cc +++ b/src/numtheory/cl_nt_sqrtmodp.cc @@ -113,12 +113,13 @@ static const sqrt_mod_p_t cantor_zassenhaus_sqrt (const cl_modint_ring& R, const }; const gcd_result gcd (const pol2& u) { - if (zerop(u.c1)) + if (zerop(u.c1)) { // constant polynomial u(X) if (zerop(u.c0)) return gcd_result(2); else return gcd_result(0); + } // u(X) = c0 + c1*X has zero -c0/c1. var cl_MI_x c1inv = R->recip(u.c1); if (c1inv.condition)