From 2ee4d3043a8babf346dc22541464971e19672c24 Mon Sep 17 00:00:00 2001
From: Relintai <relintai@protonmail.com>
Date: Wed, 27 Dec 2023 00:29:12 +0100
Subject: [PATCH] MLPPTests::test_linear_algebra partly ported.

---
 test/mlpp_tests.cpp | 67 +++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 67 insertions(+)

diff --git a/test/mlpp_tests.cpp b/test/mlpp_tests.cpp
index 8503ab6..9654921 100644
--- a/test/mlpp_tests.cpp
+++ b/test/mlpp_tests.cpp
@@ -121,6 +121,73 @@ void MLPPTests::test_statistics() {
 }
 
 void MLPPTests::test_linear_algebra() {
+	MLPPLinAlg alg;
+
+	const real_t square_arr[] = {
+		1, 1, //
+		-1, 1, //
+		1, -1, //
+		-1, -1, //
+	};
+
+	const real_t square_rot_res_arr[] = {
+		1.41421, 1.11022e-16, //
+		-1.11022e-16, 1.41421, //
+		1.11022e-16, -1.41421, //
+		-1.41421, -1.11022e-16, //
+	};
+
+	Ref<MLPPMatrix> square(memnew(MLPPMatrix(square_arr, 4, 2)));
+	Ref<MLPPMatrix> square_rot(memnew(MLPPMatrix(square_rot_res_arr, 4, 2)));
+
+	is_approx_equals_mat(square->rotaten(Math_PI / 4), square_rot, "square->rotaten(Math_PI / 4)");
+	/*
+	std::vector<std::vector<real_t>> A = {
+		{ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 },
+		{ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 },
+	};
+	std::vector<real_t> a = { 4, 3, 1, 3 };
+	std::vector<real_t> b = { 3, 5, 6, 1 };
+
+	std::vector<std::vector<real_t>> mmtr_res = {
+		{ 2, 4, 6, 8, 10, 12, 14, 16, 18, 20 },
+		{ 4, 8, 12, 16, 20, 24, 28, 32, 36, 40 },
+		{ 6, 12, 18, 24, 30, 36, 42, 48, 54, 60 },
+		{ 8, 16, 24, 32, 40, 48, 56, 64, 72, 80 },
+		{ 10, 20, 30, 40, 50, 60, 70, 80, 90, 100 },
+		{ 12, 24, 36, 48, 60, 72, 84, 96, 108, 120 },
+		{ 14, 28, 42, 56, 70, 84, 98, 112, 126, 140 },
+		{ 16, 32, 48, 64, 80, 96, 112, 128, 144, 160 },
+		{ 18, 36, 54, 72, 90, 108, 126, 144, 162, 180 },
+		{ 20, 40, 60, 80, 100, 120, 140, 160, 180, 200 }
+	};
+
+	is_approx_equals_dmat(dstd_mat_to_mat_old(alg.matmult(alg.transpose(A), A)), dstd_mat_to_mat_old(mmtr_res), "alg.matmult(alg.transpose(A), A)");
+
+	is_approx_equalsd(alg.dot(a, b), 36, "alg.dot(a, b)");
+
+	std::vector<std::vector<real_t>> had_prod_res = {
+		{ 1, 4, 9, 16, 25, 36, 49, 64, 81, 100 },
+		{ 1, 4, 9, 16, 25, 36, 49, 64, 81, 100 }
+	};
+
+	is_approx_equals_dmat(dstd_mat_to_mat_old(alg.hadamard_product(A, A)), dstd_mat_to_mat_old(had_prod_res), "alg.hadamard_product(A, A)");
+
+	std::vector<std::vector<real_t>> id_10_res = {
+		{ 1, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
+		{ 0, 1, 0, 0, 0, 0, 0, 0, 0, 0 },
+		{ 0, 0, 1, 0, 0, 0, 0, 0, 0, 0 },
+		{ 0, 0, 0, 1, 0, 0, 0, 0, 0, 0 },
+		{ 0, 0, 0, 0, 1, 0, 0, 0, 0, 0 },
+		{ 0, 0, 0, 0, 0, 1, 0, 0, 0, 0 },
+		{ 0, 0, 0, 0, 0, 0, 1, 0, 0, 0 },
+		{ 0, 0, 0, 0, 0, 0, 0, 1, 0, 0 },
+		{ 0, 0, 0, 0, 0, 0, 0, 0, 1, 0 },
+		{ 0, 0, 0, 0, 0, 0, 0, 0, 0, 1 },
+	};
+
+	is_approx_equals_dmat(dstd_mat_to_mat_old(alg.identity(10)), dstd_mat_to_mat_old(id_10_res), "alg.identity(10)");
+	*/
 }
 
 void MLPPTests::test_univariate_linear_regression() {