Method gcd

kernel/fglm/fglmvec.cc:458–503

}

number fglmVector::gcd () const
{
  int i = rep->size ();
  BOOLEAN found = FALSE;
  BOOLEAN gcdIsOne = FALSE;
  number theGcd;
  number current;
  while(i > 0 && !found)
  {
    current = rep->getconstelem (i);
    if(!nIsZero (current))
    {
      theGcd = nCopy (current);
      found = TRUE;
      if(!nGreaterZero (theGcd))
      {
        theGcd = nInpNeg (theGcd);
      }
      if(nIsOne (theGcd))
        gcdIsOne = TRUE;
    }
    i--;
  }
  if(found)
  {
    if (getCoeffType(currRing->cf)!=n_Zp)
    {
      while(i > 0 && !gcdIsOne)
      {
        current = rep->getconstelem (i);
        if(!nIsZero (current))
        {
          number temp = n_SubringGcd (theGcd, current, currRing->cf);
          nDelete (&theGcd);
          theGcd = temp;
          if(nIsOne (theGcd))
            gcdIsOne = TRUE;
        }
        i--;
      }
    }
  }
  else
    theGcd = nInit (0);
  return theGcd;
}

number fglmVector::clearDenom ()
{