include/shape/polygonizer.hpp

Go to the documentation of this file.

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

# include <shape/list.hpp>
# include <shape/mesh3d.hpp>

# define TMPL   template<class C,class V, class Input, class Output>
# define TMPL1  template<class V>
# define SELF   polygonizer<C,V,Input,Output>
//====================================================================
namespace mmx { namespace shape {
//====================================================================

template<class C, class V=default_env, class Input = Seq< cell<C,V>* >, class Output = mesh3d<C,default_env> >
struct polygonizer: public PROCESS_OF(V)
{

//  typedef Seq<Cell*>                  Input;
//  typedef tpl3d<C,V>                  Output;
  typedef typename use<mesh3d_def,C,V>::Point      Point;
  typedef typename use<mesh3d_def,C,V>::Edge       Edge;
  typedef typename use<mesh3d_def,C,V>::Face       Face;
  
  polygonizer(void);
 ~polygonizer(void) ;

  void    set_input  (Input* i)   { m_input = i; }
  Input*  get_input  (void) const { return m_input; }

  void    set_output (Output* o)  { m_output = o; }
  Output* get_output (void) const { return m_output; }

  void run();
  void run(unsigned n, C m=0, C M=1);

  void clear(void);

private:
  Input*       m_input;
  Output*      m_output;
};
//--------------------------------------------------------------------
TMPL SELF::polygonizer(void) {
  m_output = new Output;
}
//--------------------------------------------------------------------
TMPL SELF::~polygonizer(void) {
  delete m_output;
}
//--------------------------------------------------------------------
TMPL void SELF::run() {

  Point  points[8];
  C      values[8];
 
  Input L = *this->get_input();

  for(unsigned c=0; c< this->get_input()->size();c+=8){

    for(unsigned j=0; j<8; j++) {
      points[j] = L[c+j]->center();
      values[j] = L[c+j]->center_value();
    }

    marching_cube::polygonize(*this->get_output(),points, values);

//    Point* P[8];
//    if (!marching_cube::polygonize(*this->get_output(),points, values))
//    {
//      for (unsigned j=0;j<7;j++) {
//        P[j] = new Point(points[j]);
//          this->get_output()->insert(P[j]);
//      }
//    };
    
  }

  // foreach(Cell * cl, *this->input()) { 
  //     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);
  //   }
}//--------------------------------------------------------------------
TMPL void SELF::run(unsigned n, C m, C M) {

    C s=(M-m)/n;
  Point  points[8];
  C      values[8];

  Input L = *this->get_input();

  for(unsigned c=0; c< this->get_input()->size();c+=8){

    for(unsigned j=0; j<8; j++) {
      points[j] = L[c+j]->center();
      values[j] = L[c+j]->center_value()-m;
    }

    for(unsigned i=0; i<n+1; i++) {
       marching_cube::polygonize(*this->get_output(),points, values);
       for(unsigned j=0; j<8; j++) {
         values[j] += s;
       }
    }
  }

}
//--------------------------------------------------------------------
TMPL void SELF::clear(void) {  
}

//====================================================================
} ; // namespace shape
} ; // namespace mmx
//====================================================================
# undef TMPL
# undef TMPL1
# undef SELF 
# endif

---