added stdevf functions , interfaces , header, tests , everything

1 files changed, 12 insertions, 7 deletions

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);
       }
 }