Annotation of OpenXM_contrib/gmp/mpz/set_f.c, Revision 1.1.1.3
1.1 maekawa 1: /* mpz_set_f (dest_integer, src_float) -- Assign DEST_INTEGER from SRC_FLOAT.
2:
1.1.1.3 ! ohara 3: Copyright 1996, 2001 Free Software Foundation, Inc.
1.1 maekawa 4:
5: This file is part of the GNU MP Library.
6:
7: The GNU MP Library is free software; you can redistribute it and/or modify
1.1.1.2 maekawa 8: it under the terms of the GNU Lesser General Public License as published by
9: the Free Software Foundation; either version 2.1 of the License, or (at your
1.1 maekawa 10: option) any later version.
11:
12: The GNU MP Library is distributed in the hope that it will be useful, but
13: WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
1.1.1.2 maekawa 14: or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
1.1 maekawa 15: License for more details.
16:
1.1.1.2 maekawa 17: You should have received a copy of the GNU Lesser General Public License
1.1 maekawa 18: along with the GNU MP Library; see the file COPYING.LIB. If not, write to
19: the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
20: MA 02111-1307, USA. */
21:
22: #include "gmp.h"
23: #include "gmp-impl.h"
24:
1.1.1.3 ! ohara 25:
1.1 maekawa 26: void
27: mpz_set_f (mpz_ptr w, mpf_srcptr u)
28: {
1.1.1.3 ! ohara 29: mp_ptr wp, up;
! 30: mp_size_t size;
! 31: mp_exp_t exp;
1.1 maekawa 32:
1.1.1.3 ! ohara 33: /* abs(u)<1 truncates to zero */
1.1 maekawa 34: exp = EXP (u);
35: if (exp <= 0)
36: {
1.1.1.3 ! ohara 37: SIZ(w) = 0;
1.1 maekawa 38: return;
39: }
1.1.1.3 ! ohara 40:
! 41: MPZ_REALLOC (w, exp);
! 42: wp = PTR(w);
! 43: up = PTR(u);
! 44:
! 45: size = SIZ (u);
! 46: SIZ(w) = (size >= 0 ? exp : -exp);
! 47: size = ABS (size);
! 48:
! 49: if (exp > size)
1.1 maekawa 50: {
1.1.1.3 ! ohara 51: /* pad with low zeros to get a total "exp" many limbs */
! 52: mp_size_t zeros = exp - size;
! 53: MPN_ZERO (wp, zeros);
! 54: wp += zeros;
1.1 maekawa 55: }
56: else
57: {
1.1.1.3 ! ohara 58: /* exp<=size, trucate to the high "exp" many limbs */
! 59: up += (size - exp);
! 60: size = exp;
1.1 maekawa 61: }
62:
1.1.1.3 ! ohara 63: MPN_COPY (wp, up, size);
1.1 maekawa 64: }
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>