1 files changed, 213 insertions, 0 deletions
diff --git a/2.3-1/src/c/matrixOperations/multiplication/testFloatMatrixMultiplication.c b/2.3-1/src/c/matrixOperations/multiplication/testFloatMatrixMultiplication.c
new file mode 100644
index 00000000..f8ba0cde
--- /dev/null
+++ b/2.3-1/src/c/matrixOperations/multiplication/testFloatMatrixMultiplication.c
@@ -0,0 +1,213 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008-2008 - INRIA - Bruno JOFRET
+ *
+ *  This file must be used under the terms of the CeCILL.
+ *  This source file is licensed as described in the file COPYING, which
+ *  you should have received as part of this distribution.  The terms
+ *  are also available at
+ *  http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+ *
+ */
+
+#include <stdlib.h>
+#include <stdio.h>
+#include <assert.h>
+#include <math.h>
+#include "matrixMultiplication.h"
+
+int testFloatMultiplication(void);
+
+static void cmulmaTest(void) {
+  float realM1[4] = {1.0f, 2.0f, 3.0f, 4.0f};
+  float imagM1[4] = {1.0f, 2.0f, 3.0f, 4.0f};
+  float realM3[6] = {1.0f, 2.0f, 3.0f, 4.0f, 5.0f, 6.0f};
+  float imagM3[6] = {1.0f, 2.0f, 3.0f, 4.0f, 5.0f, 6.0f};
+
+  floatComplex *M1;
+  floatComplex *M2;
+  floatComplex M1_mul_M2[4];
+  floatComplex *M3;
+  floatComplex *M4;
+  floatComplex M3_mul_M4[4];
+  floatComplex miscM3_mul_M4[9];
+
+  int i = 0;
+
+  printf("\n>>>> Matrix Complex Float Multiplication Tests\n");
+
+  M1 = FloatComplexMatrix(realM1, imagM1, 4);
+  M2 = FloatComplexMatrix(realM1, imagM1, 4);
+
+  /*
+    [ 1+1.%i 3+3.%i] * [ 1+1.%i 3+3.%i] = [ 14.%i 30.%i ]
+    [ 2+2.%i 3+3.%i]   [ 2+2.%i 3+3.%i]   [ 20.%i 44.%i ]
+  */
+  cmulma(M1, 2, 2, M2, 2, 2, M1_mul_M2);
+  for (i = 0; i < 4; ++i)
+    {
+      printf("M1_mul_M2[%d] = %e + %e i\n", i, creals(M1_mul_M2[i]), cimags(M1_mul_M2[i]));
+    }
+  for (i = 0; i < 4; ++i)
+    {
+      assert(creals(M1_mul_M2[i]) == 0.0f);
+    }
+  assert(cimags(M1_mul_M2[0]) == 14.0f);
+  assert(cimags(M1_mul_M2[1]) == 20.0f);
+  assert(cimags(M1_mul_M2[2]) == 30.0f);
+  assert(cimags(M1_mul_M2[3]) == 44.0f);
+
+  M3 = FloatComplexMatrix(realM3, imagM3, 6);
+  M4 = FloatComplexMatrix(realM3, imagM3, 6);
+
+  /*
+    [ 1+1.%i 3+3.%i 5+5.%i ] * [ 1+1.%i 4+4.%i ] = [ 44.%i  98.%i ]
+    [ 2+2.%i 4+4.%i 6+6.%i ]   [ 2+2.%i 5+5.%i ]   [ 56.%i 128.%i ]
+			       [ 3+3.%i 6+6.%i ]
+   */
+  cmulma(M3, 2, 3, M4, 3, 2, M3_mul_M4);
+  for (i = 0; i < 4; ++i)
+    {
+      printf("M3_mul_M4[%d] = %e + %e i\n", i, creals(M3_mul_M4[i]), cimags(M3_mul_M4[i]));
+    }
+  for (i = 0; i < 4; ++i)
+    {
+      assert(creals(M3_mul_M4[i]) == 0.0f);
+    }
+  assert(cimags(M3_mul_M4[0]) == 44.0f);
+  assert(cimags(M3_mul_M4[1]) == 56.0f);
+  assert(cimags(M3_mul_M4[2]) == 98.0f);
+  assert(cimags(M3_mul_M4[3]) == 128.0f);
+
+  /*
+    [ 1+1.%i 4+4.%i ] * [ 1+1.%i 3+3.%i 5+5.%i ] = [ 18.%i 38.%i  58.%i ]
+    [ 2+2.%i 5+5.%i ]   [ 2+2.%i 4+4.%i 6+6.%i ]   [ 24.%i 52.%i  80.%i ]
+    [ 3+3.%i 6+6.%i ]				   [ 30.%i 66.%i 102.%i ]
+   */
+  cmulma(M3, 3, 2, M4, 2, 3, miscM3_mul_M4);
+  for (i = 0; i < 9; ++i)
+    {
+      printf("miscM3_mul_M4[%d] = %e + %e i\n", i, creals(miscM3_mul_M4[i]), cimags(miscM3_mul_M4[i]));
+    }
+  for (i = 0; i < 9; ++i)
+    {
+      assert(creals(miscM3_mul_M4[i]) == 0.0);
+    }
+  assert(cimags(miscM3_mul_M4[0]) == 18.0);
+  assert(cimags(miscM3_mul_M4[1]) == 24.0);
+  assert(cimags(miscM3_mul_M4[2]) == 30.0);
+  assert(cimags(miscM3_mul_M4[3]) == 38.0);
+  assert(cimags(miscM3_mul_M4[4]) == 52.0);
+  assert(cimags(miscM3_mul_M4[5]) == 66.0);
+  assert(cimags(miscM3_mul_M4[6]) == 58.0);
+  assert(cimags(miscM3_mul_M4[7]) == 80.0);
+  assert(cimags(miscM3_mul_M4[8]) == 102.0);
+
+}
+
+static void smulmaTest(void) {
+  float M1[4] = {1.0f, 2.0f, 3.0f, 4.0f};
+  float M2[4] = {1.0f, 2.0f, 3.0f, 4.0f};
+  float M1_by_M2[4];
+  float M3[4] = {1.0f, 0.0f, 1.0f, 0.0f};
+  float M4[4] = {0.0f, 1.0f, 0.0f, 1.0f};
+  float M3_by_M4[4];
+  float M5[4] = {1.0f, 0.0f, 0.0f, 1.0f};
+  float M6[4] = {42.0f, 51.0f, 69.0f, 1664.0f};
+  float M5_by_M6[4];
+  float M7[6] = {1.0f, 4.0f, 2.0f, 5.0f, 3.0f, 6.0f};
+  float M8[6] = {1.0f, 3.0f, 5.0f, 2.0f, 4.0f, 6.0f};
+  float M7_by_M8[4];
+  float miscM7_by_M8[9];
+
+  int i = 0;
+
+  printf("\n>>>> Matrix Float Multiplication Tests\n");
+  /*
+    [ 1 3 ]  * [ 1 3 ] = [ 7 15 ]
+    [ 2 4 ]    [ 2 4 ]   [10 22 ]
+  */
+  smulma(M1, 2, 2, M2, 2, 2, M1_by_M2);
+  for (i = 0; i < 4; ++i) {
+    printf("M1_by_M2[%d] = %e\n", i, M1_by_M2[i]);
+  }
+  assert(M1_by_M2[0] == 7.0f);
+  assert(M1_by_M2[1] == 10.0f);
+  assert(M1_by_M2[2] == 15.0f);
+  assert(M1_by_M2[3] == 22.0f);
+  /*
+    [ 1 1 ]  * [ 0 0 ] = [ 1 1 ]
+    [ 0 0 ]    [ 1 1 ]   [ 0 0 ]
+  */
+  smulma(M3, 2, 2, M4, 2, 2, M3_by_M4);
+  for (i = 0; i < 4; ++i) {
+    printf("M3_by_M4[%d] = %e\n", i, M3_by_M4[i]);
+  }
+  assert(M3_by_M4[0] == 1.0f);
+  assert(M3_by_M4[1] == 0.0f);
+  assert(M3_by_M4[2] == 1.0f);
+  assert(M3_by_M4[3] == 0.0f);
+
+  /*
+    [ 1 0 ]  * [ 42   69 ] = [ 42   69 ]
+    [ 0 1 ]    [ 51 1664 ]   [ 51 1664 ]
+  */
+  smulma(M5, 2, 2, M6, 2, 2, M5_by_M6);
+  for (i = 0; i < 4; ++i) {
+    printf("M5_by_M6[%d] = %e\n", i, M5_by_M6[i]);
+  }
+  assert(M5_by_M6[0] == 42.0f);
+  assert(M5_by_M6[1] == 51.0f);
+  assert(M5_by_M6[2] == 69.0f);
+  assert(M5_by_M6[3] == 1664.0f);
+
+  /*
+    [ 1 2 3 ]  * [ 1 2 ] = [ 22 28 ]
+    [ 4 5 6 ]    [ 3 4 ]   [ 49 64 ]
+    [ 5 6 ]
+  */
+  smulma(M7, 2, 3, M8, 3, 2, M7_by_M8);
+  for (i = 0; i < 4; ++i) {
+    printf("M7_by_M8[%d] = %e\n", i, M7_by_M8[i]);
+  }
+  assert(M7_by_M8[0] == 22.0f);
+  assert(M7_by_M8[1] == 49.0f);
+  assert(M7_by_M8[2] == 28.0f);
+  assert(M7_by_M8[3] == 64.0f);
+
+  /*
+    [ 1 5 ]  * [ 1 5 4 ] = [ 16 15 34 ]
+    [ 4 3 ]    [ 3 2 6 ]   [ 13 26 34 ]
+    [ 2 6 ]		   [ 20 22 44 ]
+  */
+  smulma(M7, 3, 2, M8, 2, 3, miscM7_by_M8);
+  for (i = 0; i < 9; ++i) {
+    printf("miscM7_by_M8[%d] = %e\n", i, miscM7_by_M8[i]);
+  }
+  assert(miscM7_by_M8[0] == 16.0f);
+  assert(miscM7_by_M8[1] == 13.0f);
+  assert(miscM7_by_M8[2] == 20.0f);
+  assert(miscM7_by_M8[3] == 15.0f);
+  assert(miscM7_by_M8[4] == 26.0f);
+  assert(miscM7_by_M8[5] == 22.0f);
+  assert(miscM7_by_M8[6] == 34.0f);
+  assert(miscM7_by_M8[7] == 34.0f);
+  assert(miscM7_by_M8[8] == 44.0f);
+
+
+
+}
+
+int testFloatMultiplication(void) {
+
+  printf("\n>>>> Matrix Multiplication Tests\n");
+  smulmaTest();
+  cmulmaTest();
+
+  return 0;
+}
+
+int main(void) {
+	assert(testFloatMultiplication() == 0);
+	return 0;
+}