Hex Artifact Content
Not logged in

Artifact 396efc44e25a3e88889041d6a727c4c7d055b849:


0000: 2f 2a 20 6d 70 66 5f 64 69 76 5f 32 65 78 70 20  /* mpf_div_2exp 
0010: 2d 2d 20 44 69 76 69 64 65 20 61 20 66 6c 6f 61  -- Divide a floa
0020: 74 20 62 79 20 32 5e 6e 2e 0a 0a 43 6f 70 79 72  t by 2^n...Copyr
0030: 69 67 68 74 20 31 39 39 33 2c 20 31 39 39 34 2c  ight 1993, 1994,
0040: 20 31 39 39 36 2c 20 32 30 30 30 2c 20 32 30 30   1996, 2000, 200
0050: 31 2c 20 32 30 30 32 20 46 72 65 65 20 53 6f 66  1, 2002 Free Sof
0060: 74 77 61 72 65 20 46 6f 75 6e 64 61 74 69 6f 6e  tware Foundation
0070: 2c 20 49 6e 63 2e 0a 0a 54 68 69 73 20 66 69 6c  , Inc...This fil
0080: 65 20 69 73 20 70 61 72 74 20 6f 66 20 74 68 65  e is part of the
0090: 20 47 4e 55 20 4d 50 20 4c 69 62 72 61 72 79 2e   GNU MP Library.
00a0: 0a 0a 54 68 65 20 47 4e 55 20 4d 50 20 4c 69 62  ..The GNU MP Lib
00b0: 72 61 72 79 20 69 73 20 66 72 65 65 20 73 6f 66  rary is free sof
00c0: 74 77 61 72 65 3b 20 79 6f 75 20 63 61 6e 20 72  tware; you can r
00d0: 65 64 69 73 74 72 69 62 75 74 65 20 69 74 20 61  edistribute it a
00e0: 6e 64 2f 6f 72 20 6d 6f 64 69 66 79 0a 69 74 20  nd/or modify.it 
00f0: 75 6e 64 65 72 20 74 68 65 20 74 65 72 6d 73 20  under the terms 
0100: 6f 66 20 74 68 65 20 47 4e 55 20 4c 65 73 73 65  of the GNU Lesse
0110: 72 20 47 65 6e 65 72 61 6c 20 50 75 62 6c 69 63  r General Public
0120: 20 4c 69 63 65 6e 73 65 20 61 73 20 70 75 62 6c   License as publ
0130: 69 73 68 65 64 20 62 79 0a 74 68 65 20 46 72 65  ished by.the Fre
0140: 65 20 53 6f 66 74 77 61 72 65 20 46 6f 75 6e 64  e Software Found
0150: 61 74 69 6f 6e 3b 20 65 69 74 68 65 72 20 76 65  ation; either ve
0160: 72 73 69 6f 6e 20 32 2e 31 20 6f 66 20 74 68 65  rsion 2.1 of the
0170: 20 4c 69 63 65 6e 73 65 2c 20 6f 72 20 28 61 74   License, or (at
0180: 20 79 6f 75 72 0a 6f 70 74 69 6f 6e 29 20 61 6e   your.option) an
0190: 79 20 6c 61 74 65 72 20 76 65 72 73 69 6f 6e 2e  y later version.
01a0: 0a 0a 54 68 65 20 47 4e 55 20 4d 50 20 4c 69 62  ..The GNU MP Lib
01b0: 72 61 72 79 20 69 73 20 64 69 73 74 72 69 62 75  rary is distribu
01c0: 74 65 64 20 69 6e 20 74 68 65 20 68 6f 70 65 20  ted in the hope 
01d0: 74 68 61 74 20 69 74 20 77 69 6c 6c 20 62 65 20  that it will be 
01e0: 75 73 65 66 75 6c 2c 20 62 75 74 0a 57 49 54 48  useful, but.WITH
01f0: 4f 55 54 20 41 4e 59 20 57 41 52 52 41 4e 54 59  OUT ANY WARRANTY
0200: 3b 20 77 69 74 68 6f 75 74 20 65 76 65 6e 20 74  ; without even t
0210: 68 65 20 69 6d 70 6c 69 65 64 20 77 61 72 72 61  he implied warra
0220: 6e 74 79 20 6f 66 20 4d 45 52 43 48 41 4e 54 41  nty of MERCHANTA
0230: 42 49 4c 49 54 59 0a 6f 72 20 46 49 54 4e 45 53  BILITY.or FITNES
0240: 53 20 46 4f 52 20 41 20 50 41 52 54 49 43 55 4c  S FOR A PARTICUL
0250: 41 52 20 50 55 52 50 4f 53 45 2e 20 20 53 65 65  AR PURPOSE.  See
0260: 20 74 68 65 20 47 4e 55 20 4c 65 73 73 65 72 20   the GNU Lesser 
0270: 47 65 6e 65 72 61 6c 20 50 75 62 6c 69 63 0a 4c  General Public.L
0280: 69 63 65 6e 73 65 20 66 6f 72 20 6d 6f 72 65 20  icense for more 
0290: 64 65 74 61 69 6c 73 2e 0a 0a 59 6f 75 20 73 68  details...You sh
02a0: 6f 75 6c 64 20 68 61 76 65 20 72 65 63 65 69 76  ould have receiv
02b0: 65 64 20 61 20 63 6f 70 79 20 6f 66 20 74 68 65  ed a copy of the
02c0: 20 47 4e 55 20 4c 65 73 73 65 72 20 47 65 6e 65   GNU Lesser Gene
02d0: 72 61 6c 20 50 75 62 6c 69 63 20 4c 69 63 65 6e  ral Public Licen
02e0: 73 65 0a 61 6c 6f 6e 67 20 77 69 74 68 20 74 68  se.along with th
02f0: 65 20 47 4e 55 20 4d 50 20 4c 69 62 72 61 72 79  e GNU MP Library
0300: 3b 20 73 65 65 20 74 68 65 20 66 69 6c 65 20 43  ; see the file C
0310: 4f 50 59 49 4e 47 2e 4c 49 42 2e 20 20 49 66 20  OPYING.LIB.  If 
0320: 6e 6f 74 2c 20 77 72 69 74 65 20 74 6f 0a 74 68  not, write to.th
0330: 65 20 46 72 65 65 20 53 6f 66 74 77 61 72 65 20  e Free Software 
0340: 46 6f 75 6e 64 61 74 69 6f 6e 2c 20 49 6e 63 2e  Foundation, Inc.
0350: 2c 20 35 39 20 54 65 6d 70 6c 65 20 50 6c 61 63  , 59 Temple Plac
0360: 65 20 2d 20 53 75 69 74 65 20 33 33 30 2c 20 42  e - Suite 330, B
0370: 6f 73 74 6f 6e 2c 0a 4d 41 20 30 32 31 31 31 2d  oston,.MA 02111-
0380: 31 33 30 37 2c 20 55 53 41 2e 20 2a 2f 0a 0a 23  1307, USA. */..#
0390: 69 6e 63 6c 75 64 65 20 22 67 6d 70 2e 68 22 0a  include "gmp.h".
03a0: 23 69 6e 63 6c 75 64 65 20 22 67 6d 70 2d 69 6d  #include "gmp-im
03b0: 70 6c 2e 68 22 0a 0a 76 6f 69 64 0a 6d 70 66 5f  pl.h"..void.mpf_
03c0: 64 69 76 5f 32 65 78 70 20 28 6d 70 66 5f 70 74  div_2exp (mpf_pt
03d0: 72 20 72 2c 20 6d 70 66 5f 73 72 63 70 74 72 20  r r, mpf_srcptr 
03e0: 75 2c 20 75 6e 73 69 67 6e 65 64 20 6c 6f 6e 67  u, unsigned long
03f0: 20 69 6e 74 20 65 78 70 29 0a 7b 0a 20 20 6d 70   int exp).{.  mp
0400: 5f 73 72 63 70 74 72 20 75 70 3b 0a 20 20 6d 70  _srcptr up;.  mp
0410: 5f 70 74 72 20 72 70 20 3d 20 72 2d 3e 5f 6d 70  _ptr rp = r->_mp
0420: 5f 64 3b 0a 20 20 6d 70 5f 73 69 7a 65 5f 74 20  _d;.  mp_size_t 
0430: 75 73 69 7a 65 3b 0a 20 20 6d 70 5f 73 69 7a 65  usize;.  mp_size
0440: 5f 74 20 61 62 73 5f 75 73 69 7a 65 3b 0a 20 20  _t abs_usize;.  
0450: 6d 70 5f 73 69 7a 65 5f 74 20 70 72 65 63 20 3d  mp_size_t prec =
0460: 20 72 2d 3e 5f 6d 70 5f 70 72 65 63 3b 0a 20 20   r->_mp_prec;.  
0470: 6d 70 5f 65 78 70 5f 74 20 75 65 78 70 20 3d 20  mp_exp_t uexp = 
0480: 75 2d 3e 5f 6d 70 5f 65 78 70 3b 0a 0a 20 20 75  u->_mp_exp;..  u
0490: 73 69 7a 65 20 3d 20 75 2d 3e 5f 6d 70 5f 73 69  size = u->_mp_si
04a0: 7a 65 3b 0a 0a 20 20 69 66 20 28 75 73 69 7a 65  ze;..  if (usize
04b0: 20 3d 3d 20 30 29 0a 20 20 20 20 7b 0a 20 20 20   == 0).    {.   
04c0: 20 20 20 72 2d 3e 5f 6d 70 5f 73 69 7a 65 20 3d     r->_mp_size =
04d0: 20 30 3b 0a 20 20 20 20 20 20 72 2d 3e 5f 6d 70   0;.      r->_mp
04e0: 5f 65 78 70 20 3d 20 30 3b 0a 20 20 20 20 20 20  _exp = 0;.      
04f0: 72 65 74 75 72 6e 3b 0a 20 20 20 20 7d 0a 0a 20  return;.    }.. 
0500: 20 61 62 73 5f 75 73 69 7a 65 20 3d 20 41 42 53   abs_usize = ABS
0510: 20 28 75 73 69 7a 65 29 3b 0a 20 20 75 70 20 3d   (usize);.  up =
0520: 20 75 2d 3e 5f 6d 70 5f 64 3b 0a 0a 20 20 69 66   u->_mp_d;..  if
0530: 20 28 65 78 70 20 25 20 47 4d 50 5f 4e 55 4d 42   (exp % GMP_NUMB
0540: 5f 42 49 54 53 20 3d 3d 20 30 29 0a 20 20 20 20  _BITS == 0).    
0550: 7b 0a 20 20 20 20 20 20 70 72 65 63 2b 2b 3b 09  {.      prec++;.
0560: 09 09 2f 2a 20 72 65 74 61 69 6e 20 6d 6f 72 65  ../* retain more
0570: 20 70 72 65 63 69 73 69 6f 6e 20 68 65 72 65 20   precision here 
0580: 61 73 20 77 65 20 64 6f 6e 27 74 20 6e 65 65 64  as we don't need
0590: 0a 09 09 09 09 20 20 20 74 6f 20 61 63 63 6f 75  .....   to accou
05a0: 6e 74 20 66 6f 72 20 63 61 72 72 79 2d 6f 75 74  nt for carry-out
05b0: 20 68 65 72 65 20 2a 2f 0a 20 20 20 20 20 20 69   here */.      i
05c0: 66 20 28 61 62 73 5f 75 73 69 7a 65 20 3e 20 70  f (abs_usize > p
05d0: 72 65 63 29 0a 09 7b 0a 09 20 20 75 70 20 2b 3d  rec)..{..  up +=
05e0: 20 61 62 73 5f 75 73 69 7a 65 20 2d 20 70 72 65   abs_usize - pre
05f0: 63 3b 0a 09 20 20 61 62 73 5f 75 73 69 7a 65 20  c;..  abs_usize 
0600: 3d 20 70 72 65 63 3b 0a 09 7d 0a 20 20 20 20 20  = prec;..}.     
0610: 20 69 66 20 28 72 70 20 21 3d 20 75 70 29 0a 09   if (rp != up)..
0620: 4d 50 4e 5f 43 4f 50 59 5f 49 4e 43 52 20 28 72  MPN_COPY_INCR (r
0630: 70 2c 20 75 70 2c 20 61 62 73 5f 75 73 69 7a 65  p, up, abs_usize
0640: 29 3b 0a 20 20 20 20 20 20 72 2d 3e 5f 6d 70 5f  );.      r->_mp_
0650: 65 78 70 20 3d 20 75 65 78 70 20 2d 20 65 78 70  exp = uexp - exp
0660: 20 2f 20 47 4d 50 5f 4e 55 4d 42 5f 42 49 54 53   / GMP_NUMB_BITS
0670: 3b 0a 20 20 20 20 7d 0a 20 20 65 6c 73 65 0a 20  ;.    }.  else. 
0680: 20 20 20 7b 0a 20 20 20 20 20 20 6d 70 5f 6c 69     {.      mp_li
0690: 6d 62 5f 74 20 63 79 5f 6c 69 6d 62 3b 0a 20 20  mb_t cy_limb;.  
06a0: 20 20 20 20 6d 70 5f 73 69 7a 65 5f 74 20 61 64      mp_size_t ad
06b0: 6a 3b 0a 20 20 20 20 20 20 69 66 20 28 61 62 73  j;.      if (abs
06c0: 5f 75 73 69 7a 65 20 3e 20 70 72 65 63 29 0a 09  _usize > prec)..
06d0: 7b 0a 09 20 20 75 70 20 2b 3d 20 61 62 73 5f 75  {..  up += abs_u
06e0: 73 69 7a 65 20 2d 20 70 72 65 63 3b 0a 09 20 20  size - prec;..  
06f0: 61 62 73 5f 75 73 69 7a 65 20 3d 20 70 72 65 63  abs_usize = prec
0700: 3b 0a 09 20 20 2f 2a 20 55 73 65 20 6d 70 6e 5f  ;..  /* Use mpn_
0710: 72 73 68 69 66 74 20 73 69 6e 63 65 20 6d 70 6e  rshift since mpn
0720: 5f 6c 73 68 69 66 74 20 6f 70 65 72 61 74 65 73  _lshift operates
0730: 20 64 6f 77 6e 77 61 72 64 73 2c 20 61 6e 64 20   downwards, and 
0740: 77 65 0a 09 20 20 20 20 20 74 68 65 72 65 66 6f  we..     therefo
0750: 72 65 20 77 6f 75 6c 64 20 63 6c 6f 62 62 65 72  re would clobber
0760: 20 70 61 72 74 20 6f 66 20 55 20 62 65 66 6f 72   part of U befor
0770: 65 20 75 73 69 6e 67 20 74 68 61 74 20 70 61 72  e using that par
0780: 74 2c 20 69 6e 20 63 61 73 65 0a 09 20 20 20 20  t, in case..    
0790: 20 52 20 69 73 20 74 68 65 20 73 61 6d 65 20 76   R is the same v
07a0: 61 72 69 61 62 6c 65 20 61 73 20 55 2e 20 20 2a  ariable as U.  *
07b0: 2f 0a 09 20 20 63 79 5f 6c 69 6d 62 20 3d 20 6d  /..  cy_limb = m
07c0: 70 6e 5f 72 73 68 69 66 74 20 28 72 70 20 2b 20  pn_rshift (rp + 
07d0: 31 2c 20 75 70 2c 20 61 62 73 5f 75 73 69 7a 65  1, up, abs_usize
07e0: 2c 20 65 78 70 20 25 20 47 4d 50 5f 4e 55 4d 42  , exp % GMP_NUMB
07f0: 5f 42 49 54 53 29 3b 0a 09 20 20 72 70 5b 30 5d  _BITS);..  rp[0]
0800: 20 3d 20 63 79 5f 6c 69 6d 62 3b 0a 09 20 20 61   = cy_limb;..  a
0810: 64 6a 20 3d 20 72 70 5b 61 62 73 5f 75 73 69 7a  dj = rp[abs_usiz
0820: 65 5d 20 21 3d 20 30 3b 0a 09 7d 0a 20 20 20 20  e] != 0;..}.    
0830: 20 20 65 6c 73 65 0a 09 7b 0a 09 20 20 63 79 5f    else..{..  cy_
0840: 6c 69 6d 62 20 3d 20 6d 70 6e 5f 6c 73 68 69 66  limb = mpn_lshif
0850: 74 20 28 72 70 2c 20 75 70 2c 20 61 62 73 5f 75  t (rp, up, abs_u
0860: 73 69 7a 65 2c 0a 09 09 09 09 47 4d 50 5f 4e 55  size,.....GMP_NU
0870: 4d 42 5f 42 49 54 53 20 2d 20 65 78 70 20 25 20  MB_BITS - exp % 
0880: 47 4d 50 5f 4e 55 4d 42 5f 42 49 54 53 29 3b 0a  GMP_NUMB_BITS);.
0890: 09 20 20 72 70 5b 61 62 73 5f 75 73 69 7a 65 5d  .  rp[abs_usize]
08a0: 20 3d 20 63 79 5f 6c 69 6d 62 3b 0a 09 20 20 61   = cy_limb;..  a
08b0: 64 6a 20 3d 20 63 79 5f 6c 69 6d 62 20 21 3d 20  dj = cy_limb != 
08c0: 30 3b 0a 09 7d 0a 0a 20 20 20 20 20 20 61 62 73  0;..}..      abs
08d0: 5f 75 73 69 7a 65 20 2b 3d 20 61 64 6a 3b 0a 20  _usize += adj;. 
08e0: 20 20 20 20 20 72 2d 3e 5f 6d 70 5f 65 78 70 20       r->_mp_exp 
08f0: 3d 20 75 65 78 70 20 2d 20 65 78 70 20 2f 20 47  = uexp - exp / G
0900: 4d 50 5f 4e 55 4d 42 5f 42 49 54 53 20 2d 20 31  MP_NUMB_BITS - 1
0910: 20 2b 20 61 64 6a 3b 0a 20 20 20 20 7d 0a 20 20   + adj;.    }.  
0920: 72 2d 3e 5f 6d 70 5f 73 69 7a 65 20 3d 20 75 73  r->_mp_size = us
0930: 69 7a 65 20 3e 3d 20 30 20 3f 20 61 62 73 5f 75  ize >= 0 ? abs_u
0940: 73 69 7a 65 20 3a 20 2d 61 62 73 5f 75 73 69 7a  size : -abs_usiz
0950: 65 3b 0a 7d 0a                                   e;.}.