diff options
Diffstat (limited to 'src/c/statisticsFunctions/meanf/crowmeanfa.c')
-rw-r--r-- | src/c/statisticsFunctions/meanf/crowmeanfa.c | 8 |
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); } } |