Function FindUnivariatePolys

kernel/fglm/fglmzero.cc:1117–1189

//<-

static ideal
FindUnivariatePolys( const idealFunctionals & l )
{
  fglmVector v;
  fglmVector p;
  ideal destIdeal = idInit( (currRing->N), 1 );

  int i;
  BOOLEAN isZero;
  int *varpermutation = (int*)omAlloc( ((currRing->N)+1)*sizeof(int) );
  ideal perm = idMaxIdeal(1);
  intvec *iv = idSort(perm,TRUE);
  idDelete(&perm);
  for(i = (currRing->N); i > 0; i--) varpermutation[(currRing->N)+1-i] = (*iv)[i-1];
  delete iv;

  for (i= 1; i <= (currRing->N); i++ )
  {
    // main loop
    STICKYPROT2( "(%i)", i /*varpermutation[i]*/);
    gaussReducer gauss( l.dimen() );
    isZero= FALSE;
    v= fglmVector( l.dimen(), 1 );
    while ( !isZero )
    {
      if ( (isZero= gauss.reduce( v )))
      {
        STICKYPROT( "+" );
        p= gauss.getDependence();
        number gcd= p.gcd();
        if ( ! nIsOne( gcd ) )
        {
          p /= gcd;
        }
        nDelete( & gcd );
        int k;
        poly temp = NULL;
        poly result=NULL;
        for ( k= p.size(); k > 0; k-- )
        {
          number n = p.getconstelem( k );
          if ( ! nIsZero( n ) )
          {
            if ( temp == NULL )
            {
              result= p_Init(currRing);
              temp= result;
            }
            else
            {
              temp->next= p_Init(currRing);
              pIter( temp );
            }
            pSetCoeff0( temp, nCopy(n) );
            pSetExp( temp, i /*varpermutation[i]*/, k-1 );
            pSetm( temp );
          }
        }