Iszero.hpp

Go to the documentation of this file.

#ifndef ALREADY_Iszero
#define ALREADY_Iszero

#define SEUIL 9.9e-20
//#define SEUIL 1e-100
#define PHTABS(x) ((x>(typcoeff)0)?x:(typcoeff)(-1)*x)

#ifdef MON_PREMIER
template<typename T>
inline int Iszero(const Z<MON_PREMIER,T> & c);
#endif

#if 1
    int Iszero(const Scl<MPQ> &c)
   {
      return c==(int)0;
    }
#endif

template<typename typcoeff>
inline int Iszero(const typcoeff & c)
{
  int res=0;
  if(PHTABS(c)<(typcoeff)SEUIL)
    res=1;
  return res;
}

template<typename pol>
int Ispolzero(const pol & p)
{
  int res=1;
  for(int i=0;i<p.size;i++)
    if (!Iszero(p.nf[i])) {res=0;break;}
  return res;
}

template<typename typP>
//cf corealgo ligne 310
void selectnozero(typP & ttmp)
{
  typedef typename typP::value_type polyalp;
  typP res;
  for(typename typP::iterator iterttmp=ttmp.begin();iterttmp!=ttmp.end();iterttmp++)
    {
      int yes=0;
      for(typename polyalp::iterator iter=iterttmp->begin();
      iter!=iterttmp->end();iter++)
    if(!Iszero(iter->GetCoeff())) yes=1;
      if (yes)
    res.push_back(*iterttmp);
    }
  ttmp.swap(res);
  return;
}
#undef SEUIL
#undef PHTABS
#endif //ALREADY_Iszero