include/shape/subdivision_with_tpl3d.hpp

Go to the documentation of this file.

/*****************************************************************************
 * M a t h e m a g i x
 *****************************************************************************
 * TopologyCurve
 * 2008-05-16
 * Julien Wintz & Bernard Mourrain
 *****************************************************************************
 *               Copyright (C) 2006 INRIA Sophia-Antipolis
 *****************************************************************************
 * Comments :
 ****************************************************************************/
# ifndef shape_subdivision_with_tpl3d_hpp
# define shape_subdivision_with_tpl3d_hpp

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

# define TMPL  template<class C,class V,class CELL>
# define TMPL1 template<class K>
# define SELF  subdivision<C, with_tpl3d<V>,CELL>

//====================================================================
namespace mmx {
namespace shape {

TMPL1 struct with_tpl3d {};
TMPL1 DECLARE_REF_OF(with_tpl3d<K>,REF_OF(K));

template<class FF, class VV> struct use<FF,with_tpl3d<VV> >: public use<FF,VV> {};

TMPL1 struct use<subdivision_def, with_tpl3d<K> >
  : public use<subdivision_def, K>
{

  template<class Self, class CELL>  static bool 
  process_singular(Self* tpl3d, CELL* cl) {
#if 0
    BoundingBox* bx=tpl3d;
    Point 
      *p0= new Point(bx->xmin(),bx->ymin(),bx->zmin()),
      *p1= new Point(bx->xmin(),bx->ymax(),bx->zmin()),
      *p2= new Point(bx->xmax(),bx->ymax(),bx->zmin()),
      *p3= new Point(bx->xmax(),bx->ymin(),bx->zmin());
    tpl3d->insert(p0);tpl3d->insert(p1); tpl3d->insert(new Edge(p0,p1));
    tpl3d->insert(p1);tpl3d->insert(p2); tpl3d->insert(new Edge(p1,p2));
    tpl3d->insert(p2);tpl3d->insert(p3); tpl3d->insert(new Edge(p2,p3));
    tpl3d->insert(p3);tpl3d->insert(p0); tpl3d->insert(new Edge(p3,p0));
    
    Point 
      *q0= new Point(bx->xmin(),bx->ymin(),bx->zmax()),
      *q1= new Point(bx->xmin(),bx->ymax(),bx->zmax()),
      *q2= new Point(bx->xmax(),bx->ymax(),bx->zmax()),
      *q3= new Point(bx->xmax(),bx->ymin(),bx->zmax());
    tpl3d->insert(q0);tpl3d->insert(q1); tpl3d->insert(new Edge(q0,q1));
    tpl3d->insert(q1);tpl3d->insert(q2); tpl3d->insert(new Edge(q1,q2));
    tpl3d->insert(q2);tpl3d->insert(q3); tpl3d->insert(new Edge(q2,q3));
    tpl3d->insert(q3);tpl3d->insert(q0); tpl3d->insert(new Edge(q3,q0));
    
    tpl3d->insert(p0);tpl3d->insert(q0);tpl3d->insert(new Edge(p0,q0));
    tpl3d->insert(p1);tpl3d->insert(q1);tpl3d->insert(new Edge(p1,q1));
    tpl3d->insert(p2);tpl3d->insert(q2);tpl3d->insert(new Edge(p2,q2));
    tpl3d->insert(p3);tpl3d->insert(q3);tpl3d->insert(new Edge(p3,q3));
#endif
    return true; 
  }

  //--------------------------------------------------------------------
  template<class Self, class CELL>  static bool 
  process_regular(Self* tpl3d, CELL* cl) {

    //    std::cout<<"process regular"<<std::endl;

    // put in the leaves and will be treated at the end
    tpl3d->m_leaves <<cl;

    marching_cube::polygonise(*tpl3d,*cl);
    
    //  if(this->m_idx == 0) 
    //  marching_cube::polygonise(*tpl3d, *this); //NOT before topology_regular
    //    tpl3d->insert(F);
  
#if 0
    process_singular(tpl3d,cl);
#endif

    //foreach(Point* p, F->points()) tpl3d->insert(p);
    return true; 
  }
};
//--------------------------------------------------------------------
TMPL1 struct adapt_tpl3d {};
template<class FF, class VV> struct use<FF,adapt_tpl3d<VV> >: public use<FF,with_tpl3d<VV> > {

  template<class Self, class CELL>  static bool 
  process_regular(Self* tpl3d, CELL* cl) {
    typedef typename Self::Scalar C;
    typedef subdivision<C,with_tpl3d<VV>,typename Self::Cell> Subdivision;
    return use<subdivision_def,with_tpl3d<VV> >::process_regular(dynamic_cast<Subdivision*>(tpl3d),cl);
  }

  template<class Self, class CELL>  static bool 
  process_singular(Self* tpl3d, CELL* cl) {
    typedef typename Self::Scalar C;
    typedef subdivision<C,with_tpl3d<VV>,typename Self::Cell> Subdivision;
    return use<subdivision_def,with_tpl3d<VV> >::process_singular(dynamic_cast<Subdivision*>(tpl3d),cl);
  }

};
//--------------------------------------------------------------------  
template<class C, class V,class CELL>
struct subdivision<C,with_tpl3d<V>,CELL>
  : public subdivision<C,adapt_tpl3d<V>,CELL>
  , public tpl3d<C,V> {
public:

  typedef subdivision<C,adapt_tpl3d<V>,CELL> Subdivision;
  typedef tpl3d<C,V>                    Topology3d;
  typedef topology<C,V>                 Topology;

  typedef typename Topology::Point      Point;  
  typedef typename Topology::Edge       Edge;
  typedef typename Topology::Face       Face;

  typedef CELL                          Cell;
  typedef typename Subdivision::Node    Node;
  typedef typename Cell::BoundingBox    BoundingBox;

  subdivision(double e1= 0.1, double e2=0.01);

  ~subdivision(void) ;
  
  void clear(void);
  //  void run(void);

  void set_precision (double e){ this->Subdivision::set_precision(e); }
  void set_smoothness(double e){ this->Subdivision::set_smoothness(e);}

//   virtual bool process_regular(Cell* cl)  { 
//     return use<subdivision_def,with_tpl3d<V> >::process_regular(this,cl);
//   }

//   virtual bool process_singular(Cell* cl) { 
//     return use<subdivision_def,with_tpl3d<V> >::process_singular(this,cl);
//   }
};
//--------------------------------------------------------------------
TMPL SELF::subdivision(double e1, double e2): Subdivision(e1, e2) { 

}

//--------------------------------------------------------------------
TMPL SELF::~subdivision(void) {

}

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

---