From 343dc58409dcc9d52d931f9ea2e4e541cfc0db4b Mon Sep 17 00:00:00 2001
From: torset
Date: Wed, 28 Jan 2009 13:25:36 +0000
Subject: Update interfaces

---
 src/operations/interface/int_OpMinus.h | 103 ++++++++++++++++++---------------
 1 file changed, 57 insertions(+), 46 deletions(-)

(limited to 'src/operations/interface/int_OpMinus.h')

diff --git a/src/operations/interface/int_OpMinus.h b/src/operations/interface/int_OpMinus.h
index a98f52af..e63f3658 100644
--- a/src/operations/interface/int_OpMinus.h
+++ b/src/operations/interface/int_OpMinus.h
@@ -37,66 +37,77 @@
 
 /* Matrix + Scalar */
 
-/*On transforme le scalaire en tableau et on se ramene à s2s2...s2*/
-#define s2s0OpMinuss2(in1,size,in2,out)		sfilla((float*)out,size[0],size[1],in2);\
-								s2s2OpMinuss2(in1,size,(float*)out,size,out)
+
+#define s2s0OpMinuss2(in1,size,in2,out)		{int i=0;\
+								for (i=0;i<size[0]*size[1];i++) out[i]=in1[i]-in2;\
+								}
+								
+
+#define d2d0OpMinusd2(in1,size,in2,out)		{int i=0;\
+								for (i=0;i<size[0]*size[1];i++) out[i]=in1[i]-in2;\
+								}
 								
-/*On transforme le scalaire en tableau et on se ramene à d2d2...d2*/
-#define d2d0OpMinusd2(in1,size,in2,out)		dfilla((double*)out,size[0],size[1],in2);\
-								d2d2OpMinusd2(in1,size,(double*)out,size,out)
 								
-/*On transforme le scalaire en tableau et on se ramene à c2c2...c2*/								
-#define c2c0OpMinusc2(in1,size,in2,out)		cfilla((floatComplex*)out,size[0],size[1],in2);\
-								c2c2OpMinusc2(in1,size,(floatComplex*)out,size,out)
+#define c2c0OpMinusc2(in1,size,in2,out)		{int i=0;\
+								for (i=0;i<size[0]*size[1];i++) out[i]=cdiffs(in1[i],in2);\
+								}
 								
 /*On transforme le scalaire float en scalaire floatComplex pour pouvoir utiliser c2c0...c2*/	
 #define c2s0OpMinusc2(in1,size,in2,out)		c2c0OpMinusc2(in1,size,FloatComplex(in2,0),out)
 
-/*On transforme le tableau de float en tableau de floatComplex pour pouvoir utiliser c2c0...c2*/	
-#define s2c0OpMinusc2(in1,size,in2,out)		sfilla((float*)out,size[0],size[1],0);\
-								c2c0OpMinusc2(FloatComplexMatrix(in1,(float*)out,size[0]*size[1]),size,in2,out)
+
+#define s2c0OpMinusc2(in1,size,in2,out)		{int i=0;\
+								for (i=0;i<size[0]*size[1];i++) out[i]=cdiffs(FloatComplex(in1[i],0),in2);\
+								}
 								
-/*On transforme le scalaire en tableau et on se ramene à z2z2...z2*/
-#define z2z0OpMinusz2(in1,size,in2,out)		zfilla((doubleComplex*)out,size[0],size[1],in2);\
-								z2z2OpMinusz2(in1,size,(doubleComplex*)out,size,out)
+
+#define z2z0OpMinusz2(in1,size,in2,out)		{int i=0;\
+								for (i=0;i<size[0]*size[1];i++) out[i]=zdiffs(in1[i],in2);\
+								}
 								
 /*On transforme le scalaire double en scalaire doubleComplex pour pouvoir utiliser z2z0...z2*/
 #define z2d0OpMinusz2(in1,size,in2,out)		z2z0OpMinusz2(in1,size,DoubleComplex(in2,0),out)
 
-/*On transforme le tableau de double en tableau de doubleComplex pour pouvoir utiliser z2z0...z2*/
-#define d2z0OpMinusz2(in1,size,in2,out)		dfilla((double*)out,size[0],size[1],0);\
-								z2z0OpMinusz2(DoubleComplexMatrix(in1,(double*)out,size[0]*size[1]),size,in2,out)
+
+#define d2z0OpMinusz2(in1,size,in2,out)		{int i=0;\
+								for (i=0;i<size[0]*size[1];i++) out[i]=zdiffs(DoubleComplex(in1[i],0),in2);\
+								}
 
 
 /* Scalar + Matrix */
 
-/*On transforme le scalaire en tableau et on se ramene à s2s2...s2*/
-#define s0s2OpMinuss2(in1,in2,size,out)		sfilla((float*)out,size[0],size[1],in1);\
-								s2s2OpMinuss2((float*)out,size,in2,size,out)
+
+#define s0s2OpMinuss2(in1,in2,size,out)		{int i=0;\
+								for (i=0;i<size[0]*size[1];i++) out[i]=in1-in2[i];\
+								}
 								
-/*On transforme le scalaire en tableau et on se ramene à d2d2...d2*/
-#define d0d2OpMinusd2(in1,in2,size,out)		dfilla((double*)out,size[0],size[1],in1);\
-								d2d2OpMinusd2((double*)out,size,in2,size,out)
+
+#define d0d2OpMinusd2(in1,in2,size,out)		{int i=0;\
+								for (i=0;i<size[0]*size[1];i++) out[i]=in1-in2[i];\
+								}
 								
-/*On transforme le scalaire en tableau et on se ramene à c2c2...c2*/
-#define c0c2OpMinusc2(in1,in2,size,out)		cfilla((floatComplex*)out,size[0],size[1],in1);\
-								c2c2OpMinusc2((floatComplex*)out,size,in2,size,out)
+
+#define c0c2OpMinusc2(in1,in2,size,out)		{int i=0;\
+								for (i=0;i<size[0]*size[1];i++) out[i]=cdiffs(in1,in2[i]);\
+								}
 								
 /*On transforme le scalaire float en scalaire floatComplex pour pouvoir utiliser c0c2...c2*/	
 #define s0c2OpMinusc2(in1,in2,size,out)		c0c2OpMinusc2(FloatComplex(in1,0),in2,size,out)	
 							
-/*On transforme le tableau de float en tableau de floatComplex pour pouvoir utiliser c0c2...c2*/
-#define c0s2OpMinusc2(in1,in2,size,out)		sfilla((float*)out,size[0],size[1],0);\
-								c0c2OpMinusc2(in1,FloatComplexMatrix(in2,(float*)out,size[0]*size[1]),size,out)
+
+#define c0s2OpMinusc2(in1,in2,size,out)		{int i=0;\
+								for (i=0;i<size[0]*size[1];i++) out[i]=cdiffs(in1,FloatComplex(in2[i],0));\
+								}
 								
-/*On transforme le scalaire en tableau et on se ramene à z2z2...z2*/
-#define z0z2OpMinusz2(in1,in2,size,out)		zfilla((doubleComplex*)out,size[0],size[1],in1);\
-								z2z2OpMinusz2((doubleComplex*)out,size,in2,size,out)
+
+#define z0z2OpMinusz2(in1,in2,size,out)		{int i=0;\
+								for (i=0;i<size[0]*size[1];i++) out[i]=zdiffs(in1,in2[i]);\
+								}
 								
 /*On transforme le scalaire double en scalaire doubleComplex pour pouvoir utiliser z0z2...z2*/
 #define d0z2OpMinusz2(in1,in2,size,out)		z0z2OpMinusz2(DoubleComplex(in1,0),in2,size,out)
 
-/*On transforme le tableau de double en tableau de doubleComplex pour pouvoir utiliser z0z2...z2*/
+
 #define z0d2OpMinusz2(in1,in2,size,out)		dfilla((double*)out,size[0],size[1],0);\
 								z0z2OpMinusz2(in1,DoubleComplexMatrix(in2,(double*)out,size[0]*size[1]),size,out)
 								
@@ -110,23 +121,23 @@
 
 #define c2c2OpMinusc2(in1,size1,in2,size2,out)		cdiffa(in1, size1[0]*size1[1], in2, size2[0]*size2[1], out)
 
-/*On transforme le tableau de float en tableau de floatComplex pour pouvoir utiliser c2c2...c2*/
-#define s2c2OpMinusc2(in1,size1,in2,size2,out)		sfilla((float*)out,size1[0],size1[1],0);\
-									c2c2OpMinusc2(FloatComplexMatrix(in1,(float*)out,size1[0]*size1[1]), size1, in2, size2, out)
+
+#define s2c2OpMinusc2(in1,size1,in2,size2,out)		{int i;\
+									for (i=0;i<size[0]*size[1];i++) out[i]=cdiffs(FloatComplex(in1[i],0),in2[i]);}
 									
-/*On transforme le tableau de float en tableau de floatComplex pour pouvoir utiliser c2c2...c2*/
-#define c2s2OpMinusc2(in1,size1,in2,size2,out)		sfilla((float*)out,size2[0],size2[1],0);\
-									c2c2OpMinusc2(in1, size1, FloatComplexMatrix(in2,(float*)out,size2[0]*size2[1]), size2, out)
+
+#define c2s2OpMinusc2(in1,size1,in2,size2,out)		{int i;\
+									for (i=0;i<size[0]*size[1];i++) out[i]=cdiffs(in1[i],FloatComplex(in2[i],0));}
 
 #define z2z2OpMinusz2(in1,size1,in2,size2,out)		zdiffa(in1, size1[0]*size1[1], in2, size2[0]*size2[1], out)
 
-/*On transforme le tableau de double en tableau de doubleComplex pour pouvoir utiliser z2z2...z2*/
-#define d2z2OpMinusz2(in1,size1,in2,size2,out)		dfilla((double*)out,size1[0],size1[1],0);\
-									z2z2OpMinusz2(DoubleComplexMatrix(in1,(double*)out,size1[0]*size1[1]), size1, in2, size2, out)
+
+#define d2z2OpMinusz2(in1,size1,in2,size2,out)		{int i;\
+									for (i=0;i<size[0]*size[1];i++) out[i]=zdiffs(DoubleComplex(in1[i],0),in2[i]);}
 									
-/*On transforme le tableau de double en tableau de doubleComplex pour pouvoir utiliser z2z2...z2*/
-#define z2d2OpMinusz2(in1,size1,in2,size2,out)		dfilla((double*)out,size2[0],size2[1],0);\
-									z2z2OpMinusz2(in1, size1, DoubleComplexMatrix(in2,(double*)out,size2[0]*size2[1]), size2, out)
+
+#define z2d2OpMinusz2(in1,size1,in2,size2,out)		{int i;\
+									for (i=0;i<size[0]*size[1];i++) out[i]=zdiffs(in1[i],DoubleComplex(in2[i],0));}
 
 
 
-- 
cgit