|
Packit |
5c3484 |
/* mpz_cmp_d -- compare absolute values of mpz and double.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
Copyright 2001-2003 Free Software Foundation, Inc.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
This file is part of the GNU MP Library.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
The GNU MP Library is free software; you can redistribute it and/or modify
|
|
Packit |
5c3484 |
it under the terms of either:
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
* the GNU Lesser General Public License as published by the Free
|
|
Packit |
5c3484 |
Software Foundation; either version 3 of the License, or (at your
|
|
Packit |
5c3484 |
option) any later version.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
or
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
* the GNU General Public License as published by the Free Software
|
|
Packit |
5c3484 |
Foundation; either version 2 of the License, or (at your option) any
|
|
Packit |
5c3484 |
later version.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
or both in parallel, as here.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
The GNU MP Library is distributed in the hope that it will be useful, but
|
|
Packit |
5c3484 |
WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
|
|
Packit |
5c3484 |
or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
|
Packit |
5c3484 |
for more details.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
You should have received copies of the GNU General Public License and the
|
|
Packit |
5c3484 |
GNU Lesser General Public License along with the GNU MP Library. If not,
|
|
Packit |
5c3484 |
see https://www.gnu.org/licenses/. */
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
#include "config.h"
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
#if HAVE_FLOAT_H
|
|
Packit |
5c3484 |
#include <float.h> /* for DBL_MAX */
|
|
Packit |
5c3484 |
#endif
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
#include "gmp.h"
|
|
Packit |
5c3484 |
#include "gmp-impl.h"
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
#define RETURN_CMP(zl, dl) \
|
|
Packit |
5c3484 |
do { \
|
|
Packit |
5c3484 |
zlimb = (zl); \
|
|
Packit |
5c3484 |
dlimb = (dl); \
|
|
Packit |
5c3484 |
if (zlimb != dlimb) \
|
|
Packit |
5c3484 |
return (zlimb >= dlimb ? ret : -ret); \
|
|
Packit |
5c3484 |
} while (0)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
#define RETURN_NONZERO(ptr, size, val) \
|
|
Packit |
5c3484 |
do { \
|
|
Packit |
5c3484 |
mp_size_t __i; \
|
|
Packit |
5c3484 |
for (__i = (size)-1; __i >= 0; __i--) \
|
|
Packit |
5c3484 |
if ((ptr)[__i] != 0) \
|
|
Packit |
5c3484 |
return val; \
|
|
Packit |
5c3484 |
return 0; \
|
|
Packit |
5c3484 |
} while (0)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
int
|
|
Packit |
5c3484 |
mpz_cmp_d (mpz_srcptr z, double d)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
mp_limb_t darray[LIMBS_PER_DOUBLE], zlimb, dlimb;
|
|
Packit |
5c3484 |
mp_srcptr zp;
|
|
Packit |
5c3484 |
mp_size_t zsize;
|
|
Packit |
5c3484 |
int dexp, ret;
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
/* d=NaN is an invalid operation, there's no sensible return value.
|
|
Packit |
5c3484 |
d=Inf or -Inf is always bigger than z. */
|
|
Packit |
5c3484 |
DOUBLE_NAN_INF_ACTION (d, __gmp_invalid_operation (), goto z_zero);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
/* 1. Either operand zero. */
|
|
Packit |
5c3484 |
zsize = SIZ(z);
|
|
Packit |
5c3484 |
if (d == 0.0)
|
|
Packit |
5c3484 |
return zsize;
|
|
Packit |
5c3484 |
if (zsize == 0)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
z_zero:
|
|
Packit |
5c3484 |
return (d < 0.0 ? 1 : -1);
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
/* 2. Opposite signs. */
|
|
Packit |
5c3484 |
if (zsize >= 0)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
if (d < 0.0)
|
|
Packit |
5c3484 |
return 1; /* >=0 cmp <0 */
|
|
Packit |
5c3484 |
ret = 1;
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
else
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
if (d >= 0.0)
|
|
Packit |
5c3484 |
return -1; /* <0 cmp >=0 */
|
|
Packit |
5c3484 |
ret = -1;
|
|
Packit |
5c3484 |
d = -d;
|
|
Packit |
5c3484 |
zsize = -zsize;
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
/* 3. Small d, knowing abs(z) >= 1. */
|
|
Packit |
5c3484 |
if (d < 1.0)
|
|
Packit |
5c3484 |
return ret;
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
dexp = __gmp_extract_double (darray, d);
|
|
Packit |
5c3484 |
ASSERT (dexp >= 1);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
/* 4. Check for different high limb positions. */
|
|
Packit |
5c3484 |
if (zsize != dexp)
|
|
Packit |
5c3484 |
return (zsize >= dexp ? ret : -ret);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
/* 5. Limb data. */
|
|
Packit |
5c3484 |
zp = PTR(z);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
#if LIMBS_PER_DOUBLE == 2
|
|
Packit |
5c3484 |
RETURN_CMP (zp[zsize-1], darray[1]);
|
|
Packit |
5c3484 |
if (zsize == 1)
|
|
Packit |
5c3484 |
return (darray[0] != 0 ? -ret : 0);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
RETURN_CMP (zp[zsize-2], darray[0]);
|
|
Packit |
5c3484 |
RETURN_NONZERO (zp, zsize-2, ret);
|
|
Packit |
5c3484 |
#endif
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
#if LIMBS_PER_DOUBLE == 3
|
|
Packit |
5c3484 |
RETURN_CMP (zp[zsize-1], darray[2]);
|
|
Packit |
5c3484 |
if (zsize == 1)
|
|
Packit |
5c3484 |
return ((darray[0] | darray[1]) != 0 ? -ret : 0);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
RETURN_CMP (zp[zsize-2], darray[1]);
|
|
Packit |
5c3484 |
if (zsize == 2)
|
|
Packit |
5c3484 |
return (darray[0] != 0 ? -ret : 0);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
RETURN_CMP (zp[zsize-3], darray[0]);
|
|
Packit |
5c3484 |
RETURN_NONZERO (zp, zsize-3, ret);
|
|
Packit |
5c3484 |
#endif
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
#if LIMBS_PER_DOUBLE >= 4
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
int i;
|
|
Packit |
5c3484 |
for (i = 1; i <= LIMBS_PER_DOUBLE; i++)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
RETURN_CMP (zp[zsize-i], darray[LIMBS_PER_DOUBLE-i]);
|
|
Packit |
5c3484 |
if (i >= zsize)
|
|
Packit |
5c3484 |
RETURN_NONZERO (darray, LIMBS_PER_DOUBLE-i, -ret);
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
RETURN_NONZERO (zp, zsize-LIMBS_PER_DOUBLE, ret);
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
#endif
|
|
Packit |
5c3484 |
}
|