clean up some files

corrected an error in ccolumnmeanfa
corrected the test in testFloatMeanfa

5 files changed, 74 insertions, 58 deletions

diff --git a/src/c/statisticsFunctions/meanf/ccolumnmeanfa.c b/src/c/statisticsFunctions/meanf/ccolumnmeanfa.c
index 76eb2eea..6e811262 100644
--- a/src/c/statisticsFunctions/meanf/ccolumnmeanfa.c
+++ b/src/c/statisticsFunctions/meanf/ccolumnmeanfa.c
@@ -14,25 +14,25 @@
 #include "sum.h"
 #include "multiplication.h"
 
-void ccolumnmeanfa(floatComplex *in1, int lines, int columns, floatComplex *out,floatComplex *in2) {
+void ccolumnmeanfa(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 < lines ; ++j)
   {
-    tempCoefSum =  FloatComplex(0.0f,0.0f);
+    tempCoefSum =  0.0f;
     out[j]=  FloatComplex(0.0f,0.0f);
     for ( i = 0 ; i < columns; ++i )
       {
 
-        tempMul = cmuls ( in1[lines*i + j] , in2[lines*i + j]); // we times by the coefficient
-        tempCoefSum = cadds ( in2[lines*i + j] ,tempCoefSum ) ;
+        tempMul = cmuls ( in1[lines*i + j] , in2[lines*i + j]); /* we times by the coefficient*/
+        tempCoefSum += creals(in2[lines*i + j]) ;
         out[j] = cadds (tempMul, out[j]) ;
       } 
 
-    out[j] = crdivs(out[j] ,tempCoefSum) ;
+    out[j] =FloatComplex(creals(out[j]) /tempCoefSum, cimags(out[j]) /tempCoefSum);
   }
 }
 
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);
   }
 }
diff --git a/src/c/statisticsFunctions/meanf/dcolumnmeanfa.c b/src/c/statisticsFunctions/meanf/dcolumnmeanfa.c
index 642807ce..e8d25a33 100644
--- a/src/c/statisticsFunctions/meanf/dcolumnmeanfa.c
+++ b/src/c/statisticsFunctions/meanf/dcolumnmeanfa.c
@@ -11,7 +11,6 @@
  */
 
 #include "meanf.h"
-#include <stdio.h>
 
 void dcolumnmeanfa(double *in1, int lines , int columns, double *in2, double *out) {
   int i = 0;
diff --git a/src/c/statisticsFunctions/meanf/srowmeanfa.c b/src/c/statisticsFunctions/meanf/srowmeanfa.c
index 6404ae10..3cdacb82 100644
--- a/src/c/statisticsFunctions/meanf/srowmeanfa.c
+++ b/src/c/statisticsFunctions/meanf/srowmeanfa.c
@@ -13,6 +13,7 @@
 #include "meanf.h"
 #include "sum.h"
 
+
 void srowmeanfa(float *in1, int lines, int columns, float *in2, float *out) {
   int i = 0;
   int j = 0;
@@ -26,7 +27,7 @@ void srowmeanfa(float *in1, int lines, int columns, float *in2, float *out) {
     for ( i = 0 ; i < lines; ++i )
       {
 
-        tempMul = in1[lines*j + i] * in2[lines*j + i]; // we times by the coefficient
+        tempMul = in1[lines*j + i] * in2[lines*j + i]; /* we times by the coefficient*/
         tempCoefSum +=in2[lines*j + i];
         out[j] +=tempMul ;
       } 
diff --git a/src/c/statisticsFunctions/meanf/testFloatMeanf.c b/src/c/statisticsFunctions/meanf/testFloatMeanf.c
index c7f1789f..ead8a86b 100644
--- a/src/c/statisticsFunctions/meanf/testFloatMeanf.c
+++ b/src/c/statisticsFunctions/meanf/testFloatMeanf.c
@@ -20,10 +20,15 @@
 #define ERROR(x) printf("diff = %e\n", x)
 
 static int smeanfsTest(void) {
-  float value1 = 3.0f;
-  float value2 = 1.123456789f;
-  float coef1  = 56.0f ;
-  float coef2  = 2.0f ;
+  float value1 ;
+  float value2 ;
+  float coef1  ;
+  float coef2  ;
+
+  value1 = 3.0f;
+  value2 = 1.123456789f;
+  coef1  = 56.0f ;
+  coef2  = 2.0f ;
 
   printf("\n>>>> Meanf Float Scalar Test\n");
   assert( ( fabs(smeanfs(value1,coef1) ) - (  3.0f ) ) / fabs ( smeanfs(value1,coef1)  ) <  1e-6 );
@@ -74,9 +79,9 @@ static int scolumnmeanfaTest(void) {
     [ 7 8 9 ]  [ 5  6  7 ]
   */
   scolumnmeanfa(table1, 3, 3, coef1 ,columnMeanmedTable1_3_3);
-  assert( ( fabs(columnMeanmedTable1_3_3[0] ) - (  68.0f /  33.0f ) ) / fabs ( columnMeanmedTable1_3_3[0]  ) <  1e-6 );
-  assert( ( fabs(columnMeanmedTable1_3_3[1] ) - (  32.0f /   6.0f ) ) / fabs ( columnMeanmedTable1_3_3[1]  ) <  1e-6 );
-  assert( ( fabs(columnMeanmedTable1_3_3[2] ) - (  146.0f / 18.0f ) ) / fabs ( columnMeanmedTable1_3_3[2]  ) <  1e-6 );
+  assert( fabs(columnMeanmedTable1_3_3[0]  - (  68.0f /  33.0f ) ) / fabs ( columnMeanmedTable1_3_3[0]  ) <  1e-6 );
+  assert( fabs(columnMeanmedTable1_3_3[1]  - (  32.0f /   6.0f ) ) / fabs ( columnMeanmedTable1_3_3[1]  ) <  1e-6 );
+  assert( fabs(columnMeanmedTable1_3_3[2]  - (  146.0f / 18.0f ) ) / fabs ( columnMeanmedTable1_3_3[2]  ) <  1e-6 );
 
   /*
     [ 1 2 3 4 5 6 7 8 9 ] => [ 45 ]
@@ -120,11 +125,11 @@ static int scolumnmeanfaTest(void) {
     [ 5 10 ]  [ 15  20 ]
   */
   scolumnmeanfa(table2, 5, 2, coef2 ,columnMeanmedTable2_5_2);
-  assert( ( fabs(columnMeanmedTable2_5_2[0] ) - (  107.0f / 27.0f ) ) / fabs ( columnMeanmedTable2_5_2[0]  ) <  1e-6 );
-  assert( ( fabs(columnMeanmedTable2_5_2[1] ) - (  143.0f / 29.0f ) ) / fabs ( columnMeanmedTable2_5_2[1]  ) <  1e-6 );
-  assert( ( fabs(columnMeanmedTable2_5_2[2] ) - (  183.0f / 31.0f ) ) / fabs ( columnMeanmedTable2_5_2[2]  ) <  1e-6 );
-  assert( ( fabs(columnMeanmedTable2_5_2[3] ) - (  227.0f / 33.0f ) ) / fabs ( columnMeanmedTable2_5_2[3]  ) <  1e-6 );
-  assert( ( fabs(columnMeanmedTable2_5_2[4] ) - (  275.0f / 35.0f ) ) / fabs ( columnMeanmedTable2_5_2[4]  ) <  1e-6 );
+  assert( fabs(columnMeanmedTable2_5_2[0]  - (  107.0f / 27.0f ) ) / fabs ( columnMeanmedTable2_5_2[0]  ) <  1e-6 );
+  assert( fabs(columnMeanmedTable2_5_2[1]  - (  143.0f / 29.0f ) ) / fabs ( columnMeanmedTable2_5_2[1]  ) <  1e-6 );
+  assert( fabs(columnMeanmedTable2_5_2[2]  - (  183.0f / 31.0f ) ) / fabs ( columnMeanmedTable2_5_2[2]  ) <  1e-6 );
+  assert( fabs(columnMeanmedTable2_5_2[3]  - (  227.0f / 33.0f ) ) / fabs ( columnMeanmedTable2_5_2[3]  ) <  1e-6 );
+  assert( fabs(columnMeanmedTable2_5_2[4]  - (  275.0f / 35.0f ) ) / fabs ( columnMeanmedTable2_5_2[4]  ) <  1e-6 );
   for ( i = 0 ; i < 5 ; ++i) {
     printf("columnMeanmedTable2_5_2[%d] = %e\n", i, columnMeanmedTable2_5_2[i]);
   }
@@ -143,7 +148,7 @@ static int srowmeanfaTest(void) {
   float rowMeanmedTable1_1_9[9] = {0};
   float rowMeanmedTable1_9_1[1] = {0};
   float rowMeanmedTable2_2_5[5] = {0};
-  float rowMeanmedTable2_5_2[2] = {0};
+  /*float rowMeanmedTable2_5_2[2] = {0};*/
 
   printf("\n>>>> Row Mean Float Array Test\n");
   /*
@@ -155,9 +160,9 @@ static int srowmeanfaTest(void) {
   for ( i = 0 ; i < 3 ; ++i) {
     printf("rowMeanmedTable1_3_3[%d] = %e\n", i, rowMeanmedTable1_3_3[i]);
   }
-  assert( ( fabs(rowMeanmedTable1_3_3[0] ) - (  49.0f / 16.0f ) ) / fabs ( rowMeanmedTable1_3_3[0]  ) <  1e-6 );
-  assert( ( fabs(rowMeanmedTable1_3_3[1] ) - (  80.0f / 19.0f ) ) / fabs ( rowMeanmedTable1_3_3[1]  ) <  1e-6 );
-  assert( ( fabs(rowMeanmedTable1_3_3[2] ) - ( 117.0f / 22.0f ) ) / fabs ( rowMeanmedTable1_3_3[2]  ) <  1e-6 );
+  assert( fabs(rowMeanmedTable1_3_3[0]  - (  49.0f / 16.0f ) ) / fabs ( rowMeanmedTable1_3_3[0]  ) <  1e-6 );
+  assert( fabs(rowMeanmedTable1_3_3[1]  - (  80.0f / 19.0f ) ) / fabs ( rowMeanmedTable1_3_3[1]  ) <  1e-6 );
+  assert( fabs(rowMeanmedTable1_3_3[2]  - ( 117.0f / 22.0f ) ) / fabs ( rowMeanmedTable1_3_3[2]  ) <  1e-6 );
 
   /*
     [ 1 ]
@@ -190,27 +195,32 @@ static int srowmeanfaTest(void) {
   for ( i = 0 ; i < 5 ; ++i) {
     printf("rowMeanmedTable2_2_5[%d] = %e\n", i, rowMeanmedTable2_2_5[i]);
   }
-  assert( ( fabs(rowMeanmedTable2_2_5[0] ) - (  35.0f / 23.0f ) ) / fabs ( rowMeanmedTable2_2_5[0]  ) <  1e-6 );
-  assert( ( fabs(rowMeanmedTable2_2_5[1] ) - (  95.0f / 27.0f ) ) / fabs ( rowMeanmedTable2_2_5[1]  ) <  1e-6 );
-  assert( ( fabs(rowMeanmedTable2_2_5[2] ) - ( 171.0f / 31.0f ) ) / fabs ( rowMeanmedTable2_2_5[2]  ) <  1e-6 );
-  assert( ( fabs(rowMeanmedTable2_2_5[3] ) - ( 263.0f / 35.0f ) ) / fabs ( rowMeanmedTable2_2_5[3]  ) <  1e-6 );
-  assert( ( fabs(rowMeanmedTable2_2_5[4] ) - ( 371.0f / 39.0f ) ) / fabs ( rowMeanmedTable2_2_5[4]  ) <  1e-6 );
+  assert( fabs(rowMeanmedTable2_2_5[0]  - (  35.0f / 23.0f ) ) / fabs ( rowMeanmedTable2_2_5[0]  ) <  1e-6 );
+  assert( fabs(rowMeanmedTable2_2_5[1]  - (  95.0f / 27.0f ) ) / fabs ( rowMeanmedTable2_2_5[1]  ) <  1e-6 );
+  assert( fabs(rowMeanmedTable2_2_5[2]  - ( 171.0f / 31.0f ) ) / fabs ( rowMeanmedTable2_2_5[2]  ) <  1e-6 );
+  assert( fabs(rowMeanmedTable2_2_5[3]  - ( 263.0f / 35.0f ) ) / fabs ( rowMeanmedTable2_2_5[3]  ) <  1e-6 );
+  assert( fabs(rowMeanmedTable2_2_5[4]  - ( 371.0f / 39.0f ) ) / fabs ( rowMeanmedTable2_2_5[4]  ) <  1e-6 );
 
 
   return 0;
 }
 
 static int cmeanfsTest(void) {
-  floatComplex value1 = FloatComplex(3.0f, 3.0f);
-  floatComplex coef1 = FloatComplex(3.0f, 0.0f);
-  floatComplex value2 = FloatComplex(1.123456789, 1.123456789);
-  floatComplex coef2 = FloatComplex(9.0f, 0.0f);
+  floatComplex value1 ;
+  floatComplex coef1 ;
+  floatComplex value2 ;
+  floatComplex coef2 ;
+
+  value1 = FloatComplex(3.0f, 3.0f);
+  coef1  = FloatComplex(3.0f, 0.0f);
+  value2 = FloatComplex(1.123456789f, 1.123456789f);
+  coef2  = FloatComplex(9.0f, 0.0f);
 
   printf("\n>>>> Mean Float Complex Scalar Test\n");
-  assert( ( fabs(creals(cmeanfs(value1,coef1)) ) - (  3.0f ) ) / fabs ( creals(cmeanfs(value1,coef1))  ) <  1e-6 );
-  assert( ( fabs(cimags(cmeanfs(value1,coef1)) ) - (  3.0f ) ) / fabs ( cimags(cmeanfs(value1,coef1))  ) <  1e-6 );
-  assert( ( fabs(creals(cmeanfs(value2,coef1)) ) - (  1.123456789 ) ) / fabs ( creals(cmeanfs(value2,coef1))  ) <  1e-6 );
-  assert( ( fabs(cimags(cmeanfs(value2,coef1)) ) - (  1.123456789 ) ) / fabs ( cimags(cmeanfs(value2,coef1))  ) <  1e-6 );
+  assert(  fabs(creals(cmeanfs(value1,coef1))  - (  3.0f ) ) / fabs ( creals(cmeanfs(value1,coef1))  ) <  1e-6 );
+  assert(  fabs(cimags(cmeanfs(value1,coef1))  - (  3.0f ) ) / fabs ( cimags(cmeanfs(value1,coef1))  ) <  1e-6 );
+  assert(  fabs(creals(cmeanfs(value2,coef1))  - (  1.123456789f ) ) / fabs ( creals(cmeanfs(value2,coef1))  ) <  1e-6 );
+  assert(  fabs(cimags(cmeanfs(value2,coef1))  - (  1.123456789f ) ) / fabs ( cimags(cmeanfs(value2,coef1))  ) <  1e-6 );
   return 0;
 }
 
@@ -230,8 +240,8 @@ static int cmeanfaTest(void) {
   printf("\n>>>> Mean Float Complex Array Test\n");
   printf("\tMeanmedTable_1_1[0] = %e + %ei\n", creals(cmeanfa(table1, 9, coef1)), cimags(cmeanfa(table1, 9, coef1)));
 
-  assert( fabs( (creals(cmeanfa(table1, 9, coef1)) ) - (  246.0f /  57.0f ) ) / fabs ( creals(cmeanfa(table1, 9, coef1))  ) <  1e-6 );
-  assert( fabs( (cimags(cmeanfa(table1, 9, coef1)) ) - (  288.0f /  57.0f ) ) / fabs ( cimags(cmeanfa(table1, 9, coef1))  ) <  1e-6 );/*
+  assert( fabs( creals(cmeanfa(table1, 9, coef1))  - (  246.0f /  57.0f ) ) / fabs ( creals(cmeanfa(table1, 9, coef1))  ) <  1e-6 );
+  assert( fabs( cimags(cmeanfa(table1, 9, coef1))  - (  288.0f /  57.0f ) ) / fabs ( cimags(cmeanfa(table1, 9, coef1))  ) <  1e-6 );/*
 */
   return 0;
 }
@@ -256,7 +266,10 @@ static int crowmeanfaTest(void) {
   floatComplex* table2 = FloatComplexMatrix (tableR2, tableI2, 9);
   floatComplex* coef2 = FloatComplexMatrix (coefR2, coefI2, 9);
 */
-  floatComplex rowMeanmedTable1_3_3[3] = {FloatComplex(0.0f, 0.0f)};
+  floatComplex rowMeanmedTable1_3_3[3] ;
+  rowMeanmedTable1_3_3[0] = FloatComplex(0.0f, 0.0f);
+  rowMeanmedTable1_3_3[1] = FloatComplex(0.0f, 0.0f);
+  rowMeanmedTable1_3_3[2] = FloatComplex(0.0f, 0.0f);
  /*
   floatComplex rowMeanmedTable1_1_9[9] = {FloatComplex(0.0f, 0.0f)};
   floatComplex rowMeanmedTable1_9_1[1] = {FloatComplex(0.0f, 0.0f)};
@@ -275,14 +288,14 @@ static int crowmeanfaTest(void) {
 
       printf("rowMeanmedTable_3_3[%d] = %e + %ei\n", i, creals(rowMeanmedTable1_3_3[i]), cimags(rowMeanmedTable1_3_3[i]));
     }
-  assert( ( fabs(cimags(rowMeanmedTable1_3_3[0]) ) - ( 27.0f / 16.0f ) ) / fabs ( cimags(rowMeanmedTable1_3_3[0])  ) <  1e-6 );
-  assert( ( fabs(creals(rowMeanmedTable1_3_3[0]) ) - ( 49.0f / 16.0f ) ) / fabs ( creals(rowMeanmedTable1_3_3[0])  ) <  1e-6 );
+  assert( fabs(cimags(rowMeanmedTable1_3_3[0])  - ( 27.0f / 16.0f ) ) / fabs ( cimags(rowMeanmedTable1_3_3[0])  ) <  1e-6 );
+  assert( fabs(creals(rowMeanmedTable1_3_3[0])  - ( 49.0f / 16.0f ) ) / fabs ( creals(rowMeanmedTable1_3_3[0])  ) <  1e-6 );
   
-  assert( ( fabs(cimags(rowMeanmedTable1_3_3[1]) ) - ( 90.0f / 19.0f ) ) / fabs ( cimags(rowMeanmedTable1_3_3[1])  ) <  1e-6 );
-  assert( ( fabs(creals(rowMeanmedTable1_3_3[1]) ) - ( 80.0f / 19.0f ) ) / fabs ( creals(rowMeanmedTable1_3_3[1])  ) <  1e-6 );
+  assert( fabs(cimags(rowMeanmedTable1_3_3[1])  - ( 90.0f / 19.0f ) ) / fabs ( cimags(rowMeanmedTable1_3_3[1])  ) <  1e-6 );
+  assert( fabs(creals(rowMeanmedTable1_3_3[1])  - ( 80.0f / 19.0f ) ) / fabs ( creals(rowMeanmedTable1_3_3[1])  ) <  1e-6 );
 
-  assert( ( fabs(cimags(rowMeanmedTable1_3_3[2]) ) - (171.0f / 22.0f ) ) / fabs ( cimags(rowMeanmedTable1_3_3[2])  ) <  1e-6 );
-  assert( ( fabs(creals(rowMeanmedTable1_3_3[2]) ) - (117.0f / 22.0f ) ) / fabs ( creals(rowMeanmedTable1_3_3[2])  ) <  1e-6 );
+  assert( fabs(cimags(rowMeanmedTable1_3_3[2])  - (171.0f / 22.0f ) ) / fabs ( cimags(rowMeanmedTable1_3_3[2])  ) <  1e-6 );
+  assert( fabs(creals(rowMeanmedTable1_3_3[2])  - (117.0f / 22.0f ) ) / fabs ( creals(rowMeanmedTable1_3_3[2])  ) <  1e-6 );
 
 
 
@@ -311,7 +324,10 @@ static int ccolumnmeanfaTest(void) {
   floatComplex* table2 = FloatComplexMatrix (tableR2, tableI2, 9);
   floatComplex* coef2 = FloatComplexMatrix (coefR2, coefI2, 9);
 */
-  floatComplex columnMeanmedTable1_3_3[3] = {FloatComplex(0.0f, 0.0f)};
+  floatComplex columnMeanmedTable1_3_3[3] ;
+  columnMeanmedTable1_3_3[0] = FloatComplex(0.0f, 0.0f);
+  columnMeanmedTable1_3_3[1] = FloatComplex(0.0f, 0.0f);
+  columnMeanmedTable1_3_3[2] = FloatComplex(0.0f, 0.0f);
  /*
   floatComplex rowMeanmedTable1_1_9[9] = {FloatComplex(0.0f, 0.0f)};
   floatComplex rowMeanmedTable1_9_1[1] = {FloatComplex(0.0f, 0.0f)};
@@ -329,14 +345,14 @@ static int ccolumnmeanfaTest(void) {
     {
       printf("columnMeanmedTable_3_3[%d] = %e + %ei\n", i, creals(columnMeanmedTable1_3_3[i]), cimags(columnMeanmedTable1_3_3[i]));
     }
-  assert( ( fabs(cimags(columnMeanmedTable1_3_3[0]) ) - (138.0f / 33.0f ) ) / fabs ( cimags(columnMeanmedTable1_3_3[0])  ) <  1e-6 );
-  assert( ( fabs(creals(columnMeanmedTable1_3_3[0]) ) - ( 68.0f / 33.0f ) ) / fabs ( creals(columnMeanmedTable1_3_3[0])  ) <  1e-6 );
+  assert( fabs(cimags(columnMeanmedTable1_3_3[0])  - (138.0f / 33.0f ) ) / fabs ( cimags(columnMeanmedTable1_3_3[0])  ) <  1e-6 );
+  assert( fabs(creals(columnMeanmedTable1_3_3[0])  - ( 68.0f / 33.0f ) ) / fabs ( creals(columnMeanmedTable1_3_3[0])  ) <  1e-6 );
   
-  assert( ( fabs(cimags(columnMeanmedTable1_3_3[1]) ) - ( 36.0f /  6.0f ) ) / fabs ( cimags(columnMeanmedTable1_3_3[1])  ) <  1e-6 );
-  assert( ( fabs(creals(columnMeanmedTable1_3_3[1]) ) - ( 32.0f /  6.0f ) ) / fabs ( creals(columnMeanmedTable1_3_3[1])  ) <  1e-6 );
+  assert( fabs(cimags(columnMeanmedTable1_3_3[1])  - ( 36.0f /  6.0f ) ) / fabs ( cimags(columnMeanmedTable1_3_3[1])  ) <  1e-6 );
+  assert( fabs(creals(columnMeanmedTable1_3_3[1])  - ( 32.0f /  6.0f ) ) / fabs ( creals(columnMeanmedTable1_3_3[1])  ) <  1e-6 );
 
-  assert( ( fabs(cimags(columnMeanmedTable1_3_3[2]) ) - (114.0f / 18.0f ) ) / fabs ( cimags(columnMeanmedTable1_3_3[2])  ) <  1e-6 );
-  assert( ( fabs(creals(columnMeanmedTable1_3_3[2]) ) - (146.0f / 18.0f ) ) / fabs ( creals(columnMeanmedTable1_3_3[2])  ) <  1e-6 );
+  assert( fabs(cimags(columnMeanmedTable1_3_3[2])  - (114.0f / 18.0f ) ) / fabs ( cimags(columnMeanmedTable1_3_3[2])  ) <  1e-6 );
+  assert( fabs(creals(columnMeanmedTable1_3_3[2])  - (146.0f / 18.0f ) ) / fabs ( creals(columnMeanmedTable1_3_3[2])  ) <  1e-6 );