dimpos.hpp

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;
}