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

Go to the documentation of this file.

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

# ifndef TOPOLOGY_H
# define TOPOLOGY_H
#include <realroot/Seq.hpp>
#include <shape/point_with_idx.hpp>
#include <shape/face.hpp>
#include <shape/edge.hpp>
#include <shape/cell.hpp>
#include <shape/bounding_box.hpp>

# define TMPL template<class C, class V>
# define SHAPE geometric<V>
# define SELF  topology<C,V>

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

TMPL struct topology;

//--------------------------------------------------------------------
struct topology_def {};

template<>struct use<topology_def> 
  :public use<vertex_def> 
{
  typedef use<scalar_def>::Scalar  Scalar;
  typedef use<shape_def>::Shape    Shape;

  typedef with_idx<double>      W;
  typedef vertex<Scalar>       Point;
  typedef edge<Scalar,default_env,Point>        Edge;
  typedef face<Scalar,default_env,Point>        Face;

  typedef  use<cell_def>::BoundingBox BoundingBox;
  typedef  use<cell_def>::Cell   Cell;  

  typedef topology<Scalar,double> Topology;
};


template<class C, class V=default_env>
class topology : public SHAPE_OF(V) {
 
public:
  typedef typename SHAPE_OF(V)     Shape;  
  typedef with_idx<V>              W;
  typedef vertex<C,V>              Point;
  typedef edge<C,V, Point>         Edge;
  typedef face<C,V, Point>         Face;

  typedef bounding_box<C,REF_OF(V)> BoundingBox;
  typedef cell<C,REF_OF(V)>         Cell;  

  topology(){};
  topology(const BoundingBox& bx);
  virtual ~topology(void) ;
  
  //virtual void run(void) = 0 ;
  virtual void insert(Point*) = 0 ;
  virtual void insert(Edge *) = 0 ;
  virtual void insert(BoundingBox*,bool);
  virtual void insert_singular(Point*)    = 0 ;

    //virtual Point* vertices(int) = 0 ;
  
  void set_precision (double);
  void set_smoothness(double);

  void push_back(Shape * s)  { m_objects.push_back(s); }
  void operator<<(Shape * s) { m_objects.push_back(s); }

  unsigned count(void)  { return m_objects.size(); }

public:

  double            m_maxprec ;
  double            m_minprec ;
  
  BoundingBox       m_bbx;
  Seq<Shape*>       m_objects;

};


TMPL SELF::topology(const typename SELF::BoundingBox& bx): m_maxprec(0.1), m_minprec(0.01), m_bbx(bx)
{
 
}
TMPL SELF::~topology(void) {}

TMPL void SELF::set_smoothness(double eps) {
  m_maxprec = eps;
}

TMPL void SELF::set_precision(double eps) {
  m_minprec = eps;
}

TMPL void SELF::insert(BoundingBox *bx, bool cross=false) {
  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());
  
  this->insert(p0);this->insert(p1); this->insert(new Edge(p0,p1));
  this->insert(p1);this->insert(p2); this->insert(new Edge(p1,p2));
  this->insert(p2);this->insert(p3); this->insert(new Edge(p2,p3));
  this->insert(p3);this->insert(p0); this->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());
  
  this->insert(q0);this->insert(q1); this->insert(new Edge(q0,q1));
  this->insert(q1);this->insert(q2); this->insert(new Edge(q1,q2));
  this->insert(q2);this->insert(q3); this->insert(new Edge(q2,q3));
  this->insert(q3);this->insert(q0); this->insert(new Edge(q3,q0));
  
  this->insert(p0);this->insert(q0);this->insert(new Edge(p0,q0));
  this->insert(p1);this->insert(q1);this->insert(new Edge(p1,q1));
  this->insert(p2);this->insert(q2);this->insert(new Edge(p2,q2));
  this->insert(p3);this->insert(q3);this->insert(new Edge(p3,q3));

  if(cross) {
    Point 
      *r0= new Point(bx->xmin(),bx->ymin(),bx->zmax()),
      *r1= new Point(bx->xmin(),bx->ymax(),bx->zmax()),
      *r2= new Point(bx->xmax(),bx->ymax(),bx->zmax()),
      *r3= new Point(bx->xmax(),bx->ymin(),bx->zmax());

    this->insert(r0);this->insert(r2);this->insert(new Edge(r0,r2));
    this->insert(r1);this->insert(r3);this->insert(new Edge(r1,r3));
  }
}

} ; // namespace shape
} ; // namespace mmx

# undef TMPL
# undef Shape
# undef SELF
# endif