Hex Artifact Content
Not logged in

Artifact 79f0cd54a0829fb9760e8fae1aaa48d04e3b6e04:


0000: 2f 2a 20 6d 70 66 5f 61 64 64 5f 75 69 20 2d 2d  /* mpf_add_ui --
0010: 20 41 64 64 20 61 20 66 6c 6f 61 74 20 61 6e 64   Add a float and
0020: 20 61 6e 20 75 6e 73 69 67 6e 65 64 20 69 6e 74   an unsigned int
0030: 65 67 65 72 2e 0a 0a 43 6f 70 79 72 69 67 68 74  eger...Copyright
0040: 20 31 39 39 33 2c 20 31 39 39 34 2c 20 31 39 39   1993, 1994, 199
0050: 36 2c 20 32 30 30 30 2c 20 32 30 30 31 20 46 72  6, 2000, 2001 Fr
0060: 65 65 20 53 6f 66 74 77 61 72 65 20 46 6f 75 6e  ee Software Foun
0070: 64 61 74 69 6f 6e 2c 20 49 6e 63 2e 0a 0a 54 68  dation, Inc...Th
0080: 69 73 20 66 69 6c 65 20 69 73 20 70 61 72 74 20  is file is part 
0090: 6f 66 20 74 68 65 20 47 4e 55 20 4d 50 20 4c 69  of the GNU MP Li
00a0: 62 72 61 72 79 2e 0a 0a 54 68 65 20 47 4e 55 20  brary...The GNU 
00b0: 4d 50 20 4c 69 62 72 61 72 79 20 69 73 20 66 72  MP Library is fr
00c0: 65 65 20 73 6f 66 74 77 61 72 65 3b 20 79 6f 75  ee software; you
00d0: 20 63 61 6e 20 72 65 64 69 73 74 72 69 62 75 74   can redistribut
00e0: 65 20 69 74 20 61 6e 64 2f 6f 72 20 6d 6f 64 69  e it and/or modi
00f0: 66 79 0a 69 74 20 75 6e 64 65 72 20 74 68 65 20  fy.it under the 
0100: 74 65 72 6d 73 20 6f 66 20 74 68 65 20 47 4e 55  terms of the GNU
0110: 20 4c 65 73 73 65 72 20 47 65 6e 65 72 61 6c 20   Lesser General 
0120: 50 75 62 6c 69 63 20 4c 69 63 65 6e 73 65 20 61  Public License a
0130: 73 20 70 75 62 6c 69 73 68 65 64 20 62 79 0a 74  s published by.t
0140: 68 65 20 46 72 65 65 20 53 6f 66 74 77 61 72 65  he Free Software
0150: 20 46 6f 75 6e 64 61 74 69 6f 6e 3b 20 65 69 74   Foundation; eit
0160: 68 65 72 20 76 65 72 73 69 6f 6e 20 32 2e 31 20  her version 2.1 
0170: 6f 66 20 74 68 65 20 4c 69 63 65 6e 73 65 2c 20  of the License, 
0180: 6f 72 20 28 61 74 20 79 6f 75 72 0a 6f 70 74 69  or (at your.opti
0190: 6f 6e 29 20 61 6e 79 20 6c 61 74 65 72 20 76 65  on) any later ve
01a0: 72 73 69 6f 6e 2e 0a 0a 54 68 65 20 47 4e 55 20  rsion...The GNU 
01b0: 4d 50 20 4c 69 62 72 61 72 79 20 69 73 20 64 69  MP Library is di
01c0: 73 74 72 69 62 75 74 65 64 20 69 6e 20 74 68 65  stributed in the
01d0: 20 68 6f 70 65 20 74 68 61 74 20 69 74 20 77 69   hope that it wi
01e0: 6c 6c 20 62 65 20 75 73 65 66 75 6c 2c 20 62 75  ll be useful, bu
01f0: 74 0a 57 49 54 48 4f 55 54 20 41 4e 59 20 57 41  t.WITHOUT ANY WA
0200: 52 52 41 4e 54 59 3b 20 77 69 74 68 6f 75 74 20  RRANTY; without 
0210: 65 76 65 6e 20 74 68 65 20 69 6d 70 6c 69 65 64  even the implied
0220: 20 77 61 72 72 61 6e 74 79 20 6f 66 20 4d 45 52   warranty of MER
0230: 43 48 41 4e 54 41 42 49 4c 49 54 59 0a 6f 72 20  CHANTABILITY.or 
0240: 46 49 54 4e 45 53 53 20 46 4f 52 20 41 20 50 41  FITNESS FOR A PA
0250: 52 54 49 43 55 4c 41 52 20 50 55 52 50 4f 53 45  RTICULAR PURPOSE
0260: 2e 20 20 53 65 65 20 74 68 65 20 47 4e 55 20 4c  .  See the GNU L
0270: 65 73 73 65 72 20 47 65 6e 65 72 61 6c 20 50 75  esser General Pu
0280: 62 6c 69 63 0a 4c 69 63 65 6e 73 65 20 66 6f 72  blic.License for
0290: 20 6d 6f 72 65 20 64 65 74 61 69 6c 73 2e 0a 0a   more details...
02a0: 59 6f 75 20 73 68 6f 75 6c 64 20 68 61 76 65 20  You should have 
02b0: 72 65 63 65 69 76 65 64 20 61 20 63 6f 70 79 20  received a copy 
02c0: 6f 66 20 74 68 65 20 47 4e 55 20 4c 65 73 73 65  of the GNU Lesse
02d0: 72 20 47 65 6e 65 72 61 6c 20 50 75 62 6c 69 63  r General Public
02e0: 20 4c 69 63 65 6e 73 65 0a 61 6c 6f 6e 67 20 77   License.along w
02f0: 69 74 68 20 74 68 65 20 47 4e 55 20 4d 50 20 4c  ith the GNU MP L
0300: 69 62 72 61 72 79 3b 20 73 65 65 20 74 68 65 20  ibrary; see the 
0310: 66 69 6c 65 20 43 4f 50 59 49 4e 47 2e 4c 49 42  file COPYING.LIB
0320: 2e 20 20 49 66 20 6e 6f 74 2c 20 77 72 69 74 65  .  If not, write
0330: 20 74 6f 0a 74 68 65 20 46 72 65 65 20 53 6f 66   to.the Free Sof
0340: 74 77 61 72 65 20 46 6f 75 6e 64 61 74 69 6f 6e  tware Foundation
0350: 2c 20 49 6e 63 2e 2c 20 35 39 20 54 65 6d 70 6c  , Inc., 59 Templ
0360: 65 20 50 6c 61 63 65 20 2d 20 53 75 69 74 65 20  e Place - Suite 
0370: 33 33 30 2c 20 42 6f 73 74 6f 6e 2c 0a 4d 41 20  330, Boston,.MA 
0380: 30 32 31 31 31 2d 31 33 30 37 2c 20 55 53 41 2e  02111-1307, USA.
0390: 20 2a 2f 0a 0a 23 69 6e 63 6c 75 64 65 20 22 67   */..#include "g
03a0: 6d 70 2e 68 22 0a 23 69 6e 63 6c 75 64 65 20 22  mp.h".#include "
03b0: 67 6d 70 2d 69 6d 70 6c 2e 68 22 0a 0a 76 6f 69  gmp-impl.h"..voi
03c0: 64 0a 6d 70 66 5f 61 64 64 5f 75 69 20 28 6d 70  d.mpf_add_ui (mp
03d0: 66 5f 70 74 72 20 73 75 6d 2c 20 6d 70 66 5f 73  f_ptr sum, mpf_s
03e0: 72 63 70 74 72 20 75 2c 20 75 6e 73 69 67 6e 65  rcptr u, unsigne
03f0: 64 20 6c 6f 6e 67 20 69 6e 74 20 76 29 0a 7b 0a  d long int v).{.
0400: 20 20 6d 70 5f 73 72 63 70 74 72 20 75 70 20 3d    mp_srcptr up =
0410: 20 75 2d 3e 5f 6d 70 5f 64 3b 0a 20 20 6d 70 5f   u->_mp_d;.  mp_
0420: 70 74 72 20 73 75 6d 70 20 3d 20 73 75 6d 2d 3e  ptr sump = sum->
0430: 5f 6d 70 5f 64 3b 0a 20 20 6d 70 5f 73 69 7a 65  _mp_d;.  mp_size
0440: 5f 74 20 75 73 69 7a 65 2c 20 73 75 6d 73 69 7a  _t usize, sumsiz
0450: 65 3b 0a 20 20 6d 70 5f 73 69 7a 65 5f 74 20 70  e;.  mp_size_t p
0460: 72 65 63 20 3d 20 73 75 6d 2d 3e 5f 6d 70 5f 70  rec = sum->_mp_p
0470: 72 65 63 3b 0a 20 20 6d 70 5f 65 78 70 5f 74 20  rec;.  mp_exp_t 
0480: 75 65 78 70 20 3d 20 75 2d 3e 5f 6d 70 5f 65 78  uexp = u->_mp_ex
0490: 70 3b 0a 0a 20 20 75 73 69 7a 65 20 3d 20 75 2d  p;..  usize = u-
04a0: 3e 5f 6d 70 5f 73 69 7a 65 3b 0a 20 20 69 66 20  >_mp_size;.  if 
04b0: 28 75 73 69 7a 65 20 3c 3d 20 30 29 0a 20 20 20  (usize <= 0).   
04c0: 20 7b 0a 20 20 20 20 20 20 69 66 20 28 75 73 69   {.      if (usi
04d0: 7a 65 20 3d 3d 20 30 29 0a 09 7b 0a 09 20 20 6d  ze == 0)..{..  m
04e0: 70 66 5f 73 65 74 5f 75 69 20 28 73 75 6d 2c 20  pf_set_ui (sum, 
04f0: 76 29 3b 0a 09 20 20 72 65 74 75 72 6e 3b 0a 09  v);..  return;..
0500: 7d 0a 20 20 20 20 20 20 65 6c 73 65 0a 09 7b 0a  }.      else..{.
0510: 09 20 20 5f 5f 6d 70 66 5f 73 74 72 75 63 74 20  .  __mpf_struct 
0520: 75 5f 6e 65 67 61 74 65 64 3b 0a 09 20 20 75 5f  u_negated;..  u_
0530: 6e 65 67 61 74 65 64 2e 5f 6d 70 5f 73 69 7a 65  negated._mp_size
0540: 20 3d 20 2d 75 73 69 7a 65 3b 0a 09 20 20 75 5f   = -usize;..  u_
0550: 6e 65 67 61 74 65 64 2e 5f 6d 70 5f 65 78 70 20  negated._mp_exp 
0560: 3d 20 75 2d 3e 5f 6d 70 5f 65 78 70 3b 0a 09 20  = u->_mp_exp;.. 
0570: 20 75 5f 6e 65 67 61 74 65 64 2e 5f 6d 70 5f 64   u_negated._mp_d
0580: 20 3d 20 75 2d 3e 5f 6d 70 5f 64 3b 0a 09 20 20   = u->_mp_d;..  
0590: 6d 70 66 5f 73 75 62 5f 75 69 20 28 73 75 6d 2c  mpf_sub_ui (sum,
05a0: 20 26 75 5f 6e 65 67 61 74 65 64 2c 20 76 29 3b   &u_negated, v);
05b0: 0a 09 20 20 73 75 6d 2d 3e 5f 6d 70 5f 73 69 7a  ..  sum->_mp_siz
05c0: 65 20 3d 20 2d 28 73 75 6d 2d 3e 5f 6d 70 5f 73  e = -(sum->_mp_s
05d0: 69 7a 65 29 3b 0a 09 20 20 72 65 74 75 72 6e 3b  ize);..  return;
05e0: 0a 09 7d 0a 20 20 20 20 7d 0a 0a 20 20 69 66 20  ..}.    }..  if 
05f0: 28 76 20 3d 3d 20 30 29 0a 20 20 20 20 7b 0a 20  (v == 0).    {. 
0600: 20 20 20 73 75 6d 5f 69 73 5f 75 3a 0a 20 20 20     sum_is_u:.   
0610: 20 20 20 69 66 20 28 75 20 21 3d 20 73 75 6d 29     if (u != sum)
0620: 0a 09 7b 0a 09 20 20 73 75 6d 73 69 7a 65 20 3d  ..{..  sumsize =
0630: 20 4d 49 4e 20 28 75 73 69 7a 65 2c 20 70 72 65   MIN (usize, pre
0640: 63 20 2b 20 31 29 3b 0a 09 20 20 4d 50 4e 5f 43  c + 1);..  MPN_C
0650: 4f 50 59 20 28 73 75 6d 2d 3e 5f 6d 70 5f 64 2c  OPY (sum->_mp_d,
0660: 20 75 70 20 2b 20 75 73 69 7a 65 20 2d 20 73 75   up + usize - su
0670: 6d 73 69 7a 65 2c 20 73 75 6d 73 69 7a 65 29 3b  msize, sumsize);
0680: 0a 09 20 20 73 75 6d 2d 3e 5f 6d 70 5f 73 69 7a  ..  sum->_mp_siz
0690: 65 20 3d 20 73 75 6d 73 69 7a 65 3b 0a 09 20 20  e = sumsize;..  
06a0: 73 75 6d 2d 3e 5f 6d 70 5f 65 78 70 20 3d 20 75  sum->_mp_exp = u
06b0: 2d 3e 5f 6d 70 5f 65 78 70 3b 0a 09 7d 0a 20 20  ->_mp_exp;..}.  
06c0: 20 20 20 20 72 65 74 75 72 6e 3b 0a 20 20 20 20      return;.    
06d0: 7d 0a 0a 20 20 69 66 20 28 75 65 78 70 20 3e 20  }..  if (uexp > 
06e0: 30 29 0a 20 20 20 20 7b 0a 20 20 20 20 20 20 2f  0).    {.      /
06f0: 2a 20 55 20 3e 3d 20 31 2e 20 20 2a 2f 0a 20 20  * U >= 1.  */.  
0700: 20 20 20 20 69 66 20 28 75 65 78 70 20 3e 20 70      if (uexp > p
0710: 72 65 63 29 0a 09 7b 0a 09 20 20 2f 2a 20 55 20  rec)..{..  /* U 
0720: 3e 3e 20 56 2c 20 56 20 69 73 20 6e 6f 74 20 70  >> V, V is not p
0730: 61 72 74 20 6f 66 20 66 69 6e 61 6c 20 72 65 73  art of final res
0740: 75 6c 74 2e 20 20 2a 2f 0a 09 20 20 67 6f 74 6f  ult.  */..  goto
0750: 20 73 75 6d 5f 69 73 5f 75 3b 0a 09 7d 0a 20 20   sum_is_u;..}.  
0760: 20 20 20 20 65 6c 73 65 0a 09 7b 0a 09 20 20 2f      else..{..  /
0770: 2a 20 55 27 73 20 22 6c 69 6d 62 20 70 6f 69 6e  * U's "limb poin
0780: 74 22 20 69 73 20 73 6f 6d 65 77 68 65 72 65 20  t" is somewhere 
0790: 62 65 74 77 65 65 6e 20 74 68 65 20 66 69 72 73  between the firs
07a0: 74 20 6c 69 6d 62 0a 09 20 20 20 20 20 61 6e 64  t limb..     and
07b0: 20 74 68 65 20 50 52 45 43 3a 74 68 20 6c 69 6d   the PREC:th lim
07c0: 62 2e 0a 09 20 20 20 20 20 42 6f 74 68 20 55 20  b...     Both U 
07d0: 61 6e 64 20 56 20 61 72 65 20 70 61 72 74 20 6f  and V are part o
07e0: 66 20 74 68 65 20 66 69 6e 61 6c 20 72 65 73 75  f the final resu
07f0: 6c 74 2e 20 20 2a 2f 0a 09 20 20 69 66 20 28 75  lt.  */..  if (u
0800: 65 78 70 20 3e 20 75 73 69 7a 65 29 0a 09 20 20  exp > usize)..  
0810: 20 20 7b 0a 09 20 20 20 20 20 20 2f 2a 20 20 20    {..      /*   
0820: 75 75 75 75 75 75 30 30 30 30 2e 20 2a 2f 0a 09  uuuuuu0000. */..
0830: 20 20 20 20 20 20 2f 2a 20 2b 20 20 20 20 20 20        /* +      
0840: 20 20 20 20 76 2e 20 2a 2f 0a 09 20 20 20 20 20      v. */..     
0850: 20 2f 2a 20 57 65 20 62 65 67 69 6e 20 77 69 74   /* We begin wit
0860: 68 20 6d 6f 76 69 6e 67 20 55 20 74 6f 20 74 68  h moving U to th
0870: 65 20 74 6f 70 20 6f 66 20 53 55 4d 2c 20 74 6f  e top of SUM, to
0880: 20 68 61 6e 64 6c 65 0a 09 09 20 73 61 6d 65 76   handle... samev
0890: 61 72 28 55 2c 53 55 4d 29 2e 20 20 2a 2f 0a 09  ar(U,SUM).  */..
08a0: 20 20 20 20 20 20 4d 50 4e 5f 43 4f 50 59 5f 44        MPN_COPY_D
08b0: 45 43 52 20 28 73 75 6d 70 20 2b 20 75 65 78 70  ECR (sump + uexp
08c0: 20 2d 20 75 73 69 7a 65 2c 20 75 70 2c 20 75 73   - usize, up, us
08d0: 69 7a 65 29 3b 0a 09 20 20 20 20 20 20 73 75 6d  ize);..      sum
08e0: 70 5b 30 5d 20 3d 20 76 3b 0a 09 20 20 20 20 20  p[0] = v;..     
08f0: 20 4d 50 4e 5f 5a 45 52 4f 20 28 73 75 6d 70 20   MPN_ZERO (sump 
0900: 2b 20 31 2c 20 75 65 78 70 20 2d 20 75 73 69 7a  + 1, uexp - usiz
0910: 65 20 2d 20 31 29 3b 0a 23 69 66 20 30 20 2f 2a  e - 1);.#if 0 /*
0920: 20 57 68 61 74 20 69 73 20 74 68 69 73 3f 3f 3f   What is this???
0930: 20 2a 2f 0a 09 20 20 20 20 20 20 69 66 20 28 73   */..      if (s
0940: 75 6d 20 3d 3d 20 75 29 0a 09 09 4d 50 4e 5f 43  um == u)...MPN_C
0950: 4f 50 59 20 28 73 75 6d 2d 3e 5f 6d 70 5f 64 2c  OPY (sum->_mp_d,
0960: 20 73 75 6d 70 2c 20 75 65 78 70 29 3b 0a 23 65   sump, uexp);.#e
0970: 6e 64 69 66 0a 09 20 20 20 20 20 20 73 75 6d 2d  ndif..      sum-
0980: 3e 5f 6d 70 5f 73 69 7a 65 20 3d 20 75 65 78 70  >_mp_size = uexp
0990: 3b 0a 09 20 20 20 20 20 20 73 75 6d 2d 3e 5f 6d  ;..      sum->_m
09a0: 70 5f 65 78 70 20 3d 20 75 65 78 70 3b 0a 09 20  p_exp = uexp;.. 
09b0: 20 20 20 7d 0a 09 20 20 65 6c 73 65 0a 09 20 20     }..  else..  
09c0: 20 20 7b 0a 09 20 20 20 20 20 20 2f 2a 20 20 20    {..      /*   
09d0: 75 75 75 75 75 75 2e 75 75 75 75 20 2a 2f 0a 09  uuuuuu.uuuu */..
09e0: 20 20 20 20 20 20 2f 2a 20 2b 20 20 20 20 20 20        /* +      
09f0: 76 2e 20 20 20 20 20 2a 2f 0a 09 20 20 20 20 20  v.     */..     
0a00: 20 6d 70 5f 6c 69 6d 62 5f 74 20 63 79 5f 6c 69   mp_limb_t cy_li
0a10: 6d 62 3b 0a 09 20 20 20 20 20 20 69 66 20 28 75  mb;..      if (u
0a20: 73 69 7a 65 20 3e 20 70 72 65 63 29 0a 09 09 7b  size > prec)...{
0a30: 0a 09 09 20 20 2f 2a 20 49 67 6e 6f 72 65 20 65  ...  /* Ignore e
0a40: 78 63 65 73 73 20 6c 69 6d 62 73 20 69 6e 20 55  xcess limbs in U
0a50: 2e 20 20 2a 2f 0a 09 09 20 20 75 70 20 2b 3d 20  .  */...  up += 
0a60: 75 73 69 7a 65 20 2d 20 70 72 65 63 3b 0a 09 09  usize - prec;...
0a70: 20 20 75 73 69 7a 65 20 2d 3d 20 75 73 69 7a 65    usize -= usize
0a80: 20 2d 20 70 72 65 63 3b 20 2f 2a 20 45 71 2e 20   - prec; /* Eq. 
0a90: 75 73 69 7a 65 20 3d 20 70 72 65 63 20 2a 2f 0a  usize = prec */.
0aa0: 09 09 7d 0a 09 20 20 20 20 20 20 69 66 20 28 73  ..}..      if (s
0ab0: 75 6d 70 20 21 3d 20 75 70 29 0a 09 09 4d 50 4e  ump != up)...MPN
0ac0: 5f 43 4f 50 59 5f 49 4e 43 52 20 28 73 75 6d 70  _COPY_INCR (sump
0ad0: 2c 20 75 70 2c 20 75 73 69 7a 65 20 2d 20 75 65  , up, usize - ue
0ae0: 78 70 29 3b 0a 09 20 20 20 20 20 20 63 79 5f 6c  xp);..      cy_l
0af0: 69 6d 62 20 3d 20 6d 70 6e 5f 61 64 64 5f 31 20  imb = mpn_add_1 
0b00: 28 73 75 6d 70 20 2b 20 75 73 69 7a 65 20 2d 20  (sump + usize - 
0b10: 75 65 78 70 2c 20 75 70 20 2b 20 75 73 69 7a 65  uexp, up + usize
0b20: 20 2d 20 75 65 78 70 2c 0a 09 09 09 09 20 20 20   - uexp,.....   
0b30: 75 65 78 70 2c 20 28 6d 70 5f 6c 69 6d 62 5f 74  uexp, (mp_limb_t
0b40: 29 20 76 29 3b 0a 09 20 20 20 20 20 20 73 75 6d  ) v);..      sum
0b50: 70 5b 75 73 69 7a 65 5d 20 3d 20 63 79 5f 6c 69  p[usize] = cy_li
0b60: 6d 62 3b 0a 09 20 20 20 20 20 20 73 75 6d 2d 3e  mb;..      sum->
0b70: 5f 6d 70 5f 73 69 7a 65 20 3d 20 75 73 69 7a 65  _mp_size = usize
0b80: 20 2b 20 63 79 5f 6c 69 6d 62 3b 0a 09 20 20 20   + cy_limb;..   
0b90: 20 20 20 73 75 6d 2d 3e 5f 6d 70 5f 65 78 70 20     sum->_mp_exp 
0ba0: 3d 20 75 65 78 70 20 2b 20 63 79 5f 6c 69 6d 62  = uexp + cy_limb
0bb0: 3b 0a 09 20 20 20 20 7d 0a 09 7d 0a 20 20 20 20  ;..    }..}.    
0bc0: 7d 0a 20 20 65 6c 73 65 0a 20 20 20 20 7b 0a 20  }.  else.    {. 
0bd0: 20 20 20 20 20 2f 2a 20 55 20 3c 20 31 2c 20 73       /* U < 1, s
0be0: 6f 20 56 20 3e 20 55 20 66 6f 72 20 73 75 72 65  o V > U for sure
0bf0: 2e 20 20 2a 2f 0a 20 20 20 20 20 20 2f 2a 20 76  .  */.      /* v
0c00: 2e 20 20 20 20 20 20 20 20 20 2a 2f 0a 20 20 20  .         */.   
0c10: 20 20 20 2f 2a 20 20 2e 30 30 30 30 75 75 75 75     /*  .0000uuuu
0c20: 20 2a 2f 0a 20 20 20 20 20 20 69 66 20 28 28 2d   */.      if ((-
0c30: 75 65 78 70 29 20 3e 3d 20 70 72 65 63 29 0a 09  uexp) >= prec)..
0c40: 7b 0a 09 20 20 73 75 6d 70 5b 30 5d 20 3d 20 76  {..  sump[0] = v
0c50: 3b 0a 09 20 20 73 75 6d 2d 3e 5f 6d 70 5f 73 69  ;..  sum->_mp_si
0c60: 7a 65 20 3d 20 31 3b 0a 09 20 20 73 75 6d 2d 3e  ze = 1;..  sum->
0c70: 5f 6d 70 5f 65 78 70 20 3d 20 31 3b 0a 09 7d 0a  _mp_exp = 1;..}.
0c80: 20 20 20 20 20 20 65 6c 73 65 0a 09 7b 0a 09 20        else..{.. 
0c90: 20 69 66 20 28 75 73 69 7a 65 20 2b 20 28 2d 75   if (usize + (-u
0ca0: 65 78 70 29 20 2b 20 31 20 3e 20 70 72 65 63 29  exp) + 1 > prec)
0cb0: 0a 09 20 20 20 20 7b 0a 09 20 20 20 20 20 20 2f  ..    {..      /
0cc0: 2a 20 49 67 6e 6f 72 65 20 65 78 63 65 73 73 20  * Ignore excess 
0cd0: 6c 69 6d 62 73 20 69 6e 20 55 2e 20 20 2a 2f 0a  limbs in U.  */.
0ce0: 09 20 20 20 20 20 20 75 70 20 2b 3d 20 75 73 69  .      up += usi
0cf0: 7a 65 20 2b 20 28 2d 75 65 78 70 29 20 2b 20 31  ze + (-uexp) + 1
0d00: 20 2d 20 70 72 65 63 3b 0a 09 20 20 20 20 20 20   - prec;..      
0d10: 75 73 69 7a 65 20 2d 3d 20 75 73 69 7a 65 20 2b  usize -= usize +
0d20: 20 28 2d 75 65 78 70 29 20 2b 20 31 20 2d 20 70   (-uexp) + 1 - p
0d30: 72 65 63 3b 0a 09 20 20 20 20 7d 0a 09 20 20 69  rec;..    }..  i
0d40: 66 20 28 73 75 6d 70 20 21 3d 20 75 70 29 0a 09  f (sump != up)..
0d50: 20 20 20 20 4d 50 4e 5f 43 4f 50 59 5f 49 4e 43      MPN_COPY_INC
0d60: 52 20 28 73 75 6d 70 2c 20 75 70 2c 20 75 73 69  R (sump, up, usi
0d70: 7a 65 29 3b 0a 09 20 20 4d 50 4e 5f 5a 45 52 4f  ze);..  MPN_ZERO
0d80: 20 28 73 75 6d 70 20 2b 20 75 73 69 7a 65 2c 20   (sump + usize, 
0d90: 2d 75 65 78 70 29 3b 0a 09 20 20 73 75 6d 70 5b  -uexp);..  sump[
0da0: 75 73 69 7a 65 20 2b 20 28 2d 75 65 78 70 29 5d  usize + (-uexp)]
0db0: 20 3d 20 76 3b 0a 09 20 20 73 75 6d 2d 3e 5f 6d   = v;..  sum->_m
0dc0: 70 5f 73 69 7a 65 20 3d 20 75 73 69 7a 65 20 2b  p_size = usize +
0dd0: 20 28 2d 75 65 78 70 29 20 2b 20 31 3b 0a 09 20   (-uexp) + 1;.. 
0de0: 20 73 75 6d 2d 3e 5f 6d 70 5f 65 78 70 20 3d 20   sum->_mp_exp = 
0df0: 31 3b 0a 09 7d 0a 20 20 20 20 7d 0a 7d 0a        1;..}.    }.}.