xref: /dflybsd-src/contrib/gmp/mpf/set_d.c (revision 86d7f5d305c6adaa56ff4582ece9859d73106103)
186d7f5d3SJohn Marino /* mpf_set_d -- Assign a float from a double.
286d7f5d3SJohn Marino 
386d7f5d3SJohn Marino Copyright 1993, 1994, 1995, 1996, 2001, 2003, 2004 Free Software Foundation,
486d7f5d3SJohn Marino Inc.
586d7f5d3SJohn Marino 
686d7f5d3SJohn Marino This file is part of the GNU MP Library.
786d7f5d3SJohn Marino 
886d7f5d3SJohn Marino The GNU MP Library is free software; you can redistribute it and/or modify
986d7f5d3SJohn Marino it under the terms of the GNU Lesser General Public License as published by
1086d7f5d3SJohn Marino the Free Software Foundation; either version 3 of the License, or (at your
1186d7f5d3SJohn Marino option) any later version.
1286d7f5d3SJohn Marino 
1386d7f5d3SJohn Marino The GNU MP Library is distributed in the hope that it will be useful, but
1486d7f5d3SJohn Marino WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
1586d7f5d3SJohn Marino or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
1686d7f5d3SJohn Marino License for more details.
1786d7f5d3SJohn Marino 
1886d7f5d3SJohn Marino You should have received a copy of the GNU Lesser General Public License
1986d7f5d3SJohn Marino along with the GNU MP Library.  If not, see http://www.gnu.org/licenses/.  */
2086d7f5d3SJohn Marino 
2186d7f5d3SJohn Marino #include "config.h"
2286d7f5d3SJohn Marino 
2386d7f5d3SJohn Marino #if HAVE_FLOAT_H
2486d7f5d3SJohn Marino #include <float.h>  /* for DBL_MAX */
2586d7f5d3SJohn Marino #endif
2686d7f5d3SJohn Marino 
2786d7f5d3SJohn Marino #include "gmp.h"
2886d7f5d3SJohn Marino #include "gmp-impl.h"
2986d7f5d3SJohn Marino 
3086d7f5d3SJohn Marino void
mpf_set_d(mpf_ptr r,double d)3186d7f5d3SJohn Marino mpf_set_d (mpf_ptr r, double d)
3286d7f5d3SJohn Marino {
3386d7f5d3SJohn Marino   int negative;
3486d7f5d3SJohn Marino 
3586d7f5d3SJohn Marino   DOUBLE_NAN_INF_ACTION (d,
3686d7f5d3SJohn Marino                          __gmp_invalid_operation (),
3786d7f5d3SJohn Marino                          __gmp_invalid_operation ());
3886d7f5d3SJohn Marino 
3986d7f5d3SJohn Marino   if (UNLIKELY (d == 0))
4086d7f5d3SJohn Marino     {
4186d7f5d3SJohn Marino       SIZ(r) = 0;
4286d7f5d3SJohn Marino       EXP(r) = 0;
4386d7f5d3SJohn Marino       return;
4486d7f5d3SJohn Marino     }
4586d7f5d3SJohn Marino   negative = d < 0;
4686d7f5d3SJohn Marino   d = ABS (d);
4786d7f5d3SJohn Marino 
4886d7f5d3SJohn Marino   SIZ(r) = negative ? -LIMBS_PER_DOUBLE : LIMBS_PER_DOUBLE;
4986d7f5d3SJohn Marino   EXP(r) = __gmp_extract_double (PTR(r), d);
5086d7f5d3SJohn Marino }
51