include/shape/algebraic_surface_axl.hpp

Go to the documentation of this file.

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

# include <shape/algebraic_surface.hpp>
# include <shape/bcell3d_algebraic_surface.hpp>
# include <shape/mesher3d.hpp>

# define TMPL template<class C, class V> 
# define Viewer viewer<axel,V>
# define AlgebraicSurface algebraic_surface<C,V>
# define Cell3dAlgebraicSurface bcell3d_algebraic_surface<C,V>

//====================================================================
namespace mmx { 
  namespace shape {
//====================================================================
template<class C,class V> 
struct use<graphic_def,AlgebraicSurface,V>  {
  typedef mesher3d<double,default_env,Cell3dAlgebraicSurface> Mesher;
  typedef typename Mesher::Output Output;
};
//--------------------------------------------------------------------
template<class C, class V, class BBOX>
typename use<graphic_def,AlgebraicSurface,V>::Output*
as_graphic(const AlgebraicSurface& s, const BBOX& bx, 
           double e_sm=0.05, double e_pr = 0.025) {

  
  typedef typename use<graphic_def,AlgebraicSurface,V>::Mesher Mesher;
  
  Cell3dAlgebraicSurface*  c= new Cell3dAlgebraicSurface(s,bx);
  
  Mesher* msh = new Mesher;
  msh->set_smoothness(e_sm);
  msh->set_precision(e_pr);
  msh->set_input(c);
  msh->run();

  return msh->get_output();

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

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

}

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

---