circle.hh

Go to the documentation of this file.

 
#ifndef circle_hh
#define circle_hh
 
#include "toolbox/array.hh"
#include "toolbox/sequence.hh"
#include "cell1D.hh"
#include "cell2D.hh"
#include "mesh.hh"
 
namespace concepts {
 
  // **************************************************************** Circle **
 
  class Circle : public Mesh2withBoundary {
  public:
    Circle(const Real r = 1.0, const Real ratio = 0.5,
           const Real innerRadius = 0.0);
    Circle(const Array<Real> rings, const Real ratio = 0.5,
           const Real innerRadius = 0.0);
    Circle(const Array<Real> rings,
           const Array<uint> ringAttrib,
           const Real ratio = 0.5, const Real innerRadius = 0.0);
    Circle(const Array<Real> rings,
           const Array<uint> ringAttrib, const Array<uint> ringEdgeAttrib,
           const Real ratio = 0.5, const Real innerRadius = 0.0);
    Circle(const Real innerRadius, const Array<Real> rings);
    Circle(const Real innerRadius, const Array<Real> rings,
           const Array<uint> ringAttrib);
 
    Circle(const Real2d center, const Real r, const uint n,
           const Attribute attrib = 0, const Real phi0 = 0,
           const Attribute edgAttr = 0);
 
    Circle(const Real2d center, const Real innerRadius,
           const Real outerRadius, const uint n,
           const Attribute attrib = 0, const Real phi0 = 0);
 
    virtual ~Circle();
 
    inline uint ncell() const { return ncell_; }
    inline Scan2* scan() { return new S(&cell_); }
 
    virtual std::ostream& info(std::ostream& os) const;
  private:
    class S : public Scan<Cell2> {
      uint idx_;
      Array<Quad2d*>* cell_;
    public:
      inline S(Array<Quad2d*>* cell) : idx_(0), 
                 cell_(cell) {}
      inline S(const S &scan) : idx_(scan.idx_),
        cell_(scan.cell_) {}
 
      inline bool eos() const { return idx_ == cell_->size(); }
      inline Cell2& operator++(int) { return *(*cell_)[idx_++]; }
 
      inline Scan2* clone() const { return new S(*this); }
    };
    uint ncell_;
 
    Array<Vertex*> vtx_;
    Array<Edge*> edg_;
    Array<Quad*> quad_;
 
    Array<CircleMappingEdge2d*> edges_;
 
    Array<Quad2d*> cell_;
 
    Real2d point_(const uint n, const Real r) const;
    void construct_(const Array<Real> rings, const Real ratio,
                    const Real innerRadius,
                    const std::map<uint, uint>* const attrib = 0,
                    const std::map<uint, uint>* const edgeAttrib = 0,
                    const bool hole = false);
  };
 
 
} // namespace concepts
 
#endif // circle_hh