Line.h

 
#ifndef ACTIVE_GEOMETRY_LINE
#define ACTIVE_GEOMETRY_LINE
 
 
#include <memory>
 
#include "Active/Geometry/XPoint.h"
#include "Active/Utility/MathFunctions.h"
 
namespace active::geometry {
    
    class Matrix3x3;
    class XList;
 
    class Line {
    public:
        
        // MARK: - Types
 
        using Unique = std::unique_ptr<Line>;
        using Shared = std::shared_ptr<Line>;
        using Option = std::optional<Line>;
        
        // MARK: - Constructors
        
        Line(const Point& origin, const Point& end) : origin(origin), end(end) {}
        Line(const Point& origin, double rad, double azim) : origin(origin), end(origin) { end.movePolar(rad, azim); }
        
        // MARK: - Variables
        
        Point origin;
        Point end;
 
        // MARK: - Operators
        
        bool operator== (const Line& ref) const { return isEqual3D(ref); }
        bool operator!= (const Line& ref) const { return !isEqual3D(ref); }
        Line operator+ (const Point& offset) const;
        Line& operator+= (const Point& offset);
        Line operator- (const Point& offset) const;
        Line& operator-= (const Point& offset);
        Line operator* (double scale) const;
        Line& operator*= (double scale);
        Line operator* (const Matrix3x3& matrix) const;
        Line& operator*= (const Matrix3x3& matrix);
        
        // MARK: - Functions (const)
        
        bool isEqual2D(const Line& ref, double prec = math::eps) const;
        bool isEqual3D(const Line& ref, double prec = math::eps) const;
        double length2D() const;
        double length3D() const;
        double azimuthAngle() const;
        double altitudeAngle() const;
        double angleTo(const Line& ref) const;
        Point midpoint() const;
        double heightAt(const Point& ref) const;
        bool isColinearTo2D(const Line& ref, double prec = math::eps) const;
        bool isColinearTo3D(const Line& ref, double prec = math::eps) const;
        bool isParallelTo2D(const Line& ref, double prec = math::eps) const;
        bool isParallelTo3D(const Line& ref, double prec = math::eps) const;
        Point closestPointTo2D(const Point& ref, double prec = math::eps) const;
        Point closestPointTo3D(const Point& ref, double prec = math::eps) const;
        Point closestPointAlong2D(const Point& ref, double prec = math::eps) const;
        Point closestPointAlong3D(const Point& ref, double prec = math::eps) const;
        XPoint::Option intersectionWith2D(const Line& ref , double prec = math::eps) const;
        vertex_index intersectionWith2D(const Line& ref, XList& inter, double prec = math::eps) const;
        vertex_index intersectionWith3D(const Line& ref, XList& inter, double prec = math::eps) const;
        Position positionOf2D(const Point& ref, double prec = math::eps) const;
        Position positionOf3D(const Point& ref, double prec = math::eps) const;
        bool encloses2D(const Point& ref, double prec) const;
        bool encloses3D(const Point& ref, double prec) const;
        double angleTo2D(const Line& ref) const;
        
        // MARK: - Functions (mutating)
        
        void extend(double len);
        void movePolar(double len, double azim, double alt);
        void movePolar(double len, double angle);
        void flip();
    };
    
}
 
#endif  //ACTIVE_GEOMETRY_LINE