/Users/mourrain/Devel/mmx/shape/include/shape/algebraic_set.hpp

Go to the documentation of this file.

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

# include <realroot/GMP.hpp>
# include <realroot/polynomial_sparse.hpp>
# define TMPL  template<class C, class V>
# define TMPL1 template<class V>
# define SELF  algebraic_set<C,V>
//====================================================================
namespace mmx {
namespace shape {
//====================================================================
struct numeric_def {};

template<> struct use<numeric_def> 
{  
  typedef ::mmx::GMP::integer  Integer;
  typedef ::mmx::GMP::rational Rational;
  typedef ::mmx::GMP::floating Floating;
};

struct algebraic_set_def {};

template<> struct use<algebraic_set_def> 
{  
  typedef use<numeric_def>::Rational         Scalar;
  //  typedef double Scalar;
  typedef polynomial<Scalar, with<Sparse, DegRevLex> > Polynomial;
};

TMPL struct use<algebraic_set_def,join2<C,V> > 
{  
  typedef typename use<numeric_def,V>::Rational         Scalar;
  //  typedef C Scalar;
  typedef polynomial<Scalar, with<Sparse, DegRevLex> > Polynomial;
};
//--------------------------------------------------------------------
TMPL
class algebraic_set 
{
public:
  //  typedef polynomial<Rational, with<Sparse, DegRevLex> > Polynomial;
  typedef typename use<algebraic_set_def,join2<C,V> >::Polynomial   Polynomial;

  algebraic_set(void) {};
  algebraic_set(const char* s1) ;
  algebraic_set(const char* s1, const char* s2) ;
  algebraic_set(const Polynomial& eq) ;
  algebraic_set(const Polynomial& eq1, const Polynomial& eq2) ;
  ~algebraic_set(void) {};
  
  const Seq<Polynomial>& equations()  const { return this->m_polynomials ; }
  const Polynomial&      equation(int i=0) const { return this->m_polynomials[i] ; }
  Polynomial&            equation(int i=0)       { return this->m_polynomials[i] ; }
  void         set_equation(const Polynomial& p) { this->m_polynomials[0]=p ; }
  void insert(const Polynomial& p)   { this->m_polynomials<<p; }
  int  nbequation()            const { return m_polynomials.size(); }

private:
  Seq<Polynomial> m_polynomials ;
} ;

TMPL
SELF::algebraic_set(const char* s1) {
  m_polynomials<<Polynomial(s1);
}

TMPL
SELF::algebraic_set(const char* s1,const char* s2) {
  m_polynomials<<Polynomial(s1);
  m_polynomials<<Polynomial(s2);
}

TMPL
SELF::algebraic_set(const Polynomial& p1) {
  m_polynomials<<p1;
}

TMPL
SELF::algebraic_set(const Polynomial& p1,const Polynomial& p2) {
  m_polynomials<<p1;
  m_polynomials<<p2;
}

//====================================================================
} ; // namespace shape
} ; // namespace mmx
//====================================================================
# undef TMPL
# undef TMPL1
# undef SELF
# endif // IMPLICITCURVE_H