dist/mpz/n_pow_ui.c

67    The initial powering for bsize==1 into blimb or blimb:blimb_low doesn't
168   mp_limb_t      blimb, rl;
209   for (blimb = *bp; blimb == 0; blimb = *++bp)
217   count_trailing_zeros (btwos, blimb);
218   blimb >>= btwos;
242       /* Power up as far as possible within blimb.  We start here with e!=0,
244 	 Notice this code works when blimb==1 too, reaching e==0.  */
246       while (blimb <= GMP_NUMB_HALFMAX)
248 	  TRACE (printf ("small e=0x%lX blimb=0x%lX rl=0x%lX\n",
249 			 e, blimb, rl));
252 	    rl *= blimb;
256 	  blimb *= blimb;
261 		     e, blimb, rl));
263       /* Can power b once more into blimb:blimb_low */
268 	  umul_ppmm (rl_high, rl, rl, blimb << GMP_NAIL_BITS);
272       umul_ppmm (blimb, blimb_low, blimb, blimb << GMP_NAIL_BITS);
276       TRACE (printf ("double power e=0x%lX blimb=0x%lX:0x%lX rl=0x%lX:%lX\n",
277 		     e, blimb, blimb_low, rl_high, rl));
303       TRACE (printf ("small power e=0x%lX blimb=0x%lX rl=0x%lX\n",
304 		     e, blimb, rl));
325 	blimb |= (bsecond << (GMP_NUMB_BITS - btwos)) & GMP_NUMB_MASK;
334       TRACE (printf ("bsize==2 using b=0x%lX:%lX", bsecond, blimb));
336       blimb_low = blimb;
339       b_twolimbs[0] = blimb;
342       blimb = bsecond;
357       blimb = bp[bsize-1];
363   /* At this point blimb is the most significant limb of the base to use.
375      If e==0 then blimb won't be anything useful (though it will be
379   ASSERT (blimb != 0);
380   count_leading_zeros (cnt, blimb);
390   TRACE (printf ("ralloc %ld, from bsize=%ld blimb=0x%lX cnt=%d\n",
391 		 ralloc, bsize, blimb, cnt));
446 	  ASSERT (blimb != 0);
453 	  rp[1] = blimb;
457 	  mult[1] = blimb;
488 	  rp[0] = blimb;
500 		MPN_MUL_1 (rp, rsize, ralloc, blimb);