std::vector<std::vector<int>> data1, data2;

for (int i = 0; i < 3; i++) {

std::vector<int> temp1, temp2;

for (int j = 0; j < 3; j++) {

temp1.emplace_back(i+j);

temp2.emplace_back(i*j);

}

data1.emplace_back(temp1);

data2.emplace_back(temp2);

}

int data3[] = {1, 2, 3};

std::cout << "Test Matrix\n========================================================" << std::endl;

std::shared_ptr<Matrix<int>> m1 = std::make_shared<Matrix<int>>(data1);

std::cout << "Matrix1: " << std::endl << *m1 << std::endl;

std::shared_ptr<Matrix<int>> m2 = std::make_shared<Matrix<int>>(data2);

std::cout << "Matrix2: " << std::endl << *m2 << std::endl;

std::cout << "Sum: " << std::endl << *m1 + *m2 << std::endl;

std::cout << "Scalar Product: " << std::endl << *m1 * 6 << std::endl;

std::cout << "Matrix Product: " << std::endl << *m1 * *m2 << std::endl;

std::shared_ptr<Vector<int>> v1 = std::make_shared<Vector<int>>(3, data3);

std::cout << "Vector Product: " << std::endl << *m1 * *v1 << std::endl;

m1->at(0, 0) = 8;

m1->at(1, 1) = 18;

m1->at(2, 2) = 28;

std::cout << "Update Matrix1: " << std::endl << *m1;

int data1[] = {1, 2, 3};

int data2[] = {8, 9, 15};

std::cout << "Test Vector\n========================================================" << std::endl;

std::shared_ptr<Vector<int>> v1 = std::make_shared<Vector<int>>(3, data1);

std::cout << "Vector1: " << *v1 << std::endl;

std::shared_ptr<Vector<int>> v2 = std::make_shared<Vector<int>>(3, data2);

std::cout << "Vector2: " << *v2 << std::endl;

std::cout << "Sum: " << *v1 + *v2 << std::endl;

std::cout << "Scalar Product: " << *v1 * 6 << std::endl;

if (v1->getSize() < 4)

std::cout << "Cross Product: " << v1->cross(*v2) << std::endl;

std::cout << "Dot Product: " << v1->dot(*v2) << std::endl;

testMatrix();

testVector();

return 0;