Add and modify sign files

3 files changed, 36 insertions, 0 deletions

diff --git a/src/auxiliaryFunctions/interfaces/int_sign.h b/src/auxiliaryFunctions/interfaces/int_sign.h
new file mode 100644
index 00000000..4463901a
--- /dev/null
+++ b/src/auxiliaryFunctions/interfaces/int_sign.h
@@ -0,0 +1,34 @@
+/*
+ * Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ * Copyright (C) 2008-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
+ *
+ */
+
+/* THIS IS AN AUTOMATICALLY GENERATED FILE : DO NOT EDIT BY HAND. */
+
+#ifndef __INT_SIGN_H__
+#define __INT_SIGN_H__
+
+#define s0signs0(in)			ssigns(in)
+
+#define d0signd0(in)			dsigns(in)
+
+#define c0signc0(in)			csigns(in)
+
+#define z0signz0(in)			zsigns(in)
+
+#define s2signs2(in,size,out)		ssigna(in, size[0]*size[1], out)
+
+#define d2signd2(in,size,out)		dsigna(in, size[0]*size[1], out)
+
+#define c2signc2(in,size,out)		csigna(in, size[0]*size[1], out)
+
+#define z2signz2(in,size,out)		zsigna(in, size[0]*size[1], out)
+
+#endif /* !__INT_SIGN_H__ */
diff --git a/src/auxiliaryFunctions/sign/csigns.c b/src/auxiliaryFunctions/sign/csigns.c
index 218314ab..b054046f 100644
--- a/src/auxiliaryFunctions/sign/csigns.c
+++ b/src/auxiliaryFunctions/sign/csigns.c
@@ -14,5 +14,6 @@
 #include "abs.h"
 
 floatComplex csigns(floatComplex in) {
+  if ( (creals(in)==0) && (cimags(in)==0) ) return FloatComplex(0,0);
   return FloatComplex(creals(in) / cabss(in), cimags(in) / cabss(in));
 }
diff --git a/src/auxiliaryFunctions/sign/zsigns.c b/src/auxiliaryFunctions/sign/zsigns.c
index 397cc33a..a4fed8d4 100644
--- a/src/auxiliaryFunctions/sign/zsigns.c
+++ b/src/auxiliaryFunctions/sign/zsigns.c
@@ -14,5 +14,6 @@
 #include "abs.h"
 
 doubleComplex zsigns(doubleComplex in) {
+  if ( (zreals(in)==0) && (zimags(in)==0) ) return DoubleComplex(0,0);
   return DoubleComplex(zreals(in) / zabss(in), zimags(in) / zabss(in));
 }