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

Go to the documentation of this file.

/*****************************************************************************
 * M a t h e m a  g i x
 *****************************************************************************
 * ParametricCell
 * 2008-03-20
 * Julien Wintz
 *****************************************************************************
 *               Copyright (C) 2008 INRIA Sophia-Antipolis
 *****************************************************************************
 * Comments :
 ****************************************************************************/

# ifndef PARAMETRICCELL_H
# define PARAMETRICCELL_H

# include <realroot/Seq.hpp>
# include <shape/cell2d.hpp>

# define TMPL template<class C, class V>
# define REF REF_OF(V)
# define Cell2d cell2d<C,REF>
# define ParametricCurve parametric_curve<C,V>
# define SELF cell2d_parametric_curve<C,V>

namespace mmx {
namespace shape {


// TMPL struct cell2d_parametric_curve;

// TMPL struct with_cell2d_parametric_curve { 
//   typedef cell2d_parametric_curve<K> Cell2dParametricCurve; 
// };


// TMPL struct cell2d_parametric_curve_def 
//   :public cell2d_def<K>
//   ,public with_parametric_curve<K>
//   ,public with_cell2d_parametric_curve<K> {
// };

TMPL
class cell2d_parametric_curve : public Cell2d
{
public:
  typedef bounding_box<C,REF>      BoundingBox;
  //typedef cell<K>              Cell;
  typedef cell<C,REF>       Cell;
  typedef typename topology<C,V>::Point       Point; 
  typedef topology<C,V>   Topology;

  //   Cell2dParametricCurve(void) ;
  //   ~Cell2dParametricCurve(void) ;
  
  cell2d_parametric_curve(ParametricCurve *, const BoundingBox &);
    

  bool  is_active (void) { return true; }
  bool  is_regular(void) ;
  bool  is_intersected(void)   { return true; }
  
  virtual Point     * pair(Point * p, int& sgn)   { return NULL; };
  virtual Point     * starting_point(int sgn)     { return NULL; };

  virtual bool        insert_regular (Topology*) { return true; } ;
  virtual bool        insert_singular(Topology*) { return true; } ;


  virtual void subdivide     (Cell*& left, Cell*& right, int v, double s);
  
};


  TMPL
  SELF::cell2d_parametric_curve(ParametricCurve * curve, const BoundingBox & bx) 
  {
    
  }

  TMPL bool 
  SELF::is_regular(void)
  {
    return false ;
  }

  TMPL void
  SELF::subdivide(Cell*& left, Cell*& right, int v, double s) {
    //   double xc = (bx->xmin()+bx->xmax())/2 ;
    //     double yc = (bx->ymin()+bx->ymax())/2 ;
    
    //    SELF * nw = new SELF();//bx->xmin(), xc, yc, bx->ymax()) ;
    //    SELF * ne = new SELF();//xc, bx->xmax(), yc, bx->ymax()) ;
    //    SELF * sw = new SELF();//bx->xmin(), xc, bx->ymin(), yc) ;
    //    SELF * se = new SELF();//xc, bx->xmax(), bx->ymin(), yc) ;
    
    //     foreach(point * p, n_intersections) {
    //  if(p->x() <= xc) nw->northIntersectionpoints() << p ;
    //  if(p->x() >  xc) ne->northIntersectionpoints() << p ;
    //     }
    
    //     foreach(point * p, s_intersections) {
    //  if(p->x() <= xc) sw->southIntersectionpoints() << p ;
    //  if(p->x() >  xc) se->southIntersectionpoints() << p ;
    //     }
    
    //     foreach(point * p, w_intersections) {
    //  if(p->y() <= yc) sw->westIntersectionpoints() << p ;
    //  if(p->y() >  yc) nw->westIntersectionpoints() << p ;
    //     }
    
    //     foreach(point * p, e_intersections) {
    //  if(p->y() <= yc) se->eastIntersectionpoints() << p ;
    //  if(p->y() >  yc) ne->eastIntersectionpoints() << p ;
    //     }
    
    //     foreach(point * p, m_critical) {
    //  if(p->x() <= xc && p->y() <= yc) sw->criticalpoints() << p ;
    //  if(p->x() <= xc && p->y() >  yc) nw->criticalpoints() << p ;
    //  if(p->x() >  xc && p->y() <= yc) se->criticalpoints() << p ;
    //  if(p->x() >  xc && p->y() >  yc) ne->criticalpoints() << p ;
    //     }
    
    //     foreach(point * p, m_extremal) {
    //  if(p->x() <= xc && p->y() <= yc) sw->extermalpoints() << p ;
    //  if(p->x() <= xc && p->y() >  yc) nw->extermalpoints() << p ;
    //  if(p->x() >  xc && p->y() <= yc) se->extermalpoints() << p ;
    //  if(p->x() >  xc && p->y() >  yc) ne->extermalpoints() << p ;
    //     }
    
    //     foreach(point * p, m_singular) {
    //  if(p->x() <= xc && p->y() <= yc) sw->singularpoints() << p ;
    //  if(p->x() <= xc && p->y() >  yc) nw->singularpoints() << p ;
    //  if(p->x() >  xc && p->y() <= yc) se->singularpoints() << p ;
    //  if(p->x() >  xc && p->y() >  yc) ne->singularpoints() << p ;
    //     }
    

    //    cells << (Cell *)nw ; cells << (Cell *)ne ;
    //    cells << (Cell *)sw ; cells << (Cell *)se ;
  }
  
} ; // namespace shape
} ; // namespace mmx
# undef TMPL
# undef SELF
# undef Cell2d
# undef ParametricCurve
# undef Cell2dList
# undef Cell2dFactory
# undef Topology2d 
# endif // PARAMETRICCELL_H