xref: /dragonfly/contrib/gmp/mpf/set.c (revision 86d7f5d3)
1*86d7f5d3SJohn Marino /* mpf_set -- Assign a float from another float.
2*86d7f5d3SJohn Marino 
3*86d7f5d3SJohn Marino Copyright 1993, 1994, 1995, 2001, 2004 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 "gmp.h"
21*86d7f5d3SJohn Marino #include "gmp-impl.h"
22*86d7f5d3SJohn Marino 
23*86d7f5d3SJohn Marino void
mpf_set(mpf_ptr r,mpf_srcptr u)24*86d7f5d3SJohn Marino mpf_set (mpf_ptr r, mpf_srcptr u)
25*86d7f5d3SJohn Marino {
26*86d7f5d3SJohn Marino   mp_ptr rp, up;
27*86d7f5d3SJohn Marino   mp_size_t size, asize;
28*86d7f5d3SJohn Marino   mp_size_t prec;
29*86d7f5d3SJohn Marino 
30*86d7f5d3SJohn Marino   prec = r->_mp_prec + 1;		/* lie not to lose precision in assignment */
31*86d7f5d3SJohn Marino   size = u->_mp_size;
32*86d7f5d3SJohn Marino   asize = ABS (size);
33*86d7f5d3SJohn Marino   rp = r->_mp_d;
34*86d7f5d3SJohn Marino   up = u->_mp_d;
35*86d7f5d3SJohn Marino 
36*86d7f5d3SJohn Marino   if (asize > prec)
37*86d7f5d3SJohn Marino     {
38*86d7f5d3SJohn Marino       up += asize - prec;
39*86d7f5d3SJohn Marino       asize = prec;
40*86d7f5d3SJohn Marino     }
41*86d7f5d3SJohn Marino 
42*86d7f5d3SJohn Marino   r->_mp_exp = u->_mp_exp;
43*86d7f5d3SJohn Marino   r->_mp_size = size >= 0 ? asize : -asize;
44*86d7f5d3SJohn Marino   MPN_COPY_INCR (rp, up, asize);
45*86d7f5d3SJohn Marino }
46