From 8ecdd376191c82b09472dc15cf72c549b223ea9c Mon Sep 17 00:00:00 2001
From: "quentin.perret" <quentinprrt33@gmail.com>
Date: Thu, 16 Jun 2022 17:05:12 +0200
Subject: [PATCH] growth ok

---
 Code-C/average.c           |  22 ++++++++++++++--------
 Code-C/average.h           |   2 +-
 Code-C/getArray.c          |   4 ++--
 Code-C/growthRate.c        |  15 +++++++++++++--
 Code-C/growthRate.h        |   2 +-
 Code-C/initialParameters.h |   3 ++-
 Code-C/main                | Bin 41760 -> 42008 bytes
 Code-C/main.c              |  30 +++++++++++++++++++-----------
 Code-C/power.c             |  12 +++++-------
 Code-C/simulateFlux.c      |   2 +-
 10 files changed, 58 insertions(+), 34 deletions(-)

diff --git a/Code-C/average.c b/Code-C/average.c
index e9d8706..8fef657 100644
--- a/Code-C/average.c
+++ b/Code-C/average.c
@@ -13,12 +13,11 @@ void averageCalculation(long **p, double averageArray[]){
     for(int i = 1; i < nCol; i++){
         int j = 0;
         averageArray[i] = 0;
-        while(j < nRow -1){
+        while(j < nRowRawData){
             averageArray[i] += p[i][j];
             j++;
         }
-        averageArray[i] /= nRow;
-        //printf("%f\n", powerArray[i]);
+        averageArray[i] /= nRowRawData;
     }
 }
 
@@ -27,12 +26,19 @@ void averageCalculation(long **p, double averageArray[]){
  * 
  * @param rawDataFileName name of the raw data file to use to realize the calcul
  */
-void averageThreadFunction(char* rawDataFileName){
+void averageFunction(char* rawDataFileName , double **aver){
     long **p = getRawDataArray(rawDataFileName);
-    double aver[nCol -1];
+    double averN[nCol -1];
     if(p !=NULL){
-        averageCalculation(p,aver);
-        appendDataInFile("averageData.csv",aver,nCol-1);
-        freeArray(p,nRow);
+        if(p !=NULL){
+            if(aver == NULL){
+                averageCalculation(p,averN);
+                appendDataInFile("averageData.csv",averN,nCol-1);
+            }else{
+                averageCalculation(p,aver[1]);
+                appendDataInFile("averageData.csv",aver[1],nCol-1);
+            }
+        freeArray(p,nRowRawData);
+    }
     }
 }
\ No newline at end of file
diff --git a/Code-C/average.h b/Code-C/average.h
index b130865..a121de4 100644
--- a/Code-C/average.h
+++ b/Code-C/average.h
@@ -1,4 +1,4 @@
 #include <math.h>
 #include <stdbool.h>
 
-void averageThreadFunction(char* rawDataFileName);
\ No newline at end of file
+void averageFunction(char* rawDataFileName , double **aver);
\ No newline at end of file
diff --git a/Code-C/getArray.c b/Code-C/getArray.c
index 7284378..f8c05ef 100644
--- a/Code-C/getArray.c
+++ b/Code-C/getArray.c
@@ -97,8 +97,8 @@ void freeArray(long **p, int N) {
  */
 long **getRawDataArray(char* rawDataFileName){
     long **p;
-    p = getlongArray(nRow, nCol);
-    fillArrayWithRawData(rawDataFileName,p ,nRow, nCol);
+    p = getlongArray(nRowRawData, nCol);
+    fillArrayWithRawData(rawDataFileName,p ,nRowRawData, nCol);
     //if(checkArrayFullyFill(p,nRow)){
         //clearRawData(N);
         return p;
diff --git a/Code-C/growthRate.c b/Code-C/growthRate.c
index 8fcca10..64a0a2e 100644
--- a/Code-C/growthRate.c
+++ b/Code-C/growthRate.c
@@ -18,8 +18,19 @@ void growthRateCalculation(double **dataArray , double *gRateArray){
  * 
  * @param dataLign 
  */
-void growthRateFunction(double **dataLign){
+void growthRateFunction(double **dataLign , char* fileName){
     double gRateArray[nCol-1];
     growthRateCalculation(dataLign , gRateArray);
-    appendDataInFile("growthRate.csv",gRateArray , nCol-1);
+    appendDataInFile(fileName,gRateArray , nCol-1);
+}
+
+void sumColArray(long **p, double res[] , int N , int M){
+    for(int i = 1; i < M; i++){
+        int j = 0;
+        res[i] = 0;
+        while(j < N){
+            res[i] += p[i][j];
+            j++;
+        }
+    }
 }
\ No newline at end of file
diff --git a/Code-C/growthRate.h b/Code-C/growthRate.h
index 7979f21..d68b096 100644
--- a/Code-C/growthRate.h
+++ b/Code-C/growthRate.h
@@ -1,4 +1,4 @@
 #include <math.h>
 #include <stdbool.h>
 
-void growthRateFunction(double **dataLign);
\ No newline at end of file
+void growthRateFunction(double **dataLign , char* fileName);
\ No newline at end of file
diff --git a/Code-C/initialParameters.h b/Code-C/initialParameters.h
index 96201f5..5403af0 100644
--- a/Code-C/initialParameters.h
+++ b/Code-C/initialParameters.h
@@ -2,7 +2,8 @@
 #include "queue.h"
 
 extern bool rawDataWriteFlag;
-extern int nRow;
+extern int nRowRawData;
+extern int nRowGR;
 extern int nCol;
 extern double freqEch;
 
diff --git a/Code-C/main b/Code-C/main
index 3cd197e6963024a135e688b05bc98a32ecf109a8..d9efec68aa3e2c7f0b12ca0fa5ab95076d77fd20 100755
GIT binary patch
delta 10306
zcmai433yaRw!U>cxA#sjA?fU$g^(4(PNFPHFhd6e1~3>9Fj2xHyCwuwY#@My5J@aY
zQ6qvfqt1w;nNf%{3U0vYD5#_O#w9wa4G5yh8+5>a|Eap&K)!FjSA3~cXFql7RNcC_
z+5VQg^`P3Gqt-mqk@vZEmroIGmydhc085osqE)<t6|&iUCEMbz9V}@P$P?6;`}}63
z^(4zb)AaGi_rE$i_r$j2RXOi;1=T!f?xqw)@wYiL1sIi`tSF(_7e^_|Pfp)8#2NMT
zXad@g{Ay2*`H(<{z1FqBQ=+jt{){)C`MciqCNMUDe-WR}KTLMnn|}?3{4M1NTea^-
z@H5E^?e`-jfAgw?igR5HpP7;n=nsijZORBmX{!qPTg&(OTi$}J6>*9Z{Nd+N=;SC6
zxBHv-nf<LU3=jHSKJvF52%fu0dPxyIT||-U#iafM)LUcE42MVPNwqAf`CC41_qWv_
zy|HvgDgPv4s(N8KAC%aa&EiuNQ*Dih{VhHIme*UgVf@a-#+0B4aPKfZzzIDt29b9u
zh#R;gDaCT2`~oHVxPgyJ%8uHAAzzb}%g-mpcoy|hlsV3vgLC}O#ruK>ehP(ndVCIF
z5%00?AF3$D<>#H77x22|6!tI=B=>hU2;v3j=F{-n-hPOmj*q8A+7A5^PUK+zeexnE
zljjz(?Htl8dH?ia_GSP3-;?)=NZ!JT-j&JgBy|JSTVqQhaw&Psk<ebrv-687mi*;s
zMRez}vgegXAByNF*yu(fp<ag7V2mZ5MrOa`gHo?ac!{|ECHFsW!BvR6h2NW6Fnuj?
zizW9($$c3hIQ~a?X<l_ydYC~{S4pXMDfKW(^(Cp+@?dNAu@+Zr`7!_2R~SiQxr%uo
z-|R6TLaN&@70<hb^NRq%%|8h9XQcUONvc9hy)UJ{CaF36C0|a$D&kI++z!cofVf%w
zye~6h5^?J!w_0+i6ZbeDnwFlBP27c&n<Tku#C?u0PRq#u?o97go&dA^4k`bxQ^flo
zK=5Y%a#~5^4&vS;xs8&02XQTD_#GL2d2xEW<EzsmBj+gvAMv^A`L`yL(nC_|^Jro6
zUjV_|zSkvNOUxv__A<@o7`GpKz1rXMwZHkhG=IygAQmGiv-JsFcK=QqY?B7>BZKMu
zLi*@24->aba)(Rq)x`Z9nQvazBW*O0)Ssl3OG@RC)K0!K<5Ana4n^4mCY4ZnW_;qq
z5&Y086uKP+_8>^X8&B~sGZOg1%x_p8pVGHy<a2gf{PLqXxc(xChjAL^D+=aUgZ)Ft
z6YAwC&G*rhguhM|Pw2gTN!GN$&<GLR{y;E`P<GV6zsEr61A<$=+uvG!6nZ$m#PRgl
z2}re8ljSh!9g=LV-UHHZU23J2I(!-F_d=?CcV){loKZbDmhy4^4zj;;SN8SnAg|9J
zmh?g{Ef}3aN4lx9UrD5#Ar#AT{^#tBKrS8yKxML49EVc~cVoBEUUW%juss+Gt>4!e
z@VCa6OM8e(w<knluM;wU2xa(NwGt9)*)O)V+)r(zKDHp=Pr|{;A`wrU-vgL;xZ5lO
z5+|gEY_ws*ZPkbUZR5={4jl0<?*&is)j8R0eb-Yt@hrL#N|C6y@ON`_@~xMXrhOX<
zHJ4vt&8se$o%h}&4!*{0MQPs0TC{~cJ#X-|Z%o2OuoR+L)iWpWZjH@?NUQc_Kb*C%
zH7-V;&rqh8d=m;C&@}L<Nxm?&4T2C7TV*<gT3*B4KjF{hWn1Rx*7x&I^HSp1LQv_w
z9a^;p?#>^cwU|y4-Bs?_Vk|7q!TC=KeZl#n=H?6Y6Re>u3c39J$%%YRe&51ZC>#oK
z3~D-T;ZMgdS@?ktTHQit7GIPfAJ`N&aZ@<(z2VbXpNNm#A|GTh;$a<~GnW(gekTd!
ze!yie7w9f`Am}*wTDPz9*rkb=1aJ_b%C#sxj8^TwzSvuDy{cj=!Ns)vEjWy|yDp>F
zNNUf4@}FqJ$J<cG;ovXQ@r{gx@mu6%jyA9Qi8;44L!3s96S0$mkU0wfV#B>Cj;`t+
zqUogLaw`Q<_C8JK^>;F_L`u;~_tEju59TgceOB&O;X$`n4F}ETKe6L)BZdRzXDITM
z^ZeY^e^3&f_UKm>TK28taA2sN@*bYoKP50S6JZBgarv3MHiZM%&320;_6p~x#HmnD
zL1sqDxSy7bL!@+_^nfIZ!dg!Sw|<{<{q5KnSf>K`!#efo>%<@GpB#8TLs6VBM!DAQ
zYjCz|9GIy1O*D}E%5?V!Z)u+ZXFob!j+GbyYR?-8&f5Q@)2;mjv{vmp^mClwANtXm
zcaS7r6_V4$xcM^a=RlXI(8bfC09IAo018fmGdK~j50CiUVpIPU3Poyi%$Jx(8`iZA
zlNapKxTr+fX<)xfX%F`MDioR{Zl$F+RJM%@q$*13)QXmmN~<ec&XrzY+H#@V-`0OA
z$WtnFFXC}!4ywF|zdYc#_YIti$M4c5P4%~3Q`b;8;FjB$D@xas0}rw6xwkE8YM3{<
zuA#oFNzgA5N|#GPa?Y0wjt%VlEEGD7s%`%w6gm&M9dJ9~`kqk8r7Fs7ggOu~AAwc@
zZUt-r)czF;Z3R3BxCd|>Hri3ZlYr*|N8z-HM+(;g76ER-8B_!K368XlfU~f3o-!x_
z#fxM2FbEZ>l{0`FQLt&KjPC$T0Q;jlY5`9Jt_GZjD&G$HTflvQ?*aAz+E9K9c0way
z8lW9jUII81uoiF^;A+6-I4-sW1_AeJ0Y$kQCD#MOAF*GfjEeF(_GtlN9uAvIz)gS)
z0gnQ128=xy3Oxh(+P_1gcK|!iBXB???vySZt#<+r1Z24ARso(wVJrpgxqxd2U<oSW
zC?K^yoDKD{^E`cMAq(<JL$edwu^#l)2{*jEC^PW4h<8A;fIpA^1Ab&^CQIh$(AV>X
zVLmpD7o&fLPeXr_uSB21A3}c(-#g4_EW{d(XZ-VFadCzsW;f~f$@GU<>WX9BIy`ma
zvsk&k<hC^0wLu;0FqStP)aN473W+CUWyfGe2SY6)XXr0%7JxJo{7G2pF~VD{&H)3e
z_*)2eE#EQR$9D4l=-=dDqMyUr2p@Z%r=$NnFGc??uN&b@H=r_ZCDm|1YQ2a?SVE;t
z;5#7u2H%hV9RCu17H361_AXCHpT<kkPvLb%nd~vXz9^ghf$sw4BmQ^var_keYq+D>
zhc~|d=uh(r#hL6OzD$gl6(<GmLQ#xGQEa#(hBc8GI>CP${1aCs<7<%&it#Vt-;SdB
z{&Ia9qg{i)Pr>)2m}U~8EZVh2t%#1l*H9jv))rMBox}65>aW#;&${aD%?5sORE@f`
zuq$Kq3DxodG{~#M9i=a_YM#C(hF7UkiL?f)62#R=u5mHlsL)0;UaH0=NcMDG9woae
zNwSxP4Y;L*5w}$-5#*I+S+-p<Vf+z(S6Nn{-O!W#AM&JyiYRi>^>*1+M$d4pjZw4L
zgX1VgltvTxRur<bu5lI18Ne&!{VeloFrt)sC<rs#0ZIEq(3%M1F|V{(t@7N;$&B;6
zD$@dFQq}A?fTfX1)g;EoRiLvJ<M;7WXFOQcSf>pmQ<>mqP#}hmtpU?yfsZ1dp`$%u
zUaGKEyH{X)0+t(y>$7J9kLq|TU{s>aZW^0cDw7GM2xJzFSy!SP0n}E8K~%?X2-|NZ
zb-OVfJ~^SrIvsFhPa`*KM~!Kw5X*D(>IvmjmBnOJ84H@8g(oniEL>VRY35a&a4J)O
zMZOGC+Wf^0N|c@lyk{$_IchdQ)W4yv+WG%Z7#^jWI~1NbF)_eFGww#mVvXX04MDNZ
z0a@oS=uM{uH-pA%hu@S7CoFbIt3XT4uM;=UZbVdujw+*~nW^+N^EmJ=)FlA5PRGE|
zQDanW?~=UDsL!B4Sq~01)@eE>SXt6hjY(Lx6_Yu&qsm^M$a^OC&2nI>dNzKAJY^$N
zr-+Q&W{nT3#EC|=G5^2Sh7opHhVg0IifChChZ^g&VkB!LBQabX=A}&#9!e}K4GY#R
z(;g=<l*ZB~aQpI$*A)ePWYu8|M<fnHVU?=MkxDAG%zcDdj!U3bj(hZtmI|DjWxm3q
z9aPS0$3h3!xgo&=_f~~Db?!W(J)u#Vd<>iEdS)KRm@*TjFB528kf1i#$bX(tz-Lw_
z@#IN~flUxonfF0(qRfw?H?@L-mlX5+K+I17VNMIhp_r{O!CW@Vh}qmk9JhmNz@=DT
zA?!n=Ka&J;6_Eg`bRq<61}T1+D(j`jwwHuIvJ*O9?U<)n?k2(;#*Z?j+HDc-O41hM
zZzP1ZCz9}+LR%>RhwvAS4j!5m$M*8L$px9Wfb<DeErx@Gq%E7!`g2&z*1_v1_Z>Nq
zRCk)m@?X=XU{%ETV3d;dy)#Ys?I6Qtpj!;P`J0pb7nPEMt;U7A@}!9J#)$IPM5#Pa
zDAR?=V%W~}rlbZ+DJvU<_kW~GO<D$OtW!Z@OYa+!FZV7DwCmGyH!MO?ZM-xsif3cI
zRH#*4?-9O2xLA-TEGy!9NiMvdw~&uwieg2Ath5<ci^e~nl3Ya-m9AFnt`wJtBo~<L
zp<JjN8!2rfUfk65co+CCpYR^gErz9h{Pn4UDHKhuMy}5Jq~=u-11gwu!nOK@WyYj|
z+DH_q$Vj!A(q7%jr6|Tlj8vyw9)&1|fT1dE#qvk8GR`<rx3C}*#)A<H<C1g>R#7|S
zawS?T)=^ehNyOen%1x264cBmL`z1hLV!@EpD!N(^smA&^c`wSAIlmq`Cv;R3lyJHV
z`*PL;b^9I<-%;(GkxOy~#!B5uX~ao=#7RMT+R{gX?ql|4M$^JU$9zg*B#Bf8dGgfs
z(U;0w2@5G9TtuT371vk&VA+Hr9W0MzO8-;!n`#@sbLy}mQ@}O!^HzXoP_+ptchv#`
zLvEi{AZhZ+q<%wEd2uBmsf+l(rsg?5Bg*!~PNIDFMB)u<AO>p2mgo;*!f2k0-jN3x
zW2>O3HgHXeG0RAP3FxMzn06XhV{DF(p@*;<q;X=Obr>7WS{zK627%<VSk2`qOEVps
z7Q6jr@D(%V(Bg_cKx2rT_(Rj81F;adT5a@DRm{nR*aX7NnS|H{qMAP<#32xaqCP}e
zw4o9%99am{>ax?^9aj<EEygMxyEf+ShK$j?551L+UB`V0!XD)~0)he0awE-pxPcc;
zPY!fL+wOFnB&$Er*zJ54V^<T#?0(ZYc*xMNTJK>@jdccK)NMcv?EW;8KzVt@!}|f-
z2LhevNg(cnHcGSFQ5#dQ+gMBn&?Nuyn1B?hfybHTajubqq_0xcc2R=RSSKB7Zeeq7
zGN~9k4&Lqg9z-_JN<6ir#Iu>u`PNAQhK}t4&n&`kb<_!bb->d_4i-c&6L`<vo}V!2
z&lWm1QCd1ay*p;N?skFh)>A<mD;$(&JQHG8^V@H%wQEGYK_g-j|Ki3tmyoKlNU4P}
zW4SiNXD=hE>nuOQ)CoRtMnC&gL>Md7Hu0Nhl(J5~b4FI?UJ7!QAdl`xb5)c|2{7wp
z&~fkJ-_6KlmE1eiXMcmFb4`?F8!wxgHE1mn^7X8jc-|%_1p+m6><M`OK<q(Q+DqbU
z=^4Ve&g^SnL}Z&V^9VmUa}ev#f0>!;R!BS6dKS`#juNBi7oIySKYutDJ4$H-QQxUu
z7}8F?bOo0Im)UVS8p&Ut70YhmM`pzbs7h6rxG<ZjN>#TQt5H#OH`MQrMylB;=6wrB
zH1j*?%~Z1{s#z_{@fisE{nJDjE=^v>O+B!!8KrHykZST`ZhD*OW@i`FhMVYGs@Vl%
zdYl*zF}9FVx<*?S6V<=wbQAzHQT^kq;(a%j`812cH3?E^HZyq4X|@>Rb&7o~|MI3j
zfxm(7D1(Gai6(0-#z<e2ErztP=KWx@SlC(}rferl94I+qN*~gUhxLtNN*YlTKxx-0
z&aH@@%@<F&YX`sr+V}W^AsXxaJ!;JLs1S*kXMd&Z1A!*r3lVG$H@)_=Cg+_3t#yA(
z=t>=}@l4lk;%7e5VsJ@Z5WXbhlUE6ek2TeUSx4__XeXHNA6pYic8SAaY!-=7;kg@0
z<r;f8aO@nnI2_7ta|P;muO;rdXwt-;)lId<{5o1i-p7lpq1rLqGn0PfLepof#H>`<
zjpAAnz`26;GLNH0Hq(MTyUEZn5$HejNlQ^SNPx}sXe7>IsT3!97GVy*GRzTrn6-+t
zj_l<KdvkPqjS+i+QxStX=~o!cz0zP_#9*J(WH3n>4DX&x24xbHd{@{?zS34o#FoM*
z-@?1j-pp0Lv3|PIsn~24G2B)^#1?&(?P^;wFg3ob-<-G2>?r?naTD|LTbE>yvb1Yu
zlgw*1*IdneR<n0&CUs_vX4|DXcWLfMbD@@J{+s4rYfjgS+5t6-IYVRS`Kn7{>=1ux
ziC-;F<ZS7@q>)u-kCrr-OlU9xO=kXBH9_oS-n?{dZTLmSU}u`+tfn<+j%iv<qvkb_
z)S}mEnbWj><~drbd4QIr_VsAd>d+X??bo8`XkJB)Qt#J%1Pm+GGN#30$i{0LfjJs8
z*xI!beB`n$%fNQcA!6jSmzA<$*YB5A8rUOUXB(<jRa5!+CT<Pq*^tP;Z5qP=-Q?l<
zD`Hf03?H*1K9KaBc^c+8U9&c7$<wv?vzqg)7W0@EJJ(z&{#9aP9`isg;XcjXpbbU@
zhC*!^foDX7hHfo`um(x&*7^}PUP@m3|49xe(L&PjYq<cJq?<td4C!8T8g4~~qPsQQ
z7R}m#sVDQpD+*Xn*Tof8hN#7pmo)M6yVBXjuDZMWMa9m>_t?t$a~D@NDdA!1a=u}8
zg^?u3P3G^e9{$3b8ivm~eEQk}iQ)g&px3kj2K<a19Mbjp+HERRc<Fj?LilEfUp%D*
zolkf>8_1V|8@{YY>)b+}Tg7*PTRTXS!@4)>+@<^oxT|$8esvV>9-Z6H*#>XI_FmjB
zox7V4+Tc|W4&viC%))c==?zV+mZ#s7%?9&|dwdDo3#Bngi#8v}qeqEWM)NiIyoWD!
zlQtHz2EJzFbu5D)*;q*RthNkD2-j=)#hSicxRiLMfLF8>=JtLS`v6=wxOV-B#uH1l
z|KaWyFFVC|w0MipL53VVG9>SiE@OZ&(ZiJ)J^^BNdL#d##p~XwPeWHyvJkR2u}GKc
zeH_iv_?We^iAtr;@BLafjkMYA)F-UhWqQBoaK3+&w}M_GD7xOCZ2qj%dw<q=U8nbc
z*E}qIoz}U%UoL(nuIds|t4bzc)9Sr>phh`OP?qXB4KE<w?<k^8{bVj`^^{cV^shaV
z@P$sFrqg@B_ty(=sme;7UzjEp!tZ%4pyT7M8pF2o(#_t4-vOt1R_O7_508rS5?{92
z8+c1_;}vtS3V2_)-uqSbE1ll^li_)t-r7e-vq+x>T{CIcy}#?lOM0I2s8>pk(-m@b
zg%Z7RX6p1JUf1UJOwzfR9^ia>TPnK~%tDsSJKDT~dm%#e4PQIL^L<Rv<CXCBCrp1*
zr-!dHSZdLp1KkBX1$t(5vhuPnAHE)i?a*%-nCW!o(%nX?9MKiRH?gn-dM~2X^nS*F
ztJ6>Plll{N&!;H^gK$j~TD@P+qjh@kk4!%9-r`LdtaHP!rAzgy97kMMLp?9qV&rjK
zJP83^wkTUV9;N%=LQ*Mo0d=eK%D`S}c@}h0c7cPse(%?Py1mhSZCTRb5Iu$;bb9ZP
zOe!)<^4IDWGFgwquG4>wlY|L6Jq>i6QG7(ZJK<_wCj5PGwk}i2>)O2uH|gAJ-NMs4
zcX3xo`w4^ol<B(T{!`3GPn^khP4p9C*PI8gQ@gsie#5%9J$$o4-I2o8Z3C>~6Allf
ztG4|HXZV_Jg_iImsh02EHo_VlA|2p?^vkvlj1A>$wvWJ%QG2(KLH{%QDZF^c7|RD`
zQtkv_vtxwiuQyBV8NMIb`~?!bmb1r4SpGzx(5Zpv?bsc~XFQ&4*|SE{j`B4SY44C&
QfbR$9d{|=RIP3KNA6bsZlK=n!

delta 9985
zcmZ`;3wTpiwm$nLCpk%zCTVHYwDbinr7g4+3k3vi!AN+h1uYMG6{V#-EesWKfWTB9
znNV%rDyUTvbkwUjXhyFVA4NwVIx~WD{g@e*&RoPOC{M?#NbkS)J}KDmyXX75*IuuG
zt-bc+?2~RksP5RWHuqPnA8jf4!v0H!BA#C|cmgZ3)!8Lm$H%c^Hj}Sr_j^_ima=Zj
zed^v```-Q28pjQ;slN<daCB(pj+W7T<9A1*YJn^NAH5YN(B%A$AfuM0DN0w;Yoiq9
zqAT+T(u_KH9TDvter-a3eV<^4v_%#rlxu7jf6<%5Ziu|^^)Xh!kEi50I|8Bi0%8Ai
zfzW~IX72Shj9xJ&5UzS1vV(!prz9SGxvMLzEhtlzrphkJUk`-#BhCZKiV}_c3A&_l
z(2;|*U&s5V)~HowoTVi5!zuAB!vC0>VWmi6?LK}iwIThqfJ_(~gM^MCCnN##Pkdoo
zdJ?!Y;ue12)pfvsiE8w99dAy{W2YWS@tOkRL40poey$!~dnsIXDeUj84Y-!>uLyMt
zphz!8d7Z&0S^|Gqj!bw2qFPt~gJjnpja-B{JEcFb?rq`gQXJ0RrHV4#-|5<F=@J^R
zSqW>K_@dqe-7|&il51-SSvNNy;@_q8k9!pZ4}|;&!&)l;XYVD!jR+QKT)RVQDCtqY
z9lucdwq6q(t5H5%BC;#`PgmEe;V9jXKx4Ze2)nN;d!6vRt=n%$AgstT{Vf5m7cj7J
zQVG=LW0j?^<sYZ#6xLo4xpg9+PNg9et0f1)+R`b-rXgG4jp~L7<TIX~af9zGblV0>
zZK2c-B5jbb%NTI`E2JGRwI@Bo_Y_EU#Cb$%TzgDL*h*F-q*c4L`hcu@kyY3q4Og8A
zxx@Ywf$eY8ywJMJ=jP0L*PF>PARYTl$ARRy=A3{(E8)L@sBN6IQl(WASxw>3X7=~(
zChZid-RBl5?gxn`@n~jW-v-jomD-I``+L&<o#$j_`=*k1vD8kG+Dg(s$!BNvDg5kg
z&r+t6{XNqDvP;D30*OxMFJzUc9)WAyW~qHfYX6J0Kc3+W`{eNa>}==9--(KJQVI6+
zso8}K$CJ||(kWg#C6dz}XAH}5xlY#0ukx%#zxmL6Re{i#fyUEWfzaA0CL=1VWh206
zB;Z@|BFA?@q8@%a`?~8t0CrnMYHOu-5oy~9zHx1b#CV#l{v@r+rPVaDdXO*Z^LX6b
zPHZqJX$g7yrlfw*&A;(2mb3#4>;Od33E%P~eSCao-;*qrkICs6dE7yh?>~kU>a^@0
z4r@mX6$Rs)hW$l*3hU(<jkoC(@pL@lC^`^b$LHiu59V|W50fy7SavM1ufs&_CZYB3
z!5%yYKkPMegx&uQtin}<9231smhzz6V^}SaR)=8~6Fn@f{Cg&bP6V22I%ZVxvVI5I
z?g#R2V|#d2-jKAN`BYv%ayimQ3;XR<Dj72QG0LCG>l2)gBOk1+Ro&N!3F9`bdD@EQ
zQW;$z?dsat-Vh9glZHq<#5eGSDe$(!CIF*8fv`4+j6(awj+Wc0X_OUM(QOEDYM98U
zDUbkKaJWqu3H?pv^p%LeC(u-NIM6g+mwDjC41E;+oY(iyV@o32`=_w@XW-PGb&$WE
z-@kCg)uM3#YV==XjcYIIuKPBN)4V26Q5xG>NSnz$1%sz2OZ#XE#SHn+o>~`9N`gsP
z+tv?9?Yj+2QRcIh)tr+!;<fJ~r6;W2kDdV9gbs|vP8koQ(7PD>A>LY$XM5Nx5{$mg
z|5A{ivYdQ-?uM{7kF&y|xy$JwF+%yi6n%kMi&p>@^cB}@8sjqyefBfC1nTcRHI4@h
zbMh&-Qmp2f^YbUJIRA%^SUCSToLl&jLSHZ#bG|*6^4qav*cd_%1w<tX?`;|?1!$rJ
z#>g&m`bxk5RlypJV1Ga`H^SW8-f-f|Xa`~M9}5J-+T|RXepuUrUTFV2Dq=+IXx)VT
z(F6XA6!Yt^D6o-HdXEenCEwAw_96?GQ`?Y@Jrad|hp|V#kY`eqa=(YRdZIMw-EHv}
zEE@e6*{Q_`{AVfTsiC}f;NbC%{<AEs?dyxljk4kXvuob$0X_99=nJwGVxlI{D`aE4
zYoiaZW|ZLHyrQTeI0rzk*UavF+gG^4nx9~4@r|U8T$$E&pX}G}2hYyC+|E37K((Nm
z=<I#xT^=nASy<E1*1FCdI`1ktNEQvka=PeGy2`o+a=DM+JJS_JQAIH9rL)**Sjk5M
zO-Uzy?&|8E+wcE{#G6pQCO_&$>YwmV?^b^b^*dCc=x6`#>Z&brW!_7iQrlcHbz;+~
z;!H)Um@+Q(X+_nz(DxO$RfH~81)2tCL451PeEk<3`+9WDzxbl!hj_hX&jW*|Oq##^
zu9Yk5CeK|_Q*`HV?k22xb63pGTCg;$er48(Eb#XF8M*TY<j$LwH7YB2p1po$#q#BI
zS3@QQi)-yl{aq_-mfrw3f{4r?bch8BNKwvS&(ciS-it%77~FH8cXb^`f7S7>u1?S`
zpgTa9baZvOvG68haDzY-F|<n1^`I+2zXjb6`uYEMb-fO{3G4hA=qI3^pt;zgDJDf(
z3_1*S<w=x!8aj12N49`|3;G=BPOR(0pne?MXF=EFq|U-@od7Kd&BY3;20aP79&{2G
z>kiN+(00%RpdFwly!8|fiGyZ=nz0JYL90NkL7zDj1OPfqu~T<|c7V2nUW0R`19TUb
zr`e(?pJMY20PTZqJ`wbHpo>8dgKh<N;==YKXxk56U4H=$bs}-l%b;!?XG?!X#X!%1
zR)VHuDb<6f*e>B!kIt$~{EKS}SO>S27PG^AcxfJ6#pjlKSruPhn#p?cM@lo*y~p_%
zrKz?oOi?@j%>47xzU(H>hGepIo{jboyaMefd@fpzZ$LYMKQ$!N;>CoVQ~2zm$wBqF
z!-~^eJowB|j$&SqLA{YA74hy(>KLa*T`!DsAi4<sD9qohSL^+aCiR71<MhV-egL~&
zu<yn(J;BIo07N69zv*O`+I20TH8iDHzEBXRLRl<^b{Y>2&19|od9-iwLuhC7uhBlw
zZDpD4eV&i@I3HJ*ss0@03(GQ?pNAlNlRuC4TYd;_I{zB&er_9<$x?Ve+6jEzu)b^?
zUotF@?c@(Z@;Ck}+IapyXh-tz(N5<{!!y}QK4f@b*21R`&%EO#roWzY>n_M5F<MG>
z4VJ}FEQ<}l(y!^(w?e-S`l4SILXMaG--P~|3;Z9$Gt~7L`A@@d)V5<0vB+G<A`?G2
zYMOenDAMP;Z&b#3QpIcPCB`e%<Wwr5s(gxjfSkMzAB9(x5kl>g>T7W6lIm1^I@0q6
zKer)~SE^<M;5hDW(xQ{T81O{2tME7BQ$ng==quGztfV&ZLqL^khD#kijXe8z0zdM}
z=n<BL989A{bJRfNtibp!7M?k-m<^5GGVX2$vVVMl_1y-AS*e4m&Yps$V>5R6%@CgO
z5~QMF$R!gevkm;-iCIB{R5iy`$TfmgEuwE(1SghZIg=t>lyj0Oto3K~tV4yi4<%x1
z*&fvE6_)8Bjj84JpzTJ|W;v{sNXv6Ui(=Lp_4V1Xu&+jw1~zUqIwQ<2dtl(8Q;In(
zu@vK=r-!vZh@KPhjx2&#Tc%kj2-AWGdBaWqnP#2fl<5So_@Obz+=4ER`N}+!8IE3@
zS}Id7W7H<IwqWTB#cYfkAF7I~o9BU;zl66M$K5v%<r8mC32p$ZJ%NTLS!(Exibw2M
zP}ce-TFYlb+lTsUOThXvvB?gqnyIDIqUic5baj0!S{pTH<YvEJsHa(!xHri@&JxRZ
zBNX(oR?EjGB^|0<i9F3BkmH6EoZ3?9IPK%ETXJ$;7?v^e(}1FpcP~*4m7WzEe`)qe
zUOp+#teVYY@)!JnlYf!Y$+5ge`w7prfTM@C+R>AfpOabxiJbhEFusOVrk3WQ4U?7Y
zNERH{Lj6i;`|$yl!-7sE$()v0V&WtdsI+{X8RSCqs|vIJZoyVMhhkIE(+Hj(+X98T
zoU~rmmc`D|hIU>|yHH_nLt97YyEGJA`8zPPi~?;zp9*wkhmxJ5%;Ba<cLY~Kt1|C4
z2+aE9XsxRu!L?rh7>wQmhPiAMUeP;=amP_1biIZck24I$t=L{9?x<y2Uoym{LlRWb
zg~)a{l43_Q$QC1xdzOqpbr9=RTj~_sViNq>a+05$l$<BDfA+<UN0RP+VJvLFCjJ#m
zIxm`>%%0``$piY*#p~~gV>9W4W#BmB^=iy3F32C8oHH_mT%Xnn`LEg1uu>8vIvDH5
zE@D5OWyE$8@My?wrnCIw<blJA31GWrjp010+j&d3^Y&EfTqm6A_F^-2@aeZ^1dFK%
zn?&?4v!o|gogUVzU<JzPn|fazT>=~~k#?RHOR#jdTp1SSvn543R4eYch+ic_{5DIv
zER;lH-qt%QMjmCc$|u*;O#2UJ?!2vcB@I-DYBxe*Qps|WUIT#Z4X}|huvWTN$$S@O
zUOnK~A-9?K@YT0v1S=?;YK=l^2Ed4J09y4DpxPL)te6DY(VfLV38YF4=@4$S^jyhe
zY&S?%`qf#8#Sk<#$53pK5XxA~Yy)AD1VY`nbt8;TGZ5@z?TpQrY_*t2xw_VO<Bg=+
z46{6*A%XUF<0U~Mhc)b4Bc*BfW{N&6PuBc4)STE+E8&E*Q{c-{*Bkf=efZ(3%$Z&a
zGQcv^2vX4<q^3K_fY`8Qi~&ZBc~^l}$C6I`fa0c!Qg-qQQ?jqWvb@V<MoJehk7kAL
zsK)WJ3SBxruDroz9vgLjdGs}<Y6bMBeqQ>$&7^8KVZp0EkSM>`T8K3H(ov6+s35r<
zB6S8&m|EcckR&@&TS@Zau5nY<AVX+NDE=rQEc$%3&J@^K!a}0{h`u#7aT(dqqJCOp
zGxf)yucsu^t5*G#`l-F{Mc<@r$$+qCLgcpDb@8rjLE&r;$6n|a{Z?Y!Ne8G8b1OeM
zJw9lHx!oQ|?@L9u6B8#GrhABS2u9WaOpH@7CPm#$T)b%_*1j_lgRr|DG<Iiyl6yp7
zrE@lpdD~!P(c5YF&^hZ|k0cys=Me}@cuiYq)Hy1jaeMFJbMSV=J5Lenaq7EWFQM<Q
zL7zQn9g7Hq3((;CND{1P`eOB*6Bb!yfd!V34^Z}S+$VyqORyqCwk?jzthQ7qHW(O7
z4GgK@nHZD?(-3h2;mx;Dk{k1+*={Nj9@bic5qJdV{NChZYB{(r;d2OCot<=QOL@XQ
z#4fN;0x`Ah2qsJ<exY-&;MWHeo}vJY;+F}&V_m|((HY1VJAY4QY58njV!IJ`kr6hb
zj67C3smypQBsTM%GpZdI2zIJQ!T^t-nS8CVnr4$$ixY1(PzF-wczWU0Iz5J=D9M8e
z^cLHB1o?^2nc2_rFbT#8z(;w@%nJ4*|72!v-)AY#C?USNpM;ty3F>hIKA!t|+N=Wn
z)8@umnT}V;I^Rm=PvgsG<(AZwpwO6)@`N`i$N<68Z8BjysY~orFiqk%nUKv7%*t`h
zAaR_4d4c~htAwTT{Mi|v??{|vzW{4fOSvT>%BRgPEYPtG%*y=`8oM}*F74v(i#94t
z&j`MI_GIRNU`|SqR<7z6w_z)-T-77`ba{fYaaFX?IyN)!e)wwoU(o8bj;*wgHM8?`
z2#o8dl`dPFylGqM&Z)_pwsjZcXjWGQrdGOmYH^}pV7-SF4#C(6s1Qze#Y+1?i+5&$
zv(i4$-2AU|MhA79$vp{HcpPOkG6+@~0^$B5gwD}WTa|c%VhNU35SuB*u!?&SO4duW
znbyZ7J4liQNmfjvlcyKhTVj%4BuSz8&4$QDpWbY)c+)*(EUq&^qKCE8-N4--EaK%W
zV50j4!S=on2G}2-J;dJS`mJEAJ+BkH%wVS_JZ(VYKVafX#In6J7`qbUm#reQA8oCH
zc0PT%;VD9r^sv@Ml-Mm!gweTVffK=VCpnC9v<bG(BMyig;^qss(j)HnH^!4AZfg`R
zATEGv%e(}Ok;_fS4?RlaYK2WPt`pUm9NkA#s?+Sc9-uI#BCnH$5a8J$3G%wDJFn6V
znb$O#*N|W5HPpzfT5;V%czFWv?w|oEKC1!BTHomg%**~2VE<nO=63`3dWQg01z>Dr
zG;}8@i<X-CE3C8_mP=NxcMObQu1V(a-5HEj-pN&!CG!8Rxt%YY?;XezapG%=DSkxU
zwQ<Aa#wDuw#I%<75a=!?4!X#O+7ETMkKc6HO8nsP{#|(k?9JNXN&0xrJzq<{pgGz!
zt2!%Dvo+|&`fiQsb*fuo>>!dhVK3J&>cz*@7pWaS{+s$bQ(@&ne!M=}`jt=*5o-6c
z{MoTD5|e{z+6B$ISJPH#&goiWgXYyoYVr4KQ>SZA^hS@={CchCJg@at^Afdw=x69f
zTAG@Zpv9}DiJE7u7GJA*71gXhq-7E@t>hb)^)fq}wPE~`Ww|DsevBVjR>3+Ww&fE|
z?1jjbhALJ4kn!Csd0edMCLiNVS7r0Xt2}Hk53cGP+^e}>)#BT<xT%_bg~m*_@mk?^
znq#l#`=e&lZ`bS%TJPJnlna{cf|fX6OPa42AJY=_L7J~l^Q_Q{kcz2TD<SftmQJKi
z%ObKu3fr_C;>K&mM8;|PAZ)zmA<{fUPb}74rusIK!TmCWWWHqe0G1tjbakbP|7VRa
zlC(D4%qMOL@W(a`i5%T<8^bRxyr8itHTK_mV<ib<Qt-c-(#Z10CsbC>J)69~*uBhc
zn3?EU!L7TVS3n!PYU00W;wd(?&3qZOJ4&Q9<~zgC?%_{C+iqx`hVN!Wdzc@AwxfqO
zVrb8Dw%M!dgL&%a*=!zPy?G^mF+8_9k7aV-mQ3H-B8diT@hrfpmY~EdH}gqbK4KX>
zJ5<b;@JXRt*kJxtsF;@E$Dtx$Y^lXQe)rHORZ^53?)!bQr{{;%C;VR_Z}22+Jc>_=
zoyvH>i03=VDP(#d$saQ0mkfDqHO01YJRHfs=bhzFN?mm+7mXnqFbaJhVwU*SkjGbI
z_}<;)n3=;%!d{=>(DwYunv4Q?*<xyg>kONoUoBe=d1o&fcZm`HF+<+-2bY%&dCwmw
z-ZbQ2W=Vi;2F@WNm!tT~u<7{;@q@_ICdSfRlQYs;X#f}`pKsXe^^Jg*YWQeBDXuoO
ze*QAFlMU@jV@8%6+SvkmxuNa(_4qsPY4ZB$F7BbT(0GX&uF4a<qRAWFLt3m9WssC(
zAmVw)z|^uNAN!#E$dLE^{q9Rc9_~@?^M<_V4{#bT7#MGX^0-$97%PdABo#`tv1n!)
z@?1XmKCkaOL)-K9TSZz{hqjpY;Vt)h=l&L!w76q;2l@M_qR_AY9&C)*ohBxK(2&P2
zB(XxYLhc4mPLJ|E$M+!$-6?6>dfrCfHqab~R}Aev!@lQhnm*eI=cj&hxHlW4{wIyA
z<Xc1A^RxPAL*Db(AP4s}dwrRB=@Mq_gR<T@Qc6jSAK|H;jy63pUJ2yMfTN7KYe<_;
z_gY#y{d<i4G03y<`*LFNRU<~v5A^p9d0MUnyvB(Ag(2_xv(OnsKEYTGlZ`xnG~}0)
zrQjw*?!>fV7LygdM}ZTd=tHTM3`CEHANx6Op5c0(aP=wGh9_9@OyYA}ysSF1y5$>F
zWZ^@nS>%pwH>>8QlkZv?`QsyBvB>tvs!i&--u$O0itMo?4R4~sPu_?FeD#yXw%Ge<
z6@U53GP_c00GCG2Ke>tFKRT<QDq}hP<)^Mkdk*acUb6Fg+lQm2-9vo!&N5r|9g=&A
zzX2|Ok>qUr$DL(%ulPw?QI1&oh+RW%(e+Z<&R6d$vpv@=xn=x~U1fIPBa$1-c-vzc
F{|CvpNt*xw

diff --git a/Code-C/main.c b/Code-C/main.c
index 32e7fd7..f48a16a 100644
--- a/Code-C/main.c
+++ b/Code-C/main.c
@@ -9,7 +9,8 @@
 #include "getArray.h"
 
 bool rawDataWriteFlag;
-int nRow = 500;
+int nRowRawData = 500;
+int nRowGR = 150;
 int nCol = 1;
 double freqEch = 250;
 
@@ -43,7 +44,7 @@ void *threadCalculAverage(void *vargp){
         while(queueGetNextE(rawDataQueue) != NULL){
             rawDataQueue = queueGetNextE(rawDataQueue);
             fileName = queueGetTabChar(rawDataQueue);
-            averageThreadFunction(fileName);
+            averageFunction(fileName, NULL);
             remove(fileName);
         }
     }
@@ -57,7 +58,7 @@ void *threadCalculBoth(void *vargp){
             rawDataQueue = queueGetNextE(rawDataQueue);
             fileName = queueGetTabChar(rawDataQueue);
             powerFunction(fileName, NULL);
-            averageThreadFunction(fileName);
+            averageFunction(fileName , NULL);
             remove(fileName);
         }
     }
@@ -67,7 +68,8 @@ void *threadCalculGrowthRate(void * vargp){
     Pqueue rawDataQueue = firstRawDataQueue;
     char* fileName;
     //double pw[nCol-1];
-    double **dataLign = getDoubleArray(2,nCol);
+    double **dataLignPw = getDoubleArray(2,nCol);
+    double **dataLignAv = getDoubleArray(2,nCol);
 
     int i = 0;
     while(rawDataWriteFlag){
@@ -76,18 +78,24 @@ void *threadCalculGrowthRate(void * vargp){
             fileName = queueGetTabChar(rawDataQueue);
             if(i < 2){
                 if(i == 1){
-                    powerFunction(fileName, dataLign);
-                    growthRateFunction(dataLign);
+                    powerFunction(fileName, dataLignPw);
+                    averageFunction(fileName , dataLignAv);
+                    growthRateFunction(dataLignPw , "growthRatePw.csv");
+                    growthRateFunction(dataLignPw,"growthRateAv.csv");
                 }else{
-                    powerFunction(fileName, dataLign);
+                    averageFunction(fileName , dataLignAv);
+                    powerFunction(fileName, dataLignPw);
                 }
                 i++;
             }else{
                 for(int y = 0; y < (nCol-1); y++){
-                    dataLign[0][y] = dataLign[1][y];
+                    dataLignPw[0][y] = dataLignPw[1][y];
+                    dataLignAv[0][y] = dataLignAv[1][y];
                 }
-                powerFunction(fileName, dataLign);
-                growthRateFunction(dataLign);
+                powerFunction(fileName, dataLignPw);
+                averageFunction(fileName , dataLignAv);
+                growthRateFunction(dataLignPw,"growthRatePw.csv");
+                growthRateFunction(dataLignPw,"growthRateAv.csv");
             }
             //remove(fileName);
         }
@@ -105,7 +113,7 @@ int main(int argc , char** argv){
     rawDataWriteFlag = true;
 
     period = 1 / freqEch;
-    invTimeBandWidth = 1 /(nRow * period);
+    invTimeBandWidth = 1 /(nRowRawData * period);
     firstRawDataQueue = queueCreateEmpty(); // change this for create empty
 
     pthread_t rawData;
diff --git a/Code-C/power.c b/Code-C/power.c
index 6fd69d1..ea62528 100644
--- a/Code-C/power.c
+++ b/Code-C/power.c
@@ -11,17 +11,15 @@
  * @param powerArray array where results are stocked
  */
 void powerCalculation(long **p, double powerArray[]){
-    printf("data in pw : \n");
     for(int i = 0; i < nCol; i++){
         int j = 0;
         powerArray[i] = 0;
 
-        while(j < nRow){
-            printf("\np[%d,%d] = %d\npwArray[%d] = %lf\n" , j,i,p[j][i+1],i,powerArray[i]);
+        while(j < nRowRawData){
             powerArray[i] += pow(p[j][i+1],2);
             j++;
         }
-        powerArray[i] /= nRow;
+        powerArray[i] /= nRowRawData;
     }
 
 }
@@ -33,7 +31,7 @@ void powerCalculation(long **p, double powerArray[]){
  * @param N number of rows in the file
  * @param M number of columns in the file
  */
-void powerFunction(char* rawDataFileName, double *pw[]){
+void powerFunction(char* rawDataFileName, double **pw){
     long **p = getRawDataArray(rawDataFileName);
     double pww[nCol-1];
     if(p !=NULL){
@@ -42,8 +40,8 @@ void powerFunction(char* rawDataFileName, double *pw[]){
             appendDataInFile("powerData.csv",pww,nCol-1);
         }else{
             powerCalculation(p,pw[1]);
-            appendDataInFile("powerData.csv",*pw,nCol-1);
+            appendDataInFile("powerData.csv",pw[1],nCol-1);
         }
-        freeArray(p,nRow);
+        freeArray(p,nRowRawData);
     }
 }
\ No newline at end of file
diff --git a/Code-C/simulateFlux.c b/Code-C/simulateFlux.c
index 9c0fbd9..c0183af 100644
--- a/Code-C/simulateFlux.c
+++ b/Code-C/simulateFlux.c
@@ -151,7 +151,7 @@ void *threadSimulateFlux(void *vargp){
     FILE *rawDataFile = fopen(fileName,"w+");
 
     while(writeOneRawData(rawDataFile)){
-        if(cptData == nRow){
+        if(cptData == nRowRawData){
             fclose(rawDataFile);
             cptData = 0;
             cptFile++;