/Users/mourrain/Devel/mmx/shape/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 Shape, class Cell> 
# define TMPL1  template<class V>
# define Viewer viewer<W>
# define SELF   mesher3d_dual<C,V,Shape,Cell> 
//====================================================================
namespace mmx {
namespace shape {
//====================================================================
template<class C,class V=default_env, class Shape=typename SHAPE_OF(V), class Cell=cell<C,V> >
struct mesher3d_dual : public mesher3d<C,V,Shape,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,Shape,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,Shape,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,Shape,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()/3<<"\n";

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

  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);
  }

  for(unsigned i=0; i< dl->output()->m_faces.size();i+=3){
    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);
  }

//  foreach(Cell * cl, this->m_leaves) { 
//     marching_cube::polygonise(*cl, *cl);
//     foreach(Point* p, cl->m_points)  m_output.insert(p);
//     foreach(Face* f,  cl->m_faces)   m_output.insert(f);
//  }

}

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