dist/src/sub1sp.c

135   mpfr_prec_t cnt, INITIALIZED(sh);
168       count_leading_zeros (cnt, a0);
169       ap[0] = a0 << cnt;
170       bx -= cnt;
211           count_leading_zeros (cnt, a0);
212           if (cnt)
213             a0 = (a0 << cnt) | (sb >> (GMP_NUMB_BITS - cnt));
214           sb <<= cnt;
215           bx -= cnt;
320   mpfr_prec_t cnt;
352       count_leading_zeros (cnt, a0);
353       ap[0] = a0 << cnt;
354       bx -= cnt;
387               count_leading_zeros (cnt, a0);
388               if (cnt)
389                 a0 = (a0 << cnt) | (sb >> (GMP_NUMB_BITS - cnt));
390               sb <<= cnt;
391               bx -= cnt;
509   mpfr_prec_t cnt, INITIALIZED(sh);
551       count_leading_zeros (cnt, a1);
552       if (cnt)
554           ap[1] = (a1 << cnt) | (a0 >> (GMP_NUMB_BITS - cnt));
555           ap[0] = a0 << cnt;
556           bx -= cnt;
635           count_leading_zeros (cnt, a1);
636           if (cnt)
638               ap[1] = (a1 << cnt) | (a0 >> (GMP_NUMB_BITS - cnt));
639               a0 = (a0 << cnt) | (sb >> (GMP_NUMB_BITS - cnt));
640               sb <<= cnt;
641               bx -= cnt;
770   mpfr_prec_t cnt;
816       count_leading_zeros (cnt, a1);
817       if (cnt)
819           /* shift [a1,a0] left by cnt bit and store in result */
820           ap[1] = (a1 << cnt) | (a0 >> (GMP_NUMB_BITS - cnt));
821           ap[0] = a0 << cnt;
822           bx -= cnt;
897               count_leading_zeros (cnt, a1);
898               if (cnt)
900                   /* shift [a1,a0,sb] left by cnt bits and adjust exponent */
901                   a1 = (a1 << cnt) | (a0 >> (GMP_NUMB_BITS - cnt));
902                   a0 = (a0 << cnt) | (sb >> (GMP_NUMB_BITS - cnt));
903                   sb <<= cnt;
904                   bx -= cnt;
1052   mpfr_prec_t cnt, INITIALIZED(sh);
1104       count_leading_zeros (cnt, a2);
1105       if (cnt)
1107           ap[2] = (a2 << cnt) | (a1 >> (GMP_NUMB_BITS - cnt));
1108           ap[1] = (a1 << cnt) | (a0 >> (GMP_NUMB_BITS - cnt));
1109           ap[0] = a0 << cnt;
1110           bx -= cnt;
1182           count_leading_zeros (cnt, a2);
1183           if (cnt)
1185               ap[2] = (a2 << cnt) | (a1 >> (GMP_NUMB_BITS - cnt));
1186               ap[1] = (a1 << cnt) | (a0 >> (GMP_NUMB_BITS - cnt));
1187               a0 = (a0 << cnt) | (sb >> (GMP_NUMB_BITS - cnt));
1188               sb <<= cnt;
1189               bx -= cnt;
1428   mpfr_prec_t p, sh, cnt;
1550               count_leading_zeros (cnt, limb);
1551               /* Warning: cnt can be 0 when we come from the case SubD1Lose
1553               if (MPFR_LIKELY(cnt))
1555                   mpn_lshift (ap, ap, n, cnt); /* Normalize number */
1556                   bx -= cnt; /* Update final expo */
1575               count_leading_zeros(cnt, limb);
1579               if (cnt)
1580                 mpn_lshift (ap + len, ap, k, cnt); /* Normalize the High Limb*/
1585               bx -= cnt + len * GMP_NUMB_BITS; /* update exponent */