Reorganize and separate operations

1 files changed, 0 insertions, 206 deletions

diff --git a/src/matrixOperations/testFloatMatrixMultiplication.c b/src/matrixOperations/testFloatMatrixMultiplication.c
deleted file mode 100644
index 6091fc66..00000000
--- a/src/matrixOperations/testFloatMatrixMultiplication.c
+++ /dev/null
@@ -1,206 +0,0 @@
-/*
- *  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;
-}