cuboid.hh

Go to the documentation of this file.

 
#ifndef Cuboid_hh
#define Cuboid_hh
 
#include "basics/typedefs.hh"
#include "geometry/mesh.hh"
 
namespace concepts {
 
// ******************************************************************** Cuboid **
 
  class Cuboid: public concepts::Mesh3 {
  public:
 
    static const uint ncell_ = 1;
    Cuboid(const Real sizex = 1.0, const Real sizey = 1.0, const Real sizez =
        1.0);
 
    virtual ~Cuboid();
 
    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_];
  };
}
#endif // Cuboid_hh