cheng/mpir
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Hopefully fixed a bug in my mpn_addmod_2expp1_1 macro.

Browse Source
This commit is contained in:
wbhart 2012-10-12 23:40:25 +00:00
parent 3b68dc2fb8
commit 0695935782
1 changed files with 5 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
gmp-impl.h
View File

@ -1091,14 +1091,16 @@ __GMP_DECLSPEC mp_limb_t mpn_divrem_hensel_r_1 __GMP_PROTO ((mp_srcptr, mp_size_
#define mpn_addmod_2expp1_1(r, limbs, c) \
do { \
mp_limb_t __sum = (r)[0] + (c); \
mp_limb_t __sum = (r)[0] + (mp_limb_signed_t)(c); \
/* check if adding c causes carry propagation */ \
if ((mp_limb_signed_t)(__sum ^ (r)[0]) >= 0) \
(r)[0] = __sum; \
else \
{ \
if ((c) >= 0) mpn_add_1((r), (r), (limbs) + 1, (c)); \
else mpn_sub_1((r), (r), (limbs) + 1, -(c)); \
if ((mp_limb_signed_t) (c) >= 0) mpn_add_1((r), \
(r), (limbs) + 1, (mp_limb_signed_t) (c)); \
else mpn_sub_1((r), (r), (limbs) + 1, \
-(mp_limb_signed_t) (c)); \
} \
} while (0)