include/borderbasix/dimpos.cc

Go to the documentation of this file.

template<typename mon, typename typB>
int has_b_index(const mon m,int index,const typB b)
{
  int res=0;
  list<mon> current_mon;
  current_mon.push_back(m);
//    cout<<"monom "<<m<<" b indx cherche "<<index<<endl;
  for(int i=0;i<index-1;i++)
    {
      list<mon> tmp;
      for(typename list<mon>::iterator iter=current_mon.begin();
          iter!=current_mon.end();iter++)
        for(int j=0;j<b.nbvars;j++)
          {
            mon tmpmon=my_divmon(*iter,j,b);
            if(IsinB(tmpmon,b))
               return 0;
            if (!member(tmp,tmpmon))
              tmp.push_back(tmpmon);
          }
      current_mon.swap(tmp);
    }
   {
      list<mon> tmp;
      for(typename list<mon>::iterator iter=current_mon.begin();
          iter!=current_mon.end();iter++)
        for(int j=0;j<b.nbvars;j++)
          {
            mon tmpmon=my_divmon(*iter,j,b);
            if (!member(tmp,tmpmon))
              tmp.push_back(tmpmon);
          }
      current_mon.swap(tmp);
   }

  for(typename list<mon>::iterator iter=current_mon.begin();
      iter!=current_mon.end();iter++)
    {
      if(IsinB(*iter,b))
        {
          res=1;
          break;
        }
    }
  return res;
}

template<typename typB>
int stop_degree(const typB &b)
{
  typedef typename typB::mon mon;
  int max_deg=0;
  int comp_index=2;
  for(typename typB::const_iterator iter=b.begin();
      iter!=b.end();iter++)
    {
      list<typename typB::mon> listppcm;
      for(int i=0;i<iter->taille2;i++)
        listppcm.push_back(iter->refuse[i]);
      for(;!listppcm.empty();comp_index++)
        {
          list<mon> tmp; 
          for(typename list<mon>::iterator iter1=listppcm.begin();
              iter1!=listppcm.end();
              iter1++)
            for(typename list<mon>::iterator iter2=(++iter1)--;
                iter2!=listppcm.end();
                iter2++)
              {
                mon tmpmon=lcm(*iter1,*iter2);
                if(!member(tmp,tmpmon))
                  tmp.push_back(tmpmon);
              }
          listppcm.erase(listppcm.begin(),listppcm.end());
          for(typename list<mon>::iterator iter=tmp.begin();
              iter!=tmp.end();
              iter++)
            {
              if(has_b_index(*iter,comp_index,b))
                {
                  listppcm.push_back(*iter);
                  if(max_deg<iter->GetDegree()-comp_index+2)
                    max_deg=iter->GetDegree()-comp_index+2;
                }
            }
        }
    }
  return max_deg;
}

---