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

Go to the documentation of this file.

/*****************************************************************************
 * M a t h e m a g i x
 *****************************************************************************
 * 2011-04-26
 * Bernard Mourrain
 *****************************************************************************
 *               Copyright (C) 2011 INRIA Sophia-Antipolis
 *****************************************************************************
 * Comments :
 ****************************************************************************/
# ifndef shape_mesher3d_hpp
# define shape_mesher3d_hpp

# include <shape/subdivision.hpp>
# include <shape/tpl3d.hpp>

# define TMPL   template<class C,class V, class Shape, class Cell> 
# define TMPL1  template<class V>
# define Viewer viewer<W>
# define SELF   mesher3d<C,V,Shape,Cell> 
//====================================================================
namespace mmx {
namespace shape {
//====================================================================
struct mesher3d_def {};

template<> struct use<mesher3d_def, default_env>
{
  template<class Topo, class Cell> static
  void polygonise(Topo* t, Cell* c) {
    c->polygonise(t);
  }
};
//--------------------------------------------------------------------

template<class C,class V=default_env, class Shape=typename SHAPE_OF(V), class Cell=cell3d<C,V> >
struct mesher3d : public PROCESS_OF(V)
{
  typedef subdivision<C,V,Shape,Cell> Subdivisor;

  typedef Cell                        Input;
  typedef tpl3d<C,V>                  Output;
  typedef typename Output::Point      Point;
  typedef typename Output::Edge       Edge;
  typedef typename Output::Face       Face;
  typedef typename Cell::BoundingBox  BoundingBox;

  mesher3d(double e1= 0.1, double e2=0.01);
 ~mesher3d(void) ;

  void set_input (Cell* cl) { m_input= cl; }
  void set_input (Shape* s, const BoundingBox& bx);

  void set_smoothness(double e) { m_smooth = e; }
  void set_precision (double e) { m_prec   = e; }

  double get_smoothness(void) { return m_smooth; }
  double get_precision (void) { return m_prec;   }


  Input*  input  (void)     { return m_input; }
  Output* output (void)     { return m_output; }

  void run(void);
  void clear(void);
 
private:
  double          m_smooth ;
  double          m_prec ;

  Cell*           m_input;
  Output*         m_output;
 
};
//--------------------------------------------------------------------
TMPL SELF::mesher3d(double e1, double e2): m_smooth(e1), m_prec(e2)
{
  m_output = new Output;
}

TMPL SELF::~mesher3d(void) {
  delete m_output;
}
//--------------------------------------------------------------------
TMPL void SELF::set_input(Shape* s, const BoundingBox& bx)
{
  this->set_input(cell3d_factory<C,V>::instance()->create(s,bx));
}
//--------------------------------------------------------------------
TMPL void SELF::run(void) {

  Subdivisor* sbd = new Subdivisor(m_smooth,m_prec);
  sbd->set_input(this->input());
  sbd->run();

  std::cout<< "leaves    = "<< sbd->output()->m_leaves.size()<<"\n";

  foreach(Cell * cl, sbd->output()->m_leaves) { 
    use<mesher3d_def,V>::polygonise(this->output(), cl);
    //    marching_cube::polygonise(*this->output(), *cl);
    //foreach(Point* p, cl->m_points) this->output()->insert(p); 
  }
}

//--------------------------------------------------------------------
TMPL void SELF::clear(void) {
  this->output()->clear();
}
//====================================================================
} ; // namespace shape
} ; // namespace mmx
//====================================================================
# undef TMPL
# undef TMPL1
# undef Shape
# undef Viewer
# undef Graphic
# undef SELF 
# endif