adjusting some headers

1 files changed, 84 insertions, 41 deletions

diff --git a/src/operations/interface/int_OpMinus.h b/src/operations/interface/int_OpMinus.h
index ccdf32df..a98f52af 100644
--- a/src/operations/interface/int_OpMinus.h
+++ b/src/operations/interface/int_OpMinus.h
@@ -34,43 +34,75 @@
 #define z0d0OpMinusz0(in1,in2)				zdiffs(in1,DoubleComplex(in2,0))
 
 
-/* Scalar - Matrix */
 
-#define s0s2OpMinuss2(in1,in2,size,out)			for(int i=0;i<size[0]*size[1];i++) out[i]= sdiffs(in1,in2[i]);
-
-#define d0d2OpMinusd2(in1,in2,size,out)			for(int i=0;i<size[0]*size[1];i++) out[i]= ddiffs(in1,in2[i]);
-
-#define c0c2OpMinusc2(in1,in2,size,out)			for(int i=0;i<size[0]*size[1];i++) out[i]= cdiffs(in1,in2[i]);
-
-#define z0z2OpMinusz2(in1,in2,size,out)			for(int i=0;i<size[0]*size[1];i++) out[i]= zdiffs(in1,in2[i]);
-
-#define s0c2OpMinusc2(in1,in2,size,out)			for(int i=0;i<size[0]*size[1];i++) out[i]= cdiffs(FloatComplex(in1,0),in2[i]);
-
-#define d0z2OpMinusz2(in1,in2,size,out)			for(int i=0;i<size[0]*size[1];i++) out[i]= zdiffs(DoubleComplex(in1,0),in2[i]);
-
-#define c0s2OpMinusc2(in1,in2,size,out)			for(int i=0;i<size[0]*size[1];i++) out[i]= cdiffs(in1,FloatComplex(in2[i],0));
-
-#define z0d2OpMinusz2(in1,in2,size,out)			for(int i=0;i<size[0]*size[1];i++) out[i]= zdiffs(in1,DoubleComplex(in2[i],0));
-
-/* Matrix - Scalar */
-
-#define s2s0OpMinuss2(in1,size,in2,out)			for(int i=0;i<size[0]*size[1];i++) out[i]= sdiffs(in1[i],in2);
-
-#define d2d0OpMinusd2(in1,size,in2,out)			for(int i=0;i<size[0]*size[1];i++) out[i]= ddiffs(in1[i],in2);
-
-#define c2c0OpMinusc2(in1,size,in2,out)			for(int i=0;i<size[0]*size[1];i++) out[i]= cdiffs(in1[i],in2);
-
-#define z2z0OpMinusz2(in1,size,in2,out)			for(int i=0;i<size[0]*size[1];i++) out[i]= zdiffs(in1[i],in2);
-
-#define s2c0OpMinusc2(in1,size,in2,out)			for(int i=0;i<size[0]*size[1];i++) out[i]= cdiffs(FloatComplex(in1[i],0),in2);
-
-#define d2z0OpMinusz2(in1,size,in2,out)			for(int i=0;i<size[0]*size[1];i++) out[i]= zdiffs(DoubleComplex(in1[i],0),in2);
-
-#define c2s0OpMinusc2(in1,size,in2,out)			for(int i=0;i<size[0]*size[1];i++) out[i]= cdiffs(in1[i],FloatComplex(in2,0));
-
-#define z2d0OpMinusz2(in1,size,in2,out)			for(int i=0;i<size[0]*size[1];i++) out[i]= zdiffs(in1[i],DoubleComplex(in2,0));
-
-/* Matrix - Matrix */
+/* 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)
+								
+/*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)
+								
+/*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)
+								
+/*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)
+								
+/*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)
+
+
+/* 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)
+								
+/*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)
+								
+/*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)
+								
+/*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)
+								
+/*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)
+								
+/*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)
+								
+
+
+/* Matrix + Matrix */
 
 #define s2s2OpMinuss2(in1,size1,in2,size2,out)		sdiffa(in1, size1[0]*size1[1], in2, size2[0]*size2[1], out)
 
@@ -78,14 +110,25 @@
 
 #define c2c2OpMinusc2(in1,size1,in2,size2,out)		cdiffa(in1, size1[0]*size1[1], in2, size2[0]*size2[1], out)
 
-#define c2s2OpMinusc2(in1,size1,in2,size2,out)		cdiffa(in1, size1[0]*size1[1], FloatComplexMatrix(in2,0,size2[0]*size2[1]), size2[0]*size2[1], out)
-
-#define s2c2OpMinusc2(in1,size1,in2,size2,out)		cdiffa(FloatComplexMatrix(in1,0,size1[0]*size1[1]), 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)
+									
+/*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 z2z2OpMinusz2(in1,size1,in2,size2,out)		zdiffa(in1, size1[0]*size1[1], in2, size2[0]*size2[1], out)
 
-#define z2d2OpMinusz2(in1,size1,in2,size2,out)		zdiffa(in1, size1[0]*size1[1], DoubleComplexMatrix(in2,0,size2[0]*size2[1]), 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)
+									
+/*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 d2z2OpMinusz2(in1,size1,in2,size2,out)		zdiffa(DoubleComplexMatrix(in1,0,size1[0]*size1[1]), size1[0]*size1[1], in2, size2[0]*size2[1], out)
 
 #endif /* !__INT_OPMINUS_H__ */