test_matrix.cpp

#include <iostream>
#include "Matrix2.h"

int main() {
    // Test default constructor
    Matrix2<double> mat1;
    std::cout << "Default Matrix: " << mat1.GetElement(0, 0) << std::endl;

    // Test constructor with dimensions
    Matrix2<double> mat2(3, 3);
    mat2.SetElement(1, 1, 42);
    std::cout << "Matrix (1,1): " << mat2.GetElement(1, 1) << std::endl;

    // Test constructor with input data
    double inputData[4] = {1, 2, 3, 4};
    Matrix2<double> mat3(2, 2, inputData);
    std::cout << "Matrix (0,1): " << mat3.GetElement(0, 1) << std::endl;

    // Test copy constructor
    Matrix2<double> mat4(mat3);
    std::cout << "Copied Matrix (0,1): " << mat4.GetElement(0, 1) << std::endl;

    std::cout << mat3 << std::endl;

    std::cout << mat4 << std::endl;

    mat4.SwapRows(0, 1);

    std::cout << mat4 << std::endl;

    mat4.JoinMatrix(mat3);

    std::cout << mat4 << std::endl;

    mat4.SeparateMatrix(2);
    std::cout << mat4 << std::endl;

    Matrix2<double> I(3, 3);
    I.GetIdentityMatrix();

    std::cout << I << std::endl;

    mat3.invertMatrix();

    double input[9] = {1,2,3,3,2,1,2,1,3};
    Matrix2<double> mat5(3, 3, input);

    // mat5.invertMatrix();

    // std::cout << mat5 << std::endl;

    std::pair<Matrix2<double>, Matrix2<double>> res = mat5.LUdecomposition();

    // Test equality operator
    std::cout << res.first << std::endl;
    std::cout << res.second << std::endl;

    std::cout << mat5.determinant() << std::endl;
    if (mat3 == mat4) {
        std::cout << "Matrix 3 and Matrix 4 are equal." << std::endl;
    } else {
        std::cout << "Matrix 3 and Matrix 4 are not equal." << std::endl;
    }

    return 0;
}