Polyline2d.h

// -*-c++-*-
//----------------------------------------------------------------------------
//  xc utils library; general purpose classes and functions.
//
//  Copyright (C)  Luis C. Pérez Tato
//
//  XC utils is free software: you can redistribute it and/or modify
//  it under the terms of the GNU General Public License as published by
//  the Free Software Foundation, either version 3 of the License, or 
//  (at your option) any later version.
//
//  This software is distributed in the hope that it will be useful, but 
//  WITHOUT ANY WARRANTY; without even the implied warranty of
//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU General Public License for more details.  
//
// You should have received a copy of the GNU General Public License 
// along with this program.
// If not, see <http://www.gnu.org/licenses/>.
//----------------------------------------------------------------------------
//Polyline2d.h

#ifndef POLYLINE2D_H
#define POLYLINE2D_H

#include "../pos_vec/Pos2d.h"
#include "Linear2d.h"

class Segment2d;
class Ray2d;
class Line2d;
class HalfPlane2d;
class Pos2dList;
class Polygon2d;
class Circle2d;
class VectorPos2d;

//
class Polyline2d : public Linear2d, public GeomObj::list_Pos2d
  {
  public:
    Polyline2d(void);
    explicit Polyline2d(const GeomObj::list_Pos2d &);
    explicit Polyline2d(const Pos2dList &);
    explicit Polyline2d(const boost::python::list &);
    explicit Polyline2d(const Polygon2d &);

    virtual GeomObj *getCopy(void) const
      { return new Polyline2d(*this); }
    const GeomObj::list_Pos2d &getVertices(void) const;
    const GeomObj::list_Pos2d &getVertexList(void) const;
    boost::python::list getVertexListPy(void) const;
    inline size_t getNumVertices(void) const
      { return GeomObj::list_Pos2d::size(); }
    size_t getNumSegments(void) const;
    void setVertices(const boost::python::list &);
    inline const_iterator vertices_begin(void) const
      { return GeomObj::list_Pos2d::begin(); }
    const const_iterator vertices_end(void) const
      { return GeomObj::list_Pos2d::end(); }

    void Transform(const Trf2d &trf2d);
    
/*     inline virtual void GiraX(double ang) */
/*       { Ref2d::GiraX(ang); TrfPoints(); } */
/*     inline virtual void GiraY(double ang) */
/*       { Ref2d::GiraY(ang); TrfPoints();  } */
/*     inline virtual void GiraZ(double ang) */
/*       { Ref2d::GiraZ(ang); TrfPoints();  } */
    const Pos2d *appendVertex(const Pos2d &);
    const Pos2d *appendVertexLeft(const Pos2d &);
    template <class InputIterator>
    inline void assign(InputIterator first, InputIterator last)
      { GeomObj::list_Pos2d::assign(first,last); }
    virtual bool In(const Pos2d &p, const double &tol= 0.0) const;
    virtual GEOM_FT dist2(const Pos2d &p) const;
    GEOM_FT dist(const Pos2d &p) const;

    virtual GEOM_FT Ix(void) const;
    virtual GEOM_FT Iy(void) const;
    virtual GEOM_FT Iz(void) const;

    GEOM_FT getCurvatureAtVertex(const_iterator) const;
    GEOM_FT getCurvatureAtLength(const GEOM_FT &) const;
    Vector2d getCurvatureVectorAtVertex(const_iterator) const;
    std::vector<Vector2d> getCurvatureVectorAtVertices(void) const;
    Vector2d getCurvatureVectorAtLength(const GEOM_FT &) const;
    Pos2d getPointAtLength(const GEOM_FT &) const;
    Vector2d getIVectorAtLength(const GEOM_FT &) const;
    Vector2d getJVectorAtLength(const GEOM_FT &) const;
    
    VectorPos2d Divide(int num_partes) const;
    boost::python::list DividePy(int num_partes) const;
    VectorPos2d Divide(const std::vector<double> &) const;
    boost::python::list DividePy(const boost::python::list &) const;
    
    inline virtual GEOM_FT Pxy(void) const
      { return NAN; }
    virtual GEOM_FT getLength(void) const
      { return GeomObj::list_Pos2d::getLength(); }
    virtual GEOM_FT getLengthUpTo(const Pos2d &) const;
    inline GEOM_FT getLambda(const Pos2d &p) const
      { return getLengthUpTo(p); }
    virtual GEOM_FT GetMax(unsigned short int i) const;
    //Return the maximum value of the i coordinate.
    virtual GEOM_FT GetMin(unsigned short int i) const;
    Polyline2d GetMayores(unsigned short int i,const GEOM_FT &d) const;
    Polyline2d GetMenores(unsigned short int i,const GEOM_FT &d) const;
    const Pos2d &Vertice(const size_t &) const;
    const Pos2d &Vertice0(const size_t &) const;
    Segment2d getSegment(const const_iterator &i) const;
    Segment2d getSegment(const size_t &i) const;
    Segment2d getSegment0(const size_t &i) const;
    std::vector<Segment2d> getSegments(void) const;
    const_iterator getNearestSegmentIter(const Pos2d &) const;    
    int getNearestSegmentIndex(const Pos2d &) const;
    Segment2d getNearestSegment(const Pos2d &) const;
    Pos2d Projection(const Pos2d &) const;
    Polyline2d offset(const GEOM_FT &d) const;
    Polygon2d getBufferPolygon(const GEOM_FT &d, const size_t &numVertices= 8) const;
    
    std::deque<const_iterator> getIntersectionIters(const Line2d &) const;
    list_Pos2d getIntersection(const Line2d &r) const;
    list_Pos2d getIntersection(const Ray2d &sr) const;
    list_Pos2d getIntersection(const Segment2d &sr) const;
    std::deque<Polyline2d> clip(const HalfPlane2d &, const GEOM_FT &tol= 0.0) const;
    boost::python::list clipPy(const HalfPlane2d &, const GEOM_FT &tol= 0.0) const;
    void insertVertex(const Pos2d &, const GEOM_FT &tol= 1e-6);
    Polyline2d getChunk(const Pos2d &,const short int &, const GEOM_FT &tol) const;
    Polyline2d getLeftChunk(const Pos2d &, const GEOM_FT &tol) const;
    Polyline2d getRightChunk(const Pos2d &, const GEOM_FT &tol) const;
    boost::python::list split(const Pos2d &) const;
    Pos2d getCenterOfMass(void) const
      { return GeomObj::list_Pos2d::getCenterOfMass(); }

    iterator getFarthestPointFromSegment(iterator it1, iterator it2, GEOM_FT &pMaxDist);
    
    boost::python::dict getPyDict(void) const;
    void setPyDict(const boost::python::dict &);

    void Print(std::ostream &stream) const;
    void Plot(Plotter &) const;
  };

GeomObj::list_Pos2d intersection(const Polyline2d &p,const Line2d &r);
GeomObj::list_Pos2d intersection(const Polyline2d &p,const Ray2d &);
GeomObj::list_Pos2d intersection(const Polyline2d &p,const Segment2d &);
GeomObj::list_Pos2d intersection(const Line2d &,const Polyline2d &);
GeomObj::list_Pos2d intersection(const Ray2d &,const Polyline2d &p);
GeomObj::list_Pos2d intersection(const Segment2d &,const Polyline2d &p);

#endif