version 1.1.1.1, 2000/01/10 15:35:27 |
version 1.1.1.3, 2003/08/25 16:06:33 |
|
|
/* mpz_set_d(integer, val) -- Assign INTEGER with a double value VAL. |
/* mpz_set_d(integer, val) -- Assign INTEGER with a double value VAL. |
|
|
Copyright (C) 1995 Free Software Foundation, Inc. |
Copyright 1995, 1996, 2000, 2001, 2002 Free Software Foundation, Inc. |
|
|
This file is part of the GNU MP Library. |
This file is part of the GNU MP Library. |
|
|
The GNU MP Library is free software; you can redistribute it and/or modify |
The GNU MP Library is free software; you can redistribute it and/or modify |
it under the terms of the GNU Library General Public License as published by |
it under the terms of the GNU Lesser General Public License as published by |
the Free Software Foundation; either version 2 of the License, or (at your |
the Free Software Foundation; either version 2.1 of the License, or (at your |
option) any later version. |
option) any later version. |
|
|
The GNU MP Library is distributed in the hope that it will be useful, but |
The GNU MP Library is distributed in the hope that it will be useful, but |
WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY |
WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY |
or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Library General Public |
or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public |
License for more details. |
License for more details. |
|
|
You should have received a copy of the GNU Library General Public License |
You should have received a copy of the GNU Lesser General Public License |
along with the GNU MP Library; see the file COPYING.LIB. If not, write to |
along with the GNU MP Library; see the file COPYING.LIB. If not, write to |
the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston, |
the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston, |
MA 02111-1307, USA. */ |
MA 02111-1307, USA. */ |
Line 23 MA 02111-1307, USA. */ |
|
Line 23 MA 02111-1307, USA. */ |
|
#include "gmp-impl.h" |
#include "gmp-impl.h" |
|
|
void |
void |
#if __STDC__ |
|
mpz_set_d (mpz_ptr r, double d) |
mpz_set_d (mpz_ptr r, double d) |
#else |
|
mpz_set_d (r, d) |
|
mpz_ptr r; |
|
double d; |
|
#endif |
|
{ |
{ |
int negative; |
int negative; |
mp_size_t size; |
mp_limb_t tp[LIMBS_PER_DOUBLE]; |
mp_limb_t tp[3]; |
|
mp_ptr rp; |
mp_ptr rp; |
|
mp_size_t rn; |
|
|
negative = d < 0; |
negative = d < 0; |
d = ABS (d); |
d = ABS (d); |
|
|
return; |
return; |
} |
} |
|
|
size = __gmp_extract_double (tp, d); |
rn = __gmp_extract_double (tp, d); |
|
|
if (ALLOC(r) < size) |
if (ALLOC(r) < rn) |
_mpz_realloc (r, size); |
_mpz_realloc (r, rn); |
|
|
rp = PTR (r); |
rp = PTR (r); |
|
|
#if BITS_PER_MP_LIMB == 32 |
#if BITS_PER_MP_LIMB == 32 |
switch (size) |
switch (rn) |
{ |
{ |
default: |
default: |
MPN_ZERO (rp, size - 3); |
MPN_ZERO (rp, rn - 3); |
rp += size - 3; |
rp += rn - 3; |
|
/* fall through */ |
case 3: |
case 3: |
rp[2] = tp[2]; |
rp[2] = tp[2]; |
rp[1] = tp[1]; |
rp[1] = tp[1]; |
|
|
rp[0] = tp[1]; |
rp[0] = tp[1]; |
break; |
break; |
case 1: |
case 1: |
|
/* handled in "small aguments" case above */ |
abort (); |
abort (); |
} |
} |
#else |
#else |
switch (size) |
switch (rn) |
{ |
{ |
default: |
default: |
MPN_ZERO (rp, size - 2); |
MPN_ZERO (rp, rn - 2); |
rp += size - 2; |
rp += rn - 2; |
|
/* fall through */ |
case 2: |
case 2: |
rp[1] = tp[1]; |
rp[1] = tp[1], rp[0] = tp[0]; |
rp[0] = tp[0]; |
|
break; |
break; |
case 1: |
case 1: |
|
/* handled in "small aguments" case above */ |
abort (); |
abort (); |
} |
} |
#endif |
#endif |
|
|
SIZ(r) = negative ? -size : size; |
SIZ(r) = negative ? -rn : rn; |
} |
} |