mpn/generic/sqrtrem.c

163   mp_limb_t qhl, q, u, np0, sp0, rp0, q2;
170   qhl = 0;
174       qhl++;
182   q += (qhl & 1) << (Prec - 1);
183   qhl >>= 1; /* if qhl=1, necessary q=0 as qhl*2^Prec + q <= 2^Prec */
184   /* now we have (initial rp0)<<Prec + np0>>Prec = (qhl<<Prec + q) * (2sp0) + u */
185   sp0 = ((sp0 + qhl) << Prec) + q;
188   /* subtract q * q or qhl*2^(2*Prec) from rp */
190   cc -= (rp0 < q2) + qhl;
192   /* now subtract 2*q*2^Prec + 2^(2*Prec) if qhl is set */