From 44a1fb77d9199ce76fec21ee8aaa6df38e141bf0 Mon Sep 17 00:00:00 2001
From: Mikko Rasa <tdb@tdb.fi>
Date: Sun, 2 Jun 2019 16:22:24 +0300
Subject: [PATCH] Copy Matrix test cases for DynamicMatrix

---
 tests/dynamicmatrix.cpp | 67 +++++++++++++++++++++++++++++++++++++++++
 1 file changed, 67 insertions(+)
 create mode 100644 tests/dynamicmatrix.cpp

diff --git a/tests/dynamicmatrix.cpp b/tests/dynamicmatrix.cpp
new file mode 100644
index 0000000..769945d
--- /dev/null
+++ b/tests/dynamicmatrix.cpp
@@ -0,0 +1,67 @@
+#include <limits>
+#include <msp/linal/dynamicmatrix.h>
+#include <msp/test/test.h>
+
+using namespace std;
+using namespace Msp;
+
+class DynamicMatrixTests: public Test::RegisteredTest<DynamicMatrixTests>
+{
+private:
+	typedef LinAl::DynamicMatrix<double> Matrixd;
+
+public:
+	DynamicMatrixTests();
+
+	static const char *get_name() { return "DynamicMatrix"; }
+
+private:
+	void multiply();
+	void invert();
+};
+
+
+DynamicMatrixTests::DynamicMatrixTests()
+{
+	add(&DynamicMatrixTests::multiply, "Multiplication");
+	add(&DynamicMatrixTests::invert, "Inversion");
+}
+
+void DynamicMatrixTests::multiply()
+{
+	static double data[] =
+	{
+		1, 2, 1, 2, 1, 2,
+		3, 0, 3, 0, 0, 3,
+		3, 6, 3, 3, 6, 3, 6, 3, 6,
+		9, 3, 9, 6
+	};
+
+	Matrixd a(3, 2, data);
+	Matrixd b(2, 3, data+6);
+	EXPECT_EQUAL(a*b, Matrixd(3, 3, data+12));
+	EXPECT_EQUAL(b*a, Matrixd(2, 2, data+21));
+}
+
+template<typename T>
+bool is_identity(const LinAl::DynamicMatrix<T> &m)
+{
+	static const T limit = numeric_limits<T>::epsilon()*4;
+	for(unsigned i=0; i<m.rows(); ++i)
+		for(unsigned j=0; j<m.columns(); ++j)
+			if(abs(T(i==j)-m(i, j))>limit)
+				return false;
+	return true;
+}
+
+void DynamicMatrixTests::invert()
+{
+	static double data[] =
+	{
+		1, 4, 2, 4, 8, 6, 2, 2, 4
+	};
+
+	Matrixd m(3, 3, data);
+	Matrixd i = LinAl::invert(m);
+	EXPECT(is_identity(i*m));
+}
-- 
2.43.0