added stdevf functions , interfaces , header, tests , everything

19 files changed, 1384 insertions, 58 deletions

diff --git a/src/c/statisticsFunctions/includes/stdevf.h b/src/c/statisticsFunctions/includes/stdevf.h
index 5ccec182..c9e48599 100644
--- a/src/c/statisticsFunctions/includes/stdevf.h
+++ b/src/c/statisticsFunctions/includes/stdevf.h
@@ -111,6 +111,16 @@ EXTERN_STATFUNC doubleComplex	zstdevfa(doubleComplex *in1, int lines, int column
 EXTERN_STATFUNC void		zrowstdevfa(doubleComplex *in1, int lines, int columns, doubleComplex *in2, doubleComplex* out);
 EXTERN_STATFUNC void		zcolumnstdevfa(doubleComplex *in1, int lines, int columns, doubleComplex *in2, doubleComplex* out);
 
+/*
+**  convenience functions
+*/
+
+EXTERN_STATFUNC floatComplex	cstdevfcs(floatComplex *in1, int lines, int columns, float *in2);
+EXTERN_STATFUNC floatComplex	cstdevfsc(float *in1, int lines, int columns, floatComplex *in2);
+
+EXTERN_STATFUNC doubleComplex	zstdevfzd(doubleComplex *in1, int lines, int columns, double *in2);
+EXTERN_STATFUNC doubleComplex	zstdevfdz(double *in1, int lines, int columns, doubleComplex *in2);
+
 #ifdef  __cplusplus
 } /* extern "C" */
 #endif
diff --git a/src/c/statisticsFunctions/interfaces/int_stdevf.h b/src/c/statisticsFunctions/interfaces/int_stdevf.h
new file mode 100644
index 00000000..940b415d
--- /dev/null
+++ b/src/c/statisticsFunctions/interfaces/int_stdevf.h
@@ -0,0 +1,188 @@
+/*
+ * Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ * Copyright (C) 2008-2008 - INRIA - Allan SIMON
+ * 
+ * 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_STDEVF_H__
+#define __INT_STDEVF_H__
+
+#define s0s0stdevfs0(in1,in2)					sstdevfs(in1,in2)
+
+#define d0d0stdevfd0(in1,in2)					dstdevfs(in1,in2)
+
+#define c0s0stdevfc0(in1,in2)					cstdevfs(in1,FloatComplex(in2 , 0.0f))
+
+#define z0d0stdevfz0(in1,in2)					zstdevfs(in1,DoubleComplex(in2, 0.0))
+
+#define s2s2stdevfs0(in1, size1, in2, size2)				sstdevfa(in1, size1[0],size1[1], in2)
+
+#define d2d2stdevfd0(in1, size1, in2, size2)				dstdevfa(in1, size1[0],size1[1], in2)
+
+#define c2s2stdevfc0(in1, size1, in2, size2)				cstdevfcs(in1, size1[0], size1[1], in2)
+
+#define z2d2stdevfz0(in1, size1, in2, size2)				zstdevfzd(in1, size1[0], size1[1], in2)
+
+
+/*by rows or columns */
+
+#define s0s0s0stdevfs0(in1,in2,in3)				(in3==1) ? srowstdevfs(in1,in2) : scolumnstdevfs(in1,in2)	
+
+#define d0d0d0stdevfd0(in1,in2,in3)				(in3==1) ? drowstdevfs(in1,in2) : dcolumnstdevfs(in1,in2)	
+
+#define c0s0s0stdevfc0(in1,in2,in3)				(in3==1) ? crowstdevfs(in1,FloatComplex(in2 , 0.0f))  : ccolumnstdevfs(in1,FloatComplex(in2 , 0.0f))	
+
+#define z0d0d0stdevfz0(in1,in2,in3)				(in3==1) ? zrowstdevfs(in1,DoubleComplex(in2, 0.0))   : zcolumnstdevfs(in1,DoubleComplex(in2, 0.0))
+	
+
+#define s2s2s0stdevfs2(in1, size1, in2, size2,in3,out)		(in3==1) ? srowstdevfa(in1,size1[0],size1[1],in2,out) : scolumnstdevfa(in1,size1[0],size1[1],in2,out)	
+
+#define d2d2d0stdevfd2(in1, size1, in2, size2,in3,out)		(in3==1) ? drowstdevfa(in1,size1[0],size1[1],in2,out) : dcolumnstdevfa(in1,size1[0],size1[1],in2,out)
+
+
+
+
+
+#define c2s2s0stdevfc2(in1, size1, in2, size2,in3,out)		if (in3 == 1) {float* temp;\
+								 temp=malloc((uint)(size1[0]*size1[1])*sizeof(float));\
+								 sfilla(temp,size1[0],size1[1],0);\
+								 crowstdevfa(in1,size1[0],size1[1],FloatComplexMatrix(in2,temp,size1[0]*size1[1]),out);\
+								 free(temp);\
+								}\
+								else\
+								{float* temp;\
+								 temp=malloc((uint)(size1[0]*size1[1])*sizeof(float));\
+								 sfilla(temp,size1[0],size1[1],0);\
+								 ccolumnstdevfa(in1,size1[0],size1[1],FloatComplexMatrix(in2,temp,size1[0]*size1[1]),out);\
+								 free(temp);\
+								}	
+
+#define z2d2d0stdevfz2(in1, size1, in2, size2,in3,out)		if (in3 == 1) {double* temp;\
+								temp=malloc((uint)(size1[0]*size1[1])*sizeof(double));\
+								dfilla(temp,size1[0],size1[1],0);\
+								zrowstdevfa(in1,size1[0],size1[1],DoubleComplexMatrix(in2,temp,size1[0]*size1[1]),out);\
+								free(temp);\
+								}\
+								else\
+								{double* temp;\
+								temp=malloc((uint)(size1[0]*size1[1])*sizeof(double));\
+								dfilla(temp,size1[0],size1[1],0);\
+								zcolumnstdevfa(in1,size1[0],size1[1],DoubleComplexMatrix(in2,temp,size1[0]*size1[1]),out);\
+								free(temp);\
+								}	
+
+		
+/**/
+
+
+#define s0c0stdevfs0(in1,in2)					creals (cstdevfs(FloatComplex(in1 , 0.0f),in2))
+#define s0c0stdevfc0(in1,in2)					cstdevfs(FloatComplex(in1 , 0.0f),in2)
+
+#define d0z0stdevfd0(in1,in2)					zreals(zstdevfs(DoubleComplex(in1, 0.0 ),in2)
+#define d0z0stdevfz0(in1,in2)					zstdevfs(DoubleComplex(in1, 0.0 ),in2)
+
+
+#define c0c0stdevfc0(in1,in2)					cstdevfs(in1,in2)
+
+#define z0z0stdevfz0(in1,in2)					zstdevfs(in1,in2)
+
+
+#define s2c2stdevfs0(in1, size1, in2, size2)				creals(cstdevfsc(in1, size1[0],size1[1], in2))
+#define s2c2stdevfc0(in1, size1, in2, size2)				cstdevfsc(in1, size1[0],size1[1], in2)
+
+#define d2z2stdevfd0(in1, size1, in2, size2)				zreals(zstdevfdz(in1, size1[0],size1[1], in2))
+#define d2z2stdevfz0(in1, size1, in2, size2)				zstdevfdz(in1, size1[0],size1[1], in2)
+
+#define c2c2stdevfc0(in1, size1, in2, size2)				cstdevfa(in1, size1[0], size1[1], in2)
+
+#define z2z2stdevfz0(in1, size1, in2, size2)				zstdevfa(in1, size1[0], size1[1], in2)
+
+
+/*by rows or columns */
+
+#define s0c0s0stdevfs0(in1,in2,in3)				(in3==1) ? creals(crowstdevfs(FloatComplex(in1, 0.0f),in2)) : creals(ccolumnstdevfs(FloatComplex(in1, 0.0f),in2))	
+
+#define d0z0d0stdevfd0(in1,in2,in3)				(in3==1) ? zreals(zrowstdevfs(DoubleComplex(in1, 0.0),in2)) : zreals(zcolumnstdevfs(DoubleComplex(in1, 0.0),in2))	
+
+#define c0c0s0stdevfc0(in1,in2,in3)				(in3==1) ? crowstdevfs(in1,in2) : ccolumnstdevfs(in1,in2)	
+
+#define z0z0d0stdevfz0(in1,in2,in3)				(in3==1) ? zrowstdevfs(in1,in2) : zcolumnstdevfs(in1,in2)
+	
+
+#define s2c2s0stdevfs2(in1, size1, in2, size2,in3,out)		if (in3 == 1) {float* temp;\
+								temp=malloc((uint)(size1[0]*size1[1])*sizeof(float));\
+								sfilla(temp,size1[0],size1[1],0);\
+								crowstdevfa(FloatComplexMatrix(in1,temp,size1[0]*size1[1]),size1[0],size1[1],in2,FloatComplexMatrix(out,temp,size1[0]*size1[1]));\
+								free(temp);\
+								}\
+								else\
+								{float* temp;\
+								temp=malloc((uint)(size1[0]*size1[1])*sizeof(float));\
+								sfilla(temp,size1[0],size1[1],0);\
+								ccolumnstdevfa(FloatComplexMatrix(in1,temp,size1[0]*size1[1]),size1[0],size1[1],in2,FloatComplexMatrix(out,temp,size1[0]*size1[1]));\
+								free(temp);\
+								}	
+
+#define d2z2d0stdevfd2(in1, size1, in2, size2,in3,out)		if (in3 == 1) {double* temp;\
+								temp=malloc((uint)(size1[0]*size1[1])*sizeof(double));\
+								dfilla(temp,size1[0],size1[1],0);\
+								zrowstdevfa(DoubleComplexMatrix(in1,temp,size1[0]*size1[1]),size1[0],size1[1],in2,DoubleComplexMatrix(out,temp,size1[0]));\
+								free(temp);\
+								}\
+								else\
+								{double* temp;\
+								temp=malloc((uint)(size1[0]*size1[1])*sizeof(double));\
+								dfilla(temp,size1[0],size1[1],0);\
+								zcolumnstdevfa(DoubleComplexMatrix(in1,temp,size1[0]*size1[1]),size1[0],size1[1],in2,DoubleComplexMatrix(out,temp,size1[0]));\
+								free(temp);\
+								}	
+
+
+/*complex output*/
+
+#define s2c2s0stdevfc2(in1, size1, in2, size2,in3,out)		if (in3 == 1) {float* temp;\
+								temp=malloc((uint)(size1[0]*size1[1])*sizeof(float));\
+								sfilla(temp,size1[0],size1[1],0);\
+								crowstdevfa(FloatComplexMatrix(in1,temp,size1[0]*size1[1]),size1[0],size1[1],in2,out);\
+								free(temp);\
+								}\
+								else\
+								{float* temp;\
+								temp=malloc((uint)(size1[0]*size1[1])*sizeof(float));\
+								sfilla(temp,size1[0],size1[1],0);\
+								ccolumnstdevfa(FloatComplexMatrix(in1,temp,size1[0]*size1[1]),size1[0],size1[1],in2,out);\
+								free(temp);\
+								}	
+
+#define d2z2d0stdevfz2(in1, size1, in2, size2,in3,out)		if (in3 == 1) {double* temp;\
+								temp=malloc((uint)(size1[0]*size1[1])*sizeof(double));\
+								dfilla(temp,size1[0],size1[1],0);\
+								zrowstdevfa(DoubleComplexMatrix(in1,temp,size1[0]*size1[1]),size1[0],size1[1],in2,out);\
+								free(temp);\
+								}\
+								else\
+								{double* temp;\
+								temp=malloc((uint)(size1[0]*size1[1])*sizeof(double));\
+								dfilla(temp,size1[0],size1[1],0);\
+								zcolumnstdevfa(DoubleComplexMatrix(in1,temp,size1[0]*size1[1]),size1[0],size1[1],in2,out);\
+								free(temp);\
+								}	
+
+
+
+#define c2c2s0stdevfc2(in1, size1, in2, size2,in3,out)		(in3==1) ? crowstdevfa(in1,size1[0],size1[1],in2,out) : ccolumnstdevfa(in1,size1[0],size1[1],in2,out)
+
+#define z2z2d0stdevfz2(in1, size1, in2, size2,in3,out)		(in3==1) ? zrowstdevfa(in1,size1[0],size1[1],in2,out) : zcolumnstdevfa(in1,size1[0],size1[1],in2,out)	
+
+
+/**/
+
+#endif /* !__INT_STDEVF_H__ */
diff --git a/src/c/statisticsFunctions/stdevf/Makefile.in b/src/c/statisticsFunctions/stdevf/Makefile.in
new file mode 100644
index 00000000..37f04b70
--- /dev/null
+++ b/src/c/statisticsFunctions/stdevf/Makefile.in
@@ -0,0 +1,818 @@
+# Makefile.in generated by automake 1.10.2 from Makefile.am.
+# @configure_input@
+
+# Copyright (C) 1994, 1995, 1996, 1997, 1998, 1999, 2000, 2001, 2002,
+# 2003, 2004, 2005, 2006, 2007, 2008  Free Software Foundation, Inc.
+# This Makefile.in is free software; the Free Software Foundation
+# gives unlimited permission to copy and/or distribute it,
+# with or without modifications, as long as this notice is preserved.
+
+# This program is distributed in the hope that it will be useful,
+# but WITHOUT ANY WARRANTY, to the extent permitted by law; without
+# even the implied warranty of MERCHANTABILITY or FITNESS FOR A
+# PARTICULAR PURPOSE.
+
+@SET_MAKE@
+
+VPATH = @srcdir@
+pkgdatadir = $(datadir)/@PACKAGE@
+pkglibdir = $(libdir)/@PACKAGE@
+pkgincludedir = $(includedir)/@PACKAGE@
+am__cd = CDPATH="$${ZSH_VERSION+.}$(PATH_SEPARATOR)" && cd
+install_sh_DATA = $(install_sh) -c -m 644
+install_sh_PROGRAM = $(install_sh) -c
+install_sh_SCRIPT = $(install_sh) -c
+INSTALL_HEADER = $(INSTALL_DATA)
+transform = $(program_transform_name)
+NORMAL_INSTALL = :
+PRE_INSTALL = :
+POST_INSTALL = :
+NORMAL_UNINSTALL = :
+PRE_UNINSTALL = :
+POST_UNINSTALL = :
+build_triplet = @build@
+host_triplet = @host@
+check_PROGRAMS = testFloatStdevf$(EXEEXT) testDoubleStdevf$(EXEEXT)
+TESTS = testFloatStdevf$(EXEEXT) testDoubleStdevf$(EXEEXT)
+subdir = src/c/statisticsFunctions/stdevf
+DIST_COMMON = $(srcdir)/Makefile.am $(srcdir)/Makefile.in
+ACLOCAL_M4 = $(top_srcdir)/aclocal.m4
+am__aclocal_m4_deps = $(top_srcdir)/configure.ac
+am__configure_deps = $(am__aclocal_m4_deps) $(CONFIGURE_DEPENDENCIES) \
+	$(ACLOCAL_M4)
+mkinstalldirs = $(install_sh) -d
+CONFIG_HEADER = $(top_builddir)/includes/machine.h
+CONFIG_CLEAN_FILES =
+am__vpath_adj_setup = srcdirstrip=`echo "$(srcdir)" | sed 's|.|.|g'`;
+am__vpath_adj = case $$p in \
+    $(srcdir)/*) f=`echo "$$p" | sed "s|^$$srcdirstrip/||"`;; \
+    *) f=$$p;; \
+  esac;
+am__strip_dir = `echo $$p | sed -e 's|^.*/||'`;
+am__installdirs = "$(DESTDIR)$(pkglibdir)"
+pkglibLTLIBRARIES_INSTALL = $(INSTALL)
+LTLIBRARIES = $(pkglib_LTLIBRARIES)
+libStdevf_la_LIBADD =
+am__objects_1 =
+am__objects_2 = libStdevf_la-sstdevfa.lo libStdevf_la-srowstdevfa.lo \
+	libStdevf_la-scolumnstdevfa.lo libStdevf_la-dstdevfa.lo \
+	libStdevf_la-drowstdevfa.lo libStdevf_la-dcolumnstdevfa.lo \
+	libStdevf_la-cstdevfa.lo libStdevf_la-crowstdevfa.lo \
+	libStdevf_la-ccolumnstdevfa.lo libStdevf_la-zstdevfa.lo \
+	libStdevf_la-zrowstdevfa.lo libStdevf_la-zcolumnstdevfa.lo
+am_libStdevf_la_OBJECTS = $(am__objects_1) $(am__objects_2)
+libStdevf_la_OBJECTS = $(am_libStdevf_la_OBJECTS)
+libStdevf_la_LINK = $(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) \
+	$(LIBTOOLFLAGS) --mode=link $(CCLD) $(libStdevf_la_CFLAGS) \
+	$(CFLAGS) $(AM_LDFLAGS) $(LDFLAGS) -o $@
+am_testDoubleStdevf_OBJECTS =  \
+	testDoubleStdevf-testDoubleStdevf.$(OBJEXT)
+testDoubleStdevf_OBJECTS = $(am_testDoubleStdevf_OBJECTS)
+am__DEPENDENCIES_1 = $(top_builddir)/src/c/type/libDoubleComplex.la \
+	$(top_builddir)/src/c/type/libFloatComplex.la \
+	$(top_builddir)/src/c/operations/addition/libAddition.la \
+	$(top_builddir)/src/fortran/lapack/libscilapack.la \
+	$(top_builddir)/src/c/operations/subtraction/libSubtraction.la \
+	$(top_builddir)/src/c/statisticsFunctions/stdevf/libStdevf.la \
+	$(top_builddir)/src/c/statisticsFunctions/sum/libSum.la \
+	$(top_builddir)/src/c/statisticsFunctions/meanf/libMeanf.la \
+	$(top_builddir)/src/c/operations/division/libDivision.la \
+	$(top_builddir)/src/c/operations/multiplication/libMultiplication.la \
+	$(top_builddir)/src/c/auxiliaryFunctions/conj/libConj.la \
+	$(top_builddir)/src/c/elementaryFunctions/pow/libPow.la \
+	$(top_builddir)/src/c/elementaryFunctions/exp/libExp.la \
+	$(top_builddir)/src/c/elementaryFunctions/sin/libSin.la \
+	$(top_builddir)/src/c/elementaryFunctions/cos/libCos.la \
+	$(top_builddir)/src/c/elementaryFunctions/cosh/libCosh.la \
+	$(top_builddir)/src/c/elementaryFunctions/log/libLog.la \
+	$(top_builddir)/src/c/elementaryFunctions/log1p/libLog1p.la \
+	$(top_builddir)/src/c/elementaryFunctions/lnp1m1/libLnp1m1.la \
+	$(top_builddir)/src/c/auxiliaryFunctions/pythag/libPythag.la \
+	$(top_builddir)/src/c/auxiliaryFunctions/abs/libAbs.la \
+	$(top_builddir)/src/c/elementaryFunctions/sinh/libSinh.la \
+	$(top_builddir)/src/c/elementaryFunctions/sqrt/libSqrt.la
+testDoubleStdevf_DEPENDENCIES = $(am__DEPENDENCIES_1)
+testDoubleStdevf_LINK = $(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) \
+	$(LIBTOOLFLAGS) --mode=link $(CCLD) $(testDoubleStdevf_CFLAGS) \
+	$(CFLAGS) $(AM_LDFLAGS) $(LDFLAGS) -o $@
+am_testFloatStdevf_OBJECTS =  \
+	testFloatStdevf-testFloatStdevf.$(OBJEXT)
+testFloatStdevf_OBJECTS = $(am_testFloatStdevf_OBJECTS)
+testFloatStdevf_DEPENDENCIES = $(am__DEPENDENCIES_1)
+testFloatStdevf_LINK = $(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) \
+	$(LIBTOOLFLAGS) --mode=link $(CCLD) $(testFloatStdevf_CFLAGS) \
+	$(CFLAGS) $(AM_LDFLAGS) $(LDFLAGS) -o $@
+DEFAULT_INCLUDES = -I.@am__isrc@ -I$(top_builddir)/includes
+depcomp = $(SHELL) $(top_srcdir)/config/depcomp
+am__depfiles_maybe = depfiles
+COMPILE = $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) \
+	$(CPPFLAGS) $(AM_CFLAGS) $(CFLAGS)
+LTCOMPILE = $(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) \
+	--mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) \
+	$(AM_CPPFLAGS) $(CPPFLAGS) $(AM_CFLAGS) $(CFLAGS)
+CCLD = $(CC)
+LINK = $(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) \
+	--mode=link $(CCLD) $(AM_CFLAGS) $(CFLAGS) $(AM_LDFLAGS) \
+	$(LDFLAGS) -o $@
+SOURCES = $(libStdevf_la_SOURCES) $(testDoubleStdevf_SOURCES) \
+	$(testFloatStdevf_SOURCES)
+DIST_SOURCES = $(libStdevf_la_SOURCES) $(testDoubleStdevf_SOURCES) \
+	$(testFloatStdevf_SOURCES)
+ETAGS = etags
+CTAGS = ctags
+DISTFILES = $(DIST_COMMON) $(DIST_SOURCES) $(TEXINFOS) $(EXTRA_DIST)
+ACLOCAL = @ACLOCAL@
+AMTAR = @AMTAR@
+AR = @AR@
+AUTOCONF = @AUTOCONF@
+AUTOHEADER = @AUTOHEADER@
+AUTOMAKE = @AUTOMAKE@
+AWK = @AWK@
+CC = @CC@
+CCDEPMODE = @CCDEPMODE@
+CFLAGS = @CFLAGS@
+CPP = @CPP@
+CPPFLAGS = @CPPFLAGS@
+CXX = @CXX@
+CXXCPP = @CXXCPP@
+CXXDEPMODE = @CXXDEPMODE@
+CXXFLAGS = @CXXFLAGS@
+CYGPATH_W = @CYGPATH_W@
+DEFS = @DEFS@
+DEPDIR = @DEPDIR@
+DSYMUTIL = @DSYMUTIL@
+DUMPBIN = @DUMPBIN@
+ECHO_C = @ECHO_C@
+ECHO_N = @ECHO_N@
+ECHO_T = @ECHO_T@
+EGREP = @EGREP@
+EXEEXT = @EXEEXT@
+F77 = @F77@
+FFLAGS = @FFLAGS@
+FGREP = @FGREP@
+GREP = @GREP@
+INSTALL = @INSTALL@
+INSTALL_DATA = @INSTALL_DATA@
+INSTALL_PROGRAM = @INSTALL_PROGRAM@
+INSTALL_SCRIPT = @INSTALL_SCRIPT@
+INSTALL_STRIP_PROGRAM = @INSTALL_STRIP_PROGRAM@
+LD = @LD@
+LDFLAGS = @LDFLAGS@
+LIBMATH = @LIBMATH@
+LIBOBJS = @LIBOBJS@
+LIBS = @LIBS@
+LIBTOOL = @LIBTOOL@
+LIPO = @LIPO@
+LN_S = @LN_S@
+LTLIBOBJS = @LTLIBOBJS@
+MAINT = @MAINT@
+MAKEINFO = @MAKEINFO@
+MKDIR_P = @MKDIR_P@
+NM = @NM@
+NMEDIT = @NMEDIT@
+OBJDUMP = @OBJDUMP@
+OBJEXT = @OBJEXT@
+OTOOL = @OTOOL@
+OTOOL64 = @OTOOL64@
+PACKAGE = @PACKAGE@
+PACKAGE_BUGREPORT = @PACKAGE_BUGREPORT@
+PACKAGE_NAME = @PACKAGE_NAME@
+PACKAGE_STRING = @PACKAGE_STRING@
+PACKAGE_TARNAME = @PACKAGE_TARNAME@
+PACKAGE_VERSION = @PACKAGE_VERSION@
+PATH_SEPARATOR = @PATH_SEPARATOR@
+RANLIB = @RANLIB@
+SED = @SED@
+SET_MAKE = @SET_MAKE@
+SHELL = @SHELL@
+STRIP = @STRIP@
+VERSION = @VERSION@
+abs_builddir = @abs_builddir@
+abs_srcdir = @abs_srcdir@
+abs_top_builddir = @abs_top_builddir@
+abs_top_srcdir = @abs_top_srcdir@
+ac_ct_CC = @ac_ct_CC@
+ac_ct_CXX = @ac_ct_CXX@
+ac_ct_DUMPBIN = @ac_ct_DUMPBIN@
+ac_ct_F77 = @ac_ct_F77@
+am__include = @am__include@
+am__leading_dot = @am__leading_dot@
+am__quote = @am__quote@
+am__tar = @am__tar@
+am__untar = @am__untar@
+bindir = @bindir@
+build = @build@
+build_alias = @build_alias@
+build_cpu = @build_cpu@
+build_os = @build_os@
+build_vendor = @build_vendor@
+builddir = @builddir@
+datadir = @datadir@
+datarootdir = @datarootdir@
+docdir = @docdir@
+dvidir = @dvidir@
+exec_prefix = @exec_prefix@
+host = @host@
+host_alias = @host_alias@
+host_cpu = @host_cpu@
+host_os = @host_os@
+host_vendor = @host_vendor@
+htmldir = @htmldir@
+includedir = @includedir@
+infodir = @infodir@
+install_sh = @install_sh@
+libdir = @libdir@
+libexecdir = @libexecdir@
+localedir = @localedir@
+localstatedir = @localstatedir@
+lt_ECHO = @lt_ECHO@
+mandir = @mandir@
+mkdir_p = @mkdir_p@
+oldincludedir = @oldincludedir@
+pdfdir = @pdfdir@
+prefix = @prefix@
+program_transform_name = @program_transform_name@
+psdir = @psdir@
+sbindir = @sbindir@
+sharedstatedir = @sharedstatedir@
+srcdir = @srcdir@
+sysconfdir = @sysconfdir@
+target_alias = @target_alias@
+top_build_prefix = @top_build_prefix@
+top_builddir = @top_builddir@
+top_srcdir = @top_srcdir@
+libStdevf_la_CFLAGS = -I $(top_builddir)/src/c/type \
+			-I $(top_builddir)/src/c/operations/includes \
+			-I $(top_builddir)/src/c/auxiliaryFunctions/includes \
+			-I $(top_builddir)/src/c/elementaryFunctions/includes \
+			-I ../includes
+
+instdir = $(top_builddir)/lib
+pkglib_LTLIBRARIES = libStdevf.la
+libStdevf_la_SOURCES = $(HEAD) $(SRC)
+SRC = sstdevfa.c \
+		srowstdevfa.c \
+		scolumnstdevfa.c \
+		dstdevfa.c \
+		drowstdevfa.c \
+		dcolumnstdevfa.c \
+		cstdevfa.c \
+		crowstdevfa.c \
+		ccolumnstdevfa.c \
+		zstdevfa.c \
+		zrowstdevfa.c \
+		zcolumnstdevfa.c
+
+HEAD = ../includes/stdevf.h
+
+####
+# Checking Part
+####
+check_INCLUDES = -I $(top_builddir)/src/c/statisticsFunctions/includes \
+			-I $(top_builddir)/src/c/operations/includes \
+			-I $(top_builddir)/src/c/auxiliaryFunctions/includes \
+			-I $(top_builddir)/src/c/elementaryFunctions/includes \
+			-I $(top_builddir)/src/c/type
+
+check_LDADD = $(top_builddir)/src/c/type/libDoubleComplex.la \
+			$(top_builddir)/src/c/type/libFloatComplex.la \
+			$(top_builddir)/src/c/operations/addition/libAddition.la \
+			$(top_builddir)/src/fortran/lapack/libscilapack.la \
+			$(top_builddir)/src/c/operations/subtraction/libSubtraction.la \
+			$(top_builddir)/src/c/statisticsFunctions/stdevf/libStdevf.la \
+			$(top_builddir)/src/c/statisticsFunctions/sum/libSum.la \
+			$(top_builddir)/src/c/statisticsFunctions/meanf/libMeanf.la \
+			$(top_builddir)/src/c/operations/division/libDivision.la \
+			$(top_builddir)/src/c/operations/multiplication/libMultiplication.la \
+			$(top_builddir)/src/c/auxiliaryFunctions/conj/libConj.la \
+			$(top_builddir)/src/c/elementaryFunctions/pow/libPow.la \
+			$(top_builddir)/src/c/elementaryFunctions/exp/libExp.la \
+			$(top_builddir)/src/c/elementaryFunctions/sin/libSin.la \
+			$(top_builddir)/src/c/elementaryFunctions/cos/libCos.la \
+			$(top_builddir)/src/c/elementaryFunctions/cosh/libCosh.la \
+			$(top_builddir)/src/c/elementaryFunctions/log/libLog.la \
+			$(top_builddir)/src/c/elementaryFunctions/log1p/libLog1p.la \
+			$(top_builddir)/src/c/elementaryFunctions/lnp1m1/libLnp1m1.la \
+			$(top_builddir)/src/c/auxiliaryFunctions/pythag/libPythag.la \
+			$(top_builddir)/src/c/auxiliaryFunctions/abs/libAbs.la \
+			$(top_builddir)/src/c/elementaryFunctions/sinh/libSinh.la \
+			$(top_builddir)/src/c/elementaryFunctions/sqrt/libSqrt.la \
+			@LIBMATH@
+
+
+#
+# -*- Stdevf Tests -*-
+#
+testFloatStdevf_SOURCES = testFloatStdevf.c
+testFloatStdevf_CFLAGS = $(check_INCLUDES)
+testFloatStdevf_LDADD = $(check_LDADD)
+testDoubleStdevf_SOURCES = testDoubleStdevf.c
+testDoubleStdevf_CFLAGS = $(check_INCLUDES)
+testDoubleStdevf_LDADD = $(check_LDADD)
+all: all-am
+
+.SUFFIXES:
+.SUFFIXES: .c .lo .o .obj
+$(srcdir)/Makefile.in: @MAINTAINER_MODE_TRUE@ $(srcdir)/Makefile.am  $(am__configure_deps)
+	@for dep in $?; do \
+	  case '$(am__configure_deps)' in \
+	    *$$dep*) \
+	      ( cd $(top_builddir) && $(MAKE) $(AM_MAKEFLAGS) am--refresh ) \
+	        && { if test -f $@; then exit 0; else break; fi; }; \
+	      exit 1;; \
+	  esac; \
+	done; \
+	echo ' cd $(top_srcdir) && $(AUTOMAKE) --foreign  src/c/statisticsFunctions/stdevf/Makefile'; \
+	cd $(top_srcdir) && \
+	  $(AUTOMAKE) --foreign  src/c/statisticsFunctions/stdevf/Makefile
+.PRECIOUS: Makefile
+Makefile: $(srcdir)/Makefile.in $(top_builddir)/config.status
+	@case '$?' in \
+	  *config.status*) \
+	    cd $(top_builddir) && $(MAKE) $(AM_MAKEFLAGS) am--refresh;; \
+	  *) \
+	    echo ' cd $(top_builddir) && $(SHELL) ./config.status $(subdir)/$@ $(am__depfiles_maybe)'; \
+	    cd $(top_builddir) && $(SHELL) ./config.status $(subdir)/$@ $(am__depfiles_maybe);; \
+	esac;
+
+$(top_builddir)/config.status: $(top_srcdir)/configure $(CONFIG_STATUS_DEPENDENCIES)
+	cd $(top_builddir) && $(MAKE) $(AM_MAKEFLAGS) am--refresh
+
+$(top_srcdir)/configure: @MAINTAINER_MODE_TRUE@ $(am__configure_deps)
+	cd $(top_builddir) && $(MAKE) $(AM_MAKEFLAGS) am--refresh
+$(ACLOCAL_M4): @MAINTAINER_MODE_TRUE@ $(am__aclocal_m4_deps)
+	cd $(top_builddir) && $(MAKE) $(AM_MAKEFLAGS) am--refresh
+install-pkglibLTLIBRARIES: $(pkglib_LTLIBRARIES)
+	@$(NORMAL_INSTALL)
+	test -z "$(pkglibdir)" || $(MKDIR_P) "$(DESTDIR)$(pkglibdir)"
+	@list='$(pkglib_LTLIBRARIES)'; for p in $$list; do \
+	  if test -f $$p; then \
+	    f=$(am__strip_dir) \
+	    echo " $(LIBTOOL) $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=install $(pkglibLTLIBRARIES_INSTALL) $(INSTALL_STRIP_FLAG) '$$p' '$(DESTDIR)$(pkglibdir)/$$f'"; \
+	    $(LIBTOOL) $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=install $(pkglibLTLIBRARIES_INSTALL) $(INSTALL_STRIP_FLAG) "$$p" "$(DESTDIR)$(pkglibdir)/$$f"; \
+	  else :; fi; \
+	done
+
+uninstall-pkglibLTLIBRARIES:
+	@$(NORMAL_UNINSTALL)
+	@list='$(pkglib_LTLIBRARIES)'; for p in $$list; do \
+	  p=$(am__strip_dir) \
+	  echo " $(LIBTOOL) $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=uninstall rm -f '$(DESTDIR)$(pkglibdir)/$$p'"; \
+	  $(LIBTOOL) $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=uninstall rm -f "$(DESTDIR)$(pkglibdir)/$$p"; \
+	done
+
+clean-pkglibLTLIBRARIES:
+	-test -z "$(pkglib_LTLIBRARIES)" || rm -f $(pkglib_LTLIBRARIES)
+	@list='$(pkglib_LTLIBRARIES)'; for p in $$list; do \
+	  dir="`echo $$p | sed -e 's|/[^/]*$$||'`"; \
+	  test "$$dir" != "$$p" || dir=.; \
+	  echo "rm -f \"$${dir}/so_locations\""; \
+	  rm -f "$${dir}/so_locations"; \
+	done
+libStdevf.la: $(libStdevf_la_OBJECTS) $(libStdevf_la_DEPENDENCIES) 
+	$(libStdevf_la_LINK) -rpath $(pkglibdir) $(libStdevf_la_OBJECTS) $(libStdevf_la_LIBADD) $(LIBS)
+
+clean-checkPROGRAMS:
+	@list='$(check_PROGRAMS)'; for p in $$list; do \
+	  f=`echo $$p|sed 's/$(EXEEXT)$$//'`; \
+	  echo " rm -f $$p $$f"; \
+	  rm -f $$p $$f ; \
+	done
+testDoubleStdevf$(EXEEXT): $(testDoubleStdevf_OBJECTS) $(testDoubleStdevf_DEPENDENCIES) 
+	@rm -f testDoubleStdevf$(EXEEXT)
+	$(testDoubleStdevf_LINK) $(testDoubleStdevf_OBJECTS) $(testDoubleStdevf_LDADD) $(LIBS)
+testFloatStdevf$(EXEEXT): $(testFloatStdevf_OBJECTS) $(testFloatStdevf_DEPENDENCIES) 
+	@rm -f testFloatStdevf$(EXEEXT)
+	$(testFloatStdevf_LINK) $(testFloatStdevf_OBJECTS) $(testFloatStdevf_LDADD) $(LIBS)
+
+mostlyclean-compile:
+	-rm -f *.$(OBJEXT)
+
+distclean-compile:
+	-rm -f *.tab.c
+
+@AMDEP_TRUE@@am__include@ @am__quote@./$(DEPDIR)/libStdevf_la-ccolumnstdevfa.Plo@am__quote@
+@AMDEP_TRUE@@am__include@ @am__quote@./$(DEPDIR)/libStdevf_la-crowstdevfa.Plo@am__quote@
+@AMDEP_TRUE@@am__include@ @am__quote@./$(DEPDIR)/libStdevf_la-cstdevfa.Plo@am__quote@
+@AMDEP_TRUE@@am__include@ @am__quote@./$(DEPDIR)/libStdevf_la-dcolumnstdevfa.Plo@am__quote@
+@AMDEP_TRUE@@am__include@ @am__quote@./$(DEPDIR)/libStdevf_la-drowstdevfa.Plo@am__quote@
+@AMDEP_TRUE@@am__include@ @am__quote@./$(DEPDIR)/libStdevf_la-dstdevfa.Plo@am__quote@
+@AMDEP_TRUE@@am__include@ @am__quote@./$(DEPDIR)/libStdevf_la-scolumnstdevfa.Plo@am__quote@
+@AMDEP_TRUE@@am__include@ @am__quote@./$(DEPDIR)/libStdevf_la-srowstdevfa.Plo@am__quote@
+@AMDEP_TRUE@@am__include@ @am__quote@./$(DEPDIR)/libStdevf_la-sstdevfa.Plo@am__quote@
+@AMDEP_TRUE@@am__include@ @am__quote@./$(DEPDIR)/libStdevf_la-zcolumnstdevfa.Plo@am__quote@
+@AMDEP_TRUE@@am__include@ @am__quote@./$(DEPDIR)/libStdevf_la-zrowstdevfa.Plo@am__quote@
+@AMDEP_TRUE@@am__include@ @am__quote@./$(DEPDIR)/libStdevf_la-zstdevfa.Plo@am__quote@
+@AMDEP_TRUE@@am__include@ @am__quote@./$(DEPDIR)/testDoubleStdevf-testDoubleStdevf.Po@am__quote@
+@AMDEP_TRUE@@am__include@ @am__quote@./$(DEPDIR)/testFloatStdevf-testFloatStdevf.Po@am__quote@
+
+.c.o:
+@am__fastdepCC_TRUE@	$(COMPILE) -MT $@ -MD -MP -MF $(DEPDIR)/$*.Tpo -c -o $@ $<
+@am__fastdepCC_TRUE@	mv -f $(DEPDIR)/$*.Tpo $(DEPDIR)/$*.Po
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	source='$<' object='$@' libtool=no @AMDEPBACKSLASH@
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	DEPDIR=$(DEPDIR) $(CCDEPMODE) $(depcomp) @AMDEPBACKSLASH@
+@am__fastdepCC_FALSE@	$(COMPILE) -c $<
+
+.c.obj:
+@am__fastdepCC_TRUE@	$(COMPILE) -MT $@ -MD -MP -MF $(DEPDIR)/$*.Tpo -c -o $@ `$(CYGPATH_W) '$<'`
+@am__fastdepCC_TRUE@	mv -f $(DEPDIR)/$*.Tpo $(DEPDIR)/$*.Po
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	source='$<' object='$@' libtool=no @AMDEPBACKSLASH@
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	DEPDIR=$(DEPDIR) $(CCDEPMODE) $(depcomp) @AMDEPBACKSLASH@
+@am__fastdepCC_FALSE@	$(COMPILE) -c `$(CYGPATH_W) '$<'`
+
+.c.lo:
+@am__fastdepCC_TRUE@	$(LTCOMPILE) -MT $@ -MD -MP -MF $(DEPDIR)/$*.Tpo -c -o $@ $<
+@am__fastdepCC_TRUE@	mv -f $(DEPDIR)/$*.Tpo $(DEPDIR)/$*.Plo
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	source='$<' object='$@' libtool=yes @AMDEPBACKSLASH@
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	DEPDIR=$(DEPDIR) $(CCDEPMODE) $(depcomp) @AMDEPBACKSLASH@
+@am__fastdepCC_FALSE@	$(LTCOMPILE) -c -o $@ $<
+
+libStdevf_la-sstdevfa.lo: sstdevfa.c
+@am__fastdepCC_TRUE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -MT libStdevf_la-sstdevfa.lo -MD -MP -MF $(DEPDIR)/libStdevf_la-sstdevfa.Tpo -c -o libStdevf_la-sstdevfa.lo `test -f 'sstdevfa.c' || echo '$(srcdir)/'`sstdevfa.c
+@am__fastdepCC_TRUE@	mv -f $(DEPDIR)/libStdevf_la-sstdevfa.Tpo $(DEPDIR)/libStdevf_la-sstdevfa.Plo
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	source='sstdevfa.c' object='libStdevf_la-sstdevfa.lo' libtool=yes @AMDEPBACKSLASH@
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	DEPDIR=$(DEPDIR) $(CCDEPMODE) $(depcomp) @AMDEPBACKSLASH@
+@am__fastdepCC_FALSE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -c -o libStdevf_la-sstdevfa.lo `test -f 'sstdevfa.c' || echo '$(srcdir)/'`sstdevfa.c
+
+libStdevf_la-srowstdevfa.lo: srowstdevfa.c
+@am__fastdepCC_TRUE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -MT libStdevf_la-srowstdevfa.lo -MD -MP -MF $(DEPDIR)/libStdevf_la-srowstdevfa.Tpo -c -o libStdevf_la-srowstdevfa.lo `test -f 'srowstdevfa.c' || echo '$(srcdir)/'`srowstdevfa.c
+@am__fastdepCC_TRUE@	mv -f $(DEPDIR)/libStdevf_la-srowstdevfa.Tpo $(DEPDIR)/libStdevf_la-srowstdevfa.Plo
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	source='srowstdevfa.c' object='libStdevf_la-srowstdevfa.lo' libtool=yes @AMDEPBACKSLASH@
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	DEPDIR=$(DEPDIR) $(CCDEPMODE) $(depcomp) @AMDEPBACKSLASH@
+@am__fastdepCC_FALSE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -c -o libStdevf_la-srowstdevfa.lo `test -f 'srowstdevfa.c' || echo '$(srcdir)/'`srowstdevfa.c
+
+libStdevf_la-scolumnstdevfa.lo: scolumnstdevfa.c
+@am__fastdepCC_TRUE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -MT libStdevf_la-scolumnstdevfa.lo -MD -MP -MF $(DEPDIR)/libStdevf_la-scolumnstdevfa.Tpo -c -o libStdevf_la-scolumnstdevfa.lo `test -f 'scolumnstdevfa.c' || echo '$(srcdir)/'`scolumnstdevfa.c
+@am__fastdepCC_TRUE@	mv -f $(DEPDIR)/libStdevf_la-scolumnstdevfa.Tpo $(DEPDIR)/libStdevf_la-scolumnstdevfa.Plo
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	source='scolumnstdevfa.c' object='libStdevf_la-scolumnstdevfa.lo' libtool=yes @AMDEPBACKSLASH@
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	DEPDIR=$(DEPDIR) $(CCDEPMODE) $(depcomp) @AMDEPBACKSLASH@
+@am__fastdepCC_FALSE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -c -o libStdevf_la-scolumnstdevfa.lo `test -f 'scolumnstdevfa.c' || echo '$(srcdir)/'`scolumnstdevfa.c
+
+libStdevf_la-dstdevfa.lo: dstdevfa.c
+@am__fastdepCC_TRUE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -MT libStdevf_la-dstdevfa.lo -MD -MP -MF $(DEPDIR)/libStdevf_la-dstdevfa.Tpo -c -o libStdevf_la-dstdevfa.lo `test -f 'dstdevfa.c' || echo '$(srcdir)/'`dstdevfa.c
+@am__fastdepCC_TRUE@	mv -f $(DEPDIR)/libStdevf_la-dstdevfa.Tpo $(DEPDIR)/libStdevf_la-dstdevfa.Plo
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	source='dstdevfa.c' object='libStdevf_la-dstdevfa.lo' libtool=yes @AMDEPBACKSLASH@
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	DEPDIR=$(DEPDIR) $(CCDEPMODE) $(depcomp) @AMDEPBACKSLASH@
+@am__fastdepCC_FALSE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -c -o libStdevf_la-dstdevfa.lo `test -f 'dstdevfa.c' || echo '$(srcdir)/'`dstdevfa.c
+
+libStdevf_la-drowstdevfa.lo: drowstdevfa.c
+@am__fastdepCC_TRUE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -MT libStdevf_la-drowstdevfa.lo -MD -MP -MF $(DEPDIR)/libStdevf_la-drowstdevfa.Tpo -c -o libStdevf_la-drowstdevfa.lo `test -f 'drowstdevfa.c' || echo '$(srcdir)/'`drowstdevfa.c
+@am__fastdepCC_TRUE@	mv -f $(DEPDIR)/libStdevf_la-drowstdevfa.Tpo $(DEPDIR)/libStdevf_la-drowstdevfa.Plo
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	source='drowstdevfa.c' object='libStdevf_la-drowstdevfa.lo' libtool=yes @AMDEPBACKSLASH@
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	DEPDIR=$(DEPDIR) $(CCDEPMODE) $(depcomp) @AMDEPBACKSLASH@
+@am__fastdepCC_FALSE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -c -o libStdevf_la-drowstdevfa.lo `test -f 'drowstdevfa.c' || echo '$(srcdir)/'`drowstdevfa.c
+
+libStdevf_la-dcolumnstdevfa.lo: dcolumnstdevfa.c
+@am__fastdepCC_TRUE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -MT libStdevf_la-dcolumnstdevfa.lo -MD -MP -MF $(DEPDIR)/libStdevf_la-dcolumnstdevfa.Tpo -c -o libStdevf_la-dcolumnstdevfa.lo `test -f 'dcolumnstdevfa.c' || echo '$(srcdir)/'`dcolumnstdevfa.c
+@am__fastdepCC_TRUE@	mv -f $(DEPDIR)/libStdevf_la-dcolumnstdevfa.Tpo $(DEPDIR)/libStdevf_la-dcolumnstdevfa.Plo
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	source='dcolumnstdevfa.c' object='libStdevf_la-dcolumnstdevfa.lo' libtool=yes @AMDEPBACKSLASH@
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	DEPDIR=$(DEPDIR) $(CCDEPMODE) $(depcomp) @AMDEPBACKSLASH@
+@am__fastdepCC_FALSE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -c -o libStdevf_la-dcolumnstdevfa.lo `test -f 'dcolumnstdevfa.c' || echo '$(srcdir)/'`dcolumnstdevfa.c
+
+libStdevf_la-cstdevfa.lo: cstdevfa.c
+@am__fastdepCC_TRUE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -MT libStdevf_la-cstdevfa.lo -MD -MP -MF $(DEPDIR)/libStdevf_la-cstdevfa.Tpo -c -o libStdevf_la-cstdevfa.lo `test -f 'cstdevfa.c' || echo '$(srcdir)/'`cstdevfa.c
+@am__fastdepCC_TRUE@	mv -f $(DEPDIR)/libStdevf_la-cstdevfa.Tpo $(DEPDIR)/libStdevf_la-cstdevfa.Plo
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	source='cstdevfa.c' object='libStdevf_la-cstdevfa.lo' libtool=yes @AMDEPBACKSLASH@
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	DEPDIR=$(DEPDIR) $(CCDEPMODE) $(depcomp) @AMDEPBACKSLASH@
+@am__fastdepCC_FALSE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -c -o libStdevf_la-cstdevfa.lo `test -f 'cstdevfa.c' || echo '$(srcdir)/'`cstdevfa.c
+
+libStdevf_la-crowstdevfa.lo: crowstdevfa.c
+@am__fastdepCC_TRUE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -MT libStdevf_la-crowstdevfa.lo -MD -MP -MF $(DEPDIR)/libStdevf_la-crowstdevfa.Tpo -c -o libStdevf_la-crowstdevfa.lo `test -f 'crowstdevfa.c' || echo '$(srcdir)/'`crowstdevfa.c
+@am__fastdepCC_TRUE@	mv -f $(DEPDIR)/libStdevf_la-crowstdevfa.Tpo $(DEPDIR)/libStdevf_la-crowstdevfa.Plo
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	source='crowstdevfa.c' object='libStdevf_la-crowstdevfa.lo' libtool=yes @AMDEPBACKSLASH@
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	DEPDIR=$(DEPDIR) $(CCDEPMODE) $(depcomp) @AMDEPBACKSLASH@
+@am__fastdepCC_FALSE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -c -o libStdevf_la-crowstdevfa.lo `test -f 'crowstdevfa.c' || echo '$(srcdir)/'`crowstdevfa.c
+
+libStdevf_la-ccolumnstdevfa.lo: ccolumnstdevfa.c
+@am__fastdepCC_TRUE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -MT libStdevf_la-ccolumnstdevfa.lo -MD -MP -MF $(DEPDIR)/libStdevf_la-ccolumnstdevfa.Tpo -c -o libStdevf_la-ccolumnstdevfa.lo `test -f 'ccolumnstdevfa.c' || echo '$(srcdir)/'`ccolumnstdevfa.c
+@am__fastdepCC_TRUE@	mv -f $(DEPDIR)/libStdevf_la-ccolumnstdevfa.Tpo $(DEPDIR)/libStdevf_la-ccolumnstdevfa.Plo
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	source='ccolumnstdevfa.c' object='libStdevf_la-ccolumnstdevfa.lo' libtool=yes @AMDEPBACKSLASH@
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	DEPDIR=$(DEPDIR) $(CCDEPMODE) $(depcomp) @AMDEPBACKSLASH@
+@am__fastdepCC_FALSE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -c -o libStdevf_la-ccolumnstdevfa.lo `test -f 'ccolumnstdevfa.c' || echo '$(srcdir)/'`ccolumnstdevfa.c
+
+libStdevf_la-zstdevfa.lo: zstdevfa.c
+@am__fastdepCC_TRUE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -MT libStdevf_la-zstdevfa.lo -MD -MP -MF $(DEPDIR)/libStdevf_la-zstdevfa.Tpo -c -o libStdevf_la-zstdevfa.lo `test -f 'zstdevfa.c' || echo '$(srcdir)/'`zstdevfa.c
+@am__fastdepCC_TRUE@	mv -f $(DEPDIR)/libStdevf_la-zstdevfa.Tpo $(DEPDIR)/libStdevf_la-zstdevfa.Plo
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	source='zstdevfa.c' object='libStdevf_la-zstdevfa.lo' libtool=yes @AMDEPBACKSLASH@
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	DEPDIR=$(DEPDIR) $(CCDEPMODE) $(depcomp) @AMDEPBACKSLASH@
+@am__fastdepCC_FALSE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -c -o libStdevf_la-zstdevfa.lo `test -f 'zstdevfa.c' || echo '$(srcdir)/'`zstdevfa.c
+
+libStdevf_la-zrowstdevfa.lo: zrowstdevfa.c
+@am__fastdepCC_TRUE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -MT libStdevf_la-zrowstdevfa.lo -MD -MP -MF $(DEPDIR)/libStdevf_la-zrowstdevfa.Tpo -c -o libStdevf_la-zrowstdevfa.lo `test -f 'zrowstdevfa.c' || echo '$(srcdir)/'`zrowstdevfa.c
+@am__fastdepCC_TRUE@	mv -f $(DEPDIR)/libStdevf_la-zrowstdevfa.Tpo $(DEPDIR)/libStdevf_la-zrowstdevfa.Plo
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	source='zrowstdevfa.c' object='libStdevf_la-zrowstdevfa.lo' libtool=yes @AMDEPBACKSLASH@
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	DEPDIR=$(DEPDIR) $(CCDEPMODE) $(depcomp) @AMDEPBACKSLASH@
+@am__fastdepCC_FALSE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -c -o libStdevf_la-zrowstdevfa.lo `test -f 'zrowstdevfa.c' || echo '$(srcdir)/'`zrowstdevfa.c
+
+libStdevf_la-zcolumnstdevfa.lo: zcolumnstdevfa.c
+@am__fastdepCC_TRUE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -MT libStdevf_la-zcolumnstdevfa.lo -MD -MP -MF $(DEPDIR)/libStdevf_la-zcolumnstdevfa.Tpo -c -o libStdevf_la-zcolumnstdevfa.lo `test -f 'zcolumnstdevfa.c' || echo '$(srcdir)/'`zcolumnstdevfa.c
+@am__fastdepCC_TRUE@	mv -f $(DEPDIR)/libStdevf_la-zcolumnstdevfa.Tpo $(DEPDIR)/libStdevf_la-zcolumnstdevfa.Plo
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	source='zcolumnstdevfa.c' object='libStdevf_la-zcolumnstdevfa.lo' libtool=yes @AMDEPBACKSLASH@
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	DEPDIR=$(DEPDIR) $(CCDEPMODE) $(depcomp) @AMDEPBACKSLASH@
+@am__fastdepCC_FALSE@	$(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libStdevf_la_CFLAGS) $(CFLAGS) -c -o libStdevf_la-zcolumnstdevfa.lo `test -f 'zcolumnstdevfa.c' || echo '$(srcdir)/'`zcolumnstdevfa.c
+
+testDoubleStdevf-testDoubleStdevf.o: testDoubleStdevf.c
+@am__fastdepCC_TRUE@	$(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(testDoubleStdevf_CFLAGS) $(CFLAGS) -MT testDoubleStdevf-testDoubleStdevf.o -MD -MP -MF $(DEPDIR)/testDoubleStdevf-testDoubleStdevf.Tpo -c -o testDoubleStdevf-testDoubleStdevf.o `test -f 'testDoubleStdevf.c' || echo '$(srcdir)/'`testDoubleStdevf.c
+@am__fastdepCC_TRUE@	mv -f $(DEPDIR)/testDoubleStdevf-testDoubleStdevf.Tpo $(DEPDIR)/testDoubleStdevf-testDoubleStdevf.Po
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	source='testDoubleStdevf.c' object='testDoubleStdevf-testDoubleStdevf.o' libtool=no @AMDEPBACKSLASH@
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	DEPDIR=$(DEPDIR) $(CCDEPMODE) $(depcomp) @AMDEPBACKSLASH@
+@am__fastdepCC_FALSE@	$(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(testDoubleStdevf_CFLAGS) $(CFLAGS) -c -o testDoubleStdevf-testDoubleStdevf.o `test -f 'testDoubleStdevf.c' || echo '$(srcdir)/'`testDoubleStdevf.c
+
+testDoubleStdevf-testDoubleStdevf.obj: testDoubleStdevf.c
+@am__fastdepCC_TRUE@	$(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(testDoubleStdevf_CFLAGS) $(CFLAGS) -MT testDoubleStdevf-testDoubleStdevf.obj -MD -MP -MF $(DEPDIR)/testDoubleStdevf-testDoubleStdevf.Tpo -c -o testDoubleStdevf-testDoubleStdevf.obj `if test -f 'testDoubleStdevf.c'; then $(CYGPATH_W) 'testDoubleStdevf.c'; else $(CYGPATH_W) '$(srcdir)/testDoubleStdevf.c'; fi`
+@am__fastdepCC_TRUE@	mv -f $(DEPDIR)/testDoubleStdevf-testDoubleStdevf.Tpo $(DEPDIR)/testDoubleStdevf-testDoubleStdevf.Po
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	source='testDoubleStdevf.c' object='testDoubleStdevf-testDoubleStdevf.obj' libtool=no @AMDEPBACKSLASH@
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	DEPDIR=$(DEPDIR) $(CCDEPMODE) $(depcomp) @AMDEPBACKSLASH@
+@am__fastdepCC_FALSE@	$(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(testDoubleStdevf_CFLAGS) $(CFLAGS) -c -o testDoubleStdevf-testDoubleStdevf.obj `if test -f 'testDoubleStdevf.c'; then $(CYGPATH_W) 'testDoubleStdevf.c'; else $(CYGPATH_W) '$(srcdir)/testDoubleStdevf.c'; fi`
+
+testFloatStdevf-testFloatStdevf.o: testFloatStdevf.c
+@am__fastdepCC_TRUE@	$(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(testFloatStdevf_CFLAGS) $(CFLAGS) -MT testFloatStdevf-testFloatStdevf.o -MD -MP -MF $(DEPDIR)/testFloatStdevf-testFloatStdevf.Tpo -c -o testFloatStdevf-testFloatStdevf.o `test -f 'testFloatStdevf.c' || echo '$(srcdir)/'`testFloatStdevf.c
+@am__fastdepCC_TRUE@	mv -f $(DEPDIR)/testFloatStdevf-testFloatStdevf.Tpo $(DEPDIR)/testFloatStdevf-testFloatStdevf.Po
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	source='testFloatStdevf.c' object='testFloatStdevf-testFloatStdevf.o' libtool=no @AMDEPBACKSLASH@
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	DEPDIR=$(DEPDIR) $(CCDEPMODE) $(depcomp) @AMDEPBACKSLASH@
+@am__fastdepCC_FALSE@	$(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(testFloatStdevf_CFLAGS) $(CFLAGS) -c -o testFloatStdevf-testFloatStdevf.o `test -f 'testFloatStdevf.c' || echo '$(srcdir)/'`testFloatStdevf.c
+
+testFloatStdevf-testFloatStdevf.obj: testFloatStdevf.c
+@am__fastdepCC_TRUE@	$(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(testFloatStdevf_CFLAGS) $(CFLAGS) -MT testFloatStdevf-testFloatStdevf.obj -MD -MP -MF $(DEPDIR)/testFloatStdevf-testFloatStdevf.Tpo -c -o testFloatStdevf-testFloatStdevf.obj `if test -f 'testFloatStdevf.c'; then $(CYGPATH_W) 'testFloatStdevf.c'; else $(CYGPATH_W) '$(srcdir)/testFloatStdevf.c'; fi`
+@am__fastdepCC_TRUE@	mv -f $(DEPDIR)/testFloatStdevf-testFloatStdevf.Tpo $(DEPDIR)/testFloatStdevf-testFloatStdevf.Po
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	source='testFloatStdevf.c' object='testFloatStdevf-testFloatStdevf.obj' libtool=no @AMDEPBACKSLASH@
+@AMDEP_TRUE@@am__fastdepCC_FALSE@	DEPDIR=$(DEPDIR) $(CCDEPMODE) $(depcomp) @AMDEPBACKSLASH@
+@am__fastdepCC_FALSE@	$(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(testFloatStdevf_CFLAGS) $(CFLAGS) -c -o testFloatStdevf-testFloatStdevf.obj `if test -f 'testFloatStdevf.c'; then $(CYGPATH_W) 'testFloatStdevf.c'; else $(CYGPATH_W) '$(srcdir)/testFloatStdevf.c'; fi`
+
+mostlyclean-libtool:
+	-rm -f *.lo
+
+clean-libtool:
+	-rm -rf .libs _libs
+
+ID: $(HEADERS) $(SOURCES) $(LISP) $(TAGS_FILES)
+	list='$(SOURCES) $(HEADERS) $(LISP) $(TAGS_FILES)'; \
+	unique=`for i in $$list; do \
+	    if test -f "$$i"; then echo $$i; else echo $(srcdir)/$$i; fi; \
+	  done | \
+	  $(AWK) '{ files[$$0] = 1; nonempty = 1; } \
+	      END { if (nonempty) { for (i in files) print i; }; }'`; \
+	mkid -fID $$unique
+tags: TAGS
+
+TAGS:  $(HEADERS) $(SOURCES)  $(TAGS_DEPENDENCIES) \
+		$(TAGS_FILES) $(LISP)
+	tags=; \
+	here=`pwd`; \
+	list='$(SOURCES) $(HEADERS)  $(LISP) $(TAGS_FILES)'; \
+	unique=`for i in $$list; do \
+	    if test -f "$$i"; then echo $$i; else echo $(srcdir)/$$i; fi; \
+	  done | \
+	  $(AWK) '{ files[$$0] = 1; nonempty = 1; } \
+	      END { if (nonempty) { for (i in files) print i; }; }'`; \
+	if test -z "$(ETAGS_ARGS)$$tags$$unique"; then :; else \
+	  test -n "$$unique" || unique=$$empty_fix; \
+	  $(ETAGS) $(ETAGSFLAGS) $(AM_ETAGSFLAGS) $(ETAGS_ARGS) \
+	    $$tags $$unique; \
+	fi
+ctags: CTAGS
+CTAGS:  $(HEADERS) $(SOURCES)  $(TAGS_DEPENDENCIES) \
+		$(TAGS_FILES) $(LISP)
+	tags=; \
+	list='$(SOURCES) $(HEADERS)  $(LISP) $(TAGS_FILES)'; \
+	unique=`for i in $$list; do \
+	    if test -f "$$i"; then echo $$i; else echo $(srcdir)/$$i; fi; \
+	  done | \
+	  $(AWK) '{ files[$$0] = 1; nonempty = 1; } \
+	      END { if (nonempty) { for (i in files) print i; }; }'`; \
+	test -z "$(CTAGS_ARGS)$$tags$$unique" \
+	  || $(CTAGS) $(CTAGSFLAGS) $(AM_CTAGSFLAGS) $(CTAGS_ARGS) \
+	     $$tags $$unique
+
+GTAGS:
+	here=`$(am__cd) $(top_builddir) && pwd` \
+	  && cd $(top_srcdir) \
+	  && gtags -i $(GTAGS_ARGS) $$here
+
+distclean-tags:
+	-rm -f TAGS ID GTAGS GRTAGS GSYMS GPATH tags
+
+check-TESTS: $(TESTS)
+	@failed=0; all=0; xfail=0; xpass=0; skip=0; \
+	srcdir=$(srcdir); export srcdir; \
+	list=' $(TESTS) '; \
+	if test -n "$$list"; then \
+	  for tst in $$list; do \
+	    if test -f ./$$tst; then dir=./; \
+	    elif test -f $$tst; then dir=; \
+	    else dir="$(srcdir)/"; fi; \
+	    if $(TESTS_ENVIRONMENT) $${dir}$$tst; then \
+	      all=`expr $$all + 1`; \
+	      case " $(XFAIL_TESTS) " in \
+	      *[\ \	]$$tst[\ \	]*) \
+		xpass=`expr $$xpass + 1`; \
+		failed=`expr $$failed + 1`; \
+		echo "XPASS: $$tst"; \
+	      ;; \
+	      *) \
+		echo "PASS: $$tst"; \
+	      ;; \
+	      esac; \
+	    elif test $$? -ne 77; then \
+	      all=`expr $$all + 1`; \
+	      case " $(XFAIL_TESTS) " in \
+	      *[\ \	]$$tst[\ \	]*) \
+		xfail=`expr $$xfail + 1`; \
+		echo "XFAIL: $$tst"; \
+	      ;; \
+	      *) \
+		failed=`expr $$failed + 1`; \
+		echo "FAIL: $$tst"; \
+	      ;; \
+	      esac; \
+	    else \
+	      skip=`expr $$skip + 1`; \
+	      echo "SKIP: $$tst"; \
+	    fi; \
+	  done; \
+	  if test "$$all" -eq 1; then \
+	    tests="test"; \
+	    All=""; \
+	  else \
+	    tests="tests"; \
+	    All="All "; \
+	  fi; \
+	  if test "$$failed" -eq 0; then \
+	    if test "$$xfail" -eq 0; then \
+	      banner="$$All$$all $$tests passed"; \
+	    else \
+	      if test "$$xfail" -eq 1; then failures=failure; else failures=failures; fi; \
+	      banner="$$All$$all $$tests behaved as expected ($$xfail expected $$failures)"; \
+	    fi; \
+	  else \
+	    if test "$$xpass" -eq 0; then \
+	      banner="$$failed of $$all $$tests failed"; \
+	    else \
+	      if test "$$xpass" -eq 1; then passes=pass; else passes=passes; fi; \
+	      banner="$$failed of $$all $$tests did not behave as expected ($$xpass unexpected $$passes)"; \
+	    fi; \
+	  fi; \
+	  dashes="$$banner"; \
+	  skipped=""; \
+	  if test "$$skip" -ne 0; then \
+	    if test "$$skip" -eq 1; then \
+	      skipped="($$skip test was not run)"; \
+	    else \
+	      skipped="($$skip tests were not run)"; \
+	    fi; \
+	    test `echo "$$skipped" | wc -c` -le `echo "$$banner" | wc -c` || \
+	      dashes="$$skipped"; \
+	  fi; \
+	  report=""; \
+	  if test "$$failed" -ne 0 && test -n "$(PACKAGE_BUGREPORT)"; then \
+	    report="Please report to $(PACKAGE_BUGREPORT)"; \
+	    test `echo "$$report" | wc -c` -le `echo "$$banner" | wc -c` || \
+	      dashes="$$report"; \
+	  fi; \
+	  dashes=`echo "$$dashes" | sed s/./=/g`; \
+	  echo "$$dashes"; \
+	  echo "$$banner"; \
+	  test -z "$$skipped" || echo "$$skipped"; \
+	  test -z "$$report" || echo "$$report"; \
+	  echo "$$dashes"; \
+	  test "$$failed" -eq 0; \
+	else :; fi
+
+distdir: $(DISTFILES)
+	@srcdirstrip=`echo "$(srcdir)" | sed 's/[].[^$$\\*]/\\\\&/g'`; \
+	topsrcdirstrip=`echo "$(top_srcdir)" | sed 's/[].[^$$\\*]/\\\\&/g'`; \
+	list='$(DISTFILES)'; \
+	  dist_files=`for file in $$list; do echo $$file; done | \
+	  sed -e "s|^$$srcdirstrip/||;t" \
+	      -e "s|^$$topsrcdirstrip/|$(top_builddir)/|;t"`; \
+	case $$dist_files in \
+	  */*) $(MKDIR_P) `echo "$$dist_files" | \
+			   sed '/\//!d;s|^|$(distdir)/|;s,/[^/]*$$,,' | \
+			   sort -u` ;; \
+	esac; \
+	for file in $$dist_files; do \
+	  if test -f $$file || test -d $$file; then d=.; else d=$(srcdir); fi; \
+	  if test -d $$d/$$file; then \
+	    dir=`echo "/$$file" | sed -e 's,/[^/]*$$,,'`; \
+	    if test -d $(srcdir)/$$file && test $$d != $(srcdir); then \
+	      cp -pR $(srcdir)/$$file $(distdir)$$dir || exit 1; \
+	    fi; \
+	    cp -pR $$d/$$file $(distdir)$$dir || exit 1; \
+	  else \
+	    test -f $(distdir)/$$file \
+	    || cp -p $$d/$$file $(distdir)/$$file \
+	    || exit 1; \
+	  fi; \
+	done
+check-am: all-am
+	$(MAKE) $(AM_MAKEFLAGS) $(check_PROGRAMS)
+	$(MAKE) $(AM_MAKEFLAGS) check-TESTS
+check: check-am
+all-am: Makefile $(LTLIBRARIES)
+installdirs:
+	for dir in "$(DESTDIR)$(pkglibdir)"; do \
+	  test -z "$$dir" || $(MKDIR_P) "$$dir"; \
+	done
+install: install-am
+install-exec: install-exec-am
+install-data: install-data-am
+uninstall: uninstall-am
+
+install-am: all-am
+	@$(MAKE) $(AM_MAKEFLAGS) install-exec-am install-data-am
+
+installcheck: installcheck-am
+install-strip:
+	$(MAKE) $(AM_MAKEFLAGS) INSTALL_PROGRAM="$(INSTALL_STRIP_PROGRAM)" \
+	  install_sh_PROGRAM="$(INSTALL_STRIP_PROGRAM)" INSTALL_STRIP_FLAG=-s \
+	  `test -z '$(STRIP)' || \
+	    echo "INSTALL_PROGRAM_ENV=STRIPPROG='$(STRIP)'"` install
+mostlyclean-generic:
+
+clean-generic:
+
+distclean-generic:
+	-test -z "$(CONFIG_CLEAN_FILES)" || rm -f $(CONFIG_CLEAN_FILES)
+
+maintainer-clean-generic:
+	@echo "This command is intended for maintainers to use"
+	@echo "it deletes files that may require special tools to rebuild."
+clean: clean-am
+
+clean-am: clean-checkPROGRAMS clean-generic clean-libtool \
+	clean-pkglibLTLIBRARIES mostlyclean-am
+
+distclean: distclean-am
+	-rm -rf ./$(DEPDIR)
+	-rm -f Makefile
+distclean-am: clean-am distclean-compile distclean-generic \
+	distclean-tags
+
+dvi: dvi-am
+
+dvi-am:
+
+html: html-am
+
+info: info-am
+
+info-am:
+
+install-data-am:
+
+install-dvi: install-dvi-am
+
+install-exec-am: install-pkglibLTLIBRARIES
+
+install-html: install-html-am
+
+install-info: install-info-am
+
+install-man:
+
+install-pdf: install-pdf-am
+
+install-ps: install-ps-am
+
+installcheck-am:
+
+maintainer-clean: maintainer-clean-am
+	-rm -rf ./$(DEPDIR)
+	-rm -f Makefile
+maintainer-clean-am: distclean-am maintainer-clean-generic
+
+mostlyclean: mostlyclean-am
+
+mostlyclean-am: mostlyclean-compile mostlyclean-generic \
+	mostlyclean-libtool
+
+pdf: pdf-am
+
+pdf-am:
+
+ps: ps-am
+
+ps-am:
+
+uninstall-am: uninstall-pkglibLTLIBRARIES
+
+.MAKE: install-am install-strip
+
+.PHONY: CTAGS GTAGS all all-am check check-TESTS check-am clean \
+	clean-checkPROGRAMS clean-generic clean-libtool \
+	clean-pkglibLTLIBRARIES ctags distclean distclean-compile \
+	distclean-generic distclean-libtool distclean-tags distdir dvi \
+	dvi-am html html-am info info-am install install-am \
+	install-data install-data-am install-dvi install-dvi-am \
+	install-exec install-exec-am install-html install-html-am \
+	install-info install-info-am install-man install-pdf \
+	install-pdf-am install-pkglibLTLIBRARIES install-ps \
+	install-ps-am install-strip installcheck installcheck-am \
+	installdirs maintainer-clean maintainer-clean-generic \
+	mostlyclean mostlyclean-compile mostlyclean-generic \
+	mostlyclean-libtool pdf pdf-am ps ps-am tags uninstall \
+	uninstall-am uninstall-pkglibLTLIBRARIES
+
+# Tell versions [3.59,3.63) of GNU make to not export all variables.
+# Otherwise a system limit (for SysV at least) may be exceeded.
+.NOEXPORT:
diff --git a/src/c/statisticsFunctions/stdevf/ccolumnstdevfa.c b/src/c/statisticsFunctions/stdevf/ccolumnstdevfa.c
index 64c2deed..14d30301 100644
--- a/src/c/statisticsFunctions/stdevf/ccolumnstdevfa.c
+++ b/src/c/statisticsFunctions/stdevf/ccolumnstdevfa.c
@@ -12,14 +12,14 @@
 
 
 #include "stdevf.h"
-
+#include "abs.h"
 
 void ccolumnstdevfa(floatComplex *in1, int lines, int columns, floatComplex*in2, floatComplex* out){
   int i = 0;
   int j = 0;
   floatComplex temp = FloatComplex(0.0f,0.0f);
   floatComplex accumulate = FloatComplex(0.0f,0.0f);
-  float accumulateFre =0.0f ;
+  floatComplex accumulateFre = FloatComplex(0.0f,0.0f);
     
   ccolumnmeanfa(in1, lines, columns, in2, out );
 
@@ -28,7 +28,7 @@ void ccolumnstdevfa(floatComplex *in1, int lines, int columns, floatComplex*in2,
   for (j = 0; j < lines; ++j)
     {
       accumulate = FloatComplex(0.0f,0.0f);
-      accumulateFre =0.0f ;
+      accumulateFre = FloatComplex(0.0f,0.0f);
       temp = FloatComplex(0.0f,0.0f);
 
       for ( i = 0 ; i < columns; ++i )
@@ -38,7 +38,7 @@ void ccolumnstdevfa(floatComplex *in1, int lines, int columns, floatComplex*in2,
          temp = cmuls( in2[lines*i + j] , temp);
 
          accumulate = cadds( temp , accumulate);
-         accumulateFre +=  creals(in2[lines*i + j]);
+         accumulateFre =  cadds (in2[lines*i + j] ,accumulateFre );
 
         } 
 
@@ -48,7 +48,9 @@ void ccolumnstdevfa(floatComplex *in1, int lines, int columns, floatComplex*in2,
         }
       else
         {
-          accumulate = FloatComplex(  creals(accumulate ) / (accumulateFre - 1) , cimags(accumulate)  / (accumulateFre - 1));
+	  if( sabss (creals(accumulate)) <= 3e-6 ) accumulate = FloatComplex(sabss(creals(accumulate)) ,cimags(accumulate));
+	  if( sabss (cimags(accumulate)) <= 3e-6 ) accumulate = FloatComplex(creals(accumulate) ,sabss(cimags(accumulate)));
+          accumulate =   crdivs (accumulate , cdiffs (accumulateFre ,FloatComplex(1.0f,0.0f))   );
           out[j] =csqrts(accumulate);
         }
     }
diff --git a/src/c/statisticsFunctions/stdevf/crowstdevfa.c b/src/c/statisticsFunctions/stdevf/crowstdevfa.c
index b22f90a9..38d8cb23 100644
--- a/src/c/statisticsFunctions/stdevf/crowstdevfa.c
+++ b/src/c/statisticsFunctions/stdevf/crowstdevfa.c
@@ -13,7 +13,7 @@
 
 #include "stdevf.h"
 #include "meanf.h"
-
+#include "abs.h"
 
 
 void crowstdevfa(floatComplex *in1, int lines, int columns, floatComplex*in2, floatComplex* out){
@@ -21,7 +21,7 @@ void crowstdevfa(floatComplex *in1, int lines, int columns, floatComplex*in2, fl
   int j = 0;
   floatComplex temp = FloatComplex(0.0f,0.0f);
   floatComplex accumulate = FloatComplex(0.0f,0.0f);
-  float accumulateFre = 0.0f ;
+  floatComplex accumulateFre = FloatComplex(0.0f,0.0f);
     
   crowmeanfa(in1, lines, columns, in2, out );
 
@@ -29,7 +29,7 @@ void crowstdevfa(floatComplex *in1, int lines, int columns, floatComplex*in2, fl
   for (j = 0; j < columns; ++j)
     {
       accumulate = FloatComplex(0.0f,0.0f);
-      accumulateFre = 0.0f ;
+      accumulateFre =FloatComplex(0.0f,0.0f);
       temp = FloatComplex(0.0f,0.0f);
 
       for ( i = 0 ; i < lines; ++i )
@@ -39,7 +39,7 @@ void crowstdevfa(floatComplex *in1, int lines, int columns, floatComplex*in2, fl
          temp = cmuls( in2[lines*j + i] , temp);
 
          accumulate = cadds( temp , accumulate);
-         accumulateFre +=  creals(in2[lines*j + i]);
+         accumulateFre =  cadds (in2[lines*j + i] ,accumulateFre );
 
         } 
 
@@ -49,8 +49,10 @@ void crowstdevfa(floatComplex *in1, int lines, int columns, floatComplex*in2, fl
         }
       else
         {
-           accumulate = FloatComplex(  creals(accumulate ) / (accumulateFre - 1) , cimags(accumulate)  / (accumulateFre - 1));
-           out[j] = csqrts(accumulate);
+	  if( sabss (creals(accumulate)) <= 3e-6 ) accumulate = FloatComplex(sabss(creals(accumulate)) ,cimags(accumulate));
+	  if( sabss (cimags(accumulate)) <= 3e-6 ) accumulate = FloatComplex(creals(accumulate) ,sabss(cimags(accumulate)));
+          accumulate =   crdivs (accumulate , cdiffs (accumulateFre ,FloatComplex(1.0f,0.0f))   );
+          out[j] =csqrts(accumulate);
         }
     }
 
diff --git a/src/c/statisticsFunctions/stdevf/cstdevfa.c b/src/c/statisticsFunctions/stdevf/cstdevfa.c
index 63ce3948..5b1631f1 100644
--- a/src/c/statisticsFunctions/stdevf/cstdevfa.c
+++ b/src/c/statisticsFunctions/stdevf/cstdevfa.c
@@ -15,9 +15,10 @@
 floatComplex cstdevfa(floatComplex *in1, int lines, int columns, floatComplex* in2)
 {
   int i = 0 ;
+
   floatComplex temp = FloatComplex(0.0f,0.0f);
   floatComplex accumulate = FloatComplex(0.0f,0.0f);
-  float accumulateFre = 0.0f ;
+  floatComplex accumulateFre = FloatComplex(0.0f,0.0f);
   floatComplex meanf = cmeanfa (in1 , lines*columns , in2);
 /*equivalent to  (in1 - meanf(x , in2 )).^2 .*in2   */
   for(i = 0 ; i < lines*columns ; ++i)
@@ -26,16 +27,17 @@ floatComplex cstdevfa(floatComplex *in1, int lines, int columns, floatComplex* i
      temp = cmuls( in2[i] , temp);
 
      accumulate = cadds( temp , accumulate);
-     accumulateFre +=  creals(in2[i]);
-    }
+     accumulateFre =  cadds (in2[i] ,accumulateFre );
+
+     } 
 
-    if (lines <= 1)
+     if (lines <= 1)
       {
-        return cmuls (FloatComplex(0.0f,0.0f) , accumulate ) ; 
+        return  cmuls (FloatComplex(0.0f,0.0f) , accumulate ) ; 
       }
-    else
+     else
       {
-        accumulate = FloatComplex(  creals(accumulate ) / (accumulateFre - 1) , cimags(accumulate)  / (accumulateFre - 1));
+        accumulate =   crdivs (accumulate , cdiffs (accumulateFre ,FloatComplex(1.0f,0.0f))   );
         return csqrts(accumulate);
       }
 }
diff --git a/src/c/statisticsFunctions/stdevf/cstdevfcs.c b/src/c/statisticsFunctions/stdevf/cstdevfcs.c
new file mode 100644
index 00000000..a3b7239f
--- /dev/null
+++ b/src/c/statisticsFunctions/stdevf/cstdevfcs.c
@@ -0,0 +1,41 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008-2008 - INRIA - Allan SIMON
+ *
+ *  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 "stdevf.h" 
+
+floatComplex cstdevfcs(floatComplex *in1, int lines, int columns, float* in2)
+{
+  int i = 0 ;
+  floatComplex temp = FloatComplex(0.0f,0.0f);
+  floatComplex accumulate = FloatComplex(0.0f,0.0f);
+  float accumulateFre = 0.0f ;
+  floatComplex meanf = cmeanfcs (in1 , lines ,columns , in2);
+/*equivalent to  (in1 - meanf(x , in2 )).^2 .*in2   */
+  for(i = 0 ; i < lines*columns ; ++i)
+    {
+     temp = cpows (  cdiffs (in1[i] , meanf ) ,FloatComplex (2.0f, 0.0f ) );
+     temp = cmuls( FloatComplex (in2[i],0.0f) , temp);
+
+     accumulate = cadds( temp , accumulate);
+     accumulateFre +=  in2[i];
+    }
+
+    if (lines <= 1)
+      {
+        return cmuls (FloatComplex(0.0f,0.0f) , accumulate ) ; 
+      }
+    else
+      {
+        accumulate = FloatComplex(  creals(accumulate ) / (accumulateFre - 1) , cimags(accumulate)  / (accumulateFre - 1));
+        return csqrts(accumulate);
+      }
+}
diff --git a/src/c/statisticsFunctions/stdevf/cstdevfsc.c b/src/c/statisticsFunctions/stdevf/cstdevfsc.c
new file mode 100644
index 00000000..bb4b6210
--- /dev/null
+++ b/src/c/statisticsFunctions/stdevf/cstdevfsc.c
@@ -0,0 +1,25 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008-2008 - INRIA - Allan SIMON
+ *
+ *  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 "stdevf.h" 
+#include "zeros.h"
+#include <malloc.h>
+
+
+floatComplex cstdevfsc(float *in1, int lines, int columns, floatComplex* in2)
+{
+    float* ZEROS;
+    ZEROS=malloc((unsigned int)(lines*columns*sizeof(float)));
+    szerosa(ZEROS,lines,columns);
+    
+    return cstdevfa(FloatComplexMatrix(in1,ZEROS,lines*columns), lines , columns , in2);
+}
diff --git a/src/c/statisticsFunctions/stdevf/dcolumnstdevfa.c b/src/c/statisticsFunctions/stdevf/dcolumnstdevfa.c
index 9c8e76d3..a6220d57 100644
--- a/src/c/statisticsFunctions/stdevf/dcolumnstdevfa.c
+++ b/src/c/statisticsFunctions/stdevf/dcolumnstdevfa.c
@@ -14,7 +14,6 @@
 #include "stdevf.h"
 #include "meanf.h"
 
-#include <stdio.h>
 
 void dcolumnstdevfa(double *in1, int lines, int columns, double *in2, double* out){
   int i = 0;
diff --git a/src/c/statisticsFunctions/stdevf/scolumnstdevfa.c b/src/c/statisticsFunctions/stdevf/scolumnstdevfa.c
index d39ce88a..f3552ea7 100644
--- a/src/c/statisticsFunctions/stdevf/scolumnstdevfa.c
+++ b/src/c/statisticsFunctions/stdevf/scolumnstdevfa.c
@@ -14,7 +14,6 @@
 #include "stdevf.h"
 #include "meanf.h"
 
-#include <stdio.h>
 
 void scolumnstdevfa(float *in1, int lines, int columns, float *in2, float* out){
   int i = 0;
diff --git a/src/c/statisticsFunctions/stdevf/srowstdevfa.c b/src/c/statisticsFunctions/stdevf/srowstdevfa.c
index 486111b6..337b80e3 100644
--- a/src/c/statisticsFunctions/stdevf/srowstdevfa.c
+++ b/src/c/statisticsFunctions/stdevf/srowstdevfa.c
@@ -14,7 +14,6 @@
 #include "stdevf.h"
 #include "meanf.h"
 
-#include <stdio.h>
 
 void srowstdevfa(float *in1, int lines, int columns, float *in2, float* out){
   int i = 0;
@@ -24,7 +23,6 @@ void srowstdevfa(float *in1, int lines, int columns, float *in2, float* out){
   float accumulateFre = 0.0f ;
     
   srowmeanfa(in1, lines, columns, in2, out );
-printf ("\nmfwerfwerfwefwef \n" );
 
   /*we first multiply each cell of the input matrix by its coefficient*/
   for (j = 0; j < columns; ++j)
diff --git a/src/c/statisticsFunctions/stdevf/sstdevfa.c b/src/c/statisticsFunctions/stdevf/sstdevfa.c
index fd67544b..5ccdab84 100644
--- a/src/c/statisticsFunctions/stdevf/sstdevfa.c
+++ b/src/c/statisticsFunctions/stdevf/sstdevfa.c
@@ -11,7 +11,7 @@
  */
 
 #include "stdevf.h" 
-#include <stdio.h>
+
 float sstdevfa(float *in1, int lines, int columns, float* in2)
 {
   int i = 0 ;
diff --git a/src/c/statisticsFunctions/stdevf/testDoubleStdevf.c b/src/c/statisticsFunctions/stdevf/testDoubleStdevf.c
index 4b979e4d..3b0b61ae 100644
--- a/src/c/statisticsFunctions/stdevf/testDoubleStdevf.c
+++ b/src/c/statisticsFunctions/stdevf/testDoubleStdevf.c
@@ -188,6 +188,12 @@ static int zstdevfaTest(void) {
   double coefI1[9] = { 0.0, 0.0, 0.0, 0.0 , 0.0, 0.0, 0.0, 0.0, 0.0}; 
 
 
+  double tableR2[9] = {1.0, 4.0, 7.0, 2.0 , 5.0, 8.0, 3.0, 6.0, 9.0};
+  double tableI2[9] = {0.0, 0.0, 0.0, 0.0 , 0.0, 0.0, 0.0, 0.0, 0.0};
+  double coefR2[9] = { 0.0, 0.0, 0.0, 0.0 , 0.0, 0.0, 0.0, 0.0, 0.0}; 
+  double coefI2[9] = {10.0, 1.0, 5.0,11.0 , 2.0, 6.0,12.0, 3.0, 7.0}; 
+
+
   double tableR3[10] = {3.0, 6.0, 9.0,10.0, 5.0,
                        6.0,18.0, 7.0,14.0, 2.0};
   double tableI3[10] = {3.0,12.0,25.0, 1.0, 2.0,
@@ -203,6 +209,9 @@ static int zstdevfaTest(void) {
   doubleComplex* table1 = DoubleComplexMatrix (tableR1, tableI1, 9);
   doubleComplex* coef1 = DoubleComplexMatrix (coefR1, coefI1, 9);
 
+  doubleComplex* table2 = DoubleComplexMatrix (tableR2, tableI2, 9);
+  doubleComplex* coef2 = DoubleComplexMatrix (coefR2, coefI2, 9);
+
   doubleComplex* table3 = DoubleComplexMatrix (tableR3, tableI3, 10);
   doubleComplex* coef3 = DoubleComplexMatrix (coefR3, coefI3, 10);
 
@@ -214,6 +223,11 @@ static int zstdevfaTest(void) {
   assert(  zreals(result) == 0.0 );
   assert(  zimags(result) == 0.0 );
 
+  result = zstdevfa(table2,1, 9, coef2);
+  printf("\nresult2_1_9 : %e \t+ %e i " ,zreals(result) ,zimags(result)) ;
+  assert(  zreals(result) == 0.0 );
+  assert(  zimags(result) == 0.0 );
+
 
   result = zstdevfa(table3,2, 5, coef3);
   printf("\nresult_2_5 : %e \t+ %e i " ,zreals(result) ,zimags(result)) ;
@@ -243,9 +257,36 @@ static int zrowstdevfaTest(void) {
 
   doubleComplex* table2 = DoubleComplexMatrix (tableR2, tableI2, 10);
   doubleComplex* coef2 = DoubleComplexMatrix (coefR2, coefI2, 10);
-
+ 
   doubleComplex rowStdevfedTable1_3_3[3];
   doubleComplex rowStdevfedTable2_2_5[5];
+  /*test when the matrix is only imaginary*/
+  double tableR3[10] = { 0.0, 0.0, 0.0, 0.0 , 0.0, 0.0, 0.0, 0.0, 0.0, 0.0}; 
+  double tableI3[10] = {1.0, 2.0, 3.0, 4.0 , 5.0, 6.0, 7.0, 8.0, 9.0, 15.0};
+  double coefR3[10]  = {10.0, 1.0, 5.0,11.0 , 2.0, 6.0,12.0, 3.0, 7.0, 19.0}; 
+  double coefI3[10]  = { 0.0, 0.0, 0.0, 0.0 , 0.0, 0.0, 0.0, 0.0, 0.0, 0.0}; 
+
+
+  doubleComplex* table3 = DoubleComplexMatrix (tableR3, tableI3, 10);
+  doubleComplex* coef3  = DoubleComplexMatrix (coefR3, coefI3, 10);
+ 
+  doubleComplex rowStdevfedTable3_2_5[5];
+ /**/
+  /*test when the matrix is only real and the coef matrix is only imaginary*/
+  double tableR4[10] = {1.0, 4.0, 7.0, 2.0 , 5.0, 8.0, 3.0, 6.0, 9.0, 10.0};
+  double tableI4[10] = { 0.0, 0.0, 0.0, 0.0 , 0.0, 0.0, 0.0, 0.0, 0.0 , 0.0};
+  double coefR4[10]  = { 0.0, 0.0, 0.0, 0.0 , 0.0, 0.0, 0.0, 0.0, 0.0 , 0.0}; 
+  double coefI4[10]   ={10.0, 1.0, 5.0,11.0 , 2.0, 6.0,12.0, 3.0, 7.0, 19.0};
+
+ 
+
+
+  doubleComplex* table4 = DoubleComplexMatrix (tableR4, tableI4, 10);
+  doubleComplex* coef4  = DoubleComplexMatrix (coefR4, coefI4, 10);
+ 
+  doubleComplex rowStdevfedTable4_2_5[5];
+ /**/
+
 
   printf("\n>>>> Row stdevf Double Complex Array Test\n");
 
@@ -259,6 +300,19 @@ static int zrowstdevfaTest(void) {
   rowStdevfedTable2_2_5[3] = DoubleComplex(0.0, 0.0);
   rowStdevfedTable2_2_5[4] = DoubleComplex(0.0, 0.0);
 
+  rowStdevfedTable3_2_5[0] = DoubleComplex(0.0, 0.0);
+  rowStdevfedTable3_2_5[1] = DoubleComplex(0.0, 0.0);
+  rowStdevfedTable3_2_5[2] = DoubleComplex(0.0, 0.0);
+  rowStdevfedTable3_2_5[3] = DoubleComplex(0.0, 0.0);
+  rowStdevfedTable3_2_5[4] = DoubleComplex(0.0, 0.0);
+
+
+  rowStdevfedTable4_2_5[0] = DoubleComplex(0.0, 0.0);
+  rowStdevfedTable4_2_5[1] = DoubleComplex(0.0, 0.0);
+  rowStdevfedTable4_2_5[2] = DoubleComplex(0.0, 0.0);
+  rowStdevfedTable4_2_5[3] = DoubleComplex(0.0, 0.0);
+  rowStdevfedTable4_2_5[4] = DoubleComplex(0.0, 0.0);
+
 
   zrowstdevfa(table1 , 3 , 3 , coef1 , rowStdevfedTable1_3_3);
 
@@ -266,34 +320,75 @@ static int zrowstdevfaTest(void) {
     {
       printf("rowStdevfedTable_3_3[%d] = %e + %ei\n", i, zreals(rowStdevfedTable1_3_3[i]), zimags(rowStdevfedTable1_3_3[i]));
     }
-  assert( fabs(zreals(rowStdevfedTable1_3_3[0])  - 2.8394542 ) / fabs ( zimags(rowStdevfedTable1_3_3[0])  ) <  1e-6 );
-  assert( fabs(zimags(rowStdevfedTable1_3_3[0])  - 0.9464847) / fabs ( zreals(rowStdevfedTable1_3_3[0])  ) <  1e-6 );
+  assert( fabs(zreals(rowStdevfedTable1_3_3[0])  - 2.8394542 ) / fabs ( zreals(rowStdevfedTable1_3_3[0])  ) <  1e-6 );
+  assert( fabs(zimags(rowStdevfedTable1_3_3[0])  - 0.9464847) / fabs ( zimags(rowStdevfedTable1_3_3[0])  ) <  1e-6 );
   
-  assert( fabs(zreals(rowStdevfedTable1_3_3[1])  - 2.8003759) / fabs ( zimags(rowStdevfedTable1_3_3[1])  ) <  1e-6 );
-  assert( fabs(zimags(rowStdevfedTable1_3_3[1])  - 0.9334586) / fabs ( zreals(rowStdevfedTable1_3_3[1])  ) <  1e-6 );
+  assert( fabs(zreals(rowStdevfedTable1_3_3[1])  - 2.8003759) / fabs ( zreals(rowStdevfedTable1_3_3[1])  ) <  1e-6 );
+  assert( fabs(zimags(rowStdevfedTable1_3_3[1])  - 0.9334586) / fabs ( zimags(rowStdevfedTable1_3_3[1])  ) <  1e-6 );
 
-  assert( fabs(zreals(rowStdevfedTable1_3_3[2])  - 2.7669196) / fabs ( zimags(rowStdevfedTable1_3_3[2])  ) <  1e-6 );
-  assert( fabs(zimags(rowStdevfedTable1_3_3[2])  - 0.9223065) / fabs ( zreals(rowStdevfedTable1_3_3[2])  ) <  1e-6 );
+  assert( fabs(zreals(rowStdevfedTable1_3_3[2])  - 2.7669196) / fabs ( zreals(rowStdevfedTable1_3_3[2])  ) <  1e-6 );
+  assert( fabs(zimags(rowStdevfedTable1_3_3[2])  - 0.9223065) / fabs ( zimags(rowStdevfedTable1_3_3[2])  ) <  1e-6 );
 /**/
   zrowstdevfa(table2 , 2 , 5 , coef2 , rowStdevfedTable2_2_5);
   for (i = 0 ; i < 5 ; ++i )
     {
       printf("rowStdevfedTable_2_5[%d] = %e + %ei\n", i, zreals(rowStdevfedTable2_2_5[i]), zimags(rowStdevfedTable2_2_5[i]));
     }
-  assert( fabs(zreals(rowStdevfedTable2_2_5[0])  - 0.9045340 ) / fabs ( zimags(rowStdevfedTable2_2_5[0])  ) <  1e-6 );
-  assert( fabs(zimags(rowStdevfedTable2_2_5[0])  - 0.3015113 ) / fabs ( zreals(rowStdevfedTable2_2_5[0])  ) <  1e-6 );
+  assert( fabs(zreals(rowStdevfedTable2_2_5[0])  - 0.9045340 ) / fabs ( zreals(rowStdevfedTable2_2_5[0])  ) <  1e-6 );
+  assert( fabs(zimags(rowStdevfedTable2_2_5[0])  - 0.3015113 ) / fabs ( zimags(rowStdevfedTable2_2_5[0])  ) <  1e-6 );
+
+  assert( fabs(zreals(rowStdevfedTable2_2_5[1])  - 2.3935678 ) / fabs ( zreals(rowStdevfedTable2_2_5[1])  ) <  1e-6 );
+  assert( fabs(zimags(rowStdevfedTable2_2_5[1])  + 0.4787136) / fabs ( zimags(rowStdevfedTable2_2_5[1])  ) <  1e-6 );
+
+  assert( fabs(zreals(rowStdevfedTable2_2_5[2])  - 1.3887301 ) / fabs ( zreals(rowStdevfedTable2_2_5[2])  ) <  1e-6 );
+  assert( fabs(zimags(rowStdevfedTable2_2_5[2])  - 0.4629100 ) / fabs ( zimags(rowStdevfedTable2_2_5[2])  ) <  1e-6 );
+
+  assert( fabs(zreals(rowStdevfedTable2_2_5[3])  - 1.2421180 ) / fabs ( zreals(rowStdevfedTable2_2_5[3])  ) <  1e-6 );
+  assert( fabs(zimags(rowStdevfedTable2_2_5[3])  - 0.4140393 ) / fabs ( zimags(rowStdevfedTable2_2_5[3])  ) <  1e-6 );
+
+  assert( fabs(zreals(rowStdevfedTable2_2_5[4])  - 0.4523443 ) / fabs ( zreals(rowStdevfedTable2_2_5[4])  ) <  1e-6 );
+  assert( fabs(zimags(rowStdevfedTable2_2_5[4])  - 2.7140659) / fabs ( zimags(rowStdevfedTable2_2_5[4])  ) <  1e-6 );
+/**/
+  zrowstdevfa(table3, 2 , 5 , coef3 , rowStdevfedTable3_2_5);
+  for (i = 0 ; i < 5 ; ++i )
+    {
+      printf("rowStdevfedTable3_2_5[%d] = %e + %ei\n", i, zreals(rowStdevfedTable3_2_5[i]), zimags(rowStdevfedTable3_2_5[i]));
+    }
+  
+  assert( fabs(zimags(rowStdevfedTable3_2_5[0])  - 0.3015113 ) / fabs (zimags(rowStdevfedTable3_2_5[0])  ) <  1e-6 );
+
+
+  assert( fabs(zimags(rowStdevfedTable3_2_5[1])  - 0.4787136) / fabs  ( zimags(rowStdevfedTable3_2_5[1])  ) <  1e-6 );
+
+
+  assert( fabs(zimags(rowStdevfedTable3_2_5[2])  - 0.4629100 ) / fabs (zimags(rowStdevfedTable3_2_5[2])  ) <  1e-6 );
+
+
+  assert( fabs(zimags(rowStdevfedTable3_2_5[3])  - 0.4140393 ) / fabs ( zimags(rowStdevfedTable3_2_5[3])  ) <  1e-6 );
+
+
+  assert( fabs(zimags(rowStdevfedTable3_2_5[4])  - 2.7140659) / fabs  ( zimags(rowStdevfedTable3_2_5[4])  ) <  1e-6 );
+
+/**/
+  zrowstdevfa(table4 , 2 , 5 , coef4 , rowStdevfedTable4_2_5);
+  for (i = 0 ; i < 5 ; ++i )
+    {
+      printf("rowStdevfedTable4_2_5[%d] = %e + %ei\n", i, zreals(rowStdevfedTable4_2_5[i]), zimags(rowStdevfedTable4_2_5[i]));
+    }
+  assert( fabs(zreals(rowStdevfedTable4_2_5[0])  - 0.8597825 ) / fabs ( zreals(rowStdevfedTable4_2_5[0])  ) <  1e-6 );
+  assert( fabs(zimags(rowStdevfedTable4_2_5[0])  + 0.0390006 ) / fabs ( zimags(rowStdevfedTable4_2_5[0])  ) <  1e-6 );
 
-  assert( fabs(zreals(rowStdevfedTable2_2_5[1])  - 2.3935678 ) / fabs ( zimags(rowStdevfedTable2_2_5[1])  ) <  1e-6 );
-  assert( fabs(zimags(rowStdevfedTable2_2_5[1])  + 0.4787136) / fabs ( zreals(rowStdevfedTable2_2_5[1])  ) <  1e-6 );
+  assert( fabs(zreals(rowStdevfedTable4_2_5[1])  - 2.3141768 ) / fabs ( zreals(rowStdevfedTable4_2_5[1])  ) <  1e-6 );
+  assert( fabs(zimags(rowStdevfedTable4_2_5[1])  + 0.0722475) / fabs ( zimags(rowStdevfedTable4_2_5[1])  ) <  1e-6 );
 
-  assert( fabs(zreals(rowStdevfedTable2_2_5[2])  - 1.3887301 ) / fabs ( zimags(rowStdevfedTable2_2_5[2])  ) <  1e-6 );
-  assert( fabs(zimags(rowStdevfedTable2_2_5[2])  - 0.4629100 ) / fabs ( zreals(rowStdevfedTable2_2_5[2])  ) <  1e-6 );
+  assert( fabs(zreals(rowStdevfedTable4_2_5[2])  - 1.2915122 ) / fabs ( zreals(rowStdevfedTable4_2_5[2])  ) <  1e-6 );
+  assert( fabs(zimags(rowStdevfedTable4_2_5[2])  + 0.0804066 ) / fabs ( zimags(rowStdevfedTable4_2_5[2])  ) <  1e-6 );
 
-  assert( fabs(zreals(rowStdevfedTable2_2_5[3])  - 1.2421180 ) / fabs ( zimags(rowStdevfedTable2_2_5[3])  ) <  1e-6 );
-  assert( fabs(zimags(rowStdevfedTable2_2_5[3])  - 0.4140393 ) / fabs ( zreals(rowStdevfedTable2_2_5[3])  ) <  1e-6 );
+  assert( fabs(zreals(rowStdevfedTable4_2_5[3])  - 1.1980065 ) / fabs ( zreals(rowStdevfedTable4_2_5[3])  ) <  1e-6 );
+  assert( fabs(zimags(rowStdevfedTable4_2_5[3])  + 0.0398893 ) / fabs ( zimags(rowStdevfedTable4_2_5[3])  ) <  1e-6 );
 
-  assert( fabs(zreals(rowStdevfedTable2_2_5[4])  - 0.4523443 ) / fabs ( zimags(rowStdevfedTable2_2_5[4])  ) <  1e-6 );
-  assert( fabs(zimags(rowStdevfedTable2_2_5[4])  - 2.7140659) / fabs ( zreals(rowStdevfedTable2_2_5[4])  ) <  1e-6 );
+  assert( fabs(zreals(rowStdevfedTable4_2_5[4])  - 0.4433143 ) / fabs ( zreals(rowStdevfedTable4_2_5[4])  ) <  1e-6 );
+  assert( fabs(zimags(rowStdevfedTable4_2_5[4])  + 0.00852212) / fabs ( zimags(rowStdevfedTable4_2_5[4])  ) <  1e-6 );
 
   return 0;
 }
@@ -315,14 +410,24 @@ static int zcolumnstdevfaTest(void) {
   double coefR2[10]  = {10.0, 1.0, 5.0,11.0 , 2.0, 6.0,12.0, 3.0, 7.0, 19.0}; 
   double coefI2[10]   = { 0.0, 0.0, 0.0, 0.0 , 0.0, 0.0, 0.0, 0.0, 0.0, 0.0}; 
 
+  double tableR3[6] = {4.0, 5.0, 6.0, 7.0, 8.0, 9.0 };
+  double tableI3[6] = {0.0, 0.0, 0.0, 0.0, 0.0, 0.0 };
+
+  double coefR3[6]  = {0.0, 0.0, 0.0, 0.0, 0.0, 0.0};  
+  double coefI3[6]  = {1.0, 2.0, 3.0, 4.0, 5.0, 6.0}; 
+
   doubleComplex* table1 = DoubleComplexMatrix (tableR1, tableI1, 9);
   doubleComplex* coef1 = DoubleComplexMatrix (coefR1, coefI1, 9);
 
   doubleComplex* table2 = DoubleComplexMatrix (tableR2, tableI2, 10);
   doubleComplex* coef2 = DoubleComplexMatrix (coefR2, coefI2, 10);
 
+  doubleComplex* table3 = DoubleComplexMatrix (tableR3, tableI3, 10);
+  doubleComplex* coef3 = DoubleComplexMatrix (coefR3, coefI3, 10);
+
   doubleComplex rowStdevfedTable1_3_3[3];
   doubleComplex rowStdevfedTable2_2_5[2];
+  doubleComplex rowStdevfedTable3_3_2[3];
  
   rowStdevfedTable1_3_3[0] = DoubleComplex(0.0, 0.0);
   rowStdevfedTable1_3_3[1] = DoubleComplex(0.0, 0.0);
@@ -331,6 +436,9 @@ static int zcolumnstdevfaTest(void) {
   rowStdevfedTable2_2_5[0] = DoubleComplex(0.0, 0.0);
   rowStdevfedTable2_2_5[1] = DoubleComplex(0.0, 0.0);
 
+  rowStdevfedTable3_3_2[0] = DoubleComplex(0.0, 0.0);
+  rowStdevfedTable3_3_2[1] = DoubleComplex(0.0, 0.0);
+  rowStdevfedTable3_3_2[2] = DoubleComplex(0.0, 0.0);
 
 
   printf("\n>>>> Column stdevf Double Complex Array Test\n");
@@ -363,7 +471,19 @@ static int zcolumnstdevfaTest(void) {
   assert( fabs(zreals(rowStdevfedTable2_2_5[1])  - 3.1890422 ) / fabs ( zimags(rowStdevfedTable2_2_5[1])  ) <  1e-6 );
   assert( fabs(zimags(rowStdevfedTable2_2_5[1])  - 4.9966774 ) / fabs ( zreals(rowStdevfedTable2_2_5[1])  ) <  1e-6 );
 
+  zcolumnstdevfa(table3 , 3 , 2 , coef3 , rowStdevfedTable3_3_2);
+  for (i = 0 ; i < 3 ; ++i )
+    {
+      printf("columnStdevfedTable3_3_2[%d] = %e + %ei\n", i, zreals(rowStdevfedTable3_3_2[i]), zimags(rowStdevfedTable3_3_2[i]));
+    }
+  assert( fabs(zreals(rowStdevfedTable3_3_2[0])  - 1.1825083) / fabs ( zimags(rowStdevfedTable3_3_2[0])  ) <  1e-6 );
+  assert( fabs(zimags(rowStdevfedTable3_3_2[0])  + 0.1170914) / fabs ( zreals(rowStdevfedTable3_3_2[0])  ) <  1e-6 );
+  
+  assert( fabs(zreals(rowStdevfedTable3_3_2[1])  - 1.3450417) / fabs ( zimags(rowStdevfedTable3_3_2[1])  ) <  1e-6 );
+  assert( fabs(zimags(rowStdevfedTable3_3_2[1])  + 0.0955892) / fabs ( zreals(rowStdevfedTable3_3_2[1])  ) <  1e-6 );
 
+  assert( fabs(zreals(rowStdevfedTable3_3_2[2])  - 1.4077246 ) / fabs ( zimags(rowStdevfedTable3_3_2[2])  ) <  1e-6 );
+  assert( fabs(zimags(rowStdevfedTable3_3_2[2])  + 0.0779670 ) / fabs ( zreals(rowStdevfedTable3_3_2[2])  ) <  1e-6 );
   return 0;
 }
 
diff --git a/src/c/statisticsFunctions/stdevf/zcolumnstdevfa.c b/src/c/statisticsFunctions/stdevf/zcolumnstdevfa.c
index 3d1e8ec9..7dadda70 100644
--- a/src/c/statisticsFunctions/stdevf/zcolumnstdevfa.c
+++ b/src/c/statisticsFunctions/stdevf/zcolumnstdevfa.c
@@ -13,15 +13,14 @@
 
 #include "stdevf.h"
 #include "meanf.h"
-
-#include <stdio.h>
+#include "abs.h"
 
 void zcolumnstdevfa(doubleComplex *in1, int lines, int columns, doubleComplex*in2, doubleComplex* out){
   int i = 0;
   int j = 0;
   doubleComplex temp = DoubleComplex(0.0,0.0);
   doubleComplex accumulate = DoubleComplex(0.0,0.0);
-  double accumulateFre = 0.0;
+  doubleComplex accumulateFre = DoubleComplex(0.0,0.0);
     
   zcolumnmeanfa(in1, lines, columns, in2, out );
 
@@ -30,7 +29,7 @@ void zcolumnstdevfa(doubleComplex *in1, int lines, int columns, doubleComplex*in
   for (j = 0; j < lines; ++j)
     {
       accumulate = DoubleComplex(0.0,0.0);
-      accumulateFre = 0.0;
+      accumulateFre = DoubleComplex(0.0,0.0);
       temp = DoubleComplex(0.0,0.0);
 
       for ( i = 0 ; i < columns; ++i )
@@ -40,7 +39,7 @@ void zcolumnstdevfa(doubleComplex *in1, int lines, int columns, doubleComplex*in
          temp = zmuls( in2[lines*i + j] , temp);
 
          accumulate = zadds( temp , accumulate);
-         accumulateFre +=  zreals(in2[lines*i + j]);
+         accumulateFre = zadds(in2[lines*i + j] ,accumulateFre );
 
         } 
 
@@ -50,7 +49,9 @@ void zcolumnstdevfa(doubleComplex *in1, int lines, int columns, doubleComplex*in
         }
       else
         {
-        accumulate = DoubleComplex(  zreals(accumulate ) / (accumulateFre - 1) , zimags(accumulate)  / (accumulateFre - 1));
+	  if( dabss (zreals(accumulate)) <= 3e-14 ) accumulate = DoubleComplex(dabss(zreals(accumulate)) ,zimags(accumulate));
+	  if( dabss (zimags(accumulate)) <= 3e-14 ) accumulate = DoubleComplex(zreals(accumulate) ,dabss(zimags(accumulate)));
+       accumulate = zrdivs(accumulate, zdiffs(accumulateFre ,DoubleComplex(1.0,0.0)));
         out[j] = zsqrts(accumulate);
         }
     }
diff --git a/src/c/statisticsFunctions/stdevf/zrowstdevfa.c b/src/c/statisticsFunctions/stdevf/zrowstdevfa.c
index 8a35d4de..bfa77aa5 100644
--- a/src/c/statisticsFunctions/stdevf/zrowstdevfa.c
+++ b/src/c/statisticsFunctions/stdevf/zrowstdevfa.c
@@ -12,8 +12,8 @@
 
 #include "stdevf.h"
 #include "meanf.h"
+#include "abs.h"
 
-#include <stdio.h>
 
 
 void zrowstdevfa(doubleComplex *in1, int lines, int columns, doubleComplex*in2, doubleComplex* out){
@@ -21,15 +21,16 @@ void zrowstdevfa(doubleComplex *in1, int lines, int columns, doubleComplex*in2,
   int j = 0;
   doubleComplex temp = DoubleComplex(0.0,0.0);
   doubleComplex accumulate = DoubleComplex(0.0,0.0);
-  double accumulateFre = 0.0;
+  doubleComplex accumulateFre = DoubleComplex(0.0,0.0);
     
   zrowmeanfa(in1, lines, columns, in2, out );
 
   /*we first multiply each zell of the input matrix by its coefficient*/
   for (j = 0; j < columns; ++j)
     {
+
       accumulate = DoubleComplex(0.0,0.0);
-      accumulateFre = 0.0;
+      accumulateFre =DoubleComplex(0.0,0.0);
       temp = DoubleComplex(0.0,0.0);
 
       for ( i = 0 ; i < lines; ++i )
@@ -39,7 +40,7 @@ void zrowstdevfa(doubleComplex *in1, int lines, int columns, doubleComplex*in2,
          temp = zmuls( in2[lines*j + i] , temp);
 
          accumulate = zadds( temp , accumulate);
-         accumulateFre +=  zreals(in2[lines*j + i]);
+         accumulateFre = zadds(in2[lines*j + i] ,accumulateFre );
 
         } 
 
@@ -49,8 +50,13 @@ void zrowstdevfa(doubleComplex *in1, int lines, int columns, doubleComplex*in2,
         }
       else
         {
-        accumulate = DoubleComplex(  zreals(accumulate ) / (accumulateFre - 1) , zimags(accumulate)  / (accumulateFre - 1));
+	  if( dabss (zreals(accumulate)) <= 3e-14 ) accumulate = DoubleComplex(dabss(zreals(accumulate)) ,zimags(accumulate));
+	  if( dabss (zimags(accumulate)) <= 3e-14 ) accumulate = DoubleComplex(zreals(accumulate) ,dabss(zimags(accumulate)));
+
+
+        accumulate = zrdivs(accumulate, zdiffs(accumulateFre ,DoubleComplex(1.0,0.0)));
         out[j] = zsqrts(accumulate);
+
         }
     }
 
diff --git a/src/c/statisticsFunctions/stdevf/zstdevfa.c b/src/c/statisticsFunctions/stdevf/zstdevfa.c
index f9a2f4c3..2a4f52b7 100644
--- a/src/c/statisticsFunctions/stdevf/zstdevfa.c
+++ b/src/c/statisticsFunctions/stdevf/zstdevfa.c
@@ -12,30 +12,35 @@
 
 #include "stdevf.h" 
 
+
 doubleComplex zstdevfa(doubleComplex *in1, int lines, int columns, doubleComplex* in2)
 {
   int i = 0 ;
+
   doubleComplex temp = DoubleComplex(0.0,0.0);
   doubleComplex accumulate = DoubleComplex(0.0,0.0);
-  double accumulateFre = 0.0 ;
+  doubleComplex accumulateFre = DoubleComplex(0.0,0.0);
   doubleComplex meanf = zmeanfa (in1 , lines*columns , in2);
-
+	
   for(i = 0 ; i < lines*columns ; ++i)
     {
-     temp = zpows (  zdiffs (in1[i] , meanf ) ,DoubleComplex (2.0, 0.0 ) );
+     temp = zpows (  zdiffs (in1[i] , meanf ) ,DoubleComplex (2.0, 0.0 ) );	
      temp = zmuls( in2[i] , temp);
 
      accumulate = zadds( temp , accumulate);
-     accumulateFre +=  zreals(in2[i]);
-    }
+     accumulateFre = zadds(in2[i ] ,accumulateFre );
+
+    } 
 
     if (lines <= 1)
       {
+	
         return zmuls (DoubleComplex(0.0,0.0) , accumulate ) ; 
       }
     else
       {
-        accumulate = DoubleComplex(  zreals(accumulate ) / (accumulateFre - 1) , zimags(accumulate)  / (accumulateFre - 1));
-        return zsqrts(accumulate);
+
+       accumulate = zrdivs(accumulate, zdiffs(accumulateFre ,DoubleComplex(1.0,0.0) ));
+       return zsqrts(accumulate);
       }
 }
diff --git a/src/c/statisticsFunctions/stdevf/zstdevfdz.c b/src/c/statisticsFunctions/stdevf/zstdevfdz.c
new file mode 100644
index 00000000..5e731beb
--- /dev/null
+++ b/src/c/statisticsFunctions/stdevf/zstdevfdz.c
@@ -0,0 +1,24 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008-2008 - INRIA - Allan SIMON
+ *
+ *  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 "stdevf.h" 
+#include "zeros.h"
+#include <malloc.h>
+
+doubleComplex zstdevfdz(double *in1, int lines, int columns, doubleComplex* in2)
+{
+    double* ZEROS;
+    ZEROS=malloc((unsigned int)(lines*columns*sizeof(double)));
+    dzerosa(ZEROS,lines,columns);
+    
+    return zstdevfa(DoubleComplexMatrix(in1,ZEROS,lines*columns), lines,columns , in2);
+}
diff --git a/src/c/statisticsFunctions/stdevf/zstdevfzd.c b/src/c/statisticsFunctions/stdevf/zstdevfzd.c
new file mode 100644
index 00000000..f7a28e30
--- /dev/null
+++ b/src/c/statisticsFunctions/stdevf/zstdevfzd.c
@@ -0,0 +1,41 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008-2008 - INRIA - Allan SIMON
+ *
+ *  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 "stdevf.h" 
+
+doubleComplex zstdevfzd(doubleComplex *in1, int lines, int columns, double* in2)
+{
+  int i = 0 ;
+  doubleComplex temp = DoubleComplex(0.0,0.0);
+  doubleComplex accumulate = DoubleComplex(0.0,0.0);
+  double accumulateFre = 0.0 ;
+  doubleComplex meanf = zmeanfzd (in1 , lines , columns , in2);
+/*equivalent to  (in1 - meanf(x , in2 )).^2 .*in2   */
+  for(i = 0 ; i < lines*columns ; ++i)
+    {
+     temp = zpows (  zdiffs (in1[i] , meanf ) ,DoubleComplex (2.0, 0.0 ) );
+     temp = zmuls( DoubleComplex(in2[i],0.0f) , temp);
+
+     accumulate = zadds( temp , accumulate);
+     accumulateFre +=  in2[i];
+    }
+
+    if (lines <= 1)
+      {
+        return zmuls (DoubleComplex(0.0,0.0) , accumulate ) ; 
+      }
+    else
+      {
+        accumulate = DoubleComplex(  zreals(accumulate ) / (accumulateFre - 1) , zimags(accumulate)  / (accumulateFre - 1));
+        return zsqrts(accumulate);
+      }
+}
diff --git a/tests/unit_tests/test_double_Stdevf.sci b/tests/unit_tests/test_double_Stdevf.sci
new file mode 100644
index 00000000..2798f98d
--- /dev/null
+++ b/tests/unit_tests/test_double_Stdevf.sci
@@ -0,0 +1,45 @@
+
+function test_double_Stdevf()
+	disp (stdevf (4 , 5 )) ;	
+	disp ("*******1*******");
+	disp (stdevf ([4 5 6], [1 2 3]));
+	disp ("*******2*******");
+	disp (stdevf (4 , 5*%i )) ;	
+	disp ("*******3*******");
+	disp (stdevf ([4 5 6], [1 2 3]*%i));
+	disp ("*******4*******");
+
+
+
+	disp (stdevf ([4 5 6 ; 7 8 9], [1 2 3 ; 4 5 6],1));
+	disp ("*******5*******");
+	disp (stdevf ([4 7; 5 8; 6 9], [1 4; 2 5; 3 6],2));
+// complex
+	disp ("*******6*******");
+	disp (stdevf (4 +4*%i , 5 )) ;	
+	disp ("*******7*******");
+	disp (stdevf ([4*%i, 5+6*%i, 6], [1 2 3]));
+        disp ("*******8*******");
+	disp (stdevf (4 +4*%i , 5*%i )) ;	
+	disp ("*******9*******");
+	disp (stdevf ([4*%i, 5+6*%i, 6], [1 2 3]*%i));
+	disp ("******10*******");
+
+
+	disp (stdevf ([4 5 6 ; 7 8 9]*%i, [1 2 3 ; 4 5 6],1));
+	disp ("******11*******");
+	disp (stdevf ([4 7; 5 8; 6 9], [1 4; 2 5; 3 6]*%i,2));
+	disp (stdevf ([4 5 6 ; 7 8 9],[1 2 3 ; 4 8 9]  + [1 2 3 ; 4 5 6]*%i,1));
+	disp ("******12*******");
+	disp (stdevf ([4 5 6 ; 7 8 9], [1 2 3 ; 4 5 6]*%i,1));
+	disp ("******13*******");
+	disp (stdevf ([4 7; 5 8; 6 9]*%i, [1 4; 2 5; 3 6],2));
+	disp ("******14*******");
+
+	disp (stdevf ([4 5 6 ; 7 8 9]*%i, [1 2 3 ; 4 5 6]*%i,1));
+	disp ("******15*******");
+	disp (stdevf ([4 7; 5 8; 6 9]*%i, [1 4; 2 5; 3 6]*%i,2));
+	disp ("******16*******");
+//
+
+endfunction