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

Go to the documentation of this file.

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

# include <shape/list.hpp>
# include <shape/subdivision_with_tpl3d.hpp>
# include <shape/tpl3d.hpp>
# include <algorithm>

# define TMPL   template<class C,class V>
# define TMPL1  template<class V>
# define Viewer viewer<V>
# define SELF   mesher3d_shape<C,V> 
//====================================================================
namespace mmx { 
namespace shape {
//====================================================================

template<class C,class V=default_env>
struct mesher3d_shape 
  : public  subdivision<C,with_tpl3d<V>, typename SHAPE_OF(V), cell3d<C,V> >  {

  typedef tpl3d<C,V>                     Topology;
  typedef subdivision<C,with_tpl3d<V>, typename SHAPE_OF(V), cell3d<C,V> >  Subdivision;
  typedef use<subdivision_def,V>              VRT;

  typedef typename SHAPE_OF(V)                Shape;
  typedef cell3d_factory<C,V>                 Cell3dFactory;

  typedef typename Topology::Point            Point;  
  typedef typename Topology::Edge             Edge;
  typedef typename Topology::Face             Face;
  typedef REF_OF(V)                           Ref;
  typedef cell3d<C,Ref>                       Cell;
  typedef cell3d<C,Ref>                       Cell3d;
  typedef surface<Ref>                        Surface;
  typedef typename Cell::BoundingBox          BoundingBox;

  typedef node<Shape *, Cell *>               Node;
  typedef kdtree<Shape *, Cell *>             Kdtree;

  
  mesher3d_shape(double e1= 0.1, double e2=0.01);
  mesher3d_shape(Shape* s, const BoundingBox& bx, double e1= 0.1, double e2=0.01);
 
 ~mesher3d_shape(void) ;

  void set_input(Cell* cl);
  void set_input(Shape* s, const BoundingBox& bx);

  bool is_regular  (Cell * cell);
  bool insert_regular (Cell * cell);
  bool sing_process(Cell * cell);   
  bool subdivide   (Cell * cell, Node * node) ;

  void run(void);
  void clear(void);
};
//--------------------------------------------------------------------
TMPL SELF::mesher3d_shape(double e1, double e2): Subdivision(e1, e2) {
}
//--------------------------------------------------------------------
TMPL SELF::mesher3d_shape(Shape* s0, const BoundingBox& bx0, double e1, double e2)
: Subdivision(e1, e2) {
  this->set_input(Cell3dFactory::instance()->create(s0,bx0));
}
//--------------------------------------------------------------------
TMPL SELF::~mesher3d_shape(void) {
  //  delete m_tree ;
}
//--------------------------------------------------------------------
TMPL void SELF::set_input(Cell* cl)
{
  this->set_input(cl);
}
//--------------------------------------------------------------------
TMPL void SELF::set_input(Shape* s0, const BoundingBox& bx0)
{
  this->set_input(Cell3dFactory::instance()->create(s0,bx0));
}
//--------------------------------------------------------------------
TMPL bool 
SELF::is_regular(Cell* cl) {
  return cl->is_regular();  
}
//--------------------------------------------------------------------
TMPL bool
SELF::insert_regular(Cell* cl) {
  this->output()->m_leaves <<cl;// put in the leaves and will be treated at the end
  return cl->insert_regular(this);
//  return true;
}
//--------------------------------------------------------------------
TMPL bool
SELF::sing_process(Cell* cl) {
  return cl->insert_singular(this) ;
}
//--------------------------------------------------------------------
TMPL bool
SELF::subdivide(Cell* cl, Node * node) {
  Cell * left, * right;
  int v = cl->subdivide(left, right) ;
  node->m_left = new Node(node, left,  Node::LEFT , v) ; this->output()->m_nodes << node->m_left ;
  node->m_right= new Node(node, right, Node::RIGHT, v) ; this->output()->m_nodes << node->m_right;
  return true ;
}
//--------------------------------------------------------------------
TMPL void SELF::run(void) {
  //this->Subdivision::run();
  double maxsz = this->get_smoothness()*this->output()->root()->get_cell()->size();
  double minsz = this->get_precision()*this->output()->root()->get_cell()->size();

  std::cout<<maxsz<<" "<<minsz<<std::endl;
  while(!this->output()->m_nodes.empty()) {
    Node* node = this->output()->m_nodes.front() ;
    Cell* cl = node->get_cell() ;

    if(cl->is_active()) {
      if(cl->size() > maxsz)
        {
          this->subdivide(cl, node) ;
        }
      else if(cl->is_regular())
        {
          this->insert_regular(cl) ;
        }
      else if(cl->size() > minsz )
        {
//        std::cout<<"subdivide not regular"<<std::endl;
          this->subdivide(cl, node) ;
        }
      else {
//      std::cout<<"insert singular"<<std::endl;
        //this->sing_process(cl);
      }
    }
    this->output()->m_nodes.pop_front() ;
  }

  foreach(Cell * cl, this->output()->m_leaves) {
     dynamic_cast<Cell3d*>(cl)->topology_regular(this);
     //foreach(Point* p, cl->m_points)  this->insert(p);
     //foreach(Face* f,  cl->m_faces)   this->insert(f);
  }

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

}
//--------------------------------------------------------------------
TMPL void SELF::clear(void) {
  //  this->output()->clear();
}


//--------------------------------------------------------------------
TMPL Viewer&
operator<<(Viewer& out, const SELF& tp) {
  use<tpl3d_def,V>::print_as_graphic(out,tp);
  return out;
}
//--------------------------------------------------------------------
TMPL typename use<tpl3d_def,V>::Graphic*
as_graphic(const SELF& tp) {
  return use<tpl3d_def,V>::as_graphic(tp);
}
//====================================================================
} ; // namespace shape
} ; // namespace mmx
//====================================================================
# undef TMPL
# undef TMPL1
# undef Shape
# undef Viewer
# undef Graphic
# undef SELF 
# endif