Hex Artifact Content
Not logged in

Artifact 32e778bf88cdfc55a57e8e19881150d72e33d758:


0000: 2f 2a 20 6d 70 66 72 5f 73 65 74 5f 64 20 2d 2d  /* mpfr_set_d --
0010: 20 63 6f 6e 76 65 72 74 20 61 20 6d 61 63 68 69   convert a machi
0020: 6e 65 20 64 6f 75 62 6c 65 20 70 72 65 63 69 73  ne double precis
0030: 69 6f 6e 20 66 6c 6f 61 74 20 74 6f 0a 20 20 20  ion float to.   
0040: 20 20 20 20 20 20 20 20 20 20 20 20 20 20 61 20                a 
0050: 6d 75 6c 74 69 70 6c 65 20 70 72 65 63 69 73 69  multiple precisi
0060: 6f 6e 20 66 6c 6f 61 74 69 6e 67 2d 70 6f 69 6e  on floating-poin
0070: 74 20 6e 75 6d 62 65 72 0a 0a 43 6f 70 79 72 69  t number..Copyri
0080: 67 68 74 20 31 39 39 39 2c 20 32 30 30 30 2c 20  ght 1999, 2000, 
0090: 32 30 30 31 2c 20 32 30 30 32 20 46 72 65 65 20  2001, 2002 Free 
00a0: 53 6f 66 74 77 61 72 65 20 46 6f 75 6e 64 61 74  Software Foundat
00b0: 69 6f 6e 2c 20 49 6e 63 2e 0a 0a 54 68 69 73 20  ion, Inc...This 
00c0: 66 69 6c 65 20 69 73 20 70 61 72 74 20 6f 66 20  file is part of 
00d0: 74 68 65 20 4d 50 46 52 20 4c 69 62 72 61 72 79  the MPFR Library
00e0: 2e 0a 0a 54 68 65 20 4d 50 46 52 20 4c 69 62 72  ...The MPFR Libr
00f0: 61 72 79 20 69 73 20 66 72 65 65 20 73 6f 66 74  ary is free soft
0100: 77 61 72 65 3b 20 79 6f 75 20 63 61 6e 20 72 65  ware; you can re
0110: 64 69 73 74 72 69 62 75 74 65 20 69 74 20 61 6e  distribute it an
0120: 64 2f 6f 72 20 6d 6f 64 69 66 79 0a 69 74 20 75  d/or modify.it u
0130: 6e 64 65 72 20 74 68 65 20 74 65 72 6d 73 20 6f  nder the terms o
0140: 66 20 74 68 65 20 47 4e 55 20 4c 65 73 73 65 72  f the GNU Lesser
0150: 20 47 65 6e 65 72 61 6c 20 50 75 62 6c 69 63 20   General Public 
0160: 4c 69 63 65 6e 73 65 20 61 73 20 70 75 62 6c 69  License as publi
0170: 73 68 65 64 20 62 79 0a 74 68 65 20 46 72 65 65  shed by.the Free
0180: 20 53 6f 66 74 77 61 72 65 20 46 6f 75 6e 64 61   Software Founda
0190: 74 69 6f 6e 3b 20 65 69 74 68 65 72 20 76 65 72  tion; either ver
01a0: 73 69 6f 6e 20 32 2e 31 20 6f 66 20 74 68 65 20  sion 2.1 of the 
01b0: 4c 69 63 65 6e 73 65 2c 20 6f 72 20 28 61 74 20  License, or (at 
01c0: 79 6f 75 72 0a 6f 70 74 69 6f 6e 29 20 61 6e 79  your.option) any
01d0: 20 6c 61 74 65 72 20 76 65 72 73 69 6f 6e 2e 0a   later version..
01e0: 0a 54 68 65 20 4d 50 46 52 20 4c 69 62 72 61 72  .The MPFR Librar
01f0: 79 20 69 73 20 64 69 73 74 72 69 62 75 74 65 64  y is distributed
0200: 20 69 6e 20 74 68 65 20 68 6f 70 65 20 74 68 61   in the hope tha
0210: 74 20 69 74 20 77 69 6c 6c 20 62 65 20 75 73 65  t it will be use
0220: 66 75 6c 2c 20 62 75 74 0a 57 49 54 48 4f 55 54  ful, but.WITHOUT
0230: 20 41 4e 59 20 57 41 52 52 41 4e 54 59 3b 20 77   ANY WARRANTY; w
0240: 69 74 68 6f 75 74 20 65 76 65 6e 20 74 68 65 20  ithout even the 
0250: 69 6d 70 6c 69 65 64 20 77 61 72 72 61 6e 74 79  implied warranty
0260: 20 6f 66 20 4d 45 52 43 48 41 4e 54 41 42 49 4c   of MERCHANTABIL
0270: 49 54 59 0a 6f 72 20 46 49 54 4e 45 53 53 20 46  ITY.or FITNESS F
0280: 4f 52 20 41 20 50 41 52 54 49 43 55 4c 41 52 20  OR A PARTICULAR 
0290: 50 55 52 50 4f 53 45 2e 20 20 53 65 65 20 74 68  PURPOSE.  See th
02a0: 65 20 47 4e 55 20 4c 65 73 73 65 72 20 47 65 6e  e GNU Lesser Gen
02b0: 65 72 61 6c 20 50 75 62 6c 69 63 0a 4c 69 63 65  eral Public.Lice
02c0: 6e 73 65 20 66 6f 72 20 6d 6f 72 65 20 64 65 74  nse for more det
02d0: 61 69 6c 73 2e 0a 0a 59 6f 75 20 73 68 6f 75 6c  ails...You shoul
02e0: 64 20 68 61 76 65 20 72 65 63 65 69 76 65 64 20  d have received 
02f0: 61 20 63 6f 70 79 20 6f 66 20 74 68 65 20 47 4e  a copy of the GN
0300: 55 20 4c 65 73 73 65 72 20 47 65 6e 65 72 61 6c  U Lesser General
0310: 20 50 75 62 6c 69 63 20 4c 69 63 65 6e 73 65 0a   Public License.
0320: 61 6c 6f 6e 67 20 77 69 74 68 20 74 68 65 20 4d  along with the M
0330: 50 46 52 20 4c 69 62 72 61 72 79 3b 20 73 65 65  PFR Library; see
0340: 20 74 68 65 20 66 69 6c 65 20 43 4f 50 59 49 4e   the file COPYIN
0350: 47 2e 4c 49 42 2e 20 20 49 66 20 6e 6f 74 2c 20  G.LIB.  If not, 
0360: 77 72 69 74 65 20 74 6f 0a 74 68 65 20 46 72 65  write to.the Fre
0370: 65 20 53 6f 66 74 77 61 72 65 20 46 6f 75 6e 64  e Software Found
0380: 61 74 69 6f 6e 2c 20 49 6e 63 2e 2c 20 35 39 20  ation, Inc., 59 
0390: 54 65 6d 70 6c 65 20 50 6c 61 63 65 20 2d 20 53  Temple Place - S
03a0: 75 69 74 65 20 33 33 30 2c 20 42 6f 73 74 6f 6e  uite 330, Boston
03b0: 2c 0a 4d 41 20 30 32 31 31 31 2d 31 33 30 37 2c  ,.MA 02111-1307,
03c0: 20 55 53 41 2e 20 2a 2f 0a 0a 23 69 6e 63 6c 75   USA. */..#inclu
03d0: 64 65 20 22 67 6d 70 2e 68 22 0a 23 69 6e 63 6c  de "gmp.h".#incl
03e0: 75 64 65 20 22 67 6d 70 2d 69 6d 70 6c 2e 68 22  ude "gmp-impl.h"
03f0: 0a 23 69 6e 63 6c 75 64 65 20 22 6c 6f 6e 67 6c  .#include "longl
0400: 6f 6e 67 2e 68 22 0a 23 69 6e 63 6c 75 64 65 20  ong.h".#include 
0410: 22 6d 70 66 72 2e 68 22 0a 23 69 6e 63 6c 75 64  "mpfr.h".#includ
0420: 65 20 22 6d 70 66 72 2d 69 6d 70 6c 2e 68 22 0a  e "mpfr-impl.h".
0430: 0a 23 69 66 20 28 42 49 54 53 5f 50 45 52 5f 4d  .#if (BITS_PER_M
0440: 50 5f 4c 49 4d 42 3d 3d 33 32 29 0a 23 64 65 66  P_LIMB==32).#def
0450: 69 6e 65 20 4d 50 46 52 5f 4c 49 4d 42 53 5f 50  ine MPFR_LIMBS_P
0460: 45 52 5f 44 4f 55 42 4c 45 20 32 0a 23 65 6c 69  ER_DOUBLE 2.#eli
0470: 66 20 28 42 49 54 53 5f 50 45 52 5f 4d 50 5f 4c  f (BITS_PER_MP_L
0480: 49 4d 42 20 3e 3d 20 36 34 29 0a 23 64 65 66 69  IMB >= 64).#defi
0490: 6e 65 20 4d 50 46 52 5f 4c 49 4d 42 53 5f 50 45  ne MPFR_LIMBS_PE
04a0: 52 5f 44 4f 55 42 4c 45 20 31 0a 23 65 6c 69 66  R_DOUBLE 1.#elif
04b0: 20 28 42 49 54 53 5f 50 45 52 5f 4d 50 5f 4c 49   (BITS_PER_MP_LI
04c0: 4d 42 20 3d 3d 20 31 36 29 0a 23 64 65 66 69 6e  MB == 16).#defin
04d0: 65 20 4d 50 46 52 5f 4c 49 4d 42 53 5f 50 45 52  e MPFR_LIMBS_PER
04e0: 5f 44 4f 55 42 4c 45 20 34 0a 23 65 6e 64 69 66  _DOUBLE 4.#endif
04f0: 0a 0a 73 74 61 74 69 63 20 69 6e 74 20 5f 5f 6d  ..static int __m
0500: 70 66 72 5f 65 78 74 72 61 63 74 5f 64 6f 75 62  pfr_extract_doub
0510: 6c 65 20 5f 50 52 4f 54 4f 20 28 28 6d 70 5f 70  le _PROTO ((mp_p
0520: 74 72 2c 20 64 6f 75 62 6c 65 29 29 3b 0a 0a 2f  tr, double));../
0530: 2a 20 49 6e 63 6c 75 64 65 64 20 66 72 6f 6d 20  * Included from 
0540: 67 6d 70 2d 32 2e 30 2e 32 2c 20 70 61 74 63 68  gmp-2.0.2, patch
0550: 65 64 20 74 6f 20 73 75 70 70 6f 72 74 20 64 65  ed to support de
0560: 6e 6f 72 6d 73 20 2a 2f 0a 0a 23 69 66 64 65 66  norms */..#ifdef
0570: 20 58 44 45 42 55 47 0a 23 75 6e 64 65 66 20 5f   XDEBUG.#undef _
0580: 47 4d 50 5f 49 45 45 45 5f 46 4c 4f 41 54 53 0a  GMP_IEEE_FLOATS.
0590: 23 65 6e 64 69 66 0a 0a 23 69 66 6e 64 65 66 20  #endif..#ifndef 
05a0: 5f 47 4d 50 5f 49 45 45 45 5f 46 4c 4f 41 54 53  _GMP_IEEE_FLOATS
05b0: 0a 23 64 65 66 69 6e 65 20 5f 47 4d 50 5f 49 45  .#define _GMP_IE
05c0: 45 45 5f 46 4c 4f 41 54 53 20 30 0a 23 65 6e 64  EE_FLOATS 0.#end
05d0: 69 66 0a 0a 73 74 61 74 69 63 20 69 6e 74 0a 5f  if..static int._
05e0: 5f 6d 70 66 72 5f 65 78 74 72 61 63 74 5f 64 6f  _mpfr_extract_do
05f0: 75 62 6c 65 20 28 6d 70 5f 70 74 72 20 72 70 2c  uble (mp_ptr rp,
0600: 20 64 6f 75 62 6c 65 20 64 29 0a 20 20 20 20 20   double d).     
0610: 2f 2a 20 65 3d 30 20 69 66 66 20 42 49 54 53 5f  /* e=0 iff BITS_
0620: 50 45 52 5f 4d 50 5f 4c 49 4d 42 3d 33 32 20 61  PER_MP_LIMB=32 a
0630: 6e 64 20 72 70 20 68 61 73 20 6f 6e 6c 79 20 6f  nd rp has only o
0640: 6e 65 20 6c 69 6d 62 20 2a 2f 0a 7b 0a 20 20 6c  ne limb */.{.  l
0650: 6f 6e 67 20 65 78 70 3b 0a 20 20 6d 70 5f 6c 69  ong exp;.  mp_li
0660: 6d 62 5f 74 20 6d 61 6e 6c 3b 0a 23 69 66 20 42  mb_t manl;.#if B
0670: 49 54 53 5f 50 45 52 5f 4d 50 5f 4c 49 4d 42 20  ITS_PER_MP_LIMB 
0680: 3d 3d 20 33 32 0a 20 20 6d 70 5f 6c 69 6d 62 5f  == 32.  mp_limb_
0690: 74 20 6d 61 6e 68 3b 0a 23 65 6e 64 69 66 0a 0a  t manh;.#endif..
06a0: 20 20 2f 2a 20 42 55 47 53 0a 0a 20 20 20 20 20    /* BUGS..     
06b0: 31 2e 20 53 68 6f 75 6c 64 20 68 61 6e 64 6c 65  1. Should handle
06c0: 20 49 6e 66 20 61 6e 64 20 4e 61 4e 20 69 6e 20   Inf and NaN in 
06d0: 49 45 45 45 20 73 70 65 63 69 66 69 63 20 63 6f  IEEE specific co
06e0: 64 65 2e 0a 20 20 20 20 20 32 2e 20 48 61 6e 64  de..     2. Hand
06f0: 6c 65 20 49 6e 66 20 61 6e 64 20 4e 61 4e 20 61  le Inf and NaN a
0700: 6c 73 6f 20 69 6e 20 64 65 66 61 75 6c 74 20 63  lso in default c
0710: 6f 64 65 2c 20 74 6f 20 61 76 6f 69 64 20 68 61  ode, to avoid ha
0720: 6e 67 73 2e 0a 20 20 20 20 20 33 2e 20 47 65 6e  ngs..     3. Gen
0730: 65 72 61 6c 69 7a 65 20 74 6f 20 68 61 6e 64 6c  eralize to handl
0740: 65 20 61 6c 6c 20 42 49 54 53 5f 50 45 52 5f 4d  e all BITS_PER_M
0750: 50 5f 4c 49 4d 42 20 3e 3d 20 33 32 2e 0a 20 20  P_LIMB >= 32..  
0760: 20 20 20 34 2e 20 54 68 69 73 20 6c 69 74 73 20     4. This lits 
0770: 69 73 20 69 6e 63 6f 6d 70 6c 65 74 65 20 61 6e  is incomplete an
0780: 64 20 6d 69 73 73 70 65 6c 6c 65 64 2e 0a 20 20  d misspelled..  
0790: 20 2a 2f 0a 0a 20 20 69 66 20 28 64 20 3d 3d 20   */..  if (d == 
07a0: 30 2e 30 29 0a 20 20 20 20 7b 0a 20 20 20 20 20  0.0).    {.     
07b0: 20 72 70 5b 30 5d 20 3d 20 30 3b 0a 20 20 20 20   rp[0] = 0;.    
07c0: 20 20 72 65 74 75 72 6e 20 30 3b 0a 20 20 20 20    return 0;.    
07d0: 7d 0a 0a 23 69 66 20 5f 47 4d 50 5f 49 45 45 45  }..#if _GMP_IEEE
07e0: 5f 46 4c 4f 41 54 53 0a 20 20 7b 0a 20 20 20 20  _FLOATS.  {.    
07f0: 75 6e 69 6f 6e 20 69 65 65 65 5f 64 6f 75 62 6c  union ieee_doubl
0800: 65 5f 65 78 74 72 61 63 74 20 78 3b 0a 20 20 20  e_extract x;.   
0810: 20 78 2e 64 20 3d 20 64 3b 0a 0a 20 20 20 20 65   x.d = d;..    e
0820: 78 70 20 3d 20 78 2e 73 2e 65 78 70 3b 0a 20 20  xp = x.s.exp;.  
0830: 20 20 69 66 20 28 65 78 70 29 0a 20 20 20 20 20    if (exp).     
0840: 20 7b 0a 23 69 66 20 42 49 54 53 5f 50 45 52 5f   {.#if BITS_PER_
0850: 4d 50 5f 4c 49 4d 42 20 3d 3d 20 36 34 0a 09 6d  MP_LIMB == 64..m
0860: 61 6e 6c 20 3d 20 28 28 4d 50 5f 4c 49 4d 42 5f  anl = ((MP_LIMB_
0870: 54 5f 4f 4e 45 20 3c 3c 20 36 33 29 0a 09 09 7c  T_ONE << 63)...|
0880: 20 28 28 6d 70 5f 6c 69 6d 62 5f 74 29 20 78 2e   ((mp_limb_t) x.
0890: 73 2e 6d 61 6e 68 20 3c 3c 20 34 33 29 20 7c 20  s.manh << 43) | 
08a0: 28 28 6d 70 5f 6c 69 6d 62 5f 74 29 20 78 2e 73  ((mp_limb_t) x.s
08b0: 2e 6d 61 6e 6c 20 3c 3c 20 31 31 29 29 3b 0a 23  .manl << 11));.#
08c0: 65 6c 73 65 0a 09 6d 61 6e 68 20 3d 20 28 4d 50  else..manh = (MP
08d0: 5f 4c 49 4d 42 5f 54 5f 4f 4e 45 20 3c 3c 20 33  _LIMB_T_ONE << 3
08e0: 31 29 20 7c 20 28 78 2e 73 2e 6d 61 6e 68 20 3c  1) | (x.s.manh <
08f0: 3c 20 31 31 29 20 7c 20 28 78 2e 73 2e 6d 61 6e  < 11) | (x.s.man
0900: 6c 20 3e 3e 20 32 31 29 3b 0a 09 6d 61 6e 6c 20  l >> 21);..manl 
0910: 3d 20 78 2e 73 2e 6d 61 6e 6c 20 3c 3c 20 31 31  = x.s.manl << 11
0920: 3b 0a 23 65 6e 64 69 66 0a 20 20 20 20 20 20 7d  ;.#endif.      }
0930: 0a 20 20 20 20 65 6c 73 65 20 2f 2a 20 64 65 6e  .    else /* den
0940: 6f 72 6d 61 6c 69 7a 65 64 20 6e 75 6d 62 65 72  ormalized number
0950: 20 2a 2f 0a 20 20 20 20 20 20 7b 0a 23 69 66 20   */.      {.#if 
0960: 42 49 54 53 5f 50 45 52 5f 4d 50 5f 4c 49 4d 42  BITS_PER_MP_LIMB
0970: 20 3d 3d 20 36 34 0a 09 6d 61 6e 6c 20 3d 20 28   == 64..manl = (
0980: 28 6d 70 5f 6c 69 6d 62 5f 74 29 20 78 2e 73 2e  (mp_limb_t) x.s.
0990: 6d 61 6e 68 20 3c 3c 20 34 33 29 20 7c 20 28 28  manh << 43) | ((
09a0: 6d 70 5f 6c 69 6d 62 5f 74 29 20 78 2e 73 2e 6d  mp_limb_t) x.s.m
09b0: 61 6e 6c 20 3c 3c 20 31 31 29 3b 0a 23 65 6c 73  anl << 11);.#els
09c0: 65 0a 20 20 20 20 20 20 20 20 6d 61 6e 68 20 3d  e.        manh =
09d0: 20 28 78 2e 73 2e 6d 61 6e 68 20 3c 3c 20 31 31   (x.s.manh << 11
09e0: 29 20 2f 2a 20 68 69 67 68 20 32 31 20 62 69 74  ) /* high 21 bit
09f0: 73 20 2a 2f 0a 20 20 20 20 20 20 20 20 20 20 7c  s */.          |
0a00: 20 28 78 2e 73 2e 6d 61 6e 6c 20 3e 3e 20 32 31   (x.s.manl >> 21
0a10: 29 3b 20 2f 2a 20 6d 69 64 64 6c 65 20 31 31 20  ); /* middle 11 
0a20: 62 69 74 73 20 2a 2f 0a 09 6d 61 6e 6c 20 3d 20  bits */..manl = 
0a30: 78 2e 73 2e 6d 61 6e 6c 20 3c 3c 20 31 31 3b 20  x.s.manl << 11; 
0a40: 2f 2a 20 6c 6f 77 20 32 31 20 62 69 74 73 20 2a  /* low 21 bits *
0a50: 2f 0a 23 65 6e 64 69 66 0a 20 20 20 20 20 20 7d  /.#endif.      }
0a60: 0a 20 20 7d 0a 23 65 6c 73 65 0a 20 20 7b 0a 20  .  }.#else.  {. 
0a70: 20 20 20 2f 2a 20 55 6e 6b 6e 6f 77 6e 20 28 6f     /* Unknown (o
0a80: 72 20 6b 6e 6f 77 6e 20 74 6f 20 62 65 20 6e 6f  r known to be no
0a90: 6e 2d 49 45 45 45 29 20 64 6f 75 62 6c 65 20 66  n-IEEE) double f
0aa0: 6f 72 6d 61 74 2e 20 20 2a 2f 0a 20 20 20 20 65  ormat.  */.    e
0ab0: 78 70 20 3d 20 30 3b 0a 20 20 20 20 69 66 20 28  xp = 0;.    if (
0ac0: 64 20 3e 3d 20 31 2e 30 29 0a 20 20 20 20 20 20  d >= 1.0).      
0ad0: 7b 0a 20 20 20 20 20 20 20 20 69 66 20 28 64 20  {.        if (d 
0ae0: 2a 20 30 2e 35 20 3d 3d 20 64 29 0a 20 20 20 20  * 0.5 == d).    
0af0: 20 20 20 20 20 20 61 62 6f 72 74 20 28 29 3b 0a        abort ();.
0b00: 0a 20 20 20 20 20 20 20 20 77 68 69 6c 65 20 28  .        while (
0b10: 64 20 3e 3d 20 33 32 37 36 38 2e 30 29 0a 20 20  d >= 32768.0).  
0b20: 20 20 20 20 20 20 20 20 7b 0a 20 20 20 20 20 20          {.      
0b30: 20 20 20 20 20 20 64 20 2a 3d 20 28 31 2e 30 20        d *= (1.0 
0b40: 2f 20 36 35 35 33 36 2e 30 29 3b 0a 20 20 20 20  / 65536.0);.    
0b50: 20 20 20 20 20 20 20 20 65 78 70 20 2b 3d 20 31          exp += 1
0b60: 36 3b 0a 20 20 20 20 20 20 20 20 20 20 7d 0a 20  6;.          }. 
0b70: 20 20 20 20 20 20 20 77 68 69 6c 65 20 28 64 20         while (d 
0b80: 3e 3d 20 31 2e 30 29 0a 20 20 20 20 20 20 20 20  >= 1.0).        
0b90: 20 20 7b 0a 20 20 20 20 20 20 20 20 20 20 20 20    {.            
0ba0: 64 20 2a 3d 20 30 2e 35 3b 0a 20 20 20 20 20 20  d *= 0.5;.      
0bb0: 20 20 20 20 20 20 65 78 70 20 2b 3d 20 31 3b 0a        exp += 1;.
0bc0: 20 20 20 20 20 20 20 20 20 20 7d 0a 20 20 20 20            }.    
0bd0: 20 20 7d 0a 20 20 20 20 65 6c 73 65 20 69 66 20    }.    else if 
0be0: 28 64 20 3c 20 30 2e 35 29 0a 20 20 20 20 20 20  (d < 0.5).      
0bf0: 7b 0a 20 20 20 20 20 20 20 20 77 68 69 6c 65 20  {.        while 
0c00: 28 64 20 3c 20 28 31 2e 30 20 2f 20 36 35 35 33  (d < (1.0 / 6553
0c10: 36 2e 30 29 29 0a 20 20 20 20 20 20 20 20 20 20  6.0)).          
0c20: 7b 0a 20 20 20 20 20 20 20 20 20 20 20 20 64 20  {.            d 
0c30: 2a 3d 20 20 36 35 35 33 36 2e 30 3b 0a 20 20 20  *=  65536.0;.   
0c40: 20 20 20 20 20 20 20 20 20 65 78 70 20 2d 3d 20           exp -= 
0c50: 31 36 3b 0a 20 20 20 20 20 20 20 20 20 20 7d 0a  16;.          }.
0c60: 20 20 20 20 20 20 20 20 77 68 69 6c 65 20 28 64          while (d
0c70: 20 3c 20 30 2e 35 29 0a 20 20 20 20 20 20 20 20   < 0.5).        
0c80: 20 20 7b 0a 20 20 20 20 20 20 20 20 20 20 20 20    {.            
0c90: 64 20 2a 3d 20 32 2e 30 3b 0a 20 20 20 20 20 20  d *= 2.0;.      
0ca0: 20 20 20 20 20 20 65 78 70 20 2d 3d 20 31 3b 0a        exp -= 1;.
0cb0: 20 20 20 20 20 20 20 20 20 20 7d 0a 20 20 20 20            }.    
0cc0: 20 20 7d 0a 0a 20 20 20 20 64 20 2a 3d 20 4d 50    }..    d *= MP
0cd0: 5f 42 41 53 45 5f 41 53 5f 44 4f 55 42 4c 45 3b  _BASE_AS_DOUBLE;
0ce0: 0a 23 69 66 20 42 49 54 53 5f 50 45 52 5f 4d 50  .#if BITS_PER_MP
0cf0: 5f 4c 49 4d 42 20 3d 3d 20 36 34 0a 20 20 20 20  _LIMB == 64.    
0d00: 6d 61 6e 6c 20 3d 20 64 3b 0a 23 65 6c 73 65 0a  manl = d;.#else.
0d10: 20 20 20 20 6d 61 6e 68 20 3d 20 64 3b 0a 20 20      manh = d;.  
0d20: 20 20 6d 61 6e 6c 20 3d 20 28 64 20 2d 20 6d 61    manl = (d - ma
0d30: 6e 68 29 20 2a 20 4d 50 5f 42 41 53 45 5f 41 53  nh) * MP_BASE_AS
0d40: 5f 44 4f 55 42 4c 45 3b 0a 23 65 6e 64 69 66 0a  _DOUBLE;.#endif.
0d50: 0a 20 20 20 20 65 78 70 20 2b 3d 20 31 30 32 32  .    exp += 1022
0d60: 3b 0a 20 20 7d 0a 23 65 6e 64 69 66 0a 0a 20 20  ;.  }.#endif..  
0d70: 69 66 20 28 65 78 70 29 20 65 78 70 20 3d 20 28  if (exp) exp = (
0d80: 75 6e 73 69 67 6e 65 64 29 20 65 78 70 20 2d 20  unsigned) exp - 
0d90: 31 30 32 32 3b 20 65 6c 73 65 20 65 78 70 20 3d  1022; else exp =
0da0: 20 2d 31 30 32 31 3b 0a 0a 23 69 66 20 42 49 54   -1021;..#if BIT
0db0: 53 5f 50 45 52 5f 4d 50 5f 4c 49 4d 42 20 3d 3d  S_PER_MP_LIMB ==
0dc0: 20 36 34 0a 20 20 72 70 5b 30 5d 20 3d 20 6d 61   64.  rp[0] = ma
0dd0: 6e 6c 3b 0a 23 65 6c 73 65 0a 20 20 72 70 5b 31  nl;.#else.  rp[1
0de0: 5d 20 3d 20 6d 61 6e 68 3b 0a 20 20 72 70 5b 30  ] = manh;.  rp[0
0df0: 5d 20 3d 20 6d 61 6e 6c 3b 0a 23 65 6e 64 69 66  ] = manl;.#endif
0e00: 0a 0a 20 20 72 65 74 75 72 6e 20 65 78 70 3b 0a  ..  return exp;.
0e10: 7d 0a 0a 2f 2a 20 45 6e 64 20 6f 66 20 70 61 72  }../* End of par
0e20: 74 20 69 6e 63 6c 75 64 65 64 20 66 72 6f 6d 20  t included from 
0e30: 67 6d 70 2d 32 2e 30 2e 32 20 2a 2f 0a 0a 69 6e  gmp-2.0.2 */..in
0e40: 74 0a 6d 70 66 72 5f 73 65 74 5f 64 20 28 6d 70  t.mpfr_set_d (mp
0e50: 66 72 5f 70 74 72 20 72 2c 20 64 6f 75 62 6c 65  fr_ptr r, double
0e60: 20 64 2c 20 6d 70 5f 72 6e 64 5f 74 20 72 6e 64   d, mp_rnd_t rnd
0e70: 5f 6d 6f 64 65 29 0a 7b 0a 20 20 69 6e 74 20 73  _mode).{.  int s
0e80: 69 67 6e 64 2c 20 73 69 7a 65 74 6d 70 2c 20 69  ignd, sizetmp, i
0e90: 6e 65 78 61 63 74 3b 0a 20 20 75 6e 73 69 67 6e  nexact;.  unsign
0ea0: 65 64 20 69 6e 74 20 63 6e 74 2c 20 6b 3b 0a 20  ed int cnt, k;. 
0eb0: 20 6d 70 66 72 5f 70 74 72 20 74 6d 70 3b 0a 20   mpfr_ptr tmp;. 
0ec0: 20 54 4d 50 5f 44 45 43 4c 28 6d 61 72 6b 65 72   TMP_DECL(marker
0ed0: 29 3b 0a 0a 20 20 54 4d 50 5f 4d 41 52 4b 28 6d  );..  TMP_MARK(m
0ee0: 61 72 6b 65 72 29 3b 0a 20 20 4d 50 46 52 5f 43  arker);.  MPFR_C
0ef0: 4c 45 41 52 5f 46 4c 41 47 53 28 72 29 3b 0a 0a  LEAR_FLAGS(r);..
0f00: 20 20 69 66 20 28 64 20 3d 3d 20 30 29 0a 20 20    if (d == 0).  
0f10: 20 20 7b 0a 20 20 20 20 20 20 75 6e 69 6f 6e 20    {.      union 
0f20: 69 65 65 65 5f 64 6f 75 62 6c 65 5f 65 78 74 72  ieee_double_extr
0f30: 61 63 74 20 78 3b 0a 0a 20 20 20 20 20 20 4d 50  act x;..      MP
0f40: 46 52 5f 53 45 54 5f 5a 45 52 4f 28 72 29 3b 0a  FR_SET_ZERO(r);.
0f50: 20 20 20 20 20 20 2f 2a 20 73 65 74 20 63 6f 72        /* set cor
0f60: 72 65 63 74 20 73 69 67 6e 20 2a 2f 0a 20 20 20  rect sign */.   
0f70: 20 20 20 78 2e 64 20 3d 20 64 3b 0a 20 20 20 20     x.d = d;.    
0f80: 20 20 69 66 20 28 28 28 78 2e 73 2e 73 69 67 20    if (((x.s.sig 
0f90: 3d 3d 20 31 29 20 26 26 20 28 4d 50 46 52 5f 53  == 1) && (MPFR_S
0fa0: 49 47 4e 28 72 29 20 3e 20 30 29 29 0a 09 20 20  IGN(r) > 0))..  
0fb0: 7c 7c 20 28 28 78 2e 73 2e 73 69 67 20 3d 3d 20  || ((x.s.sig == 
0fc0: 30 29 20 26 26 20 28 4d 50 46 52 5f 53 49 47 4e  0) && (MPFR_SIGN
0fd0: 28 72 29 20 3c 20 30 29 29 29 0a 09 4d 50 46 52  (r) < 0)))..MPFR
0fe0: 5f 43 48 41 4e 47 45 5f 53 49 47 4e 28 72 29 3b  _CHANGE_SIGN(r);
0ff0: 0a 20 20 20 20 20 20 72 65 74 75 72 6e 20 30 3b  .      return 0;
1000: 20 2f 2a 20 30 20 69 73 20 65 78 61 63 74 20 2a   /* 0 is exact *
1010: 2f 0a 20 20 20 20 7d 0a 20 20 65 6c 73 65 20 69  /.    }.  else i
1020: 66 20 28 44 4f 55 42 4c 45 5f 49 53 4e 41 4e 28  f (DOUBLE_ISNAN(
1030: 64 29 29 0a 20 20 20 20 7b 0a 20 20 20 20 20 20  d)).    {.      
1040: 4d 50 46 52 5f 53 45 54 5f 4e 41 4e 28 72 29 3b  MPFR_SET_NAN(r);
1050: 0a 20 20 20 20 20 20 4d 50 46 52 5f 52 45 54 5f  .      MPFR_RET_
1060: 4e 41 4e 3b 0a 20 20 20 20 7d 0a 20 20 65 6c 73  NAN;.    }.  els
1070: 65 20 69 66 20 28 44 4f 55 42 4c 45 5f 49 53 49  e if (DOUBLE_ISI
1080: 4e 46 28 64 29 29 0a 20 20 20 20 7b 0a 20 20 20  NF(d)).    {.   
1090: 20 20 20 4d 50 46 52 5f 53 45 54 5f 49 4e 46 28     MPFR_SET_INF(
10a0: 72 29 3b 0a 20 20 20 20 20 20 69 66 20 28 28 64  r);.      if ((d
10b0: 20 3e 20 30 20 26 26 20 28 4d 50 46 52 5f 53 49   > 0 && (MPFR_SI
10c0: 47 4e 28 72 29 20 3d 3d 20 2d 31 29 29 20 7c 7c  GN(r) == -1)) ||
10d0: 20 28 64 20 3c 20 30 20 26 26 20 28 4d 50 46 52   (d < 0 && (MPFR
10e0: 5f 53 49 47 4e 28 72 29 20 3d 3d 20 31 29 29 29  _SIGN(r) == 1)))
10f0: 0a 09 4d 50 46 52 5f 43 48 41 4e 47 45 5f 53 49  ..MPFR_CHANGE_SI
1100: 47 4e 28 72 29 3b 0a 20 20 20 20 20 20 72 65 74  GN(r);.      ret
1110: 75 72 6e 20 30 3b 20 2f 2a 20 69 6e 66 69 6e 69  urn 0; /* infini
1120: 74 79 20 69 73 20 65 78 61 63 74 20 2a 2f 0a 20  ty is exact */. 
1130: 20 20 20 7d 0a 0a 20 20 2f 2a 20 77 61 72 6e 69     }..  /* warni
1140: 6e 67 3a 20 64 6f 6e 27 74 20 75 73 65 20 74 6d  ng: don't use tm
1150: 70 3d 72 20 68 65 72 65 2c 20 65 76 65 6e 20 69  p=r here, even i
1160: 66 20 53 49 5a 45 28 72 29 20 3e 3d 20 4d 50 46  f SIZE(r) >= MPF
1170: 52 5f 4c 49 4d 42 53 5f 50 45 52 5f 44 4f 55 42  R_LIMBS_PER_DOUB
1180: 4c 45 2c 0a 20 20 20 20 20 73 69 6e 63 65 20 50  LE,.     since P
1190: 52 45 43 28 72 29 20 6d 61 79 20 62 65 20 64 69  REC(r) may be di
11a0: 66 66 65 72 65 6e 74 20 66 72 6f 6d 20 50 52 45  fferent from PRE
11b0: 43 28 74 6d 70 29 2c 20 61 6e 64 20 74 68 65 6e  C(tmp), and then
11c0: 20 62 6f 74 68 20 76 61 72 69 61 62 6c 65 73 0a   both variables.
11d0: 20 20 20 20 20 77 6f 75 6c 64 20 68 61 76 65 20       would have 
11e0: 73 61 6d 65 20 70 72 65 63 69 73 69 6f 6e 20 69  same precision i
11f0: 6e 20 74 68 65 20 6d 70 66 72 5f 73 65 74 34 20  n the mpfr_set4 
1200: 63 61 6c 6c 20 62 65 6c 6f 77 2e 20 2a 2f 0a 20  call below. */. 
1210: 20 74 6d 70 20 3d 20 28 6d 70 66 72 5f 70 74 72   tmp = (mpfr_ptr
1220: 29 20 54 4d 50 5f 41 4c 4c 4f 43 28 73 69 7a 65  ) TMP_ALLOC(size
1230: 6f 66 28 6d 70 66 72 5f 74 29 29 3b 0a 20 20 4d  of(mpfr_t));.  M
1240: 50 46 52 5f 4d 41 4e 54 28 74 6d 70 29 20 3d 20  PFR_MANT(tmp) = 
1250: 54 4d 50 5f 41 4c 4c 4f 43 28 4d 50 46 52 5f 4c  TMP_ALLOC(MPFR_L
1260: 49 4d 42 53 5f 50 45 52 5f 44 4f 55 42 4c 45 20  IMBS_PER_DOUBLE 
1270: 2a 20 42 59 54 45 53 5f 50 45 52 5f 4d 50 5f 4c  * BYTES_PER_MP_L
1280: 49 4d 42 29 3b 0a 20 20 4d 50 46 52 5f 50 52 45  IMB);.  MPFR_PRE
1290: 43 28 74 6d 70 29 20 3d 20 35 33 3b 0a 20 20 4d  C(tmp) = 53;.  M
12a0: 50 46 52 5f 53 49 5a 45 28 74 6d 70 29 20 3d 20  PFR_SIZE(tmp) = 
12b0: 4d 50 46 52 5f 4c 49 4d 42 53 5f 50 45 52 5f 44  MPFR_LIMBS_PER_D
12c0: 4f 55 42 4c 45 3b 0a 20 20 73 69 7a 65 74 6d 70  OUBLE;.  sizetmp
12d0: 20 3d 20 4d 50 46 52 5f 4c 49 4d 42 53 5f 50 45   = MPFR_LIMBS_PE
12e0: 52 5f 44 4f 55 42 4c 45 3b 0a 0a 20 20 73 69 67  R_DOUBLE;..  sig
12f0: 6e 64 20 3d 20 28 64 20 3c 20 30 29 20 3f 20 2d  nd = (d < 0) ? -
1300: 31 20 3a 20 31 3b 0a 20 20 64 20 3d 20 41 42 53  1 : 1;.  d = ABS
1310: 20 28 64 29 3b 0a 0a 20 20 4d 50 46 52 5f 45 58   (d);..  MPFR_EX
1320: 50 28 74 6d 70 29 20 3d 20 5f 5f 6d 70 66 72 5f  P(tmp) = __mpfr_
1330: 65 78 74 72 61 63 74 5f 64 6f 75 62 6c 65 20 28  extract_double (
1340: 4d 50 46 52 5f 4d 41 4e 54 28 74 6d 70 29 2c 20  MPFR_MANT(tmp), 
1350: 64 29 3b 0a 0a 20 20 2f 2a 20 64 65 74 65 72 6d  d);..  /* determ
1360: 69 6e 65 20 6e 75 6d 62 65 72 20 6b 20 6f 66 20  ine number k of 
1370: 7a 65 72 6f 20 68 69 67 68 20 6c 69 6d 62 73 20  zero high limbs 
1380: 2a 2f 0a 20 20 66 6f 72 20 28 6b 20 3d 20 30 3b  */.  for (k = 0;
1390: 20 6b 20 3c 20 73 69 7a 65 74 6d 70 20 26 26 20   k < sizetmp && 
13a0: 4d 50 46 52 5f 4d 41 4e 54 28 74 6d 70 29 5b 73  MPFR_MANT(tmp)[s
13b0: 69 7a 65 74 6d 70 20 2d 20 31 20 2d 20 6b 5d 20  izetmp - 1 - k] 
13c0: 3d 3d 20 30 3b 20 6b 2b 2b 29 3b 0a 0a 20 20 63  == 0; k++);..  c
13d0: 6f 75 6e 74 5f 6c 65 61 64 69 6e 67 5f 7a 65 72  ount_leading_zer
13e0: 6f 73 20 28 63 6e 74 2c 20 4d 50 46 52 5f 4d 41  os (cnt, MPFR_MA
13f0: 4e 54 28 74 6d 70 29 5b 73 69 7a 65 74 6d 70 20  NT(tmp)[sizetmp 
1400: 2d 20 31 20 2d 20 6b 5d 29 3b 0a 0a 20 20 69 66  - 1 - k]);..  if
1410: 20 28 63 6e 74 29 0a 20 20 20 20 6d 70 6e 5f 6c   (cnt).    mpn_l
1420: 73 68 69 66 74 20 28 4d 50 46 52 5f 4d 41 4e 54  shift (MPFR_MANT
1430: 28 74 6d 70 29 20 2b 20 6b 2c 20 4d 50 46 52 5f  (tmp) + k, MPFR_
1440: 4d 41 4e 54 28 74 6d 70 29 2c 20 73 69 7a 65 74  MANT(tmp), sizet
1450: 6d 70 20 2d 20 6b 2c 20 63 6e 74 29 3b 0a 20 20  mp - k, cnt);.  
1460: 65 6c 73 65 20 69 66 20 28 6b 29 0a 20 20 20 20  else if (k).    
1470: 4d 50 4e 5f 43 4f 50 59 20 28 4d 50 46 52 5f 4d  MPN_COPY (MPFR_M
1480: 41 4e 54 28 74 6d 70 29 20 2b 20 6b 2c 20 4d 50  ANT(tmp) + k, MP
1490: 46 52 5f 4d 41 4e 54 28 74 6d 70 29 2c 20 73 69  FR_MANT(tmp), si
14a0: 7a 65 74 6d 70 20 2d 20 6b 29 3b 0a 20 20 69 66  zetmp - k);.  if
14b0: 20 28 6b 29 0a 20 20 20 20 4d 50 4e 5f 5a 45 52   (k).    MPN_ZER
14c0: 4f 20 28 4d 50 46 52 5f 4d 41 4e 54 28 74 6d 70  O (MPFR_MANT(tmp
14d0: 29 2c 20 6b 29 3b 0a 0a 20 20 4d 50 46 52 5f 45  ), k);..  MPFR_E
14e0: 58 50 28 74 6d 70 29 20 2d 3d 20 63 6e 74 20 2b  XP(tmp) -= cnt +
14f0: 20 6b 20 2a 20 42 49 54 53 5f 50 45 52 5f 4d 50   k * BITS_PER_MP
1500: 5f 4c 49 4d 42 3b 0a 0a 20 20 2f 2a 20 74 6d 70  _LIMB;..  /* tmp
1510: 20 69 73 20 65 78 61 63 74 20 73 69 6e 63 65 20   is exact since 
1520: 50 52 45 43 28 74 6d 70 29 3d 35 33 20 2a 2f 0a  PREC(tmp)=53 */.
1530: 20 20 69 6e 65 78 61 63 74 20 3d 20 6d 70 66 72    inexact = mpfr
1540: 5f 73 65 74 34 20 28 72 2c 20 74 6d 70 2c 20 72  _set4 (r, tmp, r
1550: 6e 64 5f 6d 6f 64 65 2c 20 73 69 67 6e 64 29 3b  nd_mode, signd);
1560: 0a 0a 20 20 54 4d 50 5f 46 52 45 45 28 6d 61 72  ..  TMP_FREE(mar
1570: 6b 65 72 29 3b 0a 20 20 72 65 74 75 72 6e 20 69  ker);.  return i
1580: 6e 65 78 61 63 74 3b 0a 7d 0a                    nexact;.}.