include/shape/algebraic_curve.hpp

Go to the documentation of this file.

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


# include <shape/curve.hpp>
# include <shape/bounding_box.hpp>
# include <shape/algebraic_set.hpp>

# define TMPL template<class C, class V>
# define TMPL1 template<class V>
# define SELF algebraic_curve<C,V>
# define Viewer viewer<axel,V>
//====================================================================
namespace mmx {
namespace shape {

// TMPL struct algebraic_curve;

// struct algebraic_curve_def {};

// TMPL1 struct use<algebraic_curve_def,V>
//   :public use<algebraic_set_def,V>
// {
//   typedef curve<V>           Curve;
//   typedef algebraic_curve<typename use<scalar_def,V>::Scalar,V> AlgebraicCurve;
// };
//--------------------------------------------------------------------

template<class C, class V=default_env>
class algebraic_curve 
  :public curve<REF_OF(V)>
  ,public algebraic_set<C,REF_OF(V)>
{
public:
  
  typedef REF_OF(V) Ref;
  typedef typename use<algebraic_set_def,C,V>::Polynomial Polynomial;
  typedef curve<Ref>                                      Curve;

  algebraic_curve(void): Curve() {};
  algebraic_curve(char* s) ;
  algebraic_curve(const char* s1, const char* s2); 
  algebraic_curve(const Polynomial& eq) ;
  algebraic_curve(const Polynomial& eq1, const Polynomial& eq2);

  ~algebraic_curve(void) {};
  
};

TMPL
SELF::algebraic_curve(char * s): Curve()
{
  //typedef typename use<numeric_def,C,V>::Rational Coefficients;
  //polynomial<Coefficients, with<Sparse,DegRevLex> >t(s, variables("x y"));
  Polynomial e(s, variables("x y") ); //let::assign(s,variables("x y"));
  this->insert(e);
}

TMPL
SELF::algebraic_curve(const char * s1, const char* s2) : Curve()
{
  variables var("x y z");
  //polynomial<double,with<Sparse,DegRevLex> > t1(s1, var);
  //Polynomial e1; let::assign(e1,t1);
  Polynomial e1(s1,var);        
  this->insert(e1);  
        
  //polynomial<double,with<Sparse,DegRevLex> > t2(s2, var);
  //Polynomial e2; let::assign(e2,t2);  
  Polynomial e2(s2,var);
  this->insert(e2);
}

TMPL
SELF::algebraic_curve(const Polynomial&  s) : Curve()
{
  this->insert(s);
}

TMPL
SELF::algebraic_curve(const Polynomial&  p1, const Polynomial&  p2) : Curve()
{
  this->insert(p1);
  this->insert(p2);
}

//====================================================================
} ; // namespace shape
//====================================================================

  template<class OSTREAM, class C, class V>
OSTREAM& operator<<(OSTREAM& os, const shape::SELF& cv) {
  os<<"AlgebraicCurve(";
  for(int i=0; i<cv.nbequation();i++){
    if(i>0) os<<",";
    print(os,cv.equation(i),variables("x y z"));
  }
  os<<")";
  return os;
}

//====================================================================
} ; // namespace mmx
//====================================================================
# undef TMPL
# undef TMPL1
# undef SELF
# undef Viewer
# endif // shape_algebraic_curve_hpp

---