mpn/generic/sqrtrem.c

163   mp_limb_t qhl, q, u, np0, sp0, rp0, q2;
171   rp0 = rp[0];
172   while (rp0 >= sp0)
175       rp0 -= sp0;
177   /* now rp0 < sp0 < 2^Prec */
178   rp0 = (rp0 << Prec) + (np0 >> Prec);
180   q = rp0 / u;
181   u = rp0 - q * u;
184   /* now we have (initial rp0)<<Prec + np0>>Prec = (qhl<<Prec + q) * (2sp0) + u */
187   rp0 = ((u << Prec) & GMP_NUMB_MASK) + (np0 & (((mp_limb_t) 1 << Prec) - 1));
190   cc -= (rp0 < q2) + qhl;
191   rp0 -= q2;
197 	  rp0 += sp0;
198 	  cc += rp0 < sp0;
203       rp0 += sp0;
204       cc += rp0 < sp0;
207   rp[0] = rp0;