mirrors
/
macssh
mirror of https://github.com/macssh/macssh.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

added recast

This commit is contained in:
chombier 2002-08-24 10:21:43 +00:00
parent 4afc46b950
commit 5407c23780
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
gmp/mpfr/round.c
View File

@ -176,7 +176,7 @@ mpfr_round(x, RND_MODE, prec)
nw = prec / BITS_PER_MP_LIMB; nw = prec / BITS_PER_MP_LIMB;
if (prec & (BITS_PER_MP_LIMB - 1)) nw++; if (prec & (BITS_PER_MP_LIMB - 1)) nw++;
TMP_MARK(marker); TMP_MARK(marker);
tmp = TMP_ALLOC (nw * BYTES_PER_MP_LIMB); tmp = (mp_limb_t *)TMP_ALLOC (nw * BYTES_PER_MP_LIMB);
carry = mpfr_round_raw(tmp, MANT(x), PREC(x), (SIGN(x)<0), prec, RND_MODE); carry = mpfr_round_raw(tmp, MANT(x), PREC(x), (SIGN(x)<0), prec, RND_MODE);
if (carry) if (carry)
@ -267,7 +267,7 @@ mpfr_can_round_raw(bp, bn, neg, err, rnd1, rnd2, prec)
switch (rnd1) { switch (rnd1) {
case GMP_RNDZ: /* b <= x <= b+2^(EXP(b)-err) */ case GMP_RNDZ: /* b <= x <= b+2^(EXP(b)-err) */
tmp = TMP_ALLOC(tn*BYTES_PER_MP_LIMB); tmp = (mp_limb_t *)TMP_ALLOC(tn*BYTES_PER_MP_LIMB);
cc = (bp[bn-1]>>l1) & 1; cc = (bp[bn-1]>>l1) & 1;
cc ^= mpfr_round_raw2(bp, bn, neg, rnd2, prec); cc ^= mpfr_round_raw2(bp, bn, neg, rnd2, prec);
@ -283,7 +283,7 @@ mpfr_can_round_raw(bp, bn, neg, err, rnd1, rnd2, prec)
return (cc == cc2); return (cc == cc2);
case GMP_RNDU: /* b-2^(EXP(b)-err) <= x <= b */ case GMP_RNDU: /* b-2^(EXP(b)-err) <= x <= b */
tmp = TMP_ALLOC(tn*BYTES_PER_MP_LIMB); tmp = (mp_limb_t *)TMP_ALLOC(tn*BYTES_PER_MP_LIMB);
/* first round b */ /* first round b */
cc = (bp[bn-1]>>l1) & 1; cc = (bp[bn-1]>>l1) & 1;
cc ^= mpfr_round_raw2(bp, bn, neg, rnd2, prec); cc ^= mpfr_round_raw2(bp, bn, neg, rnd2, prec);
@ -301,7 +301,7 @@ mpfr_can_round_raw(bp, bn, neg, err, rnd1, rnd2, prec)
case GMP_RNDN: /* b-2^(EXP(b)-err-1) <= x <= b+2^(EXP(b)-err-1) */ case GMP_RNDN: /* b-2^(EXP(b)-err-1) <= x <= b+2^(EXP(b)-err-1) */
if (l==0) tn++; if (l==0) tn++;
tmp = TMP_ALLOC(tn*BYTES_PER_MP_LIMB); tmp = (mp_limb_t *)TMP_ALLOC(tn*BYTES_PER_MP_LIMB);
/* this case is the same than GMP_RNDZ, except we first have to /* this case is the same than GMP_RNDZ, except we first have to
subtract 2^(EXP(b)-err-1) from b */ subtract 2^(EXP(b)-err-1) from b */