/Users/mourrain/Devel/mmx/shape/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 algebraic_set<C,Ref>::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,V>::Rational Coefficients;
  polynomial<Coefficients, with<Sparse,DegRevLex> >t(s, variables("x y"));
  Polynomial e; let::assign(e,t);
  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);
}

//--------------------------------------------------------------------
TMPL struct viewer;
TMPL Viewer&
operator<<(Viewer& os, const SELF& c) {
  typedef SELF  AlgebraicCurve;
  typedef typename  AlgebraicCurve::Polynomial Polynomial;
  os<<"<curve type=\"implicit\" color=\""<<(int)(255*os.color.r)<<" "<<(int)(255*os.color.g)<<" "<<(int)(255*os.color.b)<<"\">\n";
  os<<"  <domain>"
    <<os(0,0)<<" "<<os(0,1)<<" "
    <<os(1,0)<<" "<<os(1,1)<<" "
    <<os(2,0)<<" "<<os(2,1)
    <<"</domain>\n";
  for(unsigned i=0; i<c.equations().size();i++)
    {
      typename use<numeric_def,V>::Integer m=1;
      Polynomial p=c.equation(i);
      for(typename Polynomial::const_iterator it=p.begin();it!=p.end();it++){
        m = lcm(denominator(it->coeff()), m);
      }
      p*=(typename Polynomial::Scalar)m;
      os<<(char*)"  <polynomial>";
      print(os,p,variables("x y z"));
      os<<(char*)"</polynomial>\n";
    }
  os<<"</curve>\n";
  return os;
}


//====================================================================
} ; // 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