|
Packit |
5c3484 |
/* mpq_set_d(mpq_t q, double d) -- Set q to d without rounding.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
Copyright 2000, 2002, 2003, 2012, 2014 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 |
#include "longlong.h"
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
#if LIMBS_PER_DOUBLE > 4
|
|
Packit |
5c3484 |
choke me
|
|
Packit |
5c3484 |
#endif
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
void
|
|
Packit |
5c3484 |
mpq_set_d (mpq_ptr dest, double d)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
int negative;
|
|
Packit |
5c3484 |
mp_exp_t exp;
|
|
Packit |
5c3484 |
mp_limb_t tp[LIMBS_PER_DOUBLE];
|
|
Packit |
5c3484 |
mp_ptr np, dp;
|
|
Packit |
5c3484 |
mp_size_t nn, dn;
|
|
Packit |
5c3484 |
int c;
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
DOUBLE_NAN_INF_ACTION (d,
|
|
Packit |
5c3484 |
__gmp_invalid_operation (),
|
|
Packit |
5c3484 |
__gmp_invalid_operation ());
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
negative = d < 0;
|
|
Packit |
5c3484 |
d = ABS (d);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
exp = __gmp_extract_double (tp, d);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
/* There are two main version of the conversion. The `then' arm handles
|
|
Packit |
5c3484 |
numbers with a fractional part, while the `else' arm handles integers. */
|
|
Packit |
5c3484 |
#if LIMBS_PER_DOUBLE == 4
|
|
Packit |
5c3484 |
if (exp <= 1 || (exp == 2 && (tp[0] | tp[1]) != 0))
|
|
Packit |
5c3484 |
#endif
|
|
Packit |
5c3484 |
#if LIMBS_PER_DOUBLE == 3
|
|
Packit |
5c3484 |
if (exp <= 1 || (exp == 2 && tp[0] != 0))
|
|
Packit |
5c3484 |
#endif
|
|
Packit |
5c3484 |
#if LIMBS_PER_DOUBLE == 2
|
|
Packit |
5c3484 |
if (exp <= 1)
|
|
Packit |
5c3484 |
#endif
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
if (d == 0.0)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
SIZ(NUM(dest)) = 0;
|
|
Packit |
5c3484 |
SIZ(DEN(dest)) = 1;
|
|
Packit |
5c3484 |
PTR(DEN(dest))[0] = 1;
|
|
Packit |
5c3484 |
return;
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
#if LIMBS_PER_DOUBLE == 4
|
|
Packit |
5c3484 |
np = MPZ_NEWALLOC (NUM(dest), 4);
|
|
Packit |
5c3484 |
if ((tp[0] | tp[1] | tp[2]) == 0)
|
|
Packit |
5c3484 |
np[0] = tp[3], nn = 1;
|
|
Packit |
5c3484 |
else if ((tp[0] | tp[1]) == 0)
|
|
Packit |
5c3484 |
np[1] = tp[3], np[0] = tp[2], nn = 2;
|
|
Packit |
5c3484 |
else if (tp[0] == 0)
|
|
Packit |
5c3484 |
np[2] = tp[3], np[1] = tp[2], np[0] = tp[1], nn = 3;
|
|
Packit |
5c3484 |
else
|
|
Packit |
5c3484 |
np[3] = tp[3], np[2] = tp[2], np[1] = tp[1], np[0] = tp[0], nn = 4;
|
|
Packit |
5c3484 |
#endif
|
|
Packit |
5c3484 |
#if LIMBS_PER_DOUBLE == 3
|
|
Packit |
5c3484 |
np = MPZ_NEWALLOC (NUM(dest), 3);
|
|
Packit |
5c3484 |
if ((tp[0] | tp[1]) == 0)
|
|
Packit |
5c3484 |
np[0] = tp[2], nn = 1;
|
|
Packit |
5c3484 |
else if (tp[0] == 0)
|
|
Packit |
5c3484 |
np[1] = tp[2], np[0] = tp[1], nn = 2;
|
|
Packit |
5c3484 |
else
|
|
Packit |
5c3484 |
np[2] = tp[2], np[1] = tp[1], np[0] = tp[0], nn = 3;
|
|
Packit |
5c3484 |
#endif
|
|
Packit |
5c3484 |
#if LIMBS_PER_DOUBLE == 2
|
|
Packit |
5c3484 |
np = MPZ_NEWALLOC (NUM(dest), 2);
|
|
Packit |
5c3484 |
if (tp[0] == 0)
|
|
Packit |
5c3484 |
np[0] = tp[1], nn = 1;
|
|
Packit |
5c3484 |
else
|
|
Packit |
5c3484 |
np[1] = tp[1], np[0] = tp[0], nn = 2;
|
|
Packit |
5c3484 |
#endif
|
|
Packit |
5c3484 |
dn = nn + 1 - exp;
|
|
Packit |
5c3484 |
ASSERT (dn > 0); /* -exp >= -1; nn >= 1*/
|
|
Packit |
5c3484 |
dp = MPZ_NEWALLOC (DEN(dest), dn);
|
|
Packit |
5c3484 |
MPN_ZERO (dp, dn - 1);
|
|
Packit |
5c3484 |
dp[dn - 1] = 1;
|
|
Packit |
5c3484 |
count_trailing_zeros (c, np[0] | dp[0]);
|
|
Packit |
5c3484 |
if (c != 0)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
mpn_rshift (np, np, nn, c);
|
|
Packit |
5c3484 |
nn -= np[nn - 1] == 0;
|
|
Packit |
5c3484 |
--dn;
|
|
Packit |
5c3484 |
dp[dn - 1] = CNST_LIMB(1) << (GMP_LIMB_BITS - c);
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
SIZ(DEN(dest)) = dn;
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
else
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
nn = exp;
|
|
Packit |
5c3484 |
np = MPZ_NEWALLOC (NUM(dest), nn);
|
|
Packit |
5c3484 |
switch (nn)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
default:
|
|
Packit |
5c3484 |
MPN_ZERO (np, nn - LIMBS_PER_DOUBLE);
|
|
Packit |
5c3484 |
np += nn - LIMBS_PER_DOUBLE;
|
|
Packit |
5c3484 |
/* fall through */
|
|
Packit |
5c3484 |
#if LIMBS_PER_DOUBLE == 2
|
|
Packit |
5c3484 |
case 2:
|
|
Packit |
5c3484 |
np[1] = tp[1], np[0] = tp[0];
|
|
Packit |
5c3484 |
break;
|
|
Packit |
5c3484 |
#endif
|
|
Packit |
5c3484 |
#if LIMBS_PER_DOUBLE == 3
|
|
Packit |
5c3484 |
case 3:
|
|
Packit |
5c3484 |
np[2] = tp[2], np[1] = tp[1], np[0] = tp[0];
|
|
Packit |
5c3484 |
break;
|
|
Packit |
5c3484 |
case 2:
|
|
Packit |
5c3484 |
np[1] = tp[2], np[0] = tp[1];
|
|
Packit |
5c3484 |
break;
|
|
Packit |
5c3484 |
#endif
|
|
Packit |
5c3484 |
#if LIMBS_PER_DOUBLE == 4
|
|
Packit |
5c3484 |
case 4:
|
|
Packit |
5c3484 |
np[3] = tp[3], np[2] = tp[2], np[1] = tp[1], np[0] = tp[0];
|
|
Packit |
5c3484 |
break;
|
|
Packit |
5c3484 |
case 3:
|
|
Packit |
5c3484 |
np[2] = tp[3], np[1] = tp[2], np[0] = tp[1];
|
|
Packit |
5c3484 |
break;
|
|
Packit |
5c3484 |
case 2:
|
|
Packit |
5c3484 |
np[1] = tp[3], np[0] = tp[2];
|
|
Packit |
5c3484 |
break;
|
|
Packit |
5c3484 |
#endif
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
*PTR(DEN(dest)) = 1;
|
|
Packit |
5c3484 |
SIZ(DEN(dest)) = 1;
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
SIZ(NUM(dest)) = negative ? -nn : nn;
|
|
Packit |
5c3484 |
}
|