OpenXM_contrib/gmp/mpz/tdiv_r_2exp.c - diff

Return to tdiv_r_2exp.c CVS log

Up to [local] / OpenXM_contrib / gmp / mpz

version 1.1.1.2, 2000/09/09 14:12:59

version 1.1.1.3, 2003/08/25 16:06:34

Line 1

/* mpz_tdiv_r_2exp -- Divide a integer by 2**CNT and produce a remainder.

This file is part of the GNU MP Library.

Line 23 MA 02111-1307, USA. */

#include "gmp-impl.h"

void

#if __STDC__

mpz_tdiv_r_2exp (mpz_ptr res, mpz_srcptr in, unsigned long int cnt)

#else

mpz_tdiv_r_2exp (res, in, cnt)

mpz_ptr res;

mpz_srcptr in;

unsigned long int cnt;

#endif

{

mp_size_t in_size = ABS (in->_mp_size);

mp_size_t res_size;

mp_size_t limb_cnt = cnt / BITS_PER_MP_LIMB;

mp_size_t limb_cnt = cnt / GMP_NUMB_BITS;

mp_srcptr in_ptr = in->_mp_d;

if (in_size > limb_cnt)

Line 42 mpz_tdiv_r_2exp (res, in, cnt)

Line 35 mpz_tdiv_r_2exp (res, in, cnt)

/* The input operand is (probably) greater than 2**CNT. */

mp_limb_t x;

x = in_ptr[limb_cnt] & (((mp_limb_t) 1 << cnt % BITS_PER_MP_LIMB) - 1);

x = in_ptr[limb_cnt] & (((mp_limb_t) 1 << cnt % GMP_NUMB_BITS) - 1);

if (x != 0)

{

res_size = limb_cnt + 1;