added test for Spec function

corrected a missing ")" in ifft interface
updated INIT_FillSCI2LibCDirs

1 files changed, 26 insertions, 9 deletions

diff --git a/src/c/matrixOperations/interfaces/int_spec.h b/src/c/matrixOperations/interfaces/int_spec.h
index a2d9e840..e25d7219 100644
--- a/src/c/matrixOperations/interfaces/int_spec.h
+++ b/src/c/matrixOperations/interfaces/int_spec.h
@@ -15,19 +15,27 @@
 #ifndef __INT_SPEC_H__
 #define __INT_SPEC_H__
 
-/* Only eigenvalues */
+/* Only eigenvalues */ 
 
-#define s0specc0(in)			sspecs(in)
+#define s0specc0(in)			FloatComplex  (sspecs(in)  , 0)
 
-#define d0specz0(in)			dspecs(in)
+#define d0specz0(in)			DoubleComplex ( dspecs(in) , 0)
 
 #define c0specc0(in)			cspecs(in)
 
 #define z0specz0(in)			zspecs(in)
 
-#define s2specc2(in,size,out)		sspeca(in, size[0], out)
+#define s2specc2(in,size,out)		{float* ZEROS;\
+						ZEROS=malloc((uint)(size[0]*size[0]*sizeof(float)));\
+						szerosa(ZEROS,size[0],size[0]);\
+						sspeca(FloatComplexMatrix(in,ZEROS,size[0]*size[0]), size[0], out);\
+						}
 
-#define d2specz2(in,size,out)		dspeca(in, size[0], out)
+#define d2specz2(in,size,out)		{double* ZEROS;\
+						ZEROS=malloc((uint)(size[0]*size[0]*sizeof(double)));\
+						dzerosa(ZEROS,size[0],size[0]);\
+						zspeca(DoubleComplexMatrix(in,ZEROS,size[0]*size[0]), size[0] , out);\
+						}
 
 #define c2specc2(in,size,out)		cspeca(in, size[0], out)
 
@@ -35,17 +43,26 @@
 
 /* Eigenvalues and eigenvectors */
 
-#define s0specc0c0(in,out)			sspec2s(in,out)
+#define s0specc0c0(in,out)			cspec2s(FloatComplex( in,0) ,out)
 
-#define d0specz0z0(in,out)			dspec2s(in,out)
+#define d0specz0z0(in,out)			zspec2s(DoubleComplex(in,0) ,out)
 
 #define c0specc0c0(in,out)			cspec2s(in,out)
 
 #define z0specz0z0(in,out)			zspec2s(in,out)
 
-#define s2specc2c2(in,size,out1,out2)		sspec2a(in, size[0], out2, out1)
+#define s2specc2c2(in,size,out1,out2)		{float* ZEROS;\
+						ZEROS=malloc((uint)(size[0]*size[0]*sizeof(float)));\
+						szerosa(ZEROS,size[0],size[0]);\
+						cspec2a(FloatComplexMatrix(in,ZEROS,size[0]*size[0]), size[0] , out1 , out2 );\
+						}
+
+#define d2specz2z2(in,size,out1,out2)		{double* ZEROS;\
+						ZEROS=malloc((uint)(size[0]*size[0]*sizeof(double)));\
+						dzerosa(ZEROS,size[0],size[0]);\
+						zspec2a(DoubleComplexMatrix(in,ZEROS,size[0]*size[0]), size[0] , out1 , out2 );\
+						}
 
-#define d2specz2z2(in,size,out1,out2)		dspec2a(in, size[0], out2, out1)
 
 #define c2specc2c2(in,size,out1,out2)		cspec2a(in, size[0], out2, out1)