1 1.1 mrg /* mpn_add_n -- Add equal length limb vectors. 2 1.1 mrg 3 1.1.1.2 mrg Copyright 1992-1994, 1996, 2000, 2002, 2009 Free Software Foundation, Inc. 4 1.1 mrg 5 1.1 mrg This file is part of the GNU MP Library. 6 1.1 mrg 7 1.1 mrg The GNU MP Library is free software; you can redistribute it and/or modify 8 1.1.1.2 mrg it under the terms of either: 9 1.1.1.2 mrg 10 1.1.1.2 mrg * the GNU Lesser General Public License as published by the Free 11 1.1.1.2 mrg Software Foundation; either version 3 of the License, or (at your 12 1.1.1.2 mrg option) any later version. 13 1.1.1.2 mrg 14 1.1.1.2 mrg or 15 1.1.1.2 mrg 16 1.1.1.2 mrg * the GNU General Public License as published by the Free Software 17 1.1.1.2 mrg Foundation; either version 2 of the License, or (at your option) any 18 1.1.1.2 mrg later version. 19 1.1.1.2 mrg 20 1.1.1.2 mrg or both in parallel, as here. 21 1.1 mrg 22 1.1 mrg The GNU MP Library is distributed in the hope that it will be useful, but 23 1.1 mrg WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 24 1.1.1.2 mrg or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License 25 1.1.1.2 mrg for more details. 26 1.1 mrg 27 1.1.1.2 mrg You should have received copies of the GNU General Public License and the 28 1.1.1.2 mrg GNU Lesser General Public License along with the GNU MP Library. If not, 29 1.1.1.2 mrg see https://www.gnu.org/licenses/. */ 30 1.1 mrg 31 1.1 mrg #include "gmp-impl.h" 32 1.1 mrg 33 1.1 mrg 34 1.1 mrg #if GMP_NAIL_BITS == 0 35 1.1 mrg 36 1.1 mrg mp_limb_t 37 1.1 mrg mpn_add_n (mp_ptr rp, mp_srcptr up, mp_srcptr vp, mp_size_t n) 38 1.1 mrg { 39 1.1 mrg mp_limb_t ul, vl, sl, rl, cy, cy1, cy2; 40 1.1 mrg 41 1.1 mrg ASSERT (n >= 1); 42 1.1 mrg ASSERT (MPN_SAME_OR_INCR_P (rp, up, n)); 43 1.1 mrg ASSERT (MPN_SAME_OR_INCR_P (rp, vp, n)); 44 1.1 mrg 45 1.1 mrg cy = 0; 46 1.1 mrg do 47 1.1 mrg { 48 1.1 mrg ul = *up++; 49 1.1 mrg vl = *vp++; 50 1.1 mrg sl = ul + vl; 51 1.1 mrg cy1 = sl < ul; 52 1.1 mrg rl = sl + cy; 53 1.1 mrg cy2 = rl < sl; 54 1.1 mrg cy = cy1 | cy2; 55 1.1 mrg *rp++ = rl; 56 1.1 mrg } 57 1.1 mrg while (--n != 0); 58 1.1 mrg 59 1.1 mrg return cy; 60 1.1 mrg } 61 1.1 mrg 62 1.1 mrg #endif 63 1.1 mrg 64 1.1 mrg #if GMP_NAIL_BITS >= 1 65 1.1 mrg 66 1.1 mrg mp_limb_t 67 1.1 mrg mpn_add_n (mp_ptr rp, mp_srcptr up, mp_srcptr vp, mp_size_t n) 68 1.1 mrg { 69 1.1 mrg mp_limb_t ul, vl, rl, cy; 70 1.1 mrg 71 1.1 mrg ASSERT (n >= 1); 72 1.1 mrg ASSERT (MPN_SAME_OR_INCR_P (rp, up, n)); 73 1.1 mrg ASSERT (MPN_SAME_OR_INCR_P (rp, vp, n)); 74 1.1 mrg 75 1.1 mrg cy = 0; 76 1.1 mrg do 77 1.1 mrg { 78 1.1 mrg ul = *up++; 79 1.1 mrg vl = *vp++; 80 1.1 mrg rl = ul + vl + cy; 81 1.1 mrg cy = rl >> GMP_NUMB_BITS; 82 1.1 mrg *rp++ = rl & GMP_NUMB_MASK; 83 1.1 mrg } 84 1.1 mrg while (--n != 0); 85 1.1 mrg 86 1.1 mrg return cy; 87 1.1 mrg } 88 1.1 mrg 89 1.1 mrg #endif 90
Indexes created Sun Mar 01 05:31:48 UTC 2026