points.hpp

#ifndef POINTS_FOR_GEOMETRY_H
#define POINTS_FOR_GEOMETRY_H

#include <memory>
#include "shared_ptr.hpp"
#include "rect_fwd.hpp"
#include <cmath>
#include <vector>

class geometry_const_Iterator;

struct point{
    double x;
    double y;
    double z=0;
    int layer;

    point():x(0.), y(0.), layer(0.){};
    point(double xx, double yy, int l = 0):x(xx), y(yy), layer(l){};
};

struct region{
    double top() const { return t;};
    double bottom() const { return b;};
    double left() const { return l;};
    double right() const { return r;};

    double width() const { return (t - b);};
    double length() const { return (r - l);};

    region(const point&, const point&);

    void include(const point&);

    static region determine_extension(const Rectangle_const_ref&,
                      double rotation = 0);
private:
    double r,t,l,b;

    region();
};

class slice_Position{
public:
    slice_Position(const geometry_const_Iterator&,
           const point& origin,
           double rotation = 0.);

    bool advance();

    bool eog();

    double x() const { return p_.x;};
    double y() const { return p_.y;};

    double degrees() const { return rotation_;};
    double rad() const { return rotation_ * M_PI / 180.;};

    double dx() const;
    double dy() const;
    double dphi() const;

    double width() const;

    void transform(point&) const;

    point upper_point() const;
    point lower_point() const;

    const geometry_const_Iterator& position() const { return *gci;};

private:
    std::unique_ptr<geometry_const_Iterator> gci;

    point p_;
    
    double rotation_;

    double cos_alpha;
    double sin_alpha;

    void update_angles();
    void rotate(double& point);
};

struct point_Positions{

    point_Positions(size_t t): r(t){};

    std::vector<point> r;

    typedef STD_TR1::shared_ptr<point_Positions> ref;

    void translate_to (const slice_Position&);
};


#endif