Section.h

#ifndef Section2_h__
#define Section2_h__
 
#include <vector>
#include <Eigen/Eigen>
 
namespace Ship{
    struct Section
    {
        EIGEN_MAKE_ALIGNED_OPERATOR_NEW;
    
        Section( const std::vector<double> &  x, const std::vector<double> & y,const std::vector<double> & z);
    
        double GetStationArea( double T );
    
        double GetStationDraught( double T  );
    
        double GetStationMaxBredth(double T);
    
        double GetStationBredth(double T);
    
        double GetStationPosition();
    
        void SectionLewisData( double B, double T, double &LewisA1, double &LewisA3, double &LewisScaleFactor );
    
        int FindSkegPoint();
    
        int m_skegPoint;
        std::vector< Eigen::Vector3d > Vertices;
        double m_minZ;
    
        static bool z_less( const Eigen::Vector3d & lhs, const Eigen::Vector3d & rhs )
        {
            return lhs.z() < rhs.z();
        }
    
        static bool z_grt( const Eigen::Vector3d & lhs, const Eigen::Vector3d & rhs )
        {
            return lhs.z() > rhs.z();
        }
    
        static bool y_less( const Eigen::Vector3d & lhs, const Eigen::Vector3d & rhs )
        {
            return lhs.y() < rhs.y();
        }
    
        static bool y_grt( const Eigen::Vector3d & lhs, const Eigen::Vector3d & rhs )
        {
            return lhs.y() > rhs.y();
        }
    
        static bool x_less( const Eigen::Vector3d & lhs, const Eigen::Vector3d & rhs )
        {
            return lhs.x() < rhs.x();
        }
    
        static bool x_grt( const Eigen::Vector3d & lhs, const Eigen::Vector3d & rhs )
        {
            return lhs.x() > rhs.x();
        }
    };
}
#endif // Section2_h__