Vector3.h

 
#ifndef ACTIVE_GEOMETRY_VECTOR_3
#define ACTIVE_GEOMETRY_VECTOR_3
 
#include "Active/Geometry/Point.h"
 
#include <array>
 
namespace active::geometry {
    
    class Line;
    class Matrix3x3;
    class Matrix4x4;
 
    class Vector3 {
    public:
        
        // MARK: - Types
        
        using Unique = std::unique_ptr<Vector3>;
        using Shared = std::shared_ptr<Vector3>;
        using Option = std::optional<Vector3>;
        
        // MARK: - Constructors
        
        Vector3();
        Vector3(double x, double y, double z = 0.0);
        Vector3(const Point& source);
        Vector3(const Line& source);
        Vector3(const Vector3& source);
        virtual ~Vector3() = default;
        
        // MARK: - Operators
        
        Vector3& operator= (const Point& source);
        Vector3& operator= (const Line& source);
        Vector3& operator= (const Vector3& source);
        Vector3 operator+ (const Vector3& vect) const;
        Vector3& operator+= (const Vector3& vect);
        Vector3 operator- (const Vector3& vect) const;
        Vector3& operator-= (const Vector3& vect);
        Vector3 operator* (double scale) const;
        Vector3 operator* (const Matrix3x3& matrix) const;
        Vector3 operator* (const Matrix4x4& matrix) const;
        Vector3& operator*= (double scale);
        Vector3& operator*= (const Matrix3x3& matrix);
        Vector3& operator*= (const Matrix4x4& matrix);
        double& operator[] (int row) { return m_vector[row % 3]; }
        const double operator[] (int row) const { return m_vector[row % 3]; }
        operator Point() const;
        
        // MARK: - Functions (const)
 
        double azimuthAngle() const;
        double altitudeAngle() const;
        bool isParallelTo(const Vector3& ref, double prec = math::eps) const;
        bool isPerpendicularTo(const Vector3& ref, double prec = math::eps) const;
        bool isSameSense(const Vector3& ref, double prec = math::eps) const;
        bool isEmpty(double prec = math::eps) const;
        bool isZAxis(double prec = math::eps) const
                { return (math::isZero(m_vector[0], prec) && math::isZero(m_vector[1], prec) && !math::isZero(m_vector[2], prec)); }
        bool isXyPlane(double prec = math::eps) const
                { return ((!math::isZero(m_vector[0], prec) || !math::isZero(m_vector[1], prec)) && math::isZero(m_vector[2], prec)); }
        double dotProduct(const Vector3& ref) const;
        Vector3 vectorProduct(const Vector3& ref) const;
        double magnitude() const;
        double modulus() const { return magnitude(); }
        Vector3 normalised() const;
        double angleTo(const Vector3& ref) const;
        
    private:
        std::array<double, 3> m_vector;
    };
 
}
 
#endif //ACTIVE_GEOMETRY_VECTOR_3