include/shape/mesher3d_dual.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_dual_hpp
# define shape_mesher3d_dual_hpp

# include <shape/mesher3d.hpp>
# include <shape/dualize.hpp>
# include <shape/tpl3d.hpp>

# define TMPL   template<class C,class V,class Cell>
# define TMPL1  template<class V>
# define Viewer viewer<W>
# define SELF   mesher3d_dual<C,V,Cell>
//====================================================================
namespace mmx {
namespace shape {
//====================================================================
template<class C,class V=default_env,class Cell=bcell<C,V> >
struct mesher3d_dual : public mesher3d<C,V,Cell>
{

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

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

//   void set_input (Cell* cl) { m_input= cl; }

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

//   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 ;
  
// private:
//   Cell*           m_input;
//   Output*         m_output;
 
};
//--------------------------------------------------------------------
TMPL SELF::mesher3d_dual(double e1, double e2): mesher3d<C,V,Cell>(e1,e2)
//m_smooth(e1), m_prec(e2)
{
  //  m_output = new Output;
}

TMPL SELF::~mesher3d_dual(void) {
  //  delete m_output;
}
//--------------------------------------------------------------------
TMPL void SELF::run(void) {

  typedef subdivision<C,V,Cell>      Subdivisor;
  Subdivisor* sbd = new Subdivisor(this->get_smoothness(),this->get_precision());
  sbd->set_input(this->input());
  sbd->run();

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

  typedef dualize<C,V,Cell>  Dualize;
  Dualize* dl = new Dualize;
  dl->set_input(sbd->output());
  dl->run();

  //this->get_adjacency();
  std::cout<< "Dual edges= "<< dl->output()->m_edges.size()/2<<"\n";
  std::cout<< "Dual faces= "<< dl->output()->m_faces.size()/4<<"\n";

  bool ctr = false;
  if (ctr) foreach(Cell * cl, sbd->output()->m_leaves) { 
      marching_cube::centralise(*cl, *cl);
      foreach(Point* p, cl->m_points) this->output()->insert(p); 
      
      //marching_cube::polygonise(*cl, *cl);
  }

  if (ctr) for(unsigned i=0; i< dl->output()->m_edges.size();i+=2){
    Edge* e=new Edge(dl->output()->m_edges[i]->get_cell()->m_points[0],
                     dl->output()->m_edges[i+1]->get_cell()->m_points[0]);
    this->output()->insert(e);
  }

  if (ctr) 
  for(unsigned i=0; i< dl->output()->m_faces.size();i+=4){
    Face* f=new Face;
    f->insert(dl->output()->m_faces[i]->get_cell()->m_points[0]); 
    f->insert(dl->output()->m_faces[i+1]->get_cell()->m_points[0]);
    f->insert(dl->output()->m_faces[i+2]->get_cell()->m_points[0]);
    this->output()->insert(f);
  }
  
  if (!ctr) 
  foreach(Cell * cl, sbd->output()->m_leaves) { 
    marching_cube::polygonise(*cl, *cl);
    foreach(Point* p, cl->m_points)  this->output()->insert(p);
    foreach(Face* f,  cl->m_faces)   this->output()->insert(f);
    std::cout<<"MC"<<std::endl;
  }

}

//--------------------------------------------------------------------
// 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

---