Tested Addition and Subtraction for new data types

43 files changed, 1118 insertions, 0 deletions

diff --git a/src/c/elementaryFunctions/includes/int16.h b/src/c/elementaryFunctions/includes/int16.h
new file mode 100644
index 0000000..f327184
--- /dev/null
+++ b/src/c/elementaryFunctions/includes/int16.h
@@ -0,0 +1,42 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2006-2008 - INRIA - Arnaud Torset
+ *
+ *  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
+ *
+ */
+
+#ifndef __INT16_H__
+#define __INT16_H__
+
+#include "dynlib_elementaryfunctions.h"
+#include "floatComplex.h"
+#include "doubleComplex.h"
+#include "types.h"
+
+#ifdef  __cplusplus
+extern "C" {
+#endif
+/*
+	int16(X) returns the signed int part of the double value
+*/
+
+EXTERN_ELEMFUNCT int16	sint16s(float in);
+
+EXTERN_ELEMFUNCT int16	dint816s(double in);
+
+EXTERN_ELEMFUNCT void	sint16a(float* in, int size, int16* out);
+
+EXTERN_ELEMFUNCT void	dint16a(double* in, int size, int16* out);
+
+#ifdef  __cplusplus
+} /* extern "C" */
+#endif
+
+
+#endif /* !__INT16_H__ */
+
diff --git a/src/c/elementaryFunctions/includes/int8.h b/src/c/elementaryFunctions/includes/int8.h
new file mode 100644
index 0000000..12fc595
--- /dev/null
+++ b/src/c/elementaryFunctions/includes/int8.h
@@ -0,0 +1,42 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2006-2008 - INRIA - Arnaud Torset
+ *
+ *  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
+ *
+ */
+
+#ifndef __INT8_H__
+#define __INT8_H__
+
+#include "dynlib_elementaryfunctions.h"
+#include "floatComplex.h"
+#include "doubleComplex.h"
+#include "types.h"
+
+#ifdef  __cplusplus
+extern "C" {
+#endif
+/*
+	int8(X) returns the signed char part of the double value
+*/
+
+EXTERN_ELEMFUNCT int8	sint8s(float in);
+
+EXTERN_ELEMFUNCT int8 dint8s(double in);
+
+EXTERN_ELEMFUNCT void	sint8a(float* in, int size, int8* out);
+
+EXTERN_ELEMFUNCT void	dint8a(double* in, int size, int8* out);
+
+#ifdef  __cplusplus
+} /* extern "C" */
+#endif
+
+
+#endif /* !__INT8_H__ */
+
diff --git a/src/c/elementaryFunctions/includes/uint16.h b/src/c/elementaryFunctions/includes/uint16.h
new file mode 100644
index 0000000..a0739b6
--- /dev/null
+++ b/src/c/elementaryFunctions/includes/uint16.h
@@ -0,0 +1,42 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2006-2008 - INRIA - Arnaud Torset
+ *
+ *  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
+ *
+ */
+
+#ifndef __UINT16_H__
+#define __UINT16_H__
+
+#include "dynlib_elementaryfunctions.h"
+#include "floatComplex.h"
+#include "doubleComplex.h"
+#include "types.h"
+
+#ifdef  __cplusplus
+extern "C" {
+#endif
+/*
+	uint16(X) returns the signed int part of the double value
+*/
+
+EXTERN_ELEMFUNCT uint16	suint16s(float in);
+
+EXTERN_ELEMFUNCT uint16	duint816s(double in);
+
+EXTERN_ELEMFUNCT void	suint16a(float* in, int size, uint16* out);
+
+EXTERN_ELEMFUNCT void	duint16a(double* in, int size, uint16* out);
+
+#ifdef  __cplusplus
+} /* extern "C" */
+#endif
+
+
+#endif /* !__UINT16_H__ */
+
diff --git a/src/c/elementaryFunctions/includes/uint8.h b/src/c/elementaryFunctions/includes/uint8.h
new file mode 100644
index 0000000..9d3dd39
--- /dev/null
+++ b/src/c/elementaryFunctions/includes/uint8.h
@@ -0,0 +1,42 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2006-2008 - INRIA - Arnaud Torset
+ *
+ *  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
+ *
+ */
+
+#ifndef __UINT8_H__
+#define __UINT8_H__
+
+#include "dynlib_elementaryfunctions.h"
+#include "floatComplex.h"
+#include "doubleComplex.h"
+#include "types.h"
+
+#ifdef  __cplusplus
+extern "C" {
+#endif
+/*
+	uint8(X) returns the unsigned char part of the double value
+*/
+
+EXTERN_ELEMFUNCT uint8	suint8s(float in);
+
+EXTERN_ELEMFUNCT uint8 duint8s(double in);
+
+EXTERN_ELEMFUNCT void	suint8a(float* in, int size, uint8* out);
+
+EXTERN_ELEMFUNCT void	duint8a(double* in, int size, uint8* out);
+
+#ifdef  __cplusplus
+} /* extern "C" */
+#endif
+
+
+#endif /* !__UINT8_H__ */
+
diff --git a/src/c/elementaryFunctions/int16/dint16a.c b/src/c/elementaryFunctions/int16/dint16a.c
new file mode 100644
index 0000000..73c9ec7
--- /dev/null
+++ b/src/c/elementaryFunctions/int16/dint16a.c
@@ -0,0 +1,18 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008-2008 - INRIA - Arnaud Torset
+ *
+ *  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 "int16.h"
+
+void dint16a(double* x, int size, int16* out) {
+	int i = 0;
+	for (i=0;i<size;i++) out[i] = dint16s(x[i]);
+}
diff --git a/src/c/elementaryFunctions/int16/dint16s.c b/src/c/elementaryFunctions/int16/dint16s.c
new file mode 100644
index 0000000..76ee96b
--- /dev/null
+++ b/src/c/elementaryFunctions/int16/dint16s.c
@@ -0,0 +1,18 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008-2008 - INRIA - Arnaud Torset
+ *
+ *  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 "int16.h"
+
+int16 dint16s(double x) {
+	return (int16) x;
+}
diff --git a/src/c/elementaryFunctions/int16/sint16a.c b/src/c/elementaryFunctions/int16/sint16a.c
new file mode 100644
index 0000000..33058ea
--- /dev/null
+++ b/src/c/elementaryFunctions/int16/sint16a.c
@@ -0,0 +1,18 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008-2008 - INRIA - Arnaud Torset
+ *
+ *  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 "int16.h"
+
+void sint16a(float* x, int size, int16* out) {
+	int i = 0;
+	for (i=0;i<size;i++) out[i] = sint16s(x[i]);
+}
diff --git a/src/c/elementaryFunctions/int16/sint16s.c b/src/c/elementaryFunctions/int16/sint16s.c
new file mode 100644
index 0000000..4634a39
--- /dev/null
+++ b/src/c/elementaryFunctions/int16/sint16s.c
@@ -0,0 +1,18 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008-2008 - INRIA - Arnaud Torset
+ *
+ *  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 "int16.h"
+
+int16 sint16s(float x) {
+	return (int16) x;
+}
diff --git a/src/c/elementaryFunctions/int8/dint8a.c b/src/c/elementaryFunctions/int8/dint8a.c
new file mode 100644
index 0000000..719d6f4
--- /dev/null
+++ b/src/c/elementaryFunctions/int8/dint8a.c
@@ -0,0 +1,18 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008-2008 - INRIA - Arnaud Torset
+ *
+ *  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 "int8.h"
+
+void dint8a(double* x, int size, int8* out) {
+	int i = 0;
+	for (i=0;i<size;i++) out[i] = dint8s(x[i]);
+}
diff --git a/src/c/elementaryFunctions/int8/dint8s.c b/src/c/elementaryFunctions/int8/dint8s.c
new file mode 100644
index 0000000..f1789e6
--- /dev/null
+++ b/src/c/elementaryFunctions/int8/dint8s.c
@@ -0,0 +1,18 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008-2008 - INRIA - Arnaud Torset
+ *
+ *  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 "int8.h"
+
+int8 dint8s(double x) {
+	return (int8) x;
+}
diff --git a/src/c/elementaryFunctions/int8/sint8a.c b/src/c/elementaryFunctions/int8/sint8a.c
new file mode 100644
index 0000000..87846a7
--- /dev/null
+++ b/src/c/elementaryFunctions/int8/sint8a.c
@@ -0,0 +1,18 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008-2008 - INRIA - Arnaud Torset
+ *
+ *  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 "int8.h"
+
+void sint8a(float* x, int size, int8* out) {
+	int i = 0;
+	for (i=0;i<size;i++) out[i] = sint8s(x[i]);
+}
diff --git a/src/c/elementaryFunctions/int8/sint8s.c b/src/c/elementaryFunctions/int8/sint8s.c
new file mode 100644
index 0000000..89d4b81
--- /dev/null
+++ b/src/c/elementaryFunctions/int8/sint8s.c
@@ -0,0 +1,18 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008-2008 - INRIA - Arnaud Torset
+ *
+ *  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 "int8.h"
+
+int8 sint8s(float x) {
+	return (int8) x;
+}
diff --git a/src/c/elementaryFunctions/interfaces/int_int16.h b/src/c/elementaryFunctions/interfaces/int_int16.h
new file mode 100644
index 0000000..759cb97
--- /dev/null
+++ b/src/c/elementaryFunctions/interfaces/int_int16.h
@@ -0,0 +1,26 @@
+/*
+ * 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_INT16_H__
+#define __INT_INT16_H__
+
+#define s0int16u80(in)			sint16s(in)
+
+#define d0int16u80(in)			dint16s(in)
+
+#define s2int16u82(in,size,out)		sint16a(in, size[0]*size[1], out)
+
+#define d2int16u82(in,size,out)		dint16a(in, size[0]*size[1], out)
+
+#endif /* !__INT_INT8_H__ */
diff --git a/src/c/elementaryFunctions/interfaces/int_int8.h b/src/c/elementaryFunctions/interfaces/int_int8.h
new file mode 100644
index 0000000..5cf129e
--- /dev/null
+++ b/src/c/elementaryFunctions/interfaces/int_int8.h
@@ -0,0 +1,26 @@
+/*
+ * 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_INT8_H__
+#define __INT_INT8_H__
+
+#define s0int8u80(in)			sint8s(in)
+
+#define d0int8u80(in)			dint8s(in)
+
+#define s2int8u82(in,size,out)		sint8a(in, size[0]*size[1], out)
+
+#define d2int8u82(in,size,out)		dint8a(in, size[0]*size[1], out)
+
+#endif /* !__INT_INT8_H__ */
diff --git a/src/c/elementaryFunctions/interfaces/int_uint16.h b/src/c/elementaryFunctions/interfaces/int_uint16.h
new file mode 100644
index 0000000..0d558b3
--- /dev/null
+++ b/src/c/elementaryFunctions/interfaces/int_uint16.h
@@ -0,0 +1,26 @@
+/*
+ * 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_UINT16_H__
+#define __INT_UINT16_H__
+
+#define s0uint16u80(in)			suint16s(in)
+
+#define d0uint16u80(in)			duint16s(in)
+
+#define s2uint16u82(in,size,out)		suint16a(in, size[0]*size[1], out)
+
+#define d2uint16u82(in,size,out)		duint16a(in, size[0]*size[1], out)
+
+#endif /* !__INT_UINT8_H__ */
diff --git a/src/c/elementaryFunctions/interfaces/int_uint8.h b/src/c/elementaryFunctions/interfaces/int_uint8.h
new file mode 100644
index 0000000..d12b275
--- /dev/null
+++ b/src/c/elementaryFunctions/interfaces/int_uint8.h
@@ -0,0 +1,26 @@
+/*
+ * 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_UINT8_H__
+#define __INT_UINT8_H__
+
+#define s0uint8u80(in)			suint8s(in)
+
+#define d0uint8u80(in)			duint8s(in)
+
+#define s2uint8u82(in,size,out)		suint8a(in, size[0]*size[1], out)
+
+#define d2uint8u82(in,size,out)		duint8a(in, size[0]*size[1], out)
+
+#endif /* !__INT_UINT8_H__ */
diff --git a/src/c/elementaryFunctions/uint16/duint16a.c b/src/c/elementaryFunctions/uint16/duint16a.c
new file mode 100644
index 0000000..c94cd18
--- /dev/null
+++ b/src/c/elementaryFunctions/uint16/duint16a.c
@@ -0,0 +1,18 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008-2008 - INRIA - Arnaud Torset
+ *
+ *  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 "uint16.h"
+
+void duint16a(double* x, int size, uint16* out) {
+	int i = 0;
+	for (i=0;i<size;i++) out[i] = duint16s(x[i]);
+}
diff --git a/src/c/elementaryFunctions/uint16/duint16s.c b/src/c/elementaryFunctions/uint16/duint16s.c
new file mode 100644
index 0000000..4c4918f
--- /dev/null
+++ b/src/c/elementaryFunctions/uint16/duint16s.c
@@ -0,0 +1,18 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008-2008 - INRIA - Arnaud Torset
+ *
+ *  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 "uint16.h"
+
+uint16 duint16s(double x) {
+	return (uint16) x;
+}
diff --git a/src/c/elementaryFunctions/uint16/suint16a.c b/src/c/elementaryFunctions/uint16/suint16a.c
new file mode 100644
index 0000000..4bc3768
--- /dev/null
+++ b/src/c/elementaryFunctions/uint16/suint16a.c
@@ -0,0 +1,18 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008-2008 - INRIA - Arnaud Torset
+ *
+ *  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 "uint16.h"
+
+void suint16a(float* x, int size, uint16* out) {
+	int i = 0;
+	for (i=0;i<size;i++) out[i] = suint16s(x[i]);
+}
diff --git a/src/c/elementaryFunctions/uint16/suint16s.c b/src/c/elementaryFunctions/uint16/suint16s.c
new file mode 100644
index 0000000..708237b
--- /dev/null
+++ b/src/c/elementaryFunctions/uint16/suint16s.c
@@ -0,0 +1,18 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008-2008 - INRIA - Arnaud Torset
+ *
+ *  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 "uint16.h"
+
+uint16 suint16s(float x) {
+	return (uint16) x;
+}
diff --git a/src/c/elementaryFunctions/uint8/duint8a.c b/src/c/elementaryFunctions/uint8/duint8a.c
new file mode 100644
index 0000000..0a3f2d2
--- /dev/null
+++ b/src/c/elementaryFunctions/uint8/duint8a.c
@@ -0,0 +1,18 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008-2008 - INRIA - Arnaud Torset
+ *
+ *  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 "uint8.h"
+
+void duint8a(double* x, int size, uint8* out) {
+	int i = 0;
+	for (i=0;i<size;i++) out[i] = duint8s(x[i]);
+}
diff --git a/src/c/elementaryFunctions/uint8/duint8s.c b/src/c/elementaryFunctions/uint8/duint8s.c
new file mode 100644
index 0000000..b797934
--- /dev/null
+++ b/src/c/elementaryFunctions/uint8/duint8s.c
@@ -0,0 +1,18 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008-2008 - INRIA - Arnaud Torset
+ *
+ *  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 "uint8.h"
+
+uint8 duint8s(double x) {
+	return (uint8) x;
+}
diff --git a/src/c/elementaryFunctions/uint8/suint8a.c b/src/c/elementaryFunctions/uint8/suint8a.c
new file mode 100644
index 0000000..cf64dd4
--- /dev/null
+++ b/src/c/elementaryFunctions/uint8/suint8a.c
@@ -0,0 +1,18 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008-2008 - INRIA - Arnaud Torset
+ *
+ *  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 "uint8.h"
+
+void suint8a(float* x, int size, uint8* out) {
+	int i = 0;
+	for (i=0;i<size;i++) out[i] = suint8s(x[i]);
+}
diff --git a/src/c/elementaryFunctions/uint8/suint8s.c b/src/c/elementaryFunctions/uint8/suint8s.c
new file mode 100644
index 0000000..e95e964
--- /dev/null
+++ b/src/c/elementaryFunctions/uint8/suint8s.c
@@ -0,0 +1,18 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008-2008 - INRIA - Arnaud Torset
+ *
+ *  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 "uint8.h"
+
+uint8 suint8s(float x) {
+	return (uint8) x;
+}
diff --git a/src/c/operations/addition/i16adda.c b/src/c/operations/addition/i16adda.c
new file mode 100644
index 0000000..8c28a98
--- /dev/null
+++ b/src/c/operations/addition/i16adda.c
@@ -0,0 +1,20 @@
+/*
+ *  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
+ *
+ */
+
+#include "addition.h"
+
+void i16adda(int16* in1, int size1, int16* in2, int size2, int16* out) {
+  int i = 0;
+  for (i = 0; i < size1 && i < size2; ++i) {
+    out[i] = i16adds(in1[i], in2[i]);
+  }
+}
diff --git a/src/c/operations/addition/i16adds.c b/src/c/operations/addition/i16adds.c
new file mode 100644
index 0000000..ea36a91
--- /dev/null
+++ b/src/c/operations/addition/i16adds.c
@@ -0,0 +1,18 @@
+/*
+ *  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
+ *
+ */
+
+#include "addition.h"
+
+int16 i16adds(int16 in1, int16 in2)
+{
+  return (in1 + in2);
+}
diff --git a/src/c/operations/addition/i8adds.c b/src/c/operations/addition/i8adds.c
new file mode 100644
index 0000000..aee4045
--- /dev/null
+++ b/src/c/operations/addition/i8adds.c
@@ -0,0 +1,18 @@
+/*
+ *  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
+ *
+ */
+
+#include "addition.h"
+
+int8 i8adds(int8 in1, int8 in2)
+{
+  return (in1 + in2);
+}
diff --git a/src/c/operations/addition/u16adda.c b/src/c/operations/addition/u16adda.c
new file mode 100644
index 0000000..884e099
--- /dev/null
+++ b/src/c/operations/addition/u16adda.c
@@ -0,0 +1,20 @@
+/*
+ *  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
+ *
+ */
+
+#include "addition.h"
+
+void u16adda(uint16* in1, int size1, uint16* in2, int size2, uint16* out) {
+  int i = 0;
+  for (i = 0; i < size1 && i < size2; ++i) {
+    out[i] = u16adds(in1[i], in2[i]);
+  }
+}
diff --git a/src/c/operations/addition/u16adds.c b/src/c/operations/addition/u16adds.c
new file mode 100644
index 0000000..f52dd05
--- /dev/null
+++ b/src/c/operations/addition/u16adds.c
@@ -0,0 +1,18 @@
+/*
+ *  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
+ *
+ */
+
+#include "addition.h"
+
+uint16 u16adds(uint16 in1, uint16 in2)
+{
+  return (in1 + in2);
+}
diff --git a/src/c/operations/addition/u8adda.c b/src/c/operations/addition/u8adda.c
new file mode 100644
index 0000000..4c195b0
--- /dev/null
+++ b/src/c/operations/addition/u8adda.c
@@ -0,0 +1,20 @@
+/*
+ *  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
+ *
+ */
+
+#include "addition.h"
+
+void u8adda(uint8* in1, int size1, uint8* in2, int size2, uint8* out) {
+  int i = 0;
+  for (i = 0; i < size1 && i < size2; ++i) {
+    out[i] = u8adds(in1[i], in2[i]);
+  }
+}
diff --git a/src/c/operations/addition/u8adds.c b/src/c/operations/addition/u8adds.c
new file mode 100644
index 0000000..b729c3b
--- /dev/null
+++ b/src/c/operations/addition/u8adds.c
@@ -0,0 +1,18 @@
+/*
+ *  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
+ *
+ */
+
+#include "addition.h"
+
+uint8 u8adds(uint8 in1, uint8 in2)
+{
+  return (in1 + in2);
+}
diff --git a/src/c/operations/includes/addition.h b/src/c/operations/includes/addition.h
index d6f0338..fbf0177 100644
--- a/src/c/operations/includes/addition.h
+++ b/src/c/operations/includes/addition.h
@@ -16,6 +16,7 @@
 #include "dynlib_operations.h"
 #include "floatComplex.h"
 #include "doubleComplex.h"
+#include "types.h"
 
 #ifdef  __cplusplus
 extern "C" {
@@ -108,6 +109,86 @@ EXTERN_OPERATIONS void	zadda(doubleComplex *in1, int size1,
 	       doubleComplex *in2, int size2,
 	       doubleComplex *out);
 
+/*
+** \brief Compute an addition with uint8.
+** \param in1 : input uint8.
+** \param in2 : input uint8.
+** \return : in1 + in2
+*/
+EXTERN_OPERATIONS uint8	u8adds(uint8 in1, uint8 in2);
+
+/*
+** \brief Compute an addition element wise for uint8.
+** \param in1 : input array.
+** \param size1 : size of in1 array.
+** \param in2 : input arry.
+** \param size2 : size of in2 array.
+** \param out : array that contains the addition in1 + in2.
+*/
+EXTERN_OPERATIONS void	u8adda(uint8 *in1, int size1,
+	       uint8 *in2, int size2,
+	       uint8 *out);
+
+/*
+** \brief Compute an addition with int8.
+** \param in1 : input int8.
+** \param in2 : input int8.
+** \return : in1 + in2
+*/
+EXTERN_OPERATIONS int8	i8adds(int8 in1, int8 in2);
+
+/*
+** \brief Compute an addition element wise for int8.
+** \param in1 : input array.
+** \param size1 : size of in1 array.
+** \param in2 : input arry.
+** \param size2 : size of in2 array.
+** \param out : array that contains the addition in1 + in2.
+*/
+EXTERN_OPERATIONS void	i8adda(int8 *in1, int size1,
+	       int8 *in2, int size2,
+	       int8 *out);
+
+/*
+** \brief Compute an addition with uint16.
+** \param in1 : input uint16.
+** \param in2 : input uint16.
+** \return : in1 + in2
+*/
+EXTERN_OPERATIONS uint16 u16adds(uint16 in1, uint16 in2);
+
+/*
+** \brief Compute an addition element wise for uint16.
+** \param in1 : input array.
+** \param size1 : size of in1 array.
+** \param in2 : input arry.
+** \param size2 : size of in2 array.
+** \param out : array that contains the addition in1 + in2.
+*/
+EXTERN_OPERATIONS void	u16adda(uint16 *in1, int size1,
+	       uint16 *in2, int size2,
+	       uint16 *out);
+
+/*
+** \brief Compute an addition with int16.
+** \param in1 : input int16.
+** \param in2 : input int16.
+** \return : in1 + in2
+*/
+EXTERN_OPERATIONS int16	i16adds(int16 in1, int16 in2);
+
+/*
+** \brief Compute an addition element wise for int16.
+** \param in1 : input array.
+** \param size1 : size of in1 array.
+** \param in2 : input arry.
+** \param size2 : size of in2 array.
+** \param out : array that contains the addition in1 + in2.
+*/
+EXTERN_OPERATIONS void	i16adda(int16 *in1, int size1,
+	       int16 *in2, int size2,
+	       int16 *out);
+
 #ifdef  __cplusplus
 } /* extern "C" */
 #endif
diff --git a/src/c/operations/includes/subtraction.h b/src/c/operations/includes/subtraction.h
index 8571a5d..798df2b 100644
--- a/src/c/operations/includes/subtraction.h
+++ b/src/c/operations/includes/subtraction.h
@@ -16,6 +16,7 @@
 #include "dynlib_operations.h"
 #include "floatComplex.h"
 #include "doubleComplex.h"
+#include "types.h" 
 
 #ifdef  __cplusplus
 extern "C" {
@@ -107,6 +108,88 @@ EXTERN_OPERATIONS void	zdiffa(doubleComplex *in1, int size1,
 	       doubleComplex *in2, int size2,
 	       doubleComplex *out);
 
+
+/*
+** \brief Compute a subtraction with uint8.
+** \param in1 : input uint8.
+** \param in2 : input uint8.
+** \return : in1 + in2
+*/
+EXTERN_OPERATIONS uint8	u8diffs(uint8 in1, uint8 in2);
+
+/*
+** \brief Compute a subtraction element wise for uint8.
+** \param in1 : input array.
+** \param size1 : size of in1 array.
+** \param in2 : input arry.
+** \param size2 : size of in2 array.
+** \param out : array that contains the subtraction in1 + in2.
+*/
+EXTERN_OPERATIONS void	u8diffa(uint8 *in1, int size1,
+	       uint8 *in2, int size2,
+	       uint8 *out);
+
+/*
+** \brief Compute a subtraction with int8.
+** \param in1 : input int8.
+** \param in2 : input int8.
+** \return : in1 + in2
+*/
+EXTERN_OPERATIONS int8	i8diffs(int8 in1, int8 in2);
+
+/*
+** \brief Compute a subtraction element wise for int8.
+** \param in1 : input array.
+** \param size1 : size of in1 array.
+** \param in2 : input arry.
+** \param size2 : size of in2 array.
+** \param out : array that contains the subtraction in1 + in2.
+*/
+EXTERN_OPERATIONS void	i8diffa(int8 *in1, int size1,
+	       int8 *in2, int size2,
+	       int8 *out);
+
+
+/*
+** \brief Compute a subtraction with uint16.
+** \param in1 : input uint16.
+** \param in2 : input uint16.
+** \return : in1 + in2
+*/
+EXTERN_OPERATIONS uint16 u16diffs(uint16 in1, uint16 in2);
+
+/*
+** \brief Compute a subtraction element wise for uint16.
+** \param in1 : input array.
+** \param size1 : size of in1 array.
+** \param in2 : input arry.
+** \param size2 : size of in2 array.
+** \param out : array that contains the subtraction in1 + in2.
+*/
+EXTERN_OPERATIONS void	u16diffa(uint16 *in1, int size1,
+	       uint16 *in2, int size2,
+	       uint16 *out);
+
+/*
+** \brief Compute a subtraction with int16.
+** \param in1 : input int16.
+** \param in2 : input int16.
+** \return : in1 + in2
+*/
+EXTERN_OPERATIONS int16	i16diffs(int16 in1, int16 in2);
+
+/*
+** \brief Compute a subtraction element wise for int16.
+** \param in1 : input array.
+** \param size1 : size of in1 array.
+** \param in2 : input arry.
+** \param size2 : size of in2 array.
+** \param out : array that contains the subtraction in1 + in2.
+*/
+EXTERN_OPERATIONS void	i16diffa(int16 *in1, int size1,
+	       int16 *in2, int size2,
+	       int16 *out);
+
 #ifdef  __cplusplus
 } /* extern "C" */
 #endif
diff --git a/src/c/operations/interfaces/int_OpMinus.h b/src/c/operations/interfaces/int_OpMinus.h
index ac27ddb..b96328e 100644
--- a/src/c/operations/interfaces/int_OpMinus.h
+++ b/src/c/operations/interfaces/int_OpMinus.h
@@ -25,6 +25,14 @@
 
 #define z0OpMinusz0(in)					DoubleComplex(-zreals(in), -zimags(in))
 
+#define u80OpMinusu80(in)				-(int8)in
+
+#define i80OpMinusi80(in)				-in
+
+#define u160OpMinusu160(in)				-(int16)in
+
+#define i160OpMinusi160(in)				-in
+
 /* - Matrix */
 
 #define s2OpMinuss2(in, size, out)			{int i=0; \
@@ -43,6 +51,19 @@
     for (i=0;i<size[0]*size[1];i++) out[i] = DoubleComplex(-zreals(in[i]), -zimags(in[i])); \
   }
 
+#define u82OpMinusu82(in, size, out)			{int i=0; \
+    for (i=0;i<size[0]*size[1];i++) out[i] = -(int8)in[i];		  \
+  }
+
+
+#define u82OpMinusu82(in, size, out)			{int i=0; \
+    for (i=0;i<size[0]*size[1];i++) out[i] = -in[i];		  \
+  }
+
+
+#define u162OpMinusu162(in, size, out)			{int i=0; \
+    for (i=0;i<size[0]*size[1];i++) out[i] = -(int16)in[i];		  \
+  }
 
 /* Scalar - Scalar */
 
@@ -62,6 +83,14 @@
 
 #define z0d0OpMinusz0(in1,in2)				zdiffs(in1,DoubleComplex(in2,0))
 
+#define u80u80OpMinusu80(in1,in2)			u8diffs(in1,in2)
+
+#define i80i80OpMinusi80(in1,in2)			i8diffs(in1,in2)
+
+#define u160u160OpMinusu160(in1,in2)			u16diffs(in1,in2)
+
+#define i160i160OpMinusi160(in1,in2)			i16diffs(in1,in2)
+
 /* Matrix - Scalar */
 
 
@@ -99,6 +128,25 @@
   }
 
 
+#define u82u80OpMinusu82(in1,size,in2,out)		{int i=0;	\
+    for (i=0;i<size[0]*size[1];i++) out[i]=in1[i]-in2;		\
+  }
+
+
+#define i82i80OpMinusi82(in1,size,in2,out)		{int i=0;	\
+    for (i=0;i<size[0]*size[1];i++) out[i]=in1[i]-in2;		\
+  }
+
+
+#define u162u160OpMinusu162(in1,size,in2,out)		{int i=0;	\
+    for (i=0;i<size[0]*size[1];i++) out[i]=in1[i]-in2;		\
+  }
+
+
+#define i162i160OpMinusi162(in1,size,in2,out)		{int i=0;	\
+    for (i=0;i<size[0]*size[1];i++) out[i]=in1[i]-in2;		\
+  }
+
 /* Scalar - Matrix */
 
 
@@ -136,6 +184,22 @@
   }
 
 
+#define u80u82OpMinusu82(in1,in2,size,out)		{int i=0;	\
+    for (i=0;i<size[0]*size[1];i++) out[i]=in1-in2[i];		\
+  }
+
+
+#define i80i82OpMinusi82(in1,in2,size,out)		{int i=0;	\
+    for (i=0;i<size[0]*size[1];i++) out[i]=in1-in2[i];		\
+  }
+
+#define u160u162OpMinusu162(in1,in2,size,out)		{int i=0;	\
+    for (i=0;i<size[0]*size[1];i++) out[i]=in1-in2[i];		\
+  }
+
+#define i160i162OpMinusi162(in1,in2,size,out)		{int i=0;	\
+    for (i=0;i<size[0]*size[1];i++) out[i]=in1-in2[i];		\
+  }
 
 /* Matrix - Matrix */
 
@@ -167,7 +231,12 @@
     for (i=0;i<size1[0]*size2[1];i++) out[i]=zdiffs(in1[i],DoubleComplex(in2[i],0)); \
   }
 
+#define u82u82OpMinusu82(in1,size1,in2,size2,out)		u8diffa(in1, size1[0]*size1[1], in2, size2[0]*size2[1], out)
+
+#define i82i82OpMinusi82(in1,size1,in2,size2,out)		i8diffa(in1, size1[0]*size1[1], in2, size2[0]*size2[1], out)
 
+#define u162u162OpMinusu162(in1,size1,in2,size2,out)		u16diffa(in1, size1[0]*size1[1], in2, size2[0]*size2[1], out)
 
+#define i162i162OpMinusu162(in1,size1,in2,size2,out)		i16diffa(in1, size1[0]*size1[1], in2, size2[0]*size2[1], out)
 
 #endif /* !__INT_OPMINUS_H__ */
diff --git a/src/c/operations/interfaces/int_OpPlus.h b/src/c/operations/interfaces/int_OpPlus.h
index cb86d81..9684d25 100644
--- a/src/c/operations/interfaces/int_OpPlus.h
+++ b/src/c/operations/interfaces/int_OpPlus.h
@@ -34,7 +34,13 @@
 
 #define z0d0OpPlusz0(in1,in2)				zadds(in1,DoubleComplex(in2,0))
 
+#define u80u80OpPlusu80(in1,in2)			u8adds(in1,in2)
 
+#define i80i80OpPlusi80(in1,in2)			i8adds(in1,in2)
+
+#define u160u160OpPlusu160(in1,in2)			u16adds(in1,in2)
+
+#define i160i160OpPlusi160(in1,in2)			i16adds(in1,in2)
 
 /* Matrix + Scalar */
 
@@ -68,6 +74,21 @@
 								for (i=0;i<size[0]*size[1];i++) out[i]=zadds(DoubleComplex(in1[i],0),in2);\
 								}
 
+#define u82u80OpPlusu82(in1,size,in2,out)		{int i=0;\
+								for (i=0;i<size[0]*size[1];i++) out[i]=in1[i]+in2;\
+								}
+
+#define i82i80OpPlusi82(in1,size,in2,out)		{int i=0;\
+								for (i=0;i<size[0]*size[1];i++) out[i]=in1[i]+in2;\
+								}
+
+#define u162u160OpPlusu162(in1,size,in2,out)		{int i=0;\
+								for (i=0;i<size[0]*size[1];i++) out[i]=in1[i]+in2;\
+								}
+
+#define i162i160OpPlusi162(in1,size,in2,out)		{int i=0;\
+								for (i=0;i<size[0]*size[1];i++) out[i]=in1[i]+in2;\
+								}
 
 /* Scalar + Matrix */
 
@@ -101,7 +122,21 @@
 								for (i=0;i<size[0]*size[1];i++) out[i]=zadds(in1,DoubleComplex(in2[i],0));\
 								}
 								
+#define u80u82OpPlusu82(in1,in2,size,out)		{int i=0;\
+								for (i=0;i<size[0]*size[1];i++) out[i]=in1+in2[i];\
+								}
+
+#define i80i82OpPlusi82(in1,in2,size,out)		{int i=0;\
+								for (i=0;i<size[0]*size[1];i++) out[i]=in1+in2[i];\
+								}
+
+#define u160u162OpPlusu162(in1,in2,size,out)		{int i=0;\
+								for (i=0;i<size[0]*size[1];i++) out[i]=in1+in2[i];\
+								}
 
+#define i160i162OpPlusi162(in1,in2,size,out)		{int i=0;\
+								for (i=0;i<size[0]*size[1];i++) out[i]=in1+in2[i];\
+								}
 
 /* Matrix + Matrix */
 
@@ -132,7 +167,13 @@
 									for (i=0;i<size1[0]*size2[1];i++) out[i]=zadds(in1[i],DoubleComplex(in2[i],0));\
 									}
 
+#define u82u82OpPlusu82(in1,size1,in2,size2,out)	u8adda(in1, size1[0]*size1[1], in2, size2[0]*size2[1], out)
+
+#define i82i82OpPlusi82(in1,size1,in2,size2,out)	i8adda(in1, size1[0]*size1[1], in2, size2[0]*size2[1], out)
+
+#define u162u162OpPlusu162(in1,size1,in2,size2,out)	u16adda(in1, size1[0]*size1[1], in2, size2[0]*size2[1], out)
 
+#define i162i162OpPlusu162(in1,size1,in2,size2,out)	i16adda(in1, size1[0]*size1[1], in2, size2[0]*size2[1], out)
 /* Strings */
 /* RNU: size1[1]-1 -> "-1" because '\0' of the first string must be removed. */
 #define g2g2OpPlusg2(in1,size1,in2,size2,out)		{int i = 0, j = 0;	\
diff --git a/src/c/operations/subtraction/i16diffa.c b/src/c/operations/subtraction/i16diffa.c
new file mode 100644
index 0000000..ca722af
--- /dev/null
+++ b/src/c/operations/subtraction/i16diffa.c
@@ -0,0 +1,20 @@
+/*
+ *  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
+ *
+ */
+
+#include "subtraction.h"
+
+void i16diffa(int16* in1, int size1, int16* in2, int size2, int16* out) {
+  int i = 0;
+  for (i = 0; i < size1 && i < size2; ++i) {
+    out[i] = i16diffs(in1[i], in2[i]);
+  }
+}
diff --git a/src/c/operations/subtraction/i16diffs.c b/src/c/operations/subtraction/i16diffs.c
new file mode 100644
index 0000000..6ca2e52
--- /dev/null
+++ b/src/c/operations/subtraction/i16diffs.c
@@ -0,0 +1,18 @@
+/*
+ *  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
+ *
+ */
+
+#include "subtraction.h"
+
+int16 i16diffs(int16 in1, int16 in2)
+{
+  return (in1 - in2);
+}
diff --git a/src/c/operations/subtraction/i8diffa.c b/src/c/operations/subtraction/i8diffa.c
new file mode 100644
index 0000000..ec966bf
--- /dev/null
+++ b/src/c/operations/subtraction/i8diffa.c
@@ -0,0 +1,20 @@
+/*
+ *  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
+ *
+ */
+
+#include "subtraction.h"
+
+void i8diffa(int8* in1, int size1, int8* in2, int size2, int8* out) {
+  int i = 0;
+  for (i = 0; i < size1 && i < size2; ++i) {
+    out[i] = i8diffs(in1[i], in2[i]);
+  }
+}
diff --git a/src/c/operations/subtraction/i8diffs.c b/src/c/operations/subtraction/i8diffs.c
new file mode 100644
index 0000000..bbc9a1a
--- /dev/null
+++ b/src/c/operations/subtraction/i8diffs.c
@@ -0,0 +1,18 @@
+/*
+ *  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
+ *
+ */
+
+#include "subtraction.h"
+
+int8 i8diffs(int8 in1, int8 in2)
+{
+  return (in1 - in2);
+}
diff --git a/src/c/operations/subtraction/u16diffa.c b/src/c/operations/subtraction/u16diffa.c
new file mode 100644
index 0000000..7c216bb
--- /dev/null
+++ b/src/c/operations/subtraction/u16diffa.c
@@ -0,0 +1,20 @@
+/*
+ *  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
+ *
+ */
+
+#include "subtraction.h"
+
+void u16diffa(uint16* in1, int size1, uint16* in2, int size2, uint16* out) {
+  int i = 0;
+  for (i = 0; i < size1 && i < size2; ++i) {
+    out[i] = u16diffs(in1[i], in2[i]);
+  }
+}
diff --git a/src/c/operations/subtraction/u16diffs.c b/src/c/operations/subtraction/u16diffs.c
new file mode 100644
index 0000000..bd52d66
--- /dev/null
+++ b/src/c/operations/subtraction/u16diffs.c
@@ -0,0 +1,18 @@
+/*
+ *  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
+ *
+ */
+
+#include "subtraction.h"
+
+uint16 u16diffs(uint16 in1, uint16 in2)
+{
+  return (in1 - in2);
+}
diff --git a/src/c/operations/subtraction/u8diffa.c b/src/c/operations/subtraction/u8diffa.c
new file mode 100644
index 0000000..42c06f7
--- /dev/null
+++ b/src/c/operations/subtraction/u8diffa.c
@@ -0,0 +1,20 @@
+/*
+ *  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
+ *
+ */
+
+#include "subtraction.h"
+
+void u8diffa(uint8* in1, int size1, uint8* in2, int size2, uint8* out) {
+  int i = 0;
+  for (i = 0; i < size1 && i < size2; ++i) {
+    out[i] = u8diffs(in1[i], in2[i]);
+  }
+}
diff --git a/src/c/operations/subtraction/u8diffs.c b/src/c/operations/subtraction/u8diffs.c
new file mode 100644
index 0000000..bc9beac
--- /dev/null
+++ b/src/c/operations/subtraction/u8diffs.c
@@ -0,0 +1,18 @@
+/*
+ *  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
+ *
+ */
+
+#include "subtraction.h"
+
+uint8 u8diffs(uint8 in1, uint8 in2)
+{
+  return (in1 - in2);
+}