From 18218c8f9c86474c76c5a9cc0aea3df47ed18dd6 Mon Sep 17 00:00:00 2001
From: Julian Kent <julian@auterion.com>
Date: Mon, 16 Sep 2019 17:40:25 +0200
Subject: [PATCH] Test non-square matrix multiplication (#91)

* Test non-square matrix multiplication
---
 test/matrixMult.cpp | 25 +++++++++++++++++++++++++
 1 file changed, 25 insertions(+)

diff --git a/test/matrixMult.cpp b/test/matrixMult.cpp
index 1699d3f3c8..fa97f95106 100644
--- a/test/matrixMult.cpp
+++ b/test/matrixMult.cpp
@@ -35,6 +35,31 @@ int main()
     TEST(isEqual(B, B_check));
     Matrix3f C = B_check.edivide(C_check);
     TEST(isEqual(C, C_check));
+
+    // Test non-square matrix
+    float data_43[12] = {1,3,2,
+                         2,2,1,
+                         5,2,1,
+                         2,3,4
+                        };
+    float data_32[6] = {2,3,
+                        1,7,
+                        5,4
+                       };
+
+    Matrix<float, 4, 3> m43(data_43);
+    Matrix<float, 3, 2> m32(data_32);
+
+    Matrix<float, 4, 2> m42 = m43 * m32;
+
+    float data_42[8] = {15,32,
+                        11,24,
+                        17,33,
+                        27,43
+                       };
+    Matrix<float, 4, 2> m42_check(data_42);
+    TEST(isEqual(m42, m42_check))
+
     return 0;
 }