include/shape/cell3d.hpp

Go to the documentation of this file.

/*****************************************************************************
 * M a t h e m a  g i x
 *****************************************************************************
 * Cell
 * 2008-03-20
 * Julien Wintz & Bernard Mourrain & Angelos Mantzaflaris
 *****************************************************************************
 *               Copyright (C) 2008 INRIA Sophia-Antipolis
 *****************************************************************************
 * Comments :
 ****************************************************************************/
# ifndef shape_cell3d_hpp
# define shape_cell3d_hpp

# include <shape/vertex.hpp>
# include <shape/edge.hpp>
# include <shape/face.hpp>
# include <shape/bounding_box.hpp>
# include <shape/cell.hpp>
# include <shape/graph.hpp>
# include <shape/topology.hpp>
# include <shape/mesh3d.hpp>

# include <shape/solver_implicit.hpp>
# define TMPL template<class C, class V>
# define TMPL1 template<class K>
# define SELF cell3d<C,V>
//====================================================================
namespace mmx { namespace shape {

// TMPL struct cell3d;
// TMPL struct tpl3d;

TMPL class cell3d;
TMPL class tpl3d;

//--------------------------------------------------------------------
struct cell3d_def {};
template<> struct use<cell3d_def>
        :public use<cell_def>
{
    typedef cell3d<double,double> Cell3d;
};

//--------------------------------------------------------------------
template<class C, class V=default_env>
class cell3d : public  bounding_box<C,V>
{
public:
    typedef topology<C,V>                 Topology;
    typedef tpl3d<C,V>                    Topology3d;
    typedef typename use<mesh3d_def,C,V>::Point      Point;
    typedef typename Topology::Edge       Edge;

    typedef cell3d<C,V>                   CellBase;
    typedef cell<C,V>                     Cell;
    typedef bounding_box<C,V>             BoundingBox;

    typedef solver_implicit<C,V> Solver;

    cell3d(void) {};
    cell3d(const BoundingBox& bx): BoundingBox(bx) {} ;
    virtual ~cell3d(void) {};

    BoundingBox boundingBox() const { return (BoundingBox)*this; }
    virtual bool        is_regular(void) = 0 ;
    virtual bool        is_active (void) const = 0 ;

    virtual int         subdivide(SELF*& left, SELF*& right);
    virtual void        subdivide(SELF*& left, SELF*& right, int v, double s)=0;
    virtual void        split_position(int& v, double& t);

    virtual bool is_adjacent(SELF* c);

    Point       center() { return Point((this->xmin()+this->xmax())/2, (this->ymin()+this->ymax())/2,
                                        (this->zmin()+this->zmax())/2); }
};

//====================================================================
TMPL int
SELF::subdivide(SELF*& Left, SELF*& Right) {
    int v; double s;
    this->split_position(v,s);
    this->subdivide(Left,Right,v,s);
    return v;
}
//====================================================================
TMPL void 
SELF::split_position(int& v, double& s) {
    double sx = (this->xmax()-this->xmin());
    double sy = (this->ymax()-this->ymin());
    double sz = (this->zmax()-this->zmin());

    if(sx<sy)
        if(sy<sz) {
            v=2;
            s=(this->zmax()+this->zmin())/2;
        } else {
            v=1;
            s=(this->ymax()+this->ymin())/2;
        }
    else
        if(sx<sz) {
            v=2;
            s=(this->zmax()+this->zmin())/2;
        } else {
            v=0;
            s=(this->xmax()+this->xmin())/2;
        }
}


TMPL
bool SELF::is_adjacent(SELF* c2){
  if(this->xmax()<c2->xmin() || c2->xmax()<this->xmin())
    return false;
  if(this->ymax()<c2->ymin() || c2->ymax()<this->ymin())
    return false;
  if(this->zmax()<c2->zmin() || c2->zmax()<this->zmin())
    return false;
  if((this->xmax()==c2->xmin() || c2->xmax()==this->xmin())) {
    if((this->ymax()==c2->ymin() || c2->ymax()==this->ymin()) ||
       (this->zmax()==c2->zmin() || c2->zmax()==this->zmin()) )
      return false;
  } else if((this->ymax()==c2->ymin() || c2->ymax()==this->ymin()) &&
            (this->zmax()==c2->zmin() || c2->zmax()==this->zmin()) )
    return false;
  return true;
}

//--------------------------------------------------------------------
template<class CELL> bool 
is_adjacentpl3d(CELL* c1,CELL* c2){
    if(c1->xmax()<c2->xmin() || c2->xmax()<c1->xmin())
        return false;
    if(c1->ymax()<c2->ymin() || c2->ymax()<c1->ymin())
        return false;
    if(c1->zmax()<c2->zmin() || c2->zmax()<c1->zmin())
        return false;
    if((c1->xmax()==c2->xmin() || c2->xmax()==c1->xmin())) {
        if((c1->ymax()==c2->ymin() || c2->ymax()==c1->ymin()) ||
                (c1->zmax()==c2->zmin() || c2->zmax()==c1->zmin()) )
            return false;
    } else if((c1->ymax()==c2->ymin() || c2->ymax()==c1->ymin()) &&
              (c1->zmax()==c2->zmin() || c2->zmax()==c1->zmin()) )
        return false;
    return true;
}
//--------------------------------------------------------------------
#define xMIN (c2->xmin()<c3->xmin()?c3->xmin():c2->xmin())
#define xMAX (c2->xmax()<c3->xmax()?c2->xmax():c3->xmax())
#define yMIN (c2->ymin()<c3->ymin()?c3->ymin():c2->ymin())
#define yMAX (c2->ymax()<c3->ymax()?c2->ymax():c3->ymax())
#define zMIN (c2->zmin()<c3->zmin()?c3->zmin():c2->zmin())
#define zMAX (c2->zmax()<c3->zmax()?c2->zmax():c3->zmax())

template<class CELL> bool 
is_adjacentpl3d(CELL* c1,CELL* c2, CELL* c3){

    if(c1->xmax()< xMIN || xMAX <c1->xmin())
        return false;
    if(c1->ymax()<yMIN || yMAX<c1->ymin())
        return false;
    if(c1->zmax()<zMIN || zMAX<c1->zmin())
        return false;

    if((c1->xmax()==xMIN || xMAX==c1->xmin())) {
        if((c1->ymax()==yMIN || yMAX==c1->ymin()) ||
                (c1->zmax()==zMIN || zMAX==c1->zmin()) )
            return false;
    } else if((c1->ymax()==yMIN || yMAX==c1->ymin()) &&
              (c1->zmax()==zMIN || zMAX==c1->zmin()) )
        return false;
    return true;
}


#undef xMIN
#undef xMAX
#undef yMIN
#undef yMAX
#undef zMIN
#undef zMAX
//====================================================================
} ; // namespace shape
              } ; // namespace mmx
# undef TMPL
# undef TMPL1
# undef SELF
# endif // shape_cell3d_hpp

---