OpenXM_contrib/gmp/mpz/gcd.c - diff

Return to gcd.c CVS log

Up to [local] / OpenXM_contrib / gmp / mpz

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

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

Line 1

/* mpz/gcd.c: Calculate the greatest common divisor of two integers.

Inc.

This file is part of the GNU MP Library.

Line 27 MA 02111-1307, USA. */

Line 28 MA 02111-1307, USA. */

#endif

#ifndef BERKELEY_MP

void

#if __STDC__

#ifndef BERKELEY_MP

mpz_gcd (mpz_ptr g, mpz_srcptr u, mpz_srcptr v)

#else

mpz_gcd (g, u, v)

mpz_ptr g;

mpz_srcptr u;

mpz_srcptr v;

#endif

#else /* BERKELEY_MP */

void

#if __STDC__

gcd (mpz_srcptr u, mpz_srcptr v, mpz_ptr g)

#else

gcd (u, v, g)

mpz_ptr g;

mpz_srcptr u;

mpz_srcptr v;

#endif

#endif /* BERKELEY_MP */

{

unsigned long int g_zero_bits, u_zero_bits, v_zero_bits;

mp_size_t g_zero_limbs, u_zero_limbs, v_zero_limbs;

Line 157 gcd (u, v, g)

Line 142 gcd (u, v, g)

if (g_zero_bits != 0)

{

mp_limb_t cy_limb;

gsize += (vp[vsize - 1] >> (BITS_PER_MP_LIMB - g_zero_bits)) != 0;

gsize += (vp[vsize - 1] >> (GMP_NUMB_BITS - g_zero_bits)) != 0;

if (g->_mp_alloc < gsize)

_mpz_realloc (g, gsize);

MPN_ZERO (g->_mp_d, g_zero_limbs);