Traitement-signal-plantes/Code-C/power.c

42 lines
1.2 KiB
C
Raw Normal View History

2022-06-01 12:32:04 +02:00
#include "power.h"
#include "getArray.h"
2022-06-02 14:29:21 +02:00
#include "fileGestion.h"
2022-06-10 17:22:01 +02:00
#include "initialParameters.h"
#include "queue.h"
2022-06-01 12:32:04 +02:00
2022-06-14 10:46:12 +02:00
/**
2022-06-15 11:57:15 +02:00
* @brief realize the powerThreadFunction calcul
2022-06-14 10:46:12 +02:00
*
* @param p array with all the values that will be used for the calcul
* @param powerArray array where results are stocked
*/
2022-06-14 18:01:52 +02:00
void powerCalculation(long **p, double powerArray[]){
for(int i = 1; i < nCol; i++){
2022-06-01 12:32:04 +02:00
int j = 0;
2022-06-09 17:34:24 +02:00
powerArray[i] = 0;
2022-06-14 18:01:52 +02:00
while(j < nRow-1){
double aire = ( pow(p[j][i],2) + pow(p[j+1][i],2) ) / 2 * period;
2022-06-09 17:34:24 +02:00
powerArray[i] += aire;
2022-06-01 12:32:04 +02:00
j++;
}
2022-06-14 18:01:52 +02:00
powerArray[i] *= invTimeBandWidth;
2022-06-01 12:32:04 +02:00
}
}
2022-06-14 10:46:12 +02:00
/**
* @brief function that realize all the action to write one lign in the file powerData.csv
*
* @param rawDataFileName name of the raw data file to use to realize the calcul
* @param N number of rows in the file
* @param M number of columns in the file
*/
2022-06-15 11:57:15 +02:00
double *powerFunction(char* rawDataFileName){
2022-06-14 18:01:52 +02:00
long **p = getRawDataArray(rawDataFileName);
double pw[nCol-1];
2022-06-07 14:28:32 +02:00
if(p !=NULL){
2022-06-14 18:01:52 +02:00
powerCalculation(p,pw);
appendDataInFile("powerData.csv",pw,nCol-1);
freeArray(p,nRow);
2022-06-07 14:28:32 +02:00
}
2022-06-15 11:57:15 +02:00
return pw;
2022-06-01 12:32:04 +02:00
}