geometry/circle.hh

Go to the documentation of this file.

/* circle made out of five quadrilaterials (4 curved, 1 straight)
 */

#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