geometry/topology3D.hh

Go to the documentation of this file.

/* Topology.
 * how to subdivide, now in 3D
 */

#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)
      throw(StrategyChange);

    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