geometry/hexsubdiv.hh

Go to the documentation of this file.

/* Topology.
 * how to subdivide a hex
 */

#ifndef hexsubdiv_hh
#define hexsubdiv_hh

#include <iostream>
#include "strategyChange.hh"
#include "subdivision.hh"

namespace concepts {

  // forward declaration
  class Hexahedron;
  class Z2;
  class Z4;
  class Quad;
  class Vertex;

  // ******************************************************** HexSubdivision **

  class HexSubdivision : public Subdivision {
  public:
    virtual void createChildren(Hexahedron& q) const = 0;
    virtual void removeChildren(Hexahedron& q) const = 0;
  protected:

    void subFace0011(Quad& q, Z2 rho, Z4 tau, Quad*& A, Quad*& B) const;
    void subFace1100(Quad& q, Z2 rho, Z4 tau, Quad*& A, Quad*& B) const;
    void subFace0110(Quad& q, Z2 rho, Z4 tau, Quad*& A, Quad*& B) const;
    void subFace1001(Quad& q, Z2 rho, Z4 tau, Quad*& A, Quad*& B) const;


    void subFace0123(Quad& q, Z2 rho, Z4 tau, Quad*& A, Quad*& B,
         Quad*& C, Quad*& D) const;
    void subFace1230(Quad& q, Z2 rho, Z4 tau, Quad*& A, Quad*& B,
         Quad*& C, Quad*& D) const;
    void subFace2301(Quad& q, Z2 rho, Z4 tau, Quad*& A, Quad*& B,
         Quad*& C, Quad*& D) const;

    bool faceSubdiv4(Quad& q) const throw(StrategyChange);

    Quad* findFace(const Vertex* const v, Quad* const q[4],
       const uint tau, const uint rho, const uint vtxIdx) const;

    void subdivFace0(Hexahedron& h, Quad*& A, Quad*& B,
         Quad*& C, Quad*& D, bool normal) const;
    void subdivFace1(Hexahedron& h, Quad*& A, Quad*& B,
         Quad*& C, Quad*& D, bool normal) const;
    void subdivFace2(Hexahedron& h, Quad*& A, Quad*& B,
         Quad*& C, Quad*& D, bool normal) const;
    void subdivFace3(Hexahedron& h, Quad*& A, Quad*& B,
         Quad*& C, Quad*& D, bool normal) const;
    void subdivFace4(Hexahedron& h, Quad*& A, Quad*& B,
         Quad*& C, Quad*& D, bool normal) const;
    void subdivFace5(Hexahedron& h, Quad*& A, Quad*& B,
         Quad*& C, Quad*& D, bool normal) const;
  };

  // ************************************************************ HexSubdiv8 **

  class HexSubdiv8 : public HexSubdivision {
  public:
    virtual void createChildren(Hexahedron& q) const;
    virtual void removeChildren(Hexahedron& q) const;
    static const HexSubdiv8* instance();
    virtual ~HexSubdiv8();
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    static std::auto_ptr<HexSubdiv8> instance_;
  };

  // *********************************************************** HexSubdiv2x **

  class HexSubdiv2x : public HexSubdivision {
  public:
    virtual void createChildren(Hexahedron& q) const;
    virtual void removeChildren(Hexahedron& q) const;
    static const HexSubdiv2x* instance();
    virtual ~HexSubdiv2x();
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    static std::auto_ptr<HexSubdiv2x> instance_;
  };

  // *********************************************************** HexSubdiv2y **

  class HexSubdiv2y : public HexSubdivision {
  public:
    virtual void createChildren(Hexahedron& q) const;
    virtual void removeChildren(Hexahedron& q) const;
    static const HexSubdiv2y* instance();
    virtual ~HexSubdiv2y();
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    static std::auto_ptr<HexSubdiv2y> instance_;
  };

  // *********************************************************** HexSubdiv2z **

  class HexSubdiv2z : public HexSubdivision {
  public:
    virtual void createChildren(Hexahedron& q) const;
    virtual void removeChildren(Hexahedron& q) const;
    static const HexSubdiv2z* instance();
    virtual ~HexSubdiv2z();
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    static std::auto_ptr<HexSubdiv2z> instance_;
  };

  // *********************************************************** HexSubdiv4x **

  class HexSubdiv4x : public HexSubdivision {
  public:
    virtual void createChildren(Hexahedron& q) const;
    virtual void removeChildren(Hexahedron& q) const;
    static const HexSubdiv4x* instance();
    virtual ~HexSubdiv4x();
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    static std::auto_ptr<HexSubdiv4x> instance_;
  };

  // *********************************************************** HexSubdiv4y **

  class HexSubdiv4y : public HexSubdivision {
  public:
    virtual void createChildren(Hexahedron& q) const;
    virtual void removeChildren(Hexahedron& q) const;
    static const HexSubdiv4y* instance();
    virtual ~HexSubdiv4y();
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    static std::auto_ptr<HexSubdiv4y> instance_;
  };

  // *********************************************************** HexSubdiv4z **

  class HexSubdiv4z : public HexSubdivision {
  public:
    virtual void createChildren(Hexahedron& q) const;
    virtual void removeChildren(Hexahedron& q) const;
    static const HexSubdiv4z* instance();
    virtual ~HexSubdiv4z();
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    static std::auto_ptr<HexSubdiv4z> instance_;
  };

} // namespace concepts

#endif // hexsubdiv_hh