version 1.1.1.2, 2000/09/09 14:12:49 |
version 1.1.1.3, 2003/08/25 16:06:32 |
|
|
always fit into the return type, the negative of the true remainder is |
always fit into the return type, the negative of the true remainder is |
returned. |
returned. |
|
|
Copyright (C) 1994, 1995, 1996 Free Software Foundation, Inc. |
Copyright 1994, 1995, 1996, 2001, 2002 Free Software Foundation, Inc. |
|
|
This file is part of the GNU MP Library. |
This file is part of the GNU MP Library. |
|
|
Line 26 MA 02111-1307, USA. */ |
|
Line 26 MA 02111-1307, USA. */ |
|
#include "gmp-impl.h" |
#include "gmp-impl.h" |
|
|
unsigned long int |
unsigned long int |
#if __STDC__ |
|
mpz_cdiv_r_ui (mpz_ptr rem, mpz_srcptr dividend, unsigned long int divisor) |
mpz_cdiv_r_ui (mpz_ptr rem, mpz_srcptr dividend, unsigned long int divisor) |
#else |
|
mpz_cdiv_r_ui (rem, dividend, divisor) |
|
mpz_ptr rem; |
|
mpz_srcptr dividend; |
|
unsigned long int divisor; |
|
#endif |
|
{ |
{ |
mp_size_t dividend_size; |
mp_size_t ns, nn; |
mp_size_t size; |
mp_ptr np; |
mp_limb_t remainder_limb; |
mp_limb_t rl; |
|
|
if (divisor == 0) |
if (divisor == 0) |
DIVIDE_BY_ZERO; |
DIVIDE_BY_ZERO; |
|
|
dividend_size = dividend->_mp_size; |
ns = SIZ(dividend); |
size = ABS (dividend_size); |
if (ns == 0) |
|
{ |
|
SIZ(rem) = 0; |
|
return 0; |
|
} |
|
|
remainder_limb = mpn_mod_1 (dividend->_mp_d, size, (mp_limb_t) divisor); |
nn = ABS(ns); |
|
np = PTR(dividend); |
|
#if GMP_NAIL_BITS != 0 |
|
if (divisor > GMP_NUMB_MAX) |
|
{ |
|
mp_limb_t dp[2]; |
|
mp_ptr rp, qp; |
|
mp_size_t rn; |
|
TMP_DECL (mark); |
|
|
if (remainder_limb != 0 && dividend_size >= 0) |
MPZ_REALLOC (rem, 2); |
remainder_limb = divisor - remainder_limb; |
rp = PTR(rem); |
|
|
rem->_mp_d[0] = remainder_limb; |
if (nn == 1) /* tdiv_qr requirements; tested above for 0 */ |
rem->_mp_size = -(remainder_limb != 0); |
{ |
|
rl = np[0]; |
|
rp[0] = rl; |
|
} |
|
else |
|
{ |
|
TMP_MARK (mark); |
|
dp[0] = divisor & GMP_NUMB_MASK; |
|
dp[1] = divisor >> GMP_NUMB_BITS; |
|
qp = TMP_ALLOC_LIMBS (nn - 2 + 1); |
|
mpn_tdiv_qr (qp, rp, (mp_size_t) 0, np, nn, dp, (mp_size_t) 2); |
|
TMP_FREE (mark); |
|
rl = rp[0] + (rp[1] << GMP_NUMB_BITS); |
|
} |
|
|
return remainder_limb; |
if (rl != 0 && ns >= 0) |
|
{ |
|
rl = divisor - rl; |
|
rp[0] = rl & GMP_NUMB_MASK; |
|
rp[1] = rl >> GMP_NUMB_BITS; |
|
} |
|
|
|
rn = 1 + (rl > GMP_NUMB_MAX); rn -= (rp[rn - 1] == 0); |
|
SIZ(rem) = -rn; |
|
} |
|
else |
|
#endif |
|
{ |
|
rl = mpn_mod_1 (np, nn, (mp_limb_t) divisor); |
|
if (rl == 0) |
|
SIZ(rem) = 0; |
|
else |
|
{ |
|
if (ns >= 0) |
|
rl = divisor - rl; |
|
|
|
PTR(rem)[0] = rl; |
|
SIZ(rem) = -1; |
|
} |
|
} |
|
|
|
return rl; |
} |
} |