[BACK]Return to gcd_ui.c CVS log [TXT][DIR] Up to [local] / OpenXM_contrib / gmp / mpz

Diff for /OpenXM_contrib/gmp/mpz/Attic/gcd_ui.c between version 1.1.1.2 and 1.1.1.3

version 1.1.1.2, 2000/09/09 14:12:52 version 1.1.1.3, 2003/08/25 16:06:33
Line 1 
Line 1 
 /* mpz_gcd_ui -- Calculate the greatest common divisior of two integers.  /* mpz_gcd_ui -- Calculate the greatest common divisior of two integers.
   
 Copyright (C) 1994, 1996, 1999, 2000 Free Software Foundation, Inc.  Copyright 1994, 1996, 1999, 2000, 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 24  MA 02111-1307, USA. */
Line 24  MA 02111-1307, USA. */
 #include "gmp-impl.h"  #include "gmp-impl.h"
   
 unsigned long int  unsigned long int
 #if __STDC__  
 mpz_gcd_ui (mpz_ptr w, mpz_srcptr u, unsigned long int v)  mpz_gcd_ui (mpz_ptr w, mpz_srcptr u, unsigned long int v)
 #else  
 mpz_gcd_ui (w, u, v)  
      mpz_ptr w;  
      mpz_srcptr u;  
      unsigned long int v;  
 #endif  
 {  {
   mp_size_t size;    mp_size_t un;
   mp_limb_t res;    mp_limb_t res;
   
   size = ABS (u->_mp_size);  #if GMP_NAIL_BITS != 0
     if (v > GMP_NUMB_MAX)
       {
         mpz_t vz;
         mp_limb_t vlimbs[2];
         vlimbs[0] = v & GMP_NUMB_MASK;
         vlimbs[1] = v >> GMP_NUMB_BITS;
         PTR(vz) = vlimbs;
         SIZ(vz) = 2;
         mpz_gcd (w, u, vz);
         return;
       }
   #endif
   
   if (size == 0)    un = ABSIZ(u);
   
     if (un == 0)
     res = v;      res = v;
   else if (v == 0)    else if (v == 0)
     {      {
       if (w != NULL && u != w)        if (w != NULL)
         {          {
           if (w->_mp_alloc < size)            if (u != w)
             _mpz_realloc (w, size);              {
                 MPZ_REALLOC (w, un);
           MPN_COPY (w->_mp_d, u->_mp_d, size);                MPN_COPY (PTR(w), PTR(u), un);
               }
             SIZ(w) = un;
         }          }
       w->_mp_size = size;  
       /* We can't return any useful result for gcd(big,0).  */        /* We can't return any useful result for gcd(big,0).  */
       return size > 1 ? 0 : w->_mp_d[0];        return un > 1 ? 0 : PTR(u)[0];
     }      }
   else    else
     res = mpn_gcd_1 (u->_mp_d, size, (mp_limb_t) v);      res = mpn_gcd_1 (PTR(u), un, (mp_limb_t) v);
   
   if (w != NULL)    if (w != NULL)
     {      {
       w->_mp_d[0] = res;        PTR(w)[0] = res;
       w->_mp_size = 1;        SIZ(w) = res != 0;
     }      }
   return res;    return res;
 }  }
Legend:
Removed from v.1.1.1.2  
changed lines
  Added in v.1.1.1.3
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>