infiniteMeshes.hh

Go to the documentation of this file.

 
#ifndef infMeshes_hh
#define infMeshes_hh
 
#include "basics/typedefs.hh"
#include "mesh.hh"
#include "topology.hh"
#include "cell2D.hh"
 
namespace concepts {
 
  // ************************************************** SquareTwoInfiniteRects **
 
  class SquareTwoInfiniteRects : public Mesh2 {
  private:
    class S : public Scan<Cell2> {
      unsigned int idx_;
      Cell2 *(&cell_)[3];
    public:
      inline S(Cell2 *(&cell)[3]) : idx_(0), cell_(cell) {}
      inline S(const S& scan) : idx_(scan.idx_),  cell_(scan.cell_) {}
      inline bool eos() const { return idx_ == 3; }
      inline Cell2& operator++(int) { return *cell_[idx_++]; }
      inline Scan2* clone() const { return new S(*this); }
    };
  
    Vertex       *vtx_    [4];
    Edge         *edg_    [4];
    InfiniteEdge *infEdg_ [4];
    Quad         *quad_   [1];
    InfiniteQuad *infQuad_[2];
    Cell2        *cell_   [3];
    bool periodic_;
  public:
    SquareTwoInfiniteRects(bool periodic = false);
    virtual ~SquareTwoInfiniteRects();
    inline unsigned int ncell() const { return 3; }
    inline Scan2* scan() { return new S(cell_); }
    virtual std::ostream& info(std::ostream& os) const;
  };
 
  // *************************************************** SquareOneInfiniteRect **
 
  class SquareOneInfiniteRect : public Mesh2 {
  private:
    class S : public Scan<Cell2> {
      unsigned int idx_;
      Cell2 *(&cell_)[2];
    public:
      inline S(Cell2 *(&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 Cell2& operator++(int) { return *cell_[idx_++]; }
      inline Scan2* clone() const { return new S(*this); }
    };
  
    Vertex       *vtx_    [4];
    Edge         *edg_    [4];
    InfiniteEdge *infEdg_ [2];
    Quad         *quad_   [1];
    InfiniteQuad *infQuad_[1];
    Cell2        *cell_   [2];
    bool periodic_;
  public:
    SquareOneInfiniteRect(bool periodic = false);
    virtual ~SquareOneInfiniteRect();
    inline unsigned int ncell() const { return 2; }
    inline Scan2* scan() { return new S(cell_); }
    virtual std::ostream& info(std::ostream& os) const;
  };
 
} // namespace concepts
 
#endif // infMeshes_hh