[gcc] GCC 4.5.0=>4.5.1

gcc/gmp/dumbmp.c

Index: gcc/gmp/dumbmp.c
diff --git a/gcc/gmp/dumbmp.c b/gcc/gmp/dumbmp.c
deleted file mode 100644
index cd8c67d369403f6b9a7411d1f39a218bc7d6b3fb..0000000000000000000000000000000000000000
--- a/gcc/gmp/dumbmp.c
+++ /dev/null
@@ -1,887 +0,0 @@
-/* dumbmp mini GMP compatible library.
-
-Copyright 2001, 2002, 2004 Free Software Foundation, Inc.
-
-This file is part of the GNU MP Library.
-
-The GNU MP Library is free software; you can redistribute it and/or modify
-it under the terms of the GNU Lesser General Public License as published by
-the Free Software Foundation; either version 3 of the License, or (at your
-option) any later version.
-
-The GNU MP Library is distributed in the hope that it will be useful, but
-WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
-or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
-License for more details.
-
-You should have received a copy of the GNU Lesser General Public License
-along with the GNU MP Library.  If not, see http://www.gnu.org/licenses/.  */
-
-
-/* The code here implements a subset (a very limited subset) of the main GMP
-   functions.  It's designed for use in a few build-time calculations and
-   will be slow, but highly portable.
-
-   None of the normal GMP configure things are used, nor any of the normal
-   gmp.h or gmp-impl.h.  To use this file in a program just #include
-   "dumbmp.c".
-
-   ANSI function definitions can be used here, since ansi2knr is run if
-   necessary.  But other ANSI-isms like "const" should be avoided.
-
-   mp_limb_t here is an unsigned long, since that's a sensible type
-   everywhere we know of, with 8*sizeof(unsigned long) giving the number of
-   bits in the type (that not being true for instance with int or short on
-   Cray vector systems.)
-
-   Only the low half of each mp_limb_t is used, so as to make carry handling
-   and limb multiplies easy.  GMP_LIMB_BITS is the number of bits used.  */
-
-#include <stdio.h>
-#include <stdlib.h>
-#include <string.h>
-
-
-typedef unsigned long mp_limb_t;
-
-typedef struct {
-  int        _mp_alloc;
-  int        _mp_size;
-  mp_limb_t *_mp_d;
-} mpz_t[1];
-
-#define GMP_LIMB_BITS  (sizeof (mp_limb_t) * 8 / 2)
-
-#define ABS(x)   ((x) >= 0 ? (x) : -(x))
-#define MIN(l,o) ((l) < (o) ? (l) : (o))
-#define MAX(h,i) ((h) > (i) ? (h) : (i))
-
-#define ALLOC(x) ((x)->_mp_alloc)
-#define PTR(x)   ((x)->_mp_d)
-#define SIZ(x)   ((x)->_mp_size)
-#define ABSIZ(x) ABS (SIZ (x))
-#define LOMASK   ((1L << GMP_LIMB_BITS) - 1)
-#define LO(x)    ((x) & LOMASK)
-#define HI(x)    ((x) >> GMP_LIMB_BITS)
-
-#define ASSERT(cond)                                    \
-  do {                                                  \
-    if (! (cond))                                       \
-      {                                                 \
-        fprintf (stderr, "Assertion failure\n");        \
-        abort ();                                       \
-      }                                                 \
-  } while (0)
-
-
-char *
-xmalloc (int n)
-{
-  char  *p;
-  p = malloc (n);
-  if (p == NULL)
-    {
-      fprintf (stderr, "Out of memory (alloc %d bytes)\n", n);
-      abort ();
-    }
-  return p;
-}
-
-mp_limb_t *
-xmalloc_limbs (int n)
-{
-  return (mp_limb_t *) xmalloc (n * sizeof (mp_limb_t));
-}
-
-void
-mem_copyi (char *dst, char *src, int size)
-{
-  int  i;
-  for (i = 0; i < size; i++)
-    dst[i] = src[i];
-}
-
-int
-isprime (int n)
-{
-  int  i;
-  if (n < 2)
-    return 0;
-  for (i = 2; i < n; i++)
-    if ((n % i) == 0)
-      return 0;
-  return 1;
-}
-
-int
-log2_ceil (int n)
-{
-  int  e;
-  ASSERT (n >= 1);
-  for (e = 0; ; e++)
-    if ((1 << e) >= n)
-      break;
-  return e;
-}
-
-void
-mpz_realloc (mpz_t r, int n)
-{
-  if (n <= ALLOC(r))
-    return;
-
-  ALLOC(r) = n;
-  PTR(r) = (mp_limb_t *) realloc (PTR(r), n * sizeof (mp_limb_t));
-  if (PTR(r) == NULL)
-    {
-      fprintf (stderr, "Out of memory (realloc to %d)\n", n);
-      abort ();
-    }
-}
-
-void
-mpn_normalize (mp_limb_t *rp, int *rnp)
-{
-  int  rn = *rnp;
-  while (rn > 0 && rp[rn-1] == 0)
-    rn--;
-  *rnp = rn;
-}
-
-void
-mpn_copyi (mp_limb_t *dst, mp_limb_t *src, int n)
-{
-  int  i;
-  for (i = 0; i < n; i++)
-    dst[i] = src[i];
-}
-
-void
-mpn_zero (mp_limb_t *rp, int rn)
-{
-  int  i;
-  for (i = 0; i < rn; i++)
-    rp[i] = 0;
-}
-
-void
-mpz_init (mpz_t r)
-{
-  ALLOC(r) = 1;
-  PTR(r) = xmalloc_limbs (ALLOC(r));
-  PTR(r)[0] = 0;
-  SIZ(r) = 0;
-}
-
-void
-mpz_clear (mpz_t r)
-{
-  free (PTR (r));
-  ALLOC(r) = -1;
-  SIZ (r) = 0xbadcafeL;
-  PTR (r) = (mp_limb_t *) 0xdeadbeefL;
-}
-
-int
-mpz_sgn (mpz_t a)
-{
-  return (SIZ(a) > 0 ? 1 : SIZ(a) == 0 ? 0 : -1);
-}
-
-int
-mpz_odd_p (mpz_t a)
-{
-  if (SIZ(a) == 0)
-    return 0;
-  else
-    return (PTR(a)[0] & 1) != 0;
-}
-
-int
-mpz_even_p (mpz_t a)
-{
-  if (SIZ(a) == 0)
-    return 1;
-  else
-    return (PTR(a)[0] & 1) == 0;
-}
-
-size_t
-mpz_sizeinbase (mpz_t a, int base)
-{
-  int an = ABSIZ (a);
-  mp_limb_t *ap = PTR (a);
-  int cnt;
-  mp_limb_t hi;
-
-  if (base != 2)
-    abort ();
-
-  if (an == 0)
-    return 1;
-
-  cnt = 0;
-  for (hi = ap[an - 1]; hi != 0; hi >>= 1)
-    cnt += 1;
-  return (an - 1) * GMP_LIMB_BITS + cnt;
-}
-
-void
-mpz_set (mpz_t r, mpz_t a)
-{
-  mpz_realloc (r, ABSIZ (a));
-  SIZ(r) = SIZ(a);
-  mpn_copyi (PTR(r), PTR(a), ABSIZ (a));
-}
-
-void
-mpz_init_set (mpz_t r, mpz_t a)
-{
-  mpz_init (r);
-  mpz_set (r, a);
-}
-
-void
-mpz_set_ui (mpz_t r, unsigned long ui)
-{
-  int  rn;
-  mpz_realloc (r, 2);
-  PTR(r)[0] = LO(ui);
-  PTR(r)[1] = HI(ui);
-  rn = 2;
-  mpn_normalize (PTR(r), &rn);
-  SIZ(r) = rn;
-}
-
-void
-mpz_init_set_ui (mpz_t r, unsigned long ui)
-{
-  mpz_init (r);
-  mpz_set_ui (r, ui);
-}
-
-void
-mpz_setbit (mpz_t r, unsigned long bit)
-{
-  int        limb, rn, extend;
-  mp_limb_t  *rp;
-
-  rn = SIZ(r);
-  if (rn < 0)
-    abort ();  /* only r>=0 */
-
-  limb = bit / GMP_LIMB_BITS;
-  bit %= GMP_LIMB_BITS;
-
-  mpz_realloc (r, limb+1);
-  rp = PTR(r);
-  extend = (limb+1) - rn;
-  if (extend > 0)
-    mpn_zero (rp + rn, extend);
-
-  rp[limb] |= (mp_limb_t) 1 << bit;
-  SIZ(r) = MAX (rn, limb+1);
-}
-
-int
-mpz_tstbit (mpz_t r, unsigned long bit)
-{
-  int  limb;
-
-  if (SIZ(r) < 0)
-    abort ();  /* only r>=0 */
-
-  limb = bit / GMP_LIMB_BITS;
-  if (SIZ(r) <= limb)
-    return 0;
-
-  bit %= GMP_LIMB_BITS;
-  return (PTR(r)[limb] >> bit) & 1;
-}
-
-int
-popc_limb (mp_limb_t a)
-{
-  int  ret = 0;
-  while (a != 0)
-    {
-      ret += (a & 1);
-      a >>= 1;
-    }
-  return ret;
-}
-
-unsigned long
-mpz_popcount (mpz_t a)
-{
-  unsigned long  ret;
-  int            i;
-
-  if (SIZ(a) < 0)
-    abort ();
-
-  ret = 0;
-  for (i = 0; i < SIZ(a); i++)
-    ret += popc_limb (PTR(a)[i]);
-  return ret;
-}
-
-void
-mpz_add (mpz_t r, mpz_t a, mpz_t b)
-{
-  int an = ABSIZ (a), bn = ABSIZ (b), rn;
-  mp_limb_t *rp, *ap, *bp;
-  int i;
-  mp_limb_t t, cy;
-
-  if ((SIZ (a) ^ SIZ (b)) < 0)
-    abort ();			/* really subtraction */
-  if (SIZ (a) < 0)
-    abort ();
-
-  mpz_realloc (r, MAX (an, bn) + 1);
-  ap = PTR (a);  bp = PTR (b);  rp = PTR (r);
-  if (an < bn)
-    {
-      mp_limb_t *tp;  int tn;
-      tn = an; an = bn; bn = tn;
-      tp = ap; ap = bp; bp = tp;
-    }
-
-  cy = 0;
-  for (i = 0; i < bn; i++)
-    {
-      t = ap[i] + bp[i] + cy;
-      rp[i] = LO (t);
-      cy = HI (t);
-    }
-  for (i = bn; i < an; i++)
-    {
-      t = ap[i] + cy;
-      rp[i] = LO (t);
-      cy = HI (t);
-    }
-  rp[an] = cy;
-  rn = an + 1;
-
-  mpn_normalize (rp, &rn);
-  SIZ (r) = rn;
-}
-
-void
-mpz_add_ui (mpz_t r, mpz_t a, unsigned long int ui)
-{
-  mpz_t b;
-
-  mpz_init (b);
-  mpz_set_ui (b, ui);
-  mpz_add (r, a, b);
-  mpz_clear (b);
-}
-
-void
-mpz_sub (mpz_t r, mpz_t a, mpz_t b)
-{
-  int an = ABSIZ (a), bn = ABSIZ (b), rn;
-  mp_limb_t *rp, *ap, *bp;
-  int i;
-  mp_limb_t t, cy;
-
-  if ((SIZ (a) ^ SIZ (b)) < 0)
-    abort ();			/* really addition */
-  if (SIZ (a) < 0)
-    abort ();
-
-  mpz_realloc (r, MAX (an, bn) + 1);
-  ap = PTR (a);  bp = PTR (b);  rp = PTR (r);
-  if (an < bn)
-    {
-      mp_limb_t *tp;  int tn;
-      tn = an; an = bn; bn = tn;
-      tp = ap; ap = bp; bp = tp;
-    }
-
-  cy = 0;
-  for (i = 0; i < bn; i++)
-    {
-      t = ap[i] - bp[i] - cy;
-      rp[i] = LO (t);
-      cy = LO (-HI (t));
-    }
-  for (i = bn; i < an; i++)
-    {
-      t = ap[i] - cy;
-      rp[i] = LO (t);
-      cy = LO (-HI (t));
-    }
-  rp[an] = cy;
-  rn = an + 1;
-
-  if (cy != 0)
-    {
-      cy = 0;
-      for (i = 0; i < rn; i++)
-	{
-	  t = -rp[i] - cy;
-	  rp[i] = LO (t);
-	  cy = LO (-HI (t));
-	}
-      SIZ (r) = -rn;
-      return;
-    }
-
-  mpn_normalize (rp, &rn);
-  SIZ (r) = rn;
-}
-
-void
-mpz_sub_ui (mpz_t r, mpz_t a, unsigned long int ui)
-{
-  mpz_t b;
-
-  mpz_init (b);
-  mpz_set_ui (b, ui);
-  mpz_sub (r, a, b);
-  mpz_clear (b);
-}
-
-void
-mpz_mul (mpz_t r, mpz_t a, mpz_t b)
-{
-  int an = ABSIZ (a), bn = ABSIZ (b), rn;
-  mp_limb_t *scratch, *tmp, *ap = PTR (a), *bp = PTR (b);
-  int ai, bi;
-  mp_limb_t t, cy;
-
-  scratch = xmalloc_limbs (an + bn);
-  tmp = scratch;
-
-  for (bi = 0; bi < bn; bi++)
-    tmp[bi] = 0;
-
-  for (ai = 0; ai < an; ai++)
-    {
-      tmp = scratch + ai;
-      cy = 0;
-      for (bi = 0; bi < bn; bi++)
-	{
-	  t = ap[ai] * bp[bi] + tmp[bi] + cy;
-	  tmp[bi] = LO (t);
-	  cy = HI (t);
-	}
-      tmp[bn] = cy;
-    }
-
-  rn = an + bn;
-  mpn_normalize (scratch, &rn);
-  free (PTR (r));
-  PTR (r) = scratch;
-  SIZ (r) = (SIZ (a) ^ SIZ (b)) >= 0 ? rn : -rn;
-  ALLOC (r) = an + bn;
-}
-
-void
-mpz_mul_ui (mpz_t r, mpz_t a, unsigned long int ui)
-{
-  mpz_t b;
-
-  mpz_init (b);
-  mpz_set_ui (b, ui);
-  mpz_mul (r, a, b);
-  mpz_clear (b);
-}
-
-void
-mpz_mul_2exp (mpz_t r, mpz_t a, unsigned long int bcnt)
-{
-  mpz_set (r, a);
-  while (bcnt)
-    {
-      mpz_add (r, r, r);
-      bcnt -= 1;
-    }
-}
-
-void
-mpz_ui_pow_ui (mpz_t r, unsigned long b, unsigned long e)
-{
-  unsigned long  i;
-  mpz_t          bz;
-
-  mpz_init (bz);
-  mpz_set_ui (bz, b);
-
-  mpz_set_ui (r, 1L);
-  for (i = 0; i < e; i++)
-    mpz_mul (r, r, bz);
-
-  mpz_clear (bz);
-}
-
-void
-mpz_tdiv_q_2exp (mpz_t r, mpz_t a, unsigned long int bcnt)
-{
-  int as, rn;
-  int cnt, tnc;
-  int lcnt;
-  mp_limb_t high_limb, low_limb;
-  int i;
-
-  as = SIZ (a);
-  lcnt = bcnt / GMP_LIMB_BITS;
-  rn = ABS (as) - lcnt;
-  if (rn <= 0)
-    SIZ (r) = 0;
-  else
-    {
-      mp_limb_t *rp, *ap;
-
-      mpz_realloc (r, rn);
-
-      rp = PTR (r);
-      ap = PTR (a);
-
-      cnt = bcnt % GMP_LIMB_BITS;
-      if (cnt != 0)
-        {
-	  ap += lcnt;
-	  tnc = GMP_LIMB_BITS - cnt;
-	  high_limb = *ap++;
-	  low_limb = high_limb >> cnt;
-
-	  for (i = rn - 1; i != 0; i--)
-	    {
-	      high_limb = *ap++;
-	      *rp++ = low_limb | LO (high_limb << tnc);
-	      low_limb = high_limb >> cnt;
-	    }
-	  *rp = low_limb;
-          rn -= low_limb == 0;
-        }
-      else
-        {
-	  ap += lcnt;
-          mpn_copyi (rp, ap, rn);
-        }
-
-      SIZ (r) = as >= 0 ? rn : -rn;
-    }
-}
-
-void
-mpz_tdiv_r_2exp (mpz_t r, mpz_t a, unsigned long int bcnt)
-{
-  int    rn, bwhole;
-
-  mpz_set (r, a);
-  rn = ABSIZ(r);
-
-  bwhole = bcnt / GMP_LIMB_BITS;
-  bcnt %= GMP_LIMB_BITS;
-  if (rn > bwhole)
-    {
-      rn = bwhole+1;
-      PTR(r)[rn-1] &= ((mp_limb_t) 1 << bcnt) - 1;
-      mpn_normalize (PTR(r), &rn);
-      SIZ(r) = (SIZ(r) >= 0 ? rn : -rn);
-    }
-}
-
-int
-mpz_cmp (mpz_t a, mpz_t b)
-{
-  mp_limb_t *ap, *bp, al, bl;
-  int as = SIZ (a), bs = SIZ (b);
-  int i;
-  int sign;
-
-  if (as != bs)
-    return as > bs ? 1 : -1;
-
-  sign = as > 0 ? 1 : -1;
-
-  ap = PTR (a);
-  bp = PTR (b);
-  for (i = ABS (as) - 1; i >= 0; i--)
-    {
-      al = ap[i];
-      bl = bp[i];
-      if (al != bl)
-	return al > bl ? sign : -sign;
-    }
-  return 0;
-}
-
-int
-mpz_cmp_ui (mpz_t a, unsigned long b)
-{
-  mpz_t  bz;
-  int    ret;
-  mpz_init_set_ui (bz, b);
-  ret = mpz_cmp (a, bz);
-  mpz_clear (bz);
-  return ret;
-}
-
-void
-mpz_tdiv_qr (mpz_t q, mpz_t r, mpz_t a, mpz_t b)
-{
-  mpz_t          tmpr, tmpb;
-  unsigned long  cnt;
-
-  ASSERT (mpz_sgn(a) >= 0);
-  ASSERT (mpz_sgn(b) > 0);
-
-  mpz_init_set (tmpr, a);
-  mpz_init_set (tmpb, b);
-  mpz_set_ui (q, 0L);
-
-  if (mpz_cmp (tmpr, tmpb) > 0)
-    {
-      cnt = mpz_sizeinbase (tmpr, 2) - mpz_sizeinbase (tmpb, 2) + 1;
-      mpz_mul_2exp (tmpb, tmpb, cnt);
-
-      for ( ; cnt > 0; cnt--)
-        {
-          mpz_mul_2exp (q, q, 1);
-          mpz_tdiv_q_2exp (tmpb, tmpb, 1L);
-          if (mpz_cmp (tmpr, tmpb) >= 0)
-            {
-              mpz_sub (tmpr, tmpr, tmpb);
-              mpz_add_ui (q, q, 1L);
-              ASSERT (mpz_cmp (tmpr, tmpb) < 0);
-            }
-        }
-    }
-
-  mpz_set (r, tmpr);
-  mpz_clear (tmpr);
-  mpz_clear (tmpb);
-}
-
-void
-mpz_tdiv_qr_ui (mpz_t q, mpz_t r, mpz_t a, unsigned long b)
-{
-  mpz_t  bz;
-  mpz_init_set_ui (bz, b);
-  mpz_tdiv_qr (q, r, a, bz);
-  mpz_clear (bz);
-}
-
-void
-mpz_tdiv_q (mpz_t q, mpz_t a, mpz_t b)
-{
-  mpz_t  r;
-
-  mpz_init (r);
-  mpz_tdiv_qr (q, r, a, b);
-  mpz_clear (r);
-}
-
-void
-mpz_tdiv_q_ui (mpz_t q, mpz_t n, unsigned long d)
-{
-  mpz_t  dz;
-  mpz_init_set_ui (dz, d);
-  mpz_tdiv_q (q, n, dz);
-  mpz_clear (dz);
-}
-
-/* Set inv to the inverse of d, in the style of invert_limb, ie. for
-   udiv_qrnnd_preinv.  */
-void
-mpz_preinv_invert (mpz_t inv, mpz_t d, int numb_bits)
-{
-  mpz_t  t;
-  int    norm;
-  ASSERT (SIZ(d) > 0);
-
-  norm = numb_bits - mpz_sizeinbase (d, 2);
-  ASSERT (norm >= 0);
-  mpz_init_set_ui (t, 1L);
-  mpz_mul_2exp (t, t, 2*numb_bits - norm);
-  mpz_tdiv_q (inv, t, d);
-  mpz_set_ui (t, 1L);
-  mpz_mul_2exp (t, t, numb_bits);
-  mpz_sub (inv, inv, t);
-
-  mpz_clear (t);
-}
-
-/* Remove leading '0' characters from the start of a string, by copying the
-   remainder down. */
-void
-strstrip_leading_zeros (char *s)
-{
-  char  c, *p;
-
-  p = s;
-  while (*s == '0')
-    s++;
-
-  do
-    {
-      c = *s++;
-      *p++ = c;
-    }
-  while (c != '\0');
-}
-
-char *
-mpz_get_str (char *buf, int base, mpz_t a)
-{
-  static char  tohex[] = "0123456789abcdef";
-
-  mp_limb_t  alimb, *ap;
-  int        an, bn, i, j;
-  char       *bp;
-
-  if (base != 16)
-    abort ();
-  if (SIZ (a) < 0)
-    abort ();
-
-  if (buf == 0)
-    buf = xmalloc (ABSIZ (a) * (GMP_LIMB_BITS / 4) + 3);
-
-  an = ABSIZ (a);
-  if (an == 0)
-    {
-      buf[0] = '0';
-      buf[1] = '\0';
-      return buf;
-    }
-
-  ap = PTR (a);
-  bn = an * (GMP_LIMB_BITS / 4);
-  bp = buf + bn;
-
-  for (i = 0; i < an; i++)
-    {
-      alimb = ap[i];
-      for (j = 0; j < GMP_LIMB_BITS / 4; j++)
-        {
-          bp--;
-          *bp = tohex [alimb & 0xF];
-          alimb >>= 4;
-        }
-      ASSERT (alimb == 0);
-    }
-  ASSERT (bp == buf);
-
-  buf[bn] = '\0';
-
-  strstrip_leading_zeros (buf);
-  return buf;
-}
-
-void
-mpz_out_str (FILE *file, int base, mpz_t a)
-{
-  char *str;
-
-  if (file == 0)
-    file = stdout;
-
-  str = mpz_get_str (0, 16, a);
-  fputs (str, file);
-  free (str);
-}
-
-/* Calculate r satisfying r*d == 1 mod 2^n. */
-void
-mpz_invert_2exp (mpz_t r, mpz_t a, unsigned long n)
-{
-  unsigned long  i;
-  mpz_t  inv, prod;
-
-  ASSERT (mpz_odd_p (a));
-
-  mpz_init_set_ui (inv, 1L);
-  mpz_init (prod);
-
-  for (i = 1; i < n; i++)
-    {
-      mpz_mul (prod, inv, a);
-      if (mpz_tstbit (prod, i) != 0)
-        mpz_setbit (inv, i);
-    }
-
-  mpz_mul (prod, inv, a);
-  mpz_tdiv_r_2exp (prod, prod, n);
-  ASSERT (mpz_cmp_ui (prod, 1L) == 0);
-
-  mpz_set (r, inv);
-
-  mpz_clear (inv);
-  mpz_clear (prod);
-}
-
-/* Calculate inv satisfying r*a == 1 mod 2^n. */
-void
-mpz_invert_ui_2exp (mpz_t r, unsigned long a, unsigned long n)
-{
-  mpz_t  az;
-  mpz_init_set_ui (az, a);
-  mpz_invert_2exp (r, az, n);
-  mpz_clear (az);
-}
-
-/* x=y^z */
-void
-mpz_pow_ui (mpz_t x, mpz_t y, unsigned long z)
-{
-  mpz_t t;
-
-  mpz_init_set_ui (t, 1);
-  for (; z != 0; z--)
-    mpz_mul (t, t, y);
-  mpz_set (x, t);
-  mpz_clear (t);
-}
-
-/* x=x+y*z */
-void
-mpz_addmul_ui (mpz_t x, mpz_t y, unsigned long z)
-{
-  mpz_t t;
-
-  mpz_init (t);
-  mpz_mul_ui (t, y, z);
-  mpz_add (x, x, t);
-  mpz_clear (t);
-}
-
-/* x=floor(y^(1/z)) */
-void
-mpz_root (mpz_t x, mpz_t y, unsigned long z)
-{
-  mpz_t t, u;
-
-  if (mpz_sgn (y) < 0)
-    {
-      fprintf (stderr, "mpz_root does not accept negative values\n");
-      abort ();
-    }
-  if (mpz_cmp_ui (y, 1) <= 0)
-    {
-      mpz_set (x, y);
-      return;
-    }
-  mpz_init (t);
-  mpz_init_set (u, y);
-  do
-    {
-      mpz_pow_ui (t, u, z - 1);
-      mpz_tdiv_q (t, y, t);
-      mpz_addmul_ui (t, u, z - 1);
-      mpz_tdiv_q_ui (t, t, z);
-      if (mpz_cmp (t, u) >= 0)
-	break;
-      mpz_set (u, t);
-    }
-  while (1);
-  mpz_set (x, u);
-  mpz_clear (t);
-  mpz_clear (u);
-}