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

Go to the documentation of this file.

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

# ifndef shape_cell2d_semialgebraic_curve_hpp
# define shape_cell2d_semialgebraic_curve_hpp

# include <realroot/Seq.hpp>
# include <shape/point.hpp>
# include <shape/solver_implicit.hpp>
# include <shape/semialgebraic_curve.hpp>
# include <shape/cell2d_algebraic_curve.hpp>


# define TMPL template<class C, class V>
# define SELF cell2d_semialgebraic_curve<C,V>
# define Cell                      cell<C,V>
# define Cell2d                    cell2d<C,V>
# define Cell2dAlgebraicCurve      cell2d_algebraic_curve<C,V>
# define SemiAlgebraicCurve        semialgebraic_curve<C,V>
# define Solver                    solver_implicit<C,V>
# undef Cell_t
# undef Polynomial
//====================================================================
namespace mmx {
namespace shape {
//====================================================================
// TMPL struct cell2d_semialgebraic_curve;

// TMPL struct with_cell2d_semialgebraic_curve { 
//   typedef cell2d_semialgebraic_curve<K> Cell2dSemiAlgebraicDomain; 
// };

// TMPL struct cell2d_semialgebraic_curve_def 
//   :public semialgebraic_curve_def<K>
//   ,public with_cell2d_semialgebraic_curve<K> 
// {
//   typedef typename cell2d_algebraic_curve_def<K>::Polynomial     Polynomial;
//   typedef typename cell2d_algebraic_curve_def<K>::Topology2d     Topology2d;
// };


TMPL
struct cell2d_semialgebraic_curve : public Cell2dAlgebraicCurve {
public:
  typedef point<C,V> Point;
  typedef typename cell2d_algebraic_curve<C,V>::Polynomial Polynomial;
  typedef typename Cell2d::BoundingBox  BoundingBox;

  typedef cell2d_semialgebraic_curve<C,V> Cell_t;
  
  cell2d_semialgebraic_curve(const Polynomial&, int, const BoundingBox&);
  cell2d_semialgebraic_curve(SemiAlgebraicCurve*, const BoundingBox&);
  
  cell2d_semialgebraic_curve(const SELF& cl)
    : Cell2dAlgebraicCurve(cl.m_polynomial, cl.boundingBox()) {}
  
  
  Polynomial inequality(void) const { return  this->m_polynomial; }
  
  virtual bool  is_active (void) ;
  //virtual bool  is_regular(void) ;

  bool is_touching(void);

  Point * pair(Point * p, int & sgn) const ;
  virtual void subdivide(Cell*& left, Cell*& right, int v, double s);

//  virtual bool insert_singular (SemiAlg2d*);// = 0;

private:
//  int m_sign;
};

TMPL
SELF::cell2d_semialgebraic_curve(const Polynomial& pol, int s, const BoundingBox& bx)
  : Cell2dAlgebraicCurve(pol,bx) { if (s<0) this->m_polynomial*=-1; }
  
TMPL
SELF::cell2d_semialgebraic_curve(SemiAlgebraicCurve* cv, const BoundingBox & b)
  : Cell2dAlgebraicCurve(cv->inequality(),b) { }


TMPL bool 
SELF::is_active() {

  return ( Cell2dAlgebraicCurve::is_active() );

//
  if( !Cell2dAlgebraicCurve::is_active() )//intersecting the bound curves
    return false;

  if( ! this->is_regular())
    return true;

  if( this->m_singular.size()!=0) 
    return true;


  if (this->is_touching() )//is touching the set
    return true;
  else
    return false;
}


TMPL  bool 
SELF::is_touching(){

  const int * sz = this->m_polynomial.rep().szs();
  const int * st = this->m_polynomial.rep().str();
  if (!((this->m_polynomial[0] >0)              ||
        (this->m_polynomial[(sz[0]-1)*st[0]] >0)||
        (this->m_polynomial[sz[0]*sz[1]-1]>0)   ||
        (this->m_polynomial[(sz[1]-1)*st[1]] >0)))
    return false;      
  else
    return true;
} 


TMPL void
SELF::subdivide(Cell*& Left, Cell*& Right, int v, double s) {
//  std::cout<<"subdivide semialg"<< this<<", intsects: "<< this->nb_intersect() <<std::endl;
  Cell2dAlgebraicCurve::subdivide(Left,Right,v,s) ;
//  std::cout<<"subdivide end"<< this<<", intsects: "<< this->nb_intersect() <<std::endl;

}


//====================================================================
} ; // namespace shape
} ; // namespace mmx
//====================================================================
# undef Cell
# undef Cell2d
//# undef AlgebraicCurve
# undef Cell2dAlgebraicCurve
# undef SemiAlgebraicCurve
# undef Solver
# undef SELF
# undef TMPL
# endif // shape_cell2d_semialgebraic_curve_hpp