Matrix3x3.h

 
#ifndef ACTIVE_GEOMETRY_MATRIX_3x3
#define ACTIVE_GEOMETRY_MATRIX_3x3
 
#include <array>
#include <memory>
#include <optional>
 
namespace active::geometry {
 
    class Matrix3x3 {
    public:
        
        // MARK: Types
        
            //Matrix row
        using row_t = std::array<double, 3>;
            //Matrix container
        using base_t = std::array<row_t, 3>;
        using Unique = std::unique_ptr<Matrix3x3>;
        using Shared = std::shared_ptr<Matrix3x3>;
        using Option = std::optional<Matrix3x3>;
        
        // MARK: Factory functions
        
        static Matrix3x3 createXRotate(double angle);
        static Matrix3x3 createYRotate(double angle);
        static Matrix3x3 createZRotate(double angle);
        static Matrix3x3 createScale(double x, double y, double z);
        static Matrix3x3 createTranslate(double x, double y, double z);
        static Matrix3x3 createIdentity();
        
        // MARK: Constructors
        
        Matrix3x3();
        Matrix3x3(  double a1, double b1, double c1,
                    double a2, double b2, double c2,
                    double a3, double b3, double c3);
        
        // MARK: Operators
        
        bool operator==(const Matrix3x3& ref) const;
        bool operator!=(const Matrix3x3& ref) const;
        Matrix3x3 operator*(const Matrix3x3& ref) const;
        Matrix3x3 &operator*=(const Matrix3x3& ref);
        const row_t& operator[](unsigned short index) const { return m_matrix[index]; }
        row_t& operator[](unsigned short index) { return m_matrix[index]; }
        const double& operator()(unsigned short row, unsigned short col) const { return m_matrix[row][col]; }
        double& operator()(unsigned short row, unsigned short col) { return m_matrix[row][col]; }
        
        // MARK: Functions (const)
        
        Matrix3x3 getInverse() const;
        double getDeterminant() const;
 
    private:
        base_t m_matrix;
    };
 
}
 
#endif //ACTIVE_GEOMETRY_MATRIX_3x3