1*86d7f5d3SJohn Marino /* mpf_cmp_d -- compare mpf and double.
2*86d7f5d3SJohn Marino
3*86d7f5d3SJohn Marino Copyright 2001, 2003 Free Software Foundation, Inc.
4*86d7f5d3SJohn Marino
5*86d7f5d3SJohn Marino This file is part of the GNU MP Library.
6*86d7f5d3SJohn Marino
7*86d7f5d3SJohn Marino The GNU MP Library is free software; you can redistribute it and/or modify
8*86d7f5d3SJohn Marino it under the terms of the GNU Lesser General Public License as published by
9*86d7f5d3SJohn Marino the Free Software Foundation; either version 3 of the License, or (at your
10*86d7f5d3SJohn Marino option) any later version.
11*86d7f5d3SJohn Marino
12*86d7f5d3SJohn Marino The GNU MP Library is distributed in the hope that it will be useful, but
13*86d7f5d3SJohn Marino WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
14*86d7f5d3SJohn Marino or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
15*86d7f5d3SJohn Marino License for more details.
16*86d7f5d3SJohn Marino
17*86d7f5d3SJohn Marino You should have received a copy of the GNU Lesser General Public License
18*86d7f5d3SJohn Marino along with the GNU MP Library. If not, see http://www.gnu.org/licenses/. */
19*86d7f5d3SJohn Marino
20*86d7f5d3SJohn Marino #include "config.h"
21*86d7f5d3SJohn Marino
22*86d7f5d3SJohn Marino #if HAVE_FLOAT_H
23*86d7f5d3SJohn Marino #include <float.h> /* for DBL_MAX */
24*86d7f5d3SJohn Marino #endif
25*86d7f5d3SJohn Marino
26*86d7f5d3SJohn Marino #include "gmp.h"
27*86d7f5d3SJohn Marino #include "gmp-impl.h"
28*86d7f5d3SJohn Marino
29*86d7f5d3SJohn Marino int
mpf_cmp_d(mpf_srcptr f,double d)30*86d7f5d3SJohn Marino mpf_cmp_d (mpf_srcptr f, double d)
31*86d7f5d3SJohn Marino {
32*86d7f5d3SJohn Marino mp_limb_t darray[LIMBS_PER_DOUBLE];
33*86d7f5d3SJohn Marino mpf_t df;
34*86d7f5d3SJohn Marino
35*86d7f5d3SJohn Marino /* d=NaN has no sensible return value, so raise an exception.
36*86d7f5d3SJohn Marino d=Inf or -Inf is always bigger than z. */
37*86d7f5d3SJohn Marino DOUBLE_NAN_INF_ACTION (d,
38*86d7f5d3SJohn Marino __gmp_invalid_operation (),
39*86d7f5d3SJohn Marino return (d < 0.0 ? 1 : -1));
40*86d7f5d3SJohn Marino
41*86d7f5d3SJohn Marino if (d == 0.0)
42*86d7f5d3SJohn Marino return SIZ(f);
43*86d7f5d3SJohn Marino
44*86d7f5d3SJohn Marino PTR(df) = darray;
45*86d7f5d3SJohn Marino SIZ(df) = (d >= 0.0 ? LIMBS_PER_DOUBLE : -LIMBS_PER_DOUBLE);
46*86d7f5d3SJohn Marino EXP(df) = __gmp_extract_double (darray, ABS(d));
47*86d7f5d3SJohn Marino
48*86d7f5d3SJohn Marino return mpf_cmp (f, df);
49*86d7f5d3SJohn Marino }
50