OpenXM_contrib/gmp/mpz/remove.c - diff

Return to remove.c CVS log

Up to [local] / OpenXM_contrib / gmp / mpz

version 1.1.1.1, 2000/09/09 14:12:56

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

Line 1

/* mpz_remove -- divide out a factor and return its multiplicity.

This file is part of the GNU MP Library.

Line 23 MA 02111-1307, USA. */

#include "gmp-impl.h"

unsigned long int

#if __STDC__

mpz_remove (mpz_ptr dest, mpz_srcptr src, mpz_srcptr f)

#else

mpz_remove (dest, src, f)

mpz_ptr dest;

mpz_srcptr src;

mpz_srcptr f;

#endif

{

mpz_t fpow[40]; /* inexhaustible...until year 2020 or so */

mpz_t x, rem;

unsigned long int pwr;

int p;

if (mpz_cmp_ui (f, 1) <= 0 || mpz_sgn (src) == 0)

if (mpz_cmp_ui (f, 1) <= 0)

DIVIDE_BY_ZERO;

if (SIZ (src) == 0)

{

if (src != dest)

mpz_set (dest, src);

return 0;

}

if (mpz_cmp_ui (f, 2) == 0)

{

unsigned long int s0;