diff options
author | zimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4> | 2007-10-08 10:09:56 +0000 |
---|---|---|
committer | zimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4> | 2007-10-08 10:09:56 +0000 |
commit | 5ae8239ff164cfd433a3d5320d2bc1c7c7809d70 (patch) | |
tree | 426c86121b17480771244c8a8846aaa9bf34d8d5 | |
parent | ae59366dc468b61fbeac5bad3c26219a482e562d (diff) | |
download | mpfr-5ae8239ff164cfd433a3d5320d2bc1c7c7809d70.tar.gz |
changed mpfr_can_round() call into MPFR_CAN_ROUND() call
git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@4879 280ebfd0-de03-0410-8827-d642c229c3f4
-rw-r--r-- | sin.c | 7 |
1 files changed, 1 insertions, 6 deletions
@@ -132,12 +132,7 @@ mpfr_sin (mpfr_ptr y, mpfr_srcptr x, mp_rnd_t rnd_mode) Since EXP(c) <= 1, 3-m-EXP(c) >= 2-m, thus the error is at most 2^(3-m-EXP(c)) in case of argument reduction. */ err = 2 * MPFR_GET_EXP (c) + (mp_exp_t) m - 3 - (reduce != 0); - if (mpfr_can_round (c, err, GMP_RNDN, GMP_RNDZ, - precy + (rnd_mode == GMP_RNDN))) - /* WARNING: even if we know c <= sin(x), don't give GMP_RNDZ - as 3rd argument to mpfr_can_round, since if c is exactly - representable to the target precision (inexact = 0 below), - we would have to add one ulp when rounding away from 0. */ + if (MPFR_CAN_ROUND (c, err, precy, rnd_mode)) break; /* check for huge cancellation (Near 0) */ |