Pos3dList.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/>.
//----------------------------------------------------------------------------
//Pos3dList.

#ifndef LISTAPOS3D_H
#define LISTAPOS3D_H

#include "Pos3d.h"
#include "Vector3d.h"
#include "../d3/GeomObj3d.h"

class Trf3d;


//
class Pos3dList: public GeomObj3d
  {
  public:
    typedef GeomObj::list_Pos3d list_Pos3d;
    typedef GeomObj::list_Pos3d::iterator point_iterator;
    typedef GeomObj::list_Pos3d::const_iterator point_const_iterator;
  protected:
    list_Pos3d point_lst;
    Pos3d &operator[](const size_t &i);
  public:
    Pos3dList(void);
    explicit Pos3dList(const GeomObj::list_Pos3d &l);
    explicit Pos3dList(const boost::python::list &);
    template <typename InputIterator>
    explicit Pos3dList(InputIterator first,InputIterator last);
    virtual GeomObj *getCopy(void) const
      { return new Pos3dList(*this); }
    inline size_t getNumberOfPoints(void) const
      { return point_lst.size(); }
    inline point_const_iterator points_begin(void) const
      { return point_lst.begin(); }
    const point_const_iterator points_end(void) const
      { return point_lst.end(); }
    inline short unsigned int Dimension(void) const
      { return 0; }

    const GeomObj::list_Pos3d &getPoints(void) const
      { return point_lst; }
    const Pos3d &operator[](const size_t &i) const;

    const Pos3d *appendPoint(const Pos3d &p);
    template <class InputIterator>
    inline void assign(InputIterator first, InputIterator last)
      { point_lst.assign(first,last); }
    virtual bool In(const Pos3d &p, const double &tol= 0.0) const;

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

    inline virtual GEOM_FT Pxy(void) const
      { return NAN; }
    virtual GEOM_FT getLength(void) const
      { return 0.0; }
    virtual GEOM_FT getArea(void) const
      { return 0.0; }
    virtual GEOM_FT getVolume(void) const
      { return 0.0; }
    virtual GEOM_FT GetMax(unsigned short int i) const;
    virtual GEOM_FT GetMin(unsigned short int i) const;
    Pos3dList GetMayores(unsigned short int i,const GEOM_FT &d) const;
    Pos3dList GetMenores(unsigned short int i,const GEOM_FT &d) const;
    const Pos3d &Point(const size_t &i) const;
    Pos3d getCenterOfMass(void) const
      { return point_lst.getCenterOfMass(); }
    Vector3d getAvgKVector(void) const;
    bool clockwise(const Pos3d &) const;
    bool counterclockwise(const Pos3d &) const;
    std::string orientation(const Pos3d &) const;

    void Transform(const Trf3d &trf3d);

    std::deque<GEOM_FT> &GetSeparaciones(void) const;
    double GetSeparacionMedia(void) const;

    CGPlane_3 linearLeastSquaresFittingPlane(GEOM_FT &quality);
    void Print(std::ostream &stream) const;
    void Plot(Plotter &) const;
  };

template <typename InputIterator>
Pos3dList::Pos3dList(InputIterator begin,InputIterator end)
  : GeomObj3d()
  {
    for(InputIterator i= begin; i!= end; i++)
      point_lst.push_back(*i);
  }

#endif