include/shape/algebraic_curve_axl.hpp

Go to the documentation of this file.

/*****************************************************************************
 * M a t h e m a g i x
 *****************************************************************************
 * AlgebraicCurve  with axel
 * 2012-01-28
 * Bernard Mourrain
 *****************************************************************************
 *               Copyright (C) 2012 INRIA Sophia-Antipolis
 *****************************************************************************
 * Comments :
 ****************************************************************************/
# ifndef algebraic_curve_axl_hpp
# define algebraic_curve_axl_hpp

# include <shape/algebraic_curve.hpp>
# include <shape/mesher2d.hpp>

# define TMPL template<class C, class V>
# define AlgebraicCurve algebraic_curve<C,V>
# define Viewer viewer<axel,V>
//====================================================================
namespace mmx {
  namespace shape {
//====================================================================
template<class C, class V> 
struct use<graphic_def,AlgebraicCurve,V> {
  typedef mesher2d<C,V> Output;
};

//--------------------------------------------------------------------
//template<class MESH, class C, class V, class BBOX>
//void
//as_graphic(MESH* mesh, AlgebraicCurve* cv, const BBOX& box, 

//template<class C, class V, class BBOX>
//typename use<graphic_def,AlgebraicCurve,V>::Output*
//as_graphic(AlgebraicCurve* cv, const BBOX& box,
//         double e_sm=0.05, double e_pr = 0.05) {

//  typedef mesher2d<double> Mesher2d;

//  Mesher2d* mesher= new Mesher2d(box);
//  mesher->set_smoothness(e_sm);
//  mesher->set_precision (e_pr);
//  mesher->push_back(cv);
//  mesher->run();
  
//  return mesher;

//  // foreach(Mesher2d::Point* p, mesher->vertices()) {
//  //   Point* pt=new Point(p->x(),p->y(),0.);
//  //   mesh->push_back_vertex(pt);
//  // }

//  // foreach(Mesher2d::Edge* e, mesher->edges()) {
//  //   mesh->push_back_edge(e->source()->index(), e->destination()->index());
//  // }

//}

//--------------------------------------------------------------------
TMPL struct viewer;
TMPL Viewer&
operator<<(Viewer& os, const AlgebraicCurve& c) {
  typedef typename  AlgebraicCurve::Polynomial Polynomial;
  os<<"<curve type=\"algebraic\" 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;
}

//====================================================================
}// end shape namespace
}// end mmx namespace
//====================================================================
# undef TMPL
# undef AlgebraicCurve
# undef Viewer
#endif //shape_algebraic_curve

---