diff options
author | Torbjorn Granlund <tege@gmplib.org> | 2008-11-09 00:34:40 +0100 |
---|---|---|
committer | Torbjorn Granlund <tege@gmplib.org> | 2008-11-09 00:34:40 +0100 |
commit | b9458898ec6b21df2e34be9dc660552998e6f92e (patch) | |
tree | 561dbd275feba7db9e1fceecae77d1e3cc045e0f /mpf | |
parent | 4783e96e2c0125f5d0d80afc7c85da9bac30686f (diff) | |
download | gmp-b9458898ec6b21df2e34be9dc660552998e6f92e.tar.gz |
Partial rewrite; Compare the right number of bits.
Diffstat (limited to 'mpf')
-rw-r--r-- | mpf/eq.c | 24 |
1 files changed, 16 insertions, 8 deletions
@@ -19,6 +19,7 @@ along with the GNU MP Library. If not, see http://www.gnu.org/licenses/. */ #include "gmp.h" #include "gmp-impl.h" +#include "longlong.h" int mpf_eq (mpf_srcptr u, mpf_srcptr v, unsigned long int n_bits) @@ -27,6 +28,7 @@ mpf_eq (mpf_srcptr u, mpf_srcptr v, unsigned long int n_bits) mp_size_t usize, vsize, size, i; mp_exp_t uexp, vexp; mp_limb_t diff; + int cnt; uexp = u->_mp_exp; vexp = v->_mp_exp; @@ -54,10 +56,8 @@ mpf_eq (mpf_srcptr u, mpf_srcptr v, unsigned long int n_bits) /* U and V have the same sign and are both non-zero. */ /* 2. Are the exponents different? */ - if (uexp > vexp) - return 0; /* ??? handle (uexp = vexp + 1) */ - if (vexp > uexp) - return 0; /* ??? handle (vexp = uexp + 1) */ + if (uexp != vexp) + return 0; usize = ABS (usize); vsize = ABS (vsize); @@ -94,11 +94,19 @@ mpf_eq (mpf_srcptr u, mpf_srcptr v, unsigned long int n_bits) size = usize; } - if (size > (n_bits + GMP_NUMB_BITS - 1) / GMP_NUMB_BITS) - size = (n_bits + GMP_NUMB_BITS - 1) / GMP_NUMB_BITS; + up += usize; /* point just above most significant limb */ + vp += vsize; /* point just above most significant limb */ - up += usize - size; - vp += vsize - size; + count_leading_zeros (cnt, up[-1]); + if ((vp[-1] >> (GMP_LIMB_BITS - 1 - cnt)) != 1) + return 0; /* msb positions different */ + + n_bits += cnt - GMP_NAIL_BITS; + + size = MIN (size, (n_bits + GMP_NUMB_BITS - 1) / GMP_NUMB_BITS); + + up -= size; /* point at least significant relevant limb */ + vp -= size; /* point at least significant relevant limb */ for (i = size - 1; i > 0; i--) { |