topology3D.hh

Go to the documentation of this file.

 
#ifndef topology3D_hh
#define topology3D_hh
 
#include <memory>
#include "basics/typedefs.hh"
#include "basics/Zm.hh"
#include "connector.hh"
#include "topology.hh"
#include "hexsubdiv.hh"
 
namespace concepts {
 
  // *********************************************************** Tetrahedron **
 
  class Tetrahedron : public Connector3 {
  public:
    Tetrahedron(Triangle& tri0, Triangle& tri1, Triangle& tri2, 
    Triangle& tri3, Attribute attrib = Attribute());
    virtual ~Tetrahedron();
 
    virtual Tetrahedron* child(uint i, bool mode = 0);
    virtual const Tetrahedron* child(uint i) const;
 
    Edge* edge(uint i) const;
    Triangle* face(uint i) const { return i < 4 ? tri_[i] : 0; }
    Vertex* vertex(uint i) const;
 
    int rho(int i) const { return rho_[i]; }
 
    int tau(int i) const { return tau_[i]; }
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    Tetrahedron* chld_;
 
    Tetrahedron* lnk_;
 
    Triangle*    tri_[4];
 
    Z3               tau_[4];
 
    Z2               rho_[4];
  };
 
  // ************************************************************ Hexahedron **
 
  class Hexahedron : public Connector3 {
    friend class HexSubdiv8;
    friend class HexSubdiv2x; 
    friend class HexSubdiv2y;
    friend class HexSubdiv2z;
    friend class HexSubdiv4x;
    friend class HexSubdiv4y;
    friend class HexSubdiv4z;
 public:
    Hexahedron(Quad& quad0, Quad& quad1, Quad& quad2, 
         Quad& quad3, Quad& quad4, Quad& quad5, 
         const Attribute attrib = Attribute());
    virtual ~Hexahedron();
 
    virtual Hexahedron* child(uint i, bool mode = 0);
    virtual const Hexahedron* child(uint i) const;
 
    Edge* edge(uint i) const;
    Quad* face(uint i) const { return i < 6 ? quad_[i] : 0; }
    Vertex* vertex(uint i) const;
 
    Z2 rho(int i) const { return rho_[i]; }
 
    Z4 tau(int i) const { return tau_[i]; }
 
 
 
    void setStrategy(const HexSubdivision* strategy = 0);
 
    const HexSubdivision* getStrategy() const { return subdivStrategy_; }
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    Hexahedron* chld_;
 
    Hexahedron* lnk_;
 
    Quad* quad_[6];
 
    Z4 tau_[6];
 
    Z2 rho_[6];
    const HexSubdivision* subdivStrategy_;
  };
 
} // namespace concepts
 
#endif // topology3D_hh