summaryrefslogtreecommitdiff
path: root/src/c/statisticsFunctions/meanf/crowmeanfa.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/c/statisticsFunctions/meanf/crowmeanfa.c')
-rw-r--r--src/c/statisticsFunctions/meanf/crowmeanfa.c8
1 files changed, 4 insertions, 4 deletions
diff --git a/src/c/statisticsFunctions/meanf/crowmeanfa.c b/src/c/statisticsFunctions/meanf/crowmeanfa.c
index 4c2b4866..14dc3271 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;
- float tempCoefSum = 0.0f;
+ floatComplex tempCoefSum = FloatComplex(0.0f,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 = 0.0f;
+ tempCoefSum = FloatComplex(0.0f,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 += creals(in2[lines*j + i]) ;
+ tempCoefSum = cadds ( in2[lines*j + i] ,tempCoefSum ) ;
out[j] = cadds (tempMul, out[j]) ;
}
- out[j] = FloatComplex(creals(out[j]) /tempCoefSum, cimags(out[j]) /tempCoefSum);
+ out[j] = crdivs(out[j] ,tempCoefSum);
}
}