include/shape/semialgebraic_set.hpp

Go to the documentation of this file.

/*****************************************************************************
 * M a t h e m a g i x
 *****************************************************************************
 * AlgebraicCurve
 * 2008-03-28
 * Julien Wintz & Bernard Mourrain
 *****************************************************************************
 *               Copyright (C) 2006 INRIA Sophia-Antipolis
 *****************************************************************************
 * Comments :
 ****************************************************************************/

# ifndef shape_semialgebraic_set_hpp
# define shape_semialgebraic_set_hpp

# include <shape/algebraic_set.hpp>
# define TMPL template<class K>
# define SemialgebraicSet semialgebraic_set<K>
//====================================================================
namespace mmx {

namespace shape {
TMPL struct semialgebraic_set_def { 
  typedef typename algebraic_set_def<K>::Polynomial Polynomial;
};

TMPL
class semialgebraic_set
{
public:
  typedef typename semialgebraic_set_def<K>::Polynomial Polynomial;

  semialgebraic_set(void) {};
  semialgebraic_set(const char* str1, int s1) ;
  semialgebraic_set(const char* str1, int s1, const char* str2, int s2) ;
  semialgebraic_set(const Polynomial& eq, int s) ;
  semialgebraic_set(const Polynomial& eq1, int s1, const Polynomial& eq2, int s2) ;
  ~semialgebraic_set(void) {};
  
  Seq<Polynomial> inequations()  const { return this->m_polynomials ; }

  Polynomial inequation(int i=0) const { return this->m_polynomials[i] ; }
  int        sign(int i=0)       const { return this->m_signs[i] ; }

  void insert(const Polynomial& p, int s)   { 
    this->m_polynomials<<p; this->m_signs<<s; 
  }
  int  nbineq()                  const { return m_polynomials.size(); }

private:
  Seq<Polynomial> m_polynomials ;
  Seq<int>        m_signs;
} ;

TMPL
SemialgebraicSet::semialgebraic_set(const char* s1,int i1) {
  m_polynomials<<Polynomial(s1);
  m_signs<<i1;
}

TMPL
SemialgebraicSet::semialgebraic_set(const char* s1, int i1, const char* s2, int i2) {
    m_polynomials<< (i1<0? (-1)*Polynomial(s1) : s1) ;
    m_polynomials<< (i2<0? (-1)*Polynomial(s2) : s2) ;
}

TMPL
SemialgebraicSet::semialgebraic_set(const Polynomial& p1, int i1) {
    m_polynomials<< (i1<0? (-1)*p1 : p1) ;
  //m_signs<<i1;
}

TMPL
SemialgebraicSet::semialgebraic_set(const Polynomial& p1, int i1, const Polynomial& p2, int i2) {
  m_polynomials<<p1;
  m_polynomials<<p2;
}


} ; // namespace shape
} ; // namespace mmx

# undef TMPL
# undef SemialgebraicSet
# endif // shape_semialgebraic_set_hpp

---