include/borderbasix/proj.cc

Go to the documentation of this file.

template<typename typPk,typename typdump,typename Base> 
void rewrite(typPk &Pk,const typdump &dump,const Base &b)
{
  typedef typename typPk::value_type polyalp;
  typedef typename polyalp::monom_t mon;
  typedef typename typdump::value_type dumpstruct;
  typedef typename dumpstruct::pol_t pol;
  typename typPk::iterator iter=Pk.begin();
  //cout<<"Pk.size dans rewrite "<<Pk.size()<<endl;
  for(;iter!=Pk.end();iter++)
    //pour chaque poly de Pk
    { 
      polyalp tmp(mon(0));
      for(typename polyalp::iterator itermon=iter->begin();
          itermon!=iter->end();itermon++)
      //pour chaque mon du polynome *iter
      {
          if (!IsinB(*itermon,b))
            {
              pol ttmp;
              ttmp.size=-1;
              ttmp.nf=NULL;
              //cout<<"findcorresppol "<<*itermon<<endl;
              findcorresppol(*itermon,dump,ttmp);
              //cout<<"invconv"<<endl;
              //fincorresppol ne touche pas son parametre de sortie qd 
              //le monome recherche n'existe pas
              //il existe ce monome?
              if(ttmp.size!=-1)
                {
                  //  cout<<"trouve "<<endl;
                  //tmp*=mon(ttmp.ind.GetCoeff());
                  tmp+=invconv<polyalp>(ttmp)*mon(itermon->GetCoeff()/ttmp.ind.GetCoeff());
                }
                }
        
      }
      //cout<<"le pol que je vais enlever"<<endl;
      //cout<<tmp<<endl;
      //cout<<"fin du pol"<<endl;
      //cout<<"le poly avant reecriture "<<*iter<<endl;
      (*iter)-=tmp;
      //cout<<"poly apres reecriture!!\n"<<(*iter)<<endl;
    }
  return;
}
#ifdef uninf
void reinit_multiple();
#endif
template<typename T,typename typdump,typename Base> 
void proj(T &Crochet,const typdump &dump, const Base &b)
{
  typedef typename typdump::value_type::pol_t pol;
  typedef typename T::value_type polyalp;
  typedef typename polyalp::monom_t mon;
  typedef typename mon::coeff_t coeff;
  typename T::iterator iter=Crochet.begin();
#ifdef unimul 
  reinit_multiple();
#endif
#ifdef uninf
  reinit_nf<pol>();
#endif
  for(;iter!=Crochet.end();iter++)
    {
      list<mon> monared;
      //cout<<"le pol que je vais reduire"<<endl;
      //cout<<*iter<<endl;
      for(typename polyalp::iterator iterpol=iter->begin();iterpol!=iter->end();iterpol++)
        {
          if (!IsinB(*iterpol,b)) monared.push_back(*iterpol);
          //cout<<"mon dans proj "<<*iterpol<<" b? "<<IsinB(*iterpol,b)<<endl;
        }
      for(typename list<mon>::iterator itermon=monared.begin();itermon!=monared.end();itermon++)
        {
          //  cout<<"dans le for sur monared "<<b.nbvar()<<endl;
#ifdef projfast
          polyalp result(0);
          pol accumul;
          mon *exce;
          int sizeex;
          accumul.size=1;
          accumul.nf=(coeff*)MAC_REV_MALLOC<coeff>(sizeof(coeff));
          accumul.nf[0]=1;
          for(int i=0;i<b.nbvar();i++)
            for(int j=0;j<itermon->GetDegree(i);j++)
              {
                pol tmp;
                //if(result.size()!=0)
                result*=mon(i,1);
                tmp=mult(i,accumul,dump,b,&exce,sizeex);
                MAC_REV_FREE<coeff>(accumul.nf,accumul.size*sizeof(coeff));
                accumul=tmp;
                for(int k=0;k<sizeex-1;k++)
                  result+=exce[k];
                MAC_REV_FREE<mon>(exce,sizeex*sizeof(mon));
              }
          for(int k=0;k<accumul.size;k++)
            if(!Iszero(accumul.nf[k]))
              {
                mon tmpmon;
                int2mon(k,tmpmon);
                //      cout<<"tmpmon "<<tmpmon<<endl;
                result+=tmpmon*accumul.nf[k];
              }
          MAC_REV_FREE<coeff>(accumul.nf,accumul.size*sizeof(coeff));
#else
          polyalp result(1);
          for(int i=0;i<b.nbvar();i++)
            for(int j=0;j<itermon->GetDegree(i);j++)
              {
                result*=mon(i,1);
                int flag=0;
                //cout<<"nbtour dans la boucle "<<itermon->GetDegree(i)<<"J"<<j<<endl;
                for (typename polyalp::iterator iterres=result.begin();
                     iterres!=result.end();iterres++)
                 {
                   if (!IsinB(*iterres,b)) flag=1;
                   // cout<<"IsinB("<<*iterres<<") "<<IsinB(*iterres,b)<<endl;
                 }
                if (flag)
                  {
                    T tmp;
                    tmp.push_back(result);
                    //cout<<"rewrite"<<endl;
                    rewrite(tmp,dump,b);
                    //cout<<"sortie rewrite"<<endl;
                    result=tmp.front();
                  }
              }
#endif
          coeff Coeff;
          for(typename polyalp::iterator iterpol=iter->begin();
              iterpol!=iter->end();iterpol++)
            if ((itermon->rep)==(iterpol->rep))
              {
                Coeff=iterpol->GetCoeff();
                break;
              }
          //cout<<"iter avant reecriture "<<*iter<<endl;
          *iter+=result*mon(Coeff);
          *iter-=polyalp(mon((*itermon)*(coeff)(Coeff/itermon->GetCoeff())));
          //cout<<"iter apres reecriture "<<*iter<<endl<<endl;
        }
    }
  //stab(Crochet);
  return;

}

---