1 files changed, 5 insertions, 5 deletions
diff --git a/src/c/statisticsFunctions/meanf/crowmeanfa.c b/src/c/statisticsFunctions/meanf/crowmeanfa.c
index 309aa930..4c2b4866 100644
--- a/src/c/statisticsFunctions/meanf/crowmeanfa.c
+++ b/src/c/statisticsFunctions/meanf/crowmeanfa.c
@@ -16,21 +16,21 @@
 void crowmeanfa(floatComplex *in1, int lines, int columns, floatComplex *in2, floatComplex *out) {
   int i = 0;
   int j = 0;
-  floatComplex tempCoefSum = FloatComplex(0.0f,0.0f);
+  float tempCoefSum = 0.0f;
   floatComplex tempMul = FloatComplex(0.0f,0.0f);
   /*we first multiply each cell of the input matrix by its coefficient*/
   for (j = 0; j < columns; ++j)
   {
-    tempCoefSum =  FloatComplex(0.0f,0.0f);
+    tempCoefSum =  0.0f;
     out[j]=  FloatComplex(0.0f,0.0f);
     for ( i = 0 ; i < lines; ++i )
       {
 
-        tempMul = cmuls ( in1[lines*j + i] , in2[lines*j + i]); // we times by the coefficient
-        tempCoefSum = cadds ( in2[lines*j + i] ,tempCoefSum ) ;
+        tempMul = cmuls ( in1[lines*j + i] , in2[lines*j + i]); /* we times by the coefficient*/
+        tempCoefSum += creals(in2[lines*j + i]) ;
         out[j] = cadds (tempMul, out[j]) ;
       } 
 
-    out[j] = crdivs(out[j] ,tempCoefSum) ;
+    out[j] = FloatComplex(creals(out[j]) /tempCoefSum, cimags(out[j]) /tempCoefSum);
   }
 }