1 files changed, 95 insertions, 0 deletions
diff --git a/src/operations/multiplication/testMultiplication.c b/src/operations/multiplication/testMultiplication.c
new file mode 100644
index 00000000..38c43438
--- /dev/null
+++ b/src/operations/multiplication/testMultiplication.c
@@ -0,0 +1,95 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2006-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 "multiplication.h"
+
+#define SIZE 10000000
+/*
+#define LOCAL_DEBUG 
+*/
+static void ztimessTest(void) {
+  doubleComplex M1;
+  doubleComplex M2;
+  doubleComplex M1_by_M2;
+
+  int i = 0;
+
+  srand(1);
+  for (i = 0; i < SIZE; ++i) {
+    M1 = DoubleComplex((double) rand(), (double) rand());
+    M2 = DoubleComplex((double) rand(), (double) rand());
+
+    M1_by_M2 = ztimess(M1, M2);
+#ifdef LOCAL_DEBUG
+    /*ATTENTION à modifier (c'est du copier coller de testDivision)
+    printf("erreurR = %e\n", fabs(zreals(M1_by_M2) - zreals(num)/zreals(den) )/fabs(zreals(M1_by_M2)));
+    printf("erreurI = %e\n", fabs(zimags(M1_by_M2) - zimags(num)/zreals(den) )/fabs(zimags(M1_by_M2)));*/
+#endif
+    if (zreals(M1_by_M2)<1e-14 && (zreals(M1)*zreals(M2) - zimags(M1)*zimags(M2)) < 1e-18 ) 
+	assert(1);
+   else 
+	assert( fabs(zreals(M1_by_M2) - (zreals(M1)*zreals(M2) - zimags(M1)*zimags(M2)) )/fabs(zreals(M1_by_M2)) <= 1e-16);
+
+   if (  zimags(M1_by_M2) < 1e-14 && (zreals(M1)*zimags(M2) + zimags(M1)*zreals(M2)) < 1e-18 )
+        assert ( 1 ) ;
+   else
+	assert( fabs(zimags(M1_by_M2) - (zreals(M1)*zimags(M2) + zimags(M1)*zreals(M2)) )/fabs(zimags(M1_by_M2)) <= 1e-16);
+  }
+}
+
+static void ctimessTest(void) {
+  floatComplex M1;
+  floatComplex M2;
+  floatComplex M1_by_M2;
+
+  int i = 0;
+
+  srand(1);
+  for (i = 0; i < SIZE; ++i) {
+    M1 = FloatComplex((float) rand(), (float) rand());
+    M2 = FloatComplex((float) rand(), (float) rand());
+
+    M1_by_M2 = ctimess(M1, M2);
+#ifdef LOCAL_DEBUG
+    /*ATTENTION à modifier (c'est du copier coller de testDivision)
+    printf("erreurR = %e\n", creals(M1_by_M2) - (creals(num)/creals(den))/creals(M1_by_M2));
+    printf("erreurI = %e\n", cimags(M1_by_M2) - (cimags(num)/creals(den))/cimags(M1_by_M2));*/
+#endif
+  if (creals(M1_by_M2)<1e-14 && (creals(M1)*creals(M2) - cimags(M1)*cimags(M2)) < 1e-18 ) 
+	assert(1);
+  else 
+	assert( fabs(creals(M1_by_M2) - (creals(M1)*creals(M2) - cimags(M1)*cimags(M2)) )/fabs(creals(M1_by_M2)) < 1e-7);
+
+  if (  cimags(M1_by_M2) < 1e-14 && (creals(M1)*cimags(M2) + cimags(M1)*creals(M2)) < 1e-18 )
+        assert ( 1 ) ;
+  else
+    assert( fabs(cimags(M1_by_M2) - (creals(M1)*cimags(M2) + cimags(M1)*creals(M2)) )/fabs(cimags(M1_by_M2)) < 1e-7);
+  }
+}
+
+
+static int testMult(void) {
+
+  ztimessTest();
+  ctimessTest();
+  return 0;
+}
+
+
+
+int main(void) {
+  assert(testMult() == 0);
+  return 0;
+}