added division element by element

1 files changed, 206 insertions, 0 deletions

diff --git a/src/operations/division/testLDivision.c b/src/operations/division/testLDivision.c
new file mode 100644
index 00000000..137de8e1
--- /dev/null
+++ b/src/operations/division/testLDivision.c
@@ -0,0 +1,206 @@
+/*
+ *  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 "division.h"
+/*
+#define SIZE 10000000
+
+#define LOCAL_DEBUG 
+*/
+
+/*
+static void sldivsTest()
+{
+	float M1;
+	float M2;
+	int i=0,j=0;
+
+	srand(20);
+	for (i=0;i<10;i++){
+		for (j=0;j<10;j++){
+			M1=(float)rand();
+			M2=(float)rand();
+			printf("%f \\ %f = %1.7f\n" , M1, M2, sldivs(M1,M2));
+		}
+	}		
+
+}
+
+
+static void dldivsTest()
+{
+	double M1;
+	double M2;
+	int i=0,j=0;
+
+	srand(20);
+	for (i=0;i<10;i++){
+		for (j=0;j<10;j++){
+			M1=(double)rand();
+			M2=(double)rand();
+			printf("%f \\ %f = %1.7f\n" , M1, M2, dldivs(M1,M2));
+		}
+	}		
+
+}
+
+
+static void cldivsTest()
+{
+	floatComplex M1;
+	floatComplex M2;
+	int i=0,j=0;
+
+	srand(20);
+	for (i=0;i<10;i++){
+		for (j=0;j<10;j++){
+			M1=FloatComplex((float)rand(),(float)rand());
+			M2=FloatComplex((float)rand(),(float)rand());
+			printf("(%f + %f *%%i)\\ (%f + %f *%%i) = %1.7f + %1.7f i\n" , creals(M1),cimags(M1),creals(M2),cimags(M2),creals(cldivs(M1,M2)),cimags(cldivs(M1,M2)));
+		}
+	}		
+
+}
+
+
+static void zldivsTest()
+{
+	doubleComplex M1;
+	doubleComplex M2;
+	int i=0,j=0;
+
+	srand(20);
+	for (i=0;i<10;i++){
+		for (j=0;j<10;j++){
+			M1=DoubleComplex((double)rand(),(double)rand());
+			M2=DoubleComplex((double)rand(),(double)rand());
+			printf("(%f + %f *%%i)\\ (%f + %f *%%i) = %1.7f + %1.7f i\n\n" , zreals(M1),zimags(M1),zreals(M2),zimags(M2),zreals(zldivs(M1,M2)),zimags(zldivs(M1,M2)));
+		}
+	}		
+
+}
+*/
+
+static void sldivaTest()
+{
+	float M1[4]={1, 3, 2, 4};
+	float M2[4]={23, 3, 2, 13};
+	float MRes[4];
+	int i=0;
+
+	sldiva(M1,4,M2,4,MRes);
+	
+	for (i=0;i<4;i++){
+		printf("%1.7f\n",MRes[i]);
+	}		
+
+}
+
+
+static void dldivaTest()
+{
+	double M1[4]={1, 3, 2, 4};
+	double M2[4]={23, 3, 2, 13};
+	double MRes[4];
+	int i=0;
+
+	dldiva(M1,4,M2,4,MRes);
+
+	for (i=0;i<4;i++){
+		printf("%1.7f\n",MRes[i]);
+	}		
+
+}
+
+
+static void cldivaTest()
+{	
+	int i=0;
+	floatComplex M1[4];	
+	floatComplex M2[4];
+	floatComplex MRes[4];
+
+	M1[0]=FloatComplex(1,0);
+	M1[1]=FloatComplex(3,0);
+	M1[2]=FloatComplex(2,0);
+	M1[3]=FloatComplex(4,0);
+	M2[0]=FloatComplex(23,0);
+	M2[1]=FloatComplex(3,0);
+	M2[2]=FloatComplex(2,0);
+	M2[3]=FloatComplex(13,0);
+
+	cldiva(M1,4,M2,4,MRes);
+	for (i=0;i<4;i++){
+		printf("%1.7f + %1.7f i\n",creals(MRes[i]),cimags(MRes[i]));
+	}		
+
+}
+
+
+static void zldivaTest()
+{
+	int i=0;
+	doubleComplex M1[4];
+	doubleComplex M2[4];
+	doubleComplex MRes[4];
+
+	M1[0]=DoubleComplex(1,0);
+	M1[1]=DoubleComplex(3,0);
+	M1[2]=DoubleComplex(2,0);
+	M1[3]=DoubleComplex(4,0);
+	M2[0]=DoubleComplex(23,0);
+	M2[1]=DoubleComplex(3,0);
+	M2[2]=DoubleComplex(2,0);
+	M2[3]=DoubleComplex(13,0);
+
+	zldiva(M1,4,M2,4,MRes);
+	for (i=0;i<4;i++){
+		printf("%1.7f + %1.7f i\n",zreals(MRes[i]),zimags(MRes[i]));
+	}		
+		
+
+}
+
+
+static int testLDivision(void) {
+/* Resultats a comparer avec ceux de Scilab
+   Ok chez moi
+
+	printf(">>>SCALAIRES<<<\n");printf("\n");printf("\n");
+  sldivsTest();
+  printf("\n");printf("\n");printf("\n");printf("\n");
+  dldivsTest();	
+  printf("\n");printf("\n");printf("\n");printf("\n");
+  cldivsTest();
+  printf("\n");printf("\n");printf("\n");printf("\n");
+  zldivsTest();
+  printf("\n");printf("\n");printf("\n");printf("\n")};*/
+	printf(">>>ARRAYS<<<\n");printf("\n");printf("\n");
+  sldivaTest();
+  printf("\n");printf("\n");printf("\n");printf("\n");  
+  dldivaTest();	
+  printf("\n");printf("\n");printf("\n");printf("\n");
+  cldivaTest();
+  printf("\n");printf("\n");printf("\n");printf("\n");
+  zldivaTest();
+  return 0;
+}
+
+
+
+int main(void) {
+  assert(testLDivision() == 0);
+  return 0;
+}