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;..}. }.}.