cube.hh

 
#ifndef cube_hh
#define cube_hh
 
#include "basics/typedefs.hh"
#include "geometry/mesh.hh"
 
using concepts::Real;
 
// ******************************************************************** Cube **
 
class Cube : public concepts::Mesh3 {
public:
  static const uint ncell_ = 1;
  Cube(const Real sizex = 1.0, const Real sizey = 1.0, const Real sizez = 1.0);
  virtual ~Cube();
 
  inline uint ncell() const { return ncell_; }
  inline concepts::Scan3* scan() { return new S(cell_); }
 
  virtual std::ostream& info(std::ostream& os) const;
private:
  class S : public concepts::Scan<concepts::Cell3> {
    uint idx_;
    concepts::Hexahedron3d *(&cell_)[ncell_];
  public:
    inline S(concepts::Hexahedron3d *(&cell)[ncell_]) : idx_(0), cell_(cell) {}
    inline S(const S &scan) : idx_(scan.idx_), cell_(scan.cell_) {}
 
    inline bool eos() const { return idx_ == ncell_; }
    inline concepts::Cell3& operator++(int) { return *cell_[idx_++]; }
 
    inline concepts::Scan3* clone() const { return new S(*this); }
  };
 
  concepts::Vertex* vtx_[8];
  concepts::Edge* edg_[12];
  concepts::Quad* quad_[6];
  concepts::Hexahedron* hex_[ncell_];
 
  concepts::Hexahedron3d* cell_[ncell_];
};
 
// ******************************************************************* Cube2 **
 
class Cube2 : public concepts::Mesh3 {
public:
  Cube2(const uint variant = 0);
  virtual ~Cube2();
 
  inline uint ncell() const { return 2; }
  inline concepts::Scan3* scan() { return new S(cell_); }
 
  virtual std::ostream& info(std::ostream& os) const;
private:
  class S : public concepts::Scan<concepts::Cell3> {
    uint idx_;
    concepts::Hexahedron3d *(&cell_)[2];
  public:
    inline S(concepts::Hexahedron3d *(&cell)[2]) : idx_(0), cell_(cell) {}
    inline S(const S &scan) : idx_(scan.idx_), cell_(scan.cell_) {}
 
    inline bool eos() const { return idx_ == 2; }
    inline concepts::Cell3& operator++(int) { return *cell_[idx_++]; }
 
    inline concepts::Scan3* clone() const { return new S(*this); }
  };
 
  concepts::Vertex* vtx_[12];
  concepts::Edge* edg_[20];
  concepts::Quad* quad_[11];
  concepts::Hexahedron* hex_[2];
 
  concepts::Hexahedron3d* cell_[2];
};
 
// ******************************************************************* Cube3 **
 
class Cube3 : public concepts::Mesh3 {
public:
  Cube3();
  virtual ~Cube3();
 
  inline uint ncell() const { return 2; }
  inline concepts::Scan3* scan() { return new S(cell_); }
 
  virtual std::ostream& info(std::ostream& os) const;
private:
  class S : public concepts::Scan<concepts::Cell3> {
    uint idx_;
    concepts::Hexahedron3d *(&cell_)[2];
  public:
    inline S(concepts::Hexahedron3d *(&cell)[2]) : idx_(0), cell_(cell) {}
    inline S(const S &scan) : idx_(scan.idx_), cell_(scan.cell_) {}
 
    inline bool eos() const { return idx_ == 2; }
    inline concepts::Cell3& operator++(int) { return *cell_[idx_++]; }
 
    inline concepts::Scan3* clone() const { return new S(*this); }
  };
 
  concepts::Vertex* vtx_[12];
  concepts::Edge* edg_[20];
  concepts::Quad* quad_[11];
  concepts::Hexahedron* hex_[2];
 
  concepts::Hexahedron3d* cell_[2];
};
 
#endif // cube_hh