From e23d744cdc5b66a309b20167742b132efb00e4be Mon Sep 17 00:00:00 2001
From: "quentin.perret" <quentinprrt33@gmail.com>
Date: Thu, 16 Jun 2022 14:35:02 +0200
Subject: [PATCH] beginning change power

---
 Code-C/getArray.c     |  15 +++++++++++++--
 Code-C/getArray.h     |   3 ++-
 Code-C/growthRate.c   |  15 ++++++++-------
 Code-C/growthRate.h   |   2 +-
 Code-C/main           | Bin 0 -> 41872 bytes
 Code-C/main.c         |  23 +++++++++++------------
 Code-C/power.c        |  15 ++++++++++-----
 Code-C/power.h        |   2 +-
 Code-C/simulateFlux.c |   4 ++--
 9 files changed, 48 insertions(+), 31 deletions(-)
 create mode 100755 Code-C/main

diff --git a/Code-C/getArray.c b/Code-C/getArray.c
index eeaac8e..7284378 100644
--- a/Code-C/getArray.c
+++ b/Code-C/getArray.c
@@ -3,7 +3,7 @@
 #include "initialParameters.h"
 //#include <string.h>
 
-long **get(int N, int M) /* Allocate the array */
+long **getlongArray(int N, int M) /* Allocate the array */
 {
     /* Check if allocation succeeded. (check for NULL pointer) */
     int i;
@@ -14,6 +14,17 @@ long **get(int N, int M) /* Allocate the array */
     return array;
 }
 
+double **getDoubleArray(int N, int M) /* Allocate the array */
+{
+    /* Check if allocation succeeded. (check for NULL pointer) */
+    int i;
+    double **array;
+    array = (double **) malloc(N*sizeof(double *));
+    for(i = 0 ; i < N ; i++)
+        array[i] = (double *) malloc( M*sizeof(double));
+    return array;
+}
+
 void fillArrayWithRawData(char *rawDataFileName,long** p, int N, int M) {
     int i, j;
     char *buffer;
@@ -86,7 +97,7 @@ void freeArray(long **p, int N) {
  */
 long **getRawDataArray(char* rawDataFileName){
     long **p;
-    p = get(nRow, nCol);
+    p = getlongArray(nRow, nCol);
     fillArrayWithRawData(rawDataFileName,p ,nRow, nCol);
     //if(checkArrayFullyFill(p,nRow)){
         //clearRawData(N);
diff --git a/Code-C/getArray.h b/Code-C/getArray.h
index 02810a9..a53bd0c 100644
--- a/Code-C/getArray.h
+++ b/Code-C/getArray.h
@@ -8,4 +8,5 @@
 long **getRawDataArray(char *rawDataFileName);
 void printArrayData(long** p, int N, int M);
 void freeArray(long **p, int N);
-bool checkArrayFullyFill(long **p, int N );
\ No newline at end of file
+bool checkArrayFullyFill(long **p, int N );
+double **getDoubleArray(int N, int M);
\ No newline at end of file
diff --git a/Code-C/growthRate.c b/Code-C/growthRate.c
index c381864..8fcca10 100644
--- a/Code-C/growthRate.c
+++ b/Code-C/growthRate.c
@@ -9,16 +9,17 @@
  * @param gRateArray array that contaigns results of the growth ratecalculation
  */
 void growthRateCalculation(double **dataArray , double *gRateArray){
-    double gRate = 0;
     for(int i = 0 ; i < nCol-1 ; i++){
-        gRate = (dataArray[1][i] - dataArray[2][i]) / period;
+        gRateArray[i] = (dataArray[0][i] - dataArray[1][i]) / period;
     }
 }
-
-void growthRateFunction(double *dataLign){
-    double **dataArray;
-    //initialization + fill
+/**
+ * @brief calcul then print differentiate in csv
+ * 
+ * @param dataLign 
+ */
+void growthRateFunction(double **dataLign){
     double gRateArray[nCol-1];
-    growthRateCalculation(dataArray , gRateArray);
+    growthRateCalculation(dataLign , gRateArray);
     appendDataInFile("growthRate.csv",gRateArray , nCol-1);
 }
\ No newline at end of file
diff --git a/Code-C/growthRate.h b/Code-C/growthRate.h
index 1e0b201..7979f21 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);
\ No newline at end of file
diff --git a/Code-C/main b/Code-C/main
new file mode 100755
index 0000000000000000000000000000000000000000..640008a6248823e31c5af521edcfbba984debc5d
GIT binary patch
literal 41872
zcmeHwdwf*Ywf{LY*^`rGUP%ZLp2JHXLI_VmK@%V_LBOCOpy)6plZWJGGQ&er(13Mr
zj3)hA#g<y`ZENZ8+Lo)WrHD1C71Z7qAGO@7x2Z*&2yG);O<O9z@7j-(Gbb~xxA*hw
zr+;|z$?Uxzd+oK?-h1tR)_LtJS+Ufi>zX>8+T|LDas;MG5>9?macEkRHc!jK?>KF&
z<^nnm|3#8)Ppch*pi)EfIgwJigOYAfM@T-EDh(;6a;c7Lx1gw0^9z{DG0IBVp0<q^
zd@5CRv%W$>*KFzfOo3OawH(>z#6m$ou`)5;NCv=hQl_dLnS;umlI2dxaw^rNy;Q2~
zNj9PX=Oo>pE&+`Wl@c!<W$#2QC|^C3g9Fk%((o!(<yNB{>G^37W%jFOeW~m+PwKBy
zS-m#Bwr2C(+0$z)rq$Ne$F@(~K5y={xw8wJ8wzIfb`OFgWXff$)@m2;KY8h`MgO|$
zmp^IWx@`M(V-{_#+CzGi4kV9w=pcS&PRftLKhb%+^ncZMgG{)w_^0FL$WQ7j*M|M`
z&$X&yDF~L|OCdl`9Sc^|S^y=(FHeELH3j}BDe%oH@YJ3F{_RI45XtoUZVG&F3jEJf
z;FqSr{~!ha(G+-R3VU{>;D0;?eo_iKUrE9L^%QtV3OOUeAHcu;*Z@Q_J9nnQk4YiF
zE(QNbDe!+ufzM1K|9A@io)q}0zz6VeKXM>n%@-@J6+pI@uNAEpSmn3gmiSFJycQ~5
zvocf>X^K?UG)E&%YgR6<ZK#i|32&~AXrWM5T|<4SIT~(?hC+fhDPv>9HmxyQ-4qE|
zgv#+7jhJ9$dreepjy9Dy?$9b58YA@t)R)&aiWJz!ShT!4+yuCyrXF>Mo0}u3uQFUy
zt5sG-qRm=)xVE;ToC?&Ud?gArYnA1-4b2g)&LowGqgqp>u3>9Li`LXdnj0hKA(SSu
zk)|fCa$8dk)T$&M!B-!yi)i6!Lk+|=)zn8THSogaQKQSZheMS$_2Jr@TOw3xbEFYi
zIL2~min=fqg`|iEWF0xjxs;O@Ah?j=azR<PqIA*X(9D9F^EESDFvrZ#m`gq7z@L-<
zlBXod0e{P6PQc=q{HYbDyedI{RoAY-*dQMqz<<2s0k~PFl><iy(Mv7#*sz)ml2~uj
z9ohy-C;kU``84eoi4*(|lV7d3YMgb+NqADOrwV^SO-}pxlklmls5Y@Cq>C*7?XuvN
z2qJV_@PjS<y%v1Hg7357$5`-cp3t7fI-k@$A^bQC|Iji4kdD%191*bKsXx>)$%3b~
zPaTC8yjtstHqV0R`zTqU$b#qVro;y=cx7uMm09rUVDngS!OJZbr){#}vrKO6S1fqB
z<>S;w3*KYl-)_NsE%@CQyw8Ha+k*F7@a-0Swgtb>g3q<!4_NTDM^VQ?3qDUFksh+(
zhgk53E%>1pe5VDUZ^56i;D=f8CoTBl7W^p--nzNzvfycrQ%AQ2KT08y_FC|AOKqWk
zmi^2EpIP9a$pXQavqrF$y>*4A1^0GF9T!dpTb?tXHC4PY=WUQKOgN2S-{>Ob2(Bi2
zPuB$;6MjK3Ejc|WIs6L2wAA!;a`;7pX^H7M#NlTMrlqCl0EdqdOiM~nJBJ@5n3j^B
z-5mbE1k<wC)5zg}Bbb(yo=qJ7BEhtj^ptV<^90ip(o@9Ye<7Hbj-Eme-$^hn89f0G
z-$F1g6+M0qHxW!rM32Vdn+T?*p{Mr~ASP@k*iUd5hp!`;mV%y>9KM=h+7k41a(Fqx
zv?TN#;_xDZX({MAz~M^?rX`@Kox?K;rmaKIZVpc&n3jT`Mh=f9n3jN^O&lIZFuDGo
zG7jeuOfJ8th{Kr#ldJD3<gnv%O{;7N%#RtrEsei?V=(@9u;p}b*_st?>?XMRV4KIu
zlZ7rU1{d<4Sg2_q`bIAy%EI>{AlSy<qhx+xG!Mh)UeIayWL)U-jowXypOq>2k8%Dv
zM<|*6Nig0UJodXQg2(!tLH+sQ%b!FCK|n+bFfMde@_JPH-MqLKo4z*KvT!N5v|#+5
zs5{uYa1rEnPeYG%SL5e-_7GtG231a_J-33m7G;8WoGrvCjX(7!tPpIQlRgrn+7?nX
zyFW*$cu6nN7tmgy|AeHQ><x7H@sjh@q%$Let!ul3ZF{?b{~_(Ng7M?sXZ~{GLQBbc
zy|v{0WqV2X&i_ir_8i#ALPJ4BJzPdG{yIr`2&6MFwWHQORO>I!bN#jpMal{FF%@ow
ze&D`X<hv<98u|7Yp5@JJIa?Tv-v)H7E57zka6{8MSC9}eK0<}iK-H2RttD@^bnEj=
z&fhi{iU;=v+XnrJ+K6&pLH%T~)lc<R27UF<K+K;pd?4m9&GPPHkacGErzrK3wbW^}
zN|~wTP2WA#EviZDL6gn62L5PL{I)mau~SeHdL8R>QUiWOMv0en9WOacYMv=*EjcBq
zt!vK)%U=n$xhWC>{nvH{+eTzlnGRmY*-~;AopbiGk~eq0$eY)CL!Yo&s>(>8|7m68
z|5jyaSQ$F#1UjfQzV>K*<>BD*l2a$R&S(qt{nejEwNUZ|jNXnpOz8InTJ>Ge9io{S
zL@Q4P+t6L{pAg#z#D)RZe7vNOEOF*i$l;^p1yEYoo=>z=XHbVVGF#X75&qcQn2ayN
zpm9e_NuSQew$1TN4b}?{x`QpZb!)MMW;Y)7-P4S2JfpX(p1k)@W>2EM0UD6;*y;FM
z^xJu^)a@8I@#itrpMb|YOsMW#P)6x|jeS^ki^14h(kI6PjZK(^MBHJDxCA0FaQTSt
zi)B7d-oJbEA8Ew(q3`=HE9u(#byMxrzI(Pp?RGVDgpx|bDalp9$RZCHOXo_T$?li7
z^Cx^>+@PkRYPAwLE|oxZ+PC{Gv*Ps5pQ$c;0C=UIv}0@D4f8dvm0g41`0FV9pJJ5u
zwR)BiA@5TB23uo&O2+4eqmW}=`37i3v$&@VKho-U{m21S`hMSe&X=v7sY*sY=kptO
zeDa>+a)Q1oopyGr<eYE)F{*_361!y{qNL=L`gXP^=k}xdYrqAeFP!sDIgy;}Mky06
znN0yVoSg5TY~JCbkBZk7Ut7GkcnwYed9cgDU>uX^S=bG^<7zQD(aO$V!Ziw(V^8#+
z9*nQ;M7A{kdhydRp2qX96+gxErUN}pJA)wFf%g4?@<%xR<M)GGI{m@P&rw}1o$0}2
z-A>`#N)BUMozS9-;kph(2RI16DlEfN(#mR}LR%Sm%EMyNe;zx7?kpPbd8MFH0l|Me
zckypx8Ht>_zIbEt4aGMWhvcA?>w+4QH=<NGT)eO!d06TTKMBJe@5Fx4PBPB<Ho~la
z&sx;^x6S@ILPuhG*?E6^yHJ8VTUh&)ur~G#f~nuS-Q1$D6hFx=MwUKV`~#l<Us$(`
zuL!hS2s54KX3as}!mQ|d@+DVd1J?ar%(`G4CNx(Fj1Y5GZOwag4jS)y1;1n|;dZmh
zH?}PtN#)D_K%Pv@uU})9eTw<jQv*KX-C8=~cg<DUoKCbji!S2bLVmE#t=4LAod~uX
zmW~9d!3ngkJYS5B7H`{q*lcpI)$N0VqO!kH5lyi-&SBo3!WQ86Q|NiN1yU+~lTOf{
zptJi6A6~eCsat|onl$<px|7J$B(jS{c9F=T2_hG9k?TvSw!FV!!%VfE^WFF&vA#81
z)6V%e(9|yJO+EnO9g^&1;A2X9`<+Z4<OnC@yJs&Ag``FN81H^|GkBBJxgKy^<YJD`
zk}k#zacWIT;$mV465UJya~PWpxS73>Kn~*{;bwe$_F_iC;~dapGi)9w3m&AEeLPc5
z4e4>paW2#S5Jt|K`55*Gq=P_Kd8${j5Vf+W&2k6K6?zgCy^#i8-hPm&jF^QltC^O>
zS+$(?KRpymcE5a%Tc{7(_Gz)vI49W!c}K8e#}0kYuW`zSiFX=F_cG$dNTiDxyqdn?
z@gg2|BrNjeh^czD7@`iiqP}?S!+zn5Ncc43YnyYI5Z+3{-+NyOZ(aLg+~2<_bH9Mi
z4VHyd5Q2>nY3Ly#ScpkTu@JJ9gnXTZwA^+|aQz0TBbN%U;eu-daqSXSb+@iP9fu1(
z9o+XrTC8bRx}Qk#XV8Z3V7$ywu6!a5Pz)S$aBXe}C%2XvV%X}>_@2=3#*R*$os^v3
z`LY_IaBSQH*p2q~n;=Ab1%+NMuo>)oH<Iq(^m6Tch4%A_t5k5^Ah;Td>lks7z3B`m
zUS_eCsLJd{g9KbF1dI>@#*l#Brhp>9c~)a9`z9EVghkmdoRvexcaU^1k!4%8R3#sy
zlC`4bBckNLQ_0N!S`D*oAb#6fI+YS5aT&41(K53MF%ou*5Vk=GyN-lC^?S0x+OF1>
zXUP;|jd)S+#F7rg0N<IJ186D<{+tl33Bim6S8&19`81$QdcmQc!O0{bI&tOc@$P%z
zgE$9Y;EVjUV|_#PNs&RrXyNZ<A0a<O<a;e@b(;pTnf>d~;mAWmtr`SJVD@?>-T$fN
zwF>!f{kC8J%gFc77h5|kNZ{9nz@b9mND{bH2*maT+EZKkfIr?H5$B8j6S$DuO`h$V
z=9@rb9~WYeA*cgApGVR?fyB1l)+I*p1=MupyMpUWg6nbOx*)Xw8%CbUcUoLEsAn}Y
z?+J81&RqEx{W~TnPP6#=JKx&f3WsSojBMsha#K&u=I7vK2=pB7QSzapSjU%<o-`Jt
zqxw6SSSFY(<H2&q4b)g_*b`?{)<~y9!pJ3k6ycx&b-a6VzaAVeljbMdM519}y{eoE
zZSBw-M_Wsd@>BHFu&j1wd`rn`oP2eWd#ZV*s<N^&c<jBvu|GME^^J}ni=Ui)0?fQi
zUME+F2vq#l{x$FvvZe8)o53TPB3Qz(BiBzpO?kP+{y1Ta4X;^plKW87c`pnA#EIi2
z-CVhF>)NBi@@IlIokFqhz_C6@{Pkm>jGo*J<xU*ya_A@Y&iKod0km`$wv=@1@ne|T
z5ZN8vcemSt#{RXI9c{Vos1|!j#Vc?U*nM%ySstvwdH32ALH%Wt7To81BF%BPZU!m3
zcm52)3X)6<tbX!&G?g#42<fP}1>arp)0ZObDN5tdAsTlnBARF*Rspt5WKXuTD=x;l
z;_>eF!Yg2FC;|<A?~yU+gy#g=riVcN6s<1mG$&qi08Zv0nW{AYyg1*nnmzTI4JmOL
zk`I&Q`^<_?61ZDdBx5@7b85*U*x{5Y!Dq*xdY}$jhzsGUiDx*9<Q|gy_}2s%nMt8^
z&Q)?~Mf@~Qg1YD=$UFs-?mWJ>TdkZIsjVGth{(otM%s$A#j^v?@w<1vau^ffLjUO$
z7Jey6g)ZjN%%Kw%C69LSB=7uf`N>J~lJjPK;`Gi==p*IR+T2#sYp(r6d0XM-YkIp8
zNDsEI><b=wWA+n2zw^I`FAl~}1V4T+7(a=KYEJOjAFl|;Ukpw@8Pt1)FYi~Awi0JX
zNnP@WPewrhi8TM7&ZrAGsSKP-`t5J_t8&NLLu8%!Q|q|{$n&5SbVGN(e1zova2O{a
z@2==Z7wzrz?b(9RzwlCpR0M;nbuTa6Hs`w;B<i@JzDflk=?HdA5IN^dg4liLHy19z
zx$FwI<{g;CD?|Uv$~Y?q(#PN-4>UH|%DzNQ@nd{jC@j{ta5H3eJ7Bdl^SH8Y!5m(d
z!IkCVvI-K>hOyIzN!YgB5Ng2B#(&x!1WhRUBxVj4A=r=_+tYww*oHRi=;^#u6O{>l
zd&-Ek_I&)L-l?CY(;-akCG6zgduw3@3&CL$>WENszO76bg-k*{2+!|rz;7~%z8`iT
z>oVfUlE?%}hvy>VK)ELI84!c};`i@uL=C~VFHpi!IL&lm%p$JXF2jY#-Ui~uRxJ*1
zN3BL|CMM#=txNW{x3yn@QPAhXXO8t5){fR-XtSf=dl&s6;>CZx)E97($1L?-l%kDm
z%WdZk-@O!|Z)GPh!iiF+9%sjq;X$`EjQ5kdYbnGsbV`mdaehVXZXbr%%|dwhe#nK3
z>OFINYu;{%X=M)l&VMFW55`{dQZK%Pvv{@%SUOtSpC-s!-$#Z^1Jn1wL@j2s!N(`0
z*y@EoRzv5S#QNWN@y2zJ5cOkO7qvZr#dJ@nZ!hn3_7<wp25b%3PpM{%!!xT?{a?k&
zGGP>)fs-t17K!~K=}j$Z5~mZqReMf@flL8>B^Ac>dCd0&tA*%2ow3scE4@xtTB^zw
z*ebnWz`RnD`xRPbWxEDh>ig&eQQ!Ab-}&*W>U$3q-pK1zeLtjfmPU@Z)Mqxbt+&6@
zQ5FUM&#yV2H#jC0cmxXEIzRzFXo85fNd;Qj2jj4AKL0db#O!{R90k7!^dZ@`=Qm_r
zX~>Z_1zQA+UCoF3t|I!4jlohzo!uqZt^YRn@yEO0fx2gMZ47h<(NglEex_FK^5>#>
ze_XB?u7}5iF#3R|;E#zyxWw-e!H!;pCZ?hpeAf2vJZh_Dokfish1IyFuB0<*OzFM(
z@CvZjN{GPl=;|LHg}6F{GozDg9pz`B<PBdOtGr><*UAn7;&b&Ilo1m-{q|R%XMX~$
zLr>4=k)bflFR-K`%tN6a_7QNctP#H|#!DPsoO1)Ujy-DOJ_)>-?9Wnj;*9EoL+^!S
z$Dn^1N}Wb{&KxS6eg&>Jx`apPQAGX|Hs~g~7-SyV?Oyga8ZWM67Oz{;w(u*XG_Ckr
ziu|rEjsKx|O>zAE+F;wXcA&3WG5Igp{u<qncR~2Ezc{1AXT0GXt*wZkTM_@=lK5|m
zFAVx+u;p1jIRCBKdvtC6`WsYq7S|}XXK6UmHTZ<9?ljW=F^7)3^tN5B)kn51t*MO^
zls9kHnzW|ywk6?cm_y-<wF#AhsR8_Ef;1I35Y|Tm_0JrdP!YHk*b1VBw?>-6RS~mL
zRa3*ZX!X^&-;&EfQz{@YP*Wdh+-8>IMdeSXhFdiwO{M_V_@bL*kr;%Kf?A|DQm2KY
z;mwQbw#@p-c76k=V9g};1FLZNB;-$M4rESh#$A$uhRQ&+IufWPBHda^#EQw8+A0(d
z@FK(wxFJCDI0koyHb-mfs{)hiW3{!|@K_J7&tiTnX-QpUbcdBjc}~A35+xN^;HFRo
z<Hl94wSq~362-qP61Bq#<v~v<FY7n!ySi>gxS5L+^<7=Jw8;vS@FkI2LQ)wzxE{q7
z6;yDwR!}hgYSHIpfo270qW30LOx4EOk4&wksVNd`3e-l*qp_w)Al%sW?7`*$_zDUN
z#$jM&(!df&?-lxR=VkN^!E_)FzjfilN#yVT<%J7<NcSN<fb=#@A-_Y@Hew<cA{~cm
zScbG6X%y*SknThJTe$YaNFTr&a|-FNk@g{-gGG5L{KHP9^N@ZHYwdcZ+u_3RM*0_|
z2a!I8ttUNq@g7b<dXY|m0}Q0W<s&UZ`XbUzNLz6Zy&LIANDm;Li6BBJ(psclNDm^V
zhe~qbeF8|oiL?mmpO9`sx*SW>Zlt+bhYuj#i?kEzDWqLUgIF%ourIh3>)KSL$C0i;
zIuVO`Ez%g$dyqbc^iiZCL`{E;v;bE&dy#$>sUQ2ATx>N8k=~EA4Cym)dW}dAV{g%p
zbSS*sDWr4^5{ahYa<!&!_v^#m=|;Pr?k7BM7SMXpH>%f&H95bgjhFmgZ-~Csy#A%$
zq04<)+l<}X6~iu@GHd);RelLFcSA01ztur?ZNz`%zb{<)B_WHw{<|EDJub&DppnpJ
zC(z?DHjXFIpW*ZXaQpH966keSdP$4ZQO;R3;_1ZyAdIV2Jfz2I{NDf`9#k{+IAiLO
zhOs#fWAdBe>(}#ZCVf2UPk_EMfgUvJOF;iM=vODuSDN$+(DN|1mnF~_ne;uNF91CY
zsae0%r4z$cpuPtDM9hizt>stV=e#$qjossFNq78`iy|uT<JZ974gS4We#dUkME&|9
z=qoX&zGtPAUH?5%&Iru8Ux05O(tf>WIh0-tfIo*hc$rmii(byjya`kb*Q*x%G)F54
zVfI6aS4#bK59r01yTw-ePN5Ug9|e6g=*O&d$5FHVk3o0htZ+*Lz0#zg1$`Ch3#{~Q
z{qo)L6*qx?sg)k<r%wfa&tECO0`vo*-()TCh?w=)g8m-p4<^t*Z_@7pJq(|<I)T2z
zq(2J!=Rq$^pkHaye+>FFpudnnf6k<z1$_&GC`}3UI+N~B)3m2Rzb1jciqiwYO$EId
z^hRrcW863^ML*NnDg(~~_|sWdo-1jr-osj4&T+&^eYqWc%iv#EAm2Y0mYKTj1HBpa
zO$qetP5NQb9|wI^0=?9vp91|?pf9x2w{tQdi+!Mf0{XWwCd~Ss2l?0z05uf;yA^)?
zr`Gb07li=es0|CiGZ!1`SOQO@sl!H+0eWKseTzxo3HlqL-;qEU1D@pX2Ymt7kOHeM
z<)k9|QP7uTP5HgGJmqmL!T+0}2e7u>K!`<N|9y^9@6dalCEh?=T8VekJ#2}$@Ge)-
zJFg{unYZXB@4RAfVX=48B5z=kcPRK5c@6$N$9Swo^B{-R`0V)10-ssnA8!Hme3N>f
zNzFxi=86sl(=%39n4YPk1J?}1;TMTVf_rGgjw4?rntHAZ_jd9_?dd2EOow{Diu?~9
zivEXBE;Mi$YYZK5J^a`Q-$92=G}^b*p`HQSEA`6~Xl<xWalcK%A&X&w5omJTbR3Ww
zwU3qCU~PyXYbOvYrbFRRNrjc(y@ISsHnV<gv3Ywy=Bx-%a+u$l4&((1E4^?k!0Dr;
zU~Jtvtmv|ue|Ta(S8l6$t{ewTWqPelD`mP>rgzEo0hvA`(<fy5yi8x0>AN!hNTx2_
z3qZ$UnO-E*xiY1@C+N6Vrj;_?D$~1U`hZLyk?9jMeO{)o%k*8D3aR*AEpk}y%{ihK
zH5AeFh=ifvWs4VI8kn>svbiQ)ADB^CFtcFBw7HXIE-<rj=8VD_GYa`uT#Po+$M4YP
z3LKBfAuOQ&U)s7mw1+VQ1b?d>qZ?9;3Bi9EMxx+9EypmwHqkHE%nAP2F_Q#;mz=}&
zh<dU^@XvB9{$X;CH~&?A=;FPQ--d7)9gnNIEkN^E$bULl@fVF3q-)jOGdV>@@IM1B
zg!~8Od1*@ueZB}hbtCxof5oqR5q?&ovB^5ZLB?~uJY+r+r;+|WusPGx>6%uW;na=0
z2<=G!h!Wj(H;33t^e#(dcq$h~*l<7&h9Bp?{8}&|hx|jsDep2m2+&>kAxS5@oJyNv
zj0THQkDs(z=@H;GqYOW3b2C~gZ;;DNyMzONql4JyabT!1ktp*y5HMmyxs(Hwj4KGZ
zj01(n8>IaL4$L!fhohFZkOM`=b^?kx5H!|NnPLu<838J@hy&}5c?2vOOy_f(jGL*<
z(m|I4s5V|DN|0;aXe=e5lmk)Ya|A5UY6kaqgU&?LR_0K%b{pLUtn#e`aHsJVs`M&j
z36d7${}5%hZw-LEjaC9Scy@uU-FT6J>xYCu*=O8EWj5yDgyewnF;Q;d>Krug2Xoqu
z8Ox9yGHCmn7RsE4<gh{Kc4?cqZk>jQ1Z*}Mkvyk4#)8Ds^Y<aim`2hr_GUhW)Md;9
zgUcX77WEtH<)(o_*oy&V9_MV=QMQHiM)1pdcO%1jbHPg;jFB;zC+UY!$diE|eWdF-
zjCwkbi$v4+yD?h5lfkKv97+A}=veO^#+j!R$<a~fE#}|?8Y+&C6~0PMAC*P*IXX7^
ze<uarD+T&$G<~#3uwn$dj2sAY8fg`vduCEK87?&sc)SI8eSZfjF`iS0P+4ckk8jJa
zA~-XXXwHt4w^6-Ws$O4}rn@Du-alR{QD{-Z)2~Dh8Iv35f~GyhU~%mnFmiStMm01$
z?QV~r2>)z$(vbLUc48>}|D2sp-SG)Vp@(LGi@9kAcxckOTs$wQk2f1LLEmir5nNxE
zYc6>nI`sX7_08oRqvSmHH)?vhn+T4MYBb;ii8p9n$o>Vg`ZcaAr5GL&<YI1sgdCK+
zW{4=aK4%qxQDQ=`$-Pb%m`#=9@MF5K>7QwXPLtfFuKSUEGFY(iXe6@w{vSx$OUGDb
zFP)@hXGqzXiE>=_vVPgwFOlq7uC-F`Hzk|Ce-2h#(s|Yhi`-e`q}-g9piB`!8R<Ja
zLjbjTHv!-l&b><tB{xNfAEZmN1z$DnHHHj3!WEI*{a(R^Hoc9yh1o}pA+tL>4(}pe
zM$98!oE?XDkpV|sL86=;kM1IEMl2Q3!Cj=zNTm+dIg+b`?i$R!->59nm1TpUrN)nT
z)8gXj2o9#p4Eh*JUGIl$9H(5P=KUrSvokMan${=xa&@M{ge8@-N0Jio=cG~gcmV`S
zp+Qoijf43*kmsc$SO+k6iji{}k3;^DG_r_`4@TuV`zxUeT=|v?Tmef38Hp8S${~u?
z`aul(H2&1Cbk3hcE##c-d(2-0?-1J~hu9uD#IDI9_D~M7S#mh|X2~fwOAfJFO5=+&
zIA~MkNK3neWQi3D9JFXc2rZfpF)bNQO9ykW13RN>>2UG9YYUlRl<W8MT=0F<#(zdG
zozKNlW2xQFj(xif42sciqR<e@SW0CB?sRH!$H84LZ1;4^8#jU@a|eEg1I${4Od1A}
zOLNl|B;D;&mf%5nQ1$)ygV6)&?f^;0`kj6yXPrU>TEBbc#3C{s20Gc}q2@R`YQ5(r
z^-WF#6~oXgsTlfDiL~J_JxS9_AX#<|t;^Ik_U-&Z%m`^YY4<c~G;$i!&}r|Wv?I$$
z{!4Y4?7+lcP3TA2fvW2gJ4<vVqw^R&&lQn@0DkQNmFWB5!LRF0P8~~mOr!MwAaJmo
zdd}HV=ECOp46=qnYVloN8XA~+9Vf4MX<0{!Kg*?>ycYyI^!=_^IAwwmyxzqHXU!*R
z=mGcHe9ouw2rRN5^nlyL?a@L0_*bYVP4mHW(zph)WL*2XPn1DmqJ{q9<DSXBqhY~5
zM!OwuALU(nP+Fhpx|*Jhz`-xKAl3JeMN7ndm^dr~jF=DGz&wjc&W`q785qEm+yxxM
z^l01$>SQ;aBC6>@Z_W~9G0m9G^X60`OIfGNI8wypQpV%gg}5R#2|c5z7Mdn>ok90Z
zBEje+PuXzr;P?bE9U15(znUhVjgtCHMa4A9R@20@M^foADmuJ=Qens-vJ1A!gwno>
zzDv{I13`B3!^o2A_7m1y5QO^llj5&9H?{@UrWsB6P4sN8dr`GRAMzlle;L2&_W?6h
z_ZqJNFunmm_hpca&<rp1*8N%3afWd#QL;U>EcvyJX9#-D_1kg8h?`JEf*$vwS(z76
z!PnGoB<m?+e%(z)X+M&2C&ACTZk2^eOFH!ZE2uE%elA~dug|)gn19q?_yWO?xgHxU
z3cVl_6n?=h{Mc|=*o}8ofqB|CFwtzJmt^UW&8601{qs~o)PJE=hO=*m5F={)o~kYD
z3sl+nq%NB@FR4l;e6$>UG)8s$XA{%c461Xp)QZkO>CpGnnU84F*8);2Uhg8{Gn`$r
z-s~t|`=vjk>J=3c3$MtpF%pweLKW?EWseoYSBnId+zqOzWZy7RNj0xzAC+>Rl$8jn
zgyT}(8l!Ul66&nGxcaw@5h7`H(4p_g3$Uc*yGB^m7m_;@)W2J)FU)rLP*Gr$=6{tG
z-_EP>%PROKNM1!bRk4ti*q$$XW@Dy@g9WnRvgzI%eTNcnY1$#xH2+{KvBii$<ZA)^
z;8-9Mw1V!r5cb(J+G?M}z-KrQC?ULi3CHc@(-@WWuOsC*F;eY|0U?rKe??Z`?*LIq
z-X!fKEH5OtE6M)%Nyu7mfg7cee)X1_Le|PNaUn#E&D}~!#&<}_GS>&Ag%(>x0xcdi
z1uT<8zdDo8n`M(_A6UkEqbcudve`V>Oeyb}DK8HM(W-eWG?4iSRX%T`=%F>DhvqF5
zJrv~PMVmKC@j2g@V&_S*+3>0hl~;w`r%~-wUGK}FU9nkv#H@X)GKZ*rs;s@jQeBy9
zO>P${IGmeU`PB+jrWisa5{8hNf`v+eb|IbIbRK-Lbkp}DOK#d8Cp>^~;Qu&4;8+J~
z&WnbWLMNx#ojN?I(K&h0Z9vJep5r~@nliWuD96o|pPIj)N+gE+2!*)|{6oWe{#V3$
zV3^3me$MB{P~K(y7ctqwfOLkT!hlq_83ufn@+u5SfEflPzzhQtV1@w+FvEZZm|;Ky
z%rGDUW*CqFGYm+88E7ED3<K_lo@N-3VsT~|khC_#fCQLfz^?-M6k$NRlfw)H?gwCo
z0lxzQW*G1rpqOF6#{igNKvL&ZgaJu6GYohbW%|Q_T_D-QfVYCdWe_2Y`i=B5!+^~I
zWEk+@C@aE%baw#fJ)wAWHO~+XD_7k68c;5SI@?o#;o)lK7>5tk^kKO-qXI`W@ae;I
z+bO>Rd1GiU?X4VVDL-t`9^_>>F?|kD{)|k6p9p*dLoCDXeiC$?R1o0LJ5G79^D`)K
zT4qMZzrcMO`6!o}nMGj^%?J>X#R1)zNPwFI4&xUDcsPK!IDVNRuX6<icqT*NOuw5f
z>A9Tn**x!{olKTK<Oq0NMkjtVX(!`(0NS|IJTC)bjjp&EFC`-?Sw#K_|31z0L&$V{
zJ!b%U{y=%Z@7u`xZ$)0eKYbbMJw63LIDRS;c-&p6DmxvT>GubS1LqVuJeK-E77bx+
zb8<!!^}#Ia1#AR!GcOi}>!l1q|F>L0UylmIRAH5iR2w!)aPgZzkk$9o33xUaSv7)+
zI6F?<p7Q~adUYnz934eD?F80ju14bQIIt^cGr>1`!W`bcE9d)E!4~h$9PYY3r;q5h
z9_l7%$D6n3u9daY6Et+>tfT_lJ!E>cD0jkSj@HyA(LwO&J8}e-PG<N(q}utQ%@;7$
z&f#n=>s<(@+BuxAWo4q~JNn_g461J(BZ+TK6m=~W2`umm7==dbGi05)-I9y8V|0*4
z8<f7drni&?UX^hU_1-60`TaOx)%Vk#kAv#*f>wQX-biRd+w~!jQi;V}h3T@yCX$Rp
z-%k!Xo5w~L=TRQ#5<?O2SUg7Pvu03Tc)O-AG1V)U63FM%LF$$CtLHMP+>ocxFnuAf
z^4Aw}K5d0aP~~Z=^6U7~7xHz&(XmlzuyCx<po}zFSRgc5pZmU)aEX*aRbDQuL>~=4
zL6uETr$KsXg5cqILLsZ~rwd_X7*CPIIGA%66}VV}ML9nuRi<)~N4&mG<QbV%lAo{T
zT!NfFQ!*UbB?iGwK<3X!9U`RRr$L}+aV^y1%J13YT6_RzYE)K$8bz_QL3M~sNUN6R
z%%B>45?nJVjCa`TgZee1Zq@TLNs6<h$d&U6)ip%2uP1wqkibT57`a)ln6o2@a&-2c
zrX4_0xdouRS7{4izx21(jo$yePsf=-5i>KI!9AwydQy1!@I^!s0S;zV0Z1nY=<xH%
zQ#v_7M>fwV2C(>vU79}Rd6ZzrkMV1eJ4`2c$kIG-kyH`YN~fq66H%>n3eYeS)k;5z
zdRV&eajKKf^jVh3)0sYV^Lz%$xfo?KwRG~E%<G}+cj@FencsH=)U22BgHTC^vvD<Y
zII>RTM|ZOfr=c(uM%6vQxE1C+VZ26$(=RdZOCVSWV!K4}x*rA7LncbsyrhSo<JylT
z!#Pw5x}O9Mg6Q1}bAT{;zyuV=ppru%`)-99M3|w#v`dWdE?7ek^KJi>0!YB|dy!zG
z@Drz2&fyh1;uZcw92`N{!{Cgwzeiw$?^X_O%6^K#n<cnDXTB6k7fI>J#jq^z-yHjv
z5c?v274iI3!9u|cI`sXy=ny|&d=^cJ0E5LRdj}O*?4}-otIp<Y(UPnx4whx}=+{y&
z6@<^C%7XkH-_cQ#^Cjs#SAv(0+#Q;JtvnOkgl;wRss9Zc7QP3mmf2jbTcuieP#GNh
z{svR6+48c%dZE^wM78EhwKi$K2+5ns<z-2EjY1scJ#EUHI7T#6%A1rZZ?Y+G&~qej
z7?-DJtctwp`-RDdjZTm?JW<vNQx<|9`UG!omeV`jo2}=h`-i}GgEgo3(yYs}=4X}W
zVz!Rgva&7*2b_ivS(r&;aDw0>3~{CB>P~M~mMEqJVR$T3=Al$ZzgWzyTu>L}W7_C?
zX)f`*$7rRw)4}H{(67i>3QWu=t%`;|uqdCH5UWv?(p(1u)4ti*iz!_R;m<;ga$ulr
zP|Q{kC5LLFF8t?fu~pS8^2sdTL?e{uCYZrzH-q19hU~=h38iyTI>BbS%4UNSOQ~8P
z6voF(efrB`k+o_~1|E{UG`cQ@)%ypQ8e<qye2h^?F{^xG6<1;n39U-E+r+S{OZsJ4
z1rj1tjlC>B_JBczcHr2P%}2DRzxh<$tIU>kNF{HaRRz%|maIxS++FrIrCZCB4u+b4
z8LBOrsx4W`1fwlzrBWOr`ot#v4_cxs{$F>K-K^A23B&VKb`u|G9@R9j-BN!~(=LKb
zf{*KWT%Xd@C9}Nj(oB1oX4w^W+YLg!<SoseozDlI>Sv$wm;S`s)TEX5ksC~~W{$E@
zuC<({REo)(XaqH3|6z-Q)k@Mu=Hh;t{`Cg-fcLkyX=1*mRkV6!+ijm?l}gs__k?z{
z<yj>wxzaNXPDYdVPvI10Nj%<=*gN@Z_73|it5tH>BG+a$FYTk&+EigcL#@G9-+>i*
zwze!+*p*D;7U29W0z7KaVR|Vhxe!}zFu?MDyX6P2&gQVRTRz32#HYI!!gkTAy+mVH
z|FhSFU%Jdpx0CAAW6e$UIiu-S<>k}y9i!0P+47@CRrRsyGYY5I)R))BDk2(Zl)ym7
z%{BPi$lqp}9^KIxX~qN9c(8L>1P``0)Z=kl9Qgi7MWnKznsnkk$l}xl*KA_t(a_S;
z6(vAm57?j6*a$7dq3$9Z#(c&gql$U(&!Q{?m-L;}S%4Cwnc)i|Y5H>JImEII5H}!M
z$1(~S8^=Z*Va%A${70DnBr9GFG1HZhyy+~fl3la`@B%i$F?JBojDakN(IsUbLotUF
z2Hgp5o!T$|c}$EiwG=lsg?IeD_N}R}iPnT`%fd}~QW@X8X+{Uo!E~5i-Rv=x4e5q4
zuOyp~6gKx#(k>ulJ<B2RPBx5^%T!0W_b~m+8kTRC9i32?Lw~g_E3Ia3<9ar@g}IL~
z-%99EO*(it0A0zlHZX6D4VlioM_4v`e8>rwSHQdr+0Y>Kf==2vZdTf4{G73gOh3w|
z6+zBeK7OX7+y<5rWuvFFu_)m=!tyY*e4?zro{c3D`}sh_$l?NbD}nkkHpQ8V@%Al}
zl|7vef&||JHW_75l$z%-2c>&2(>Jg&bD8HI;`V`vk;?}#hTSLv|6K|O=aa0!nYn;v
z9L1<@VSeK~Oux$E!1z0cu5oG<=@~G^)dTTjB&6)u`RG@5Ux_nIXPzfXAI}Ds7i0NQ
z#oNM0PG@6a&hco@$kWX0n2^gY+RpN+NgGAWoI4zgbT$|m7>3WhN~SwA+u6+3#&YJb
zVuL@RPD^*(kjpY+#&qLbOgE|>ekgu|H-fIH4Qzz<DSj|$Dci9G?NIk-I5Sl{#tFM#
zL`js*q2wH!N6DjXihveUvYh2oawp>unmt5ddHGU4tW~qxnSMRHxa<TQUPTRJXb_x-
zZv5Pl4*ludTuyCH(zfu|3;d1ln1g1hdl#H;ov{y$j{ZK$#qg#+*~8{u$1XD-rU5!_
z0V@<kb)4}Iy9n7)r`c@B<Qz614AW1uX^!zZY|(nQ^c)*>CmY^@9-NcQhLkXGJ8!1r
zLD_>b$AW=92&Q??jOE6i)Q5xm`!MTU%!6*kG<lJE?li`*4C9p9qhq;Sq*(F59_3H{
z3oD4MbTGrva&^_tIZS86TNwK}dO>$ZS&qy2`Bzx>lg#}jTnH<4X2OfVBZ;{nGQEQ>
zb!IFeuP)F)TIqQmEEnU@1ECnnpkTInD0!WaS@b+dN|&{O@s6KnPVK#{U^V)}E_x8J
z`bjpHW+!`6R_N#07Mk5QTooHX0N27Uq2>&SACMYT@oack(O9b#H<5<k$ue75&UEH$
z>1COBGS5cF?rdR$m32nX>95vH1&9h4GGR$So(p2C#HCJJ6$e_5bLFe*{4LDKEA(z;
zZrBkT>KFi{7PDOVj1Mri$I{dW@;_L{RPwA77NbUL2K4pAyM4g24zo<Gxf^lOT-WJF
zDW_IO`;mW&{fTuWgH1?HU91*w4_R6p+y3{iru5mi2CSe6!}2e<MW^Axje51o!6)6K
z(_&T(-}b04w27@xf*O3PSY}O1Z%C`P3agV>&p2AevSQ4Oz0H-z+l*~t@XR@mGIW@c
zOO)Z*HaVu`up#SNU^yF7z_OkMpL+{U7rDf{X}VO3>2iYk33^!2%z0w$W$qITnI2<f
z>Q1w)s&BEO-@+;~giSfkCLdv$m29jdn9Js$X7i4)VVE)>Vj(ERAaWapRH_iZ1ml&j
zo>|xmOr#l=eCZsA@kB1j!cAUYa_e#tb*xiUm=@fHm>GN6Ftm8wa&~EyEm_W{80Uaz
zgJ-aeE2Xp%Fw!|T<VwB@J1R^8TZnPeTxJ*(S-NqWWxz2zCgno+D9gPPtMX}<y#f10
ztZeyMp>tlKjl|duEP!3W33gEno7B$6V{MPJi>mtbli(AM09nPxZDkYM*_Z`v&ckfs
zdUi!SyNu_rEMRNbv(-n~l@GJwa9&Tcp+*55L;=eyWVw}W9@gr1c8OzJ9;JMv>TGgc
zj7}$eEN5<K&cp9A|2Ayl0AI&0DrAF?vizU3OzfTpp|MAp=Sr6MBpdV<cImscqzpB-
zvN1nrS<Y=(O)y^vJKcApd-#s`CBDGjmCUtU=TC}r1$>jRM(zEr^NSB``HetzB7llM
z%}}3!p{^pm&U~G_n{Trbg>2{t%<}>I8n2MThA?X#%ZxJpYBuq58_Sb+j`4O50}jsj
zJ&P19i0op{mH!>gpzBy>jE%UCYV>`;a;sQgl`;Jk%fX<|uV&d%RshRm!I?qHqhih+
z5vM6pfjlCXpXH2JFK4OPGcIS@)MM?}V}kgcjYqfw_wi;p2QO!1E@z&1Y50o)n=J<{
zW;c~u&IZ>q|5l971#C<U27Nmle1zq<uptjKw#_(YoN{Es?}kDRm6da6hobm`<y9g6
zk;+JY1)Se45s9O7FHYnii41KHH{&~y;N2Xn3=wPMNu4RB5%=ANLXF{yIa)<SjDIn+
zt^p^Mq4H|_LL~7u<CRT9UaY>krm8+t5daaz=<|_Hk!GrbK9n0+JmV7kla7IDb<tRT
zWMNgLKGIZE4lvwQUcFF!k~c7|Dln}w+`OZ{yt=8Oz9H5;Emps+roLhtUfhIt3u(>O
z4NcL2St~tF7jDG&H!JAnPo?$r0wN6;6~jCF^`xef()ifs7P1WL+k(?0t$uMsEi|qT
zS2fSnrdK!AMKBHUPNHcurY~tIkJUx$qs`OTG==ejP0CNBhK6gWHP(jfqj<~FV%TKb
zVlC8ATR}b1jE{2iDtKL@4wx*~6e@41jn&m__yB2Dfd3?DPnwlORQXnvryi+q;1Z~l
zM2~IdO*iXMCa%SMrqEBQni>_Whwc?KEwzwo>aWHXauQ0fw$jRLBjG0brY2s87#Doh
zaI;bo(r0tHd`qlRR1ZBWBh<5GCrc?lO2E6VrXmue_j7R?ABIA0-cvLVaALcJk0XA1
zOQ1T>1(k-Y#=lj{k<^R)8>E7_ygE|8g=d%IWm-GXyR}-hVGHgA4J~Rw??Dz>R4@aO
zjT121UCVRpYUpiZ{H0{-<z$c~Y(noIQ|DWv7-8@>vc=)ra@s{8;-Tnk_@`i(#_G#C
zwz9SXA7_>9*k>f1@bm*Fe-*vaOsZ<W;mleMSy9j_w*tdWYs3_3s1UC=6Qr8@t!wZh
z*hS&`igh&=(P}|M_2HT(Sd!j@X4c3dYYp5hTi@r!5$c6s^?1>jaNNuKUqKcMm0ex3
zX3d7s(zUA=uPME16+W5UDBd^53zyM%W=pg(4xlF^TWc_DrT;JIP9xz}W}+y68C!q(
zb=3{Ek=3G*dOaDZTHlq%5zx$1x%mn*^Sy40Nxku`vWCCQjH4|tR>L$~Ob(0V)eF!R
zlPbG4TpNSSiEa(mH%H2~@`id0jc6%me-%C{+`t9EWD&SP3=*0sGu8N7N1x?giYY^-
zsj!-at73JVG2&pe(t6RM<@|MUt0LPh?@;5@H56*z(Oix%0f(AtkYkop6&lwH0xBZg
z$<xiyYVi`bQUVu;8>0<PT0?m>5}hHEnIf4*n^^i@u-c(&<&FGxby^)JNlkOV40sgU
z*Um$2+&f%VA5l{$6as-SIO3bg*I2dz_C5E_F??JZ9<g3D7q`}^=0g3=SVfvK0{Ao%
zO5@ANG$RE!Z@_G5LfVpZ<0{RXi7f9<snBYAxgX5YR7D>$hHjjL@77f}&@`u}a<@`Z
zQ-zN<SHh8)vtKks%#dZ~oHW1k9M!B0#lkTt*`j*>(PXGluJ%XP5{nuzKjE)$E-9}@
zUsd7lh|&$ATj5<2#Fx?3#pdj99SB}(2XU5eKXB)r_FPiD6R$J2bJO`-nr7fU&VD#G
zX1+(mN^@#1O}%d;5udK9_fjO{4e_N#D-wAp-iD#xlVK$}abJphpGG1cSI_-cjL19j
z){#I`yhrQJPR5U`IyuSkK25!+#7cEirG-|E$UE^R^6Ht%@HtxVtYr9H^SvilsuOR=
zQ|~)T#OG<r-mBr%25ahj(uw@IxuML85qT%xHqw_ApKpFg-CD}24YPt~-l@r_WK4|6
zIQd(uY)F7kZR9{C$IF|=ZB(Zw?=`pJaRPzS`Thk9if7&Ub!c9#(AJOWS&rXi!;g`8
z^<En;1V23Ns^x3d@`fR*kB-Io*OT-MJ@7#FcG=`_<@^Kt_Y1%$)Bmd}@ZU~>Kb!*p
zqZIi6Oo9Iy@H)n`x?f&RkFyf5Zp6ok=Z6ypA%68dfWqfXeB%tkM|P)Uy2Pu;9TdJq
z;@c&^l5@SpAF$zXmiR*w?-AlPdQAlB(`mzhS>jLH@Q+J;m&9jFIrPX6$x)9CDLv`&
zuw?!9yA*h!lk)$8%g@o)P)OW%B$MytemQr59~cC@t_}19BPG9jkCLi_-gk-%!2|r#
zVlHQ(e;AEh{v6sszq$eZIT+`=rwhr{PC9BN-hPW~l;Z`D^>=Ry{6i`5-<9(1PinuA
zg8!`)_}&zFdh9u1RZ!+fF~Se@lL7GCraCWChj3sr{HzrGbiaKz{7lhmQK2F0UB~&o
z1MW?~N%E^(8x{Y}l3!gW*(>GHgVxFP`HJLM*CXb^E_56M-jDW<m-+~-rhPw!oadz+
z`%`PLrr>`Ucxtb@eG}&5$M2*Z`-6YyC0<2=)ct;5ZfEa+Tj|TGe;otY>rpBCWrmbv
zzdtdAeotmMdPtAPzx_csdLAnoe{Bl<?NZK-0nz>iBGK-Vc>4qXUy^wH^UvRwc>BY-
zM<w3=%-;(fFDtVozmfd*=e0gg!Jm=BJ|lpq@nhdtOa`9%UEMXT{KlnRj#snajb9@1
z>TUqK3{M9=)JFEVKV5nw#|s|oZ<FMIEg-PU&UXSYQC4=X1NdyZ<HI)oA0|BXSGOyx
zd3T88v1%0xF)$ZDo(7)uv_Dh$Q;D~~T;LsvSC`@Iq?~gSuP)<eNy9PtXOeG!VF*3n
zM)-jA&o@Xp6F6Q1mW1v&Cw}{F$ZJ#Z)0fkejdLI7ogenuXS3T_slO_8v8<aOzae|7
z8?u!>@8S5Nnz~P1;qRAtb(6QU^M2r|Uc0|KDEZZm<4VpCfT#A_Z~CU&-Lo-n)g9|f
z{<D(ben0eU5^sO*<PDCu7+7Y1EBWnDcm7%8?T=c!aVKvwJCEjguldlTdah!23jUH5
z_zfxWwNk$Q0iB&G`0oRr=F22ozkdn%z(6x-bl*JL-~LqUk0k!k1W`ymwesH*zj>^{
zmr0qwmw5Y|ZvH6o_P5k{FaW534XLN<-yGnj=mc%20Z)3WtJW2gdO4SquWgcY0P<tC
z<hQ@YW|PF*U(ltMH$|JHh}Ix3hqz0m3bEQqQz%*&Dn}>>F{Mzbq9IgO+psxY8>&Fi
zra2UjZAXNuuCW&P?o<>MCg(taxV|QYs84tYBKXm!9f&R<vKy+1)zu;5jqtk(As|u{
zO=Lwdr2if+L`!-6i6Vz0uwf5js+%5oSPaA@k-&l(Gm@i8{{eBsg850Y^O9mONs67D
z6gwv=mLf<4b(obDJ2NRZdE*N0!7h1aEEHOLb@9rQP|2z#h>=o6%F+S21hHhps^XQU
zxME<XG)ECA0}2;QLM1^dA-LpfEwpUKRf~#Ogsxh;^qP`2p*6*eR+Nxb-1=7@DsPVQ
zp<s<8(shBNWh+V-Ee_4ZOD5-Njj<?Qlu5*i@Ki_y9TP|hzzqzS<9@#S^2Qwrj7^cc
zhOLQJakjd~1h!D=nw6pc!M|qZVjB90ts*1`e!EntAFfCcz$1AH7@>SZ+))O|`@`+l
zaJPzfhA`MoK{MuxD7lp}6rqSa;Z3%g^Hp$Id<~-#Gq*Xh)he*aCx{t}wgv!2OwS%j
zMe7lN?ym`bWi6!Qj9jdWhKEAo=4M2nLzUs0+JvH&=;!7H3@>N~%F*_O%P(e_Gx0Kk
z3Gvb@BreArAEmf(R+FeTXSGDtD^XXXDu{F9vRi_7aDA|KNDO=#5A6?^TO!2?m#cVH
z@>)YBzg|?CAXT+l#qL9)s=9`HQ461q70nHyYFsqI=o=VB=GQcF$47mLuA+&KhXE<V
zzXBJVq=Q_Oc>M>B?Z5aWS*(K!bQ!4@fqu)}CD&j^neoM#{?SPj*E$$1QR%vYw!}zk
zjwDXI0oR9A{8>f(L!r&h&C>h>qVPPN%`c2vmEdt~d8MGT0asfGnjWLTa$SaZl)R|Z
z6roS&B$P;cB_$kfs7c@=`{HXs2{>{b;Y4fMFMm|>>6TC`6oU6rewo~oT9Ff+k{Zow
z2&KF*!x^cy>Z>^9dPN^6yI{#L!iGZS+r#9u!?n2LnV=mPC+9i;D3G{Tlfc7AJg<UZ
z${29TTPtYZQ5OwwMjCArX|>AX%1ETCQ7fo#zz3CZlX$`ASWRukw3-TuDqd8ItLs%7
zr&WiWtF?lP9rchWQrr!|DcB_N>s(eK1X@$1HcSO%uCX?%74Y$0fNVik1HdR<SS+Bu
z3z!<Xt1XCB%iT(K1%xO_6cyVRQ9|WWL%6P{9J(|_NivipmtBCtj1LVp)X^n+;3Au2
zRVWg!ufi@%f;IJ(4Jserycw4$6-47p<zX(l+xz<yij>nS98Q<`p{~=YbcDpKM4s*`
zP#cFFwzy__bv;L=<4sPHQTtqlrYBOZ<<)f`l?D`}Ozq{D0RIu*Ew1!e*N0RpPlIuV
zNE~W^EJ1A*GKeag`p-jrhYpo)k}yt-xSdpaoRRb6TKp1=s$X4iQmJ}*f#4Q0RGw0L
z{{_7}U6ohYrBoV}%u0X7r_wE;)4SLetlr12(neXH^d~)){t~U-f(+SIl~>ofRN5;U
zRsE{GYX2Ttev#x;*TYnLKnkYvi0Yd4tMc~(M`Mcg<A*;nr83sT-v4nBt>pu<yh=~n
z2CPlI#-HDY3{}6n4yRJpUZsPARr+OH{pLxsJXyZOc4<k~ujGBzR$g7#Q)z}`l&QV_
z-<0K5|Eai-N~=wR5MZzWKWydIbwZWyvDwdF|I@be>UyF|DL!L8l-(5W1!U-&hMZ>R
zb;lD61&5tq!78QeL_SM-b$#;Wr!4;~lvQ(sYE;)TPs#FXVyI@QdX*jbp$wTw>94ML
zp1Ojh;ax={F@IHgmA(gxwY+&<v`7@`GfRq$DyPy9fVP%ban8OX!7m?k7XW@$IfeM6
z#3}uYWE^39vEW#*7-br;!Jh!8`U`D%RequkPrP)flwM4#@`56MuPhenZX3}CYJ-8W
z<EZ3I(jon1w2-5I?Y#6}NIQ;rAKENRJZ&?WowqPW`Ss+Oag0|N>;yQVt_8BZ-F`cF
zhyt(E;%GTk{nper#X0x-J4B&QnjARB;J*a_syWL4l7CTskF0qN?GV|URRKv+^(w|}
SWRjJSKO{Jo*$UV|?f(I4zMJO&

literal 0
HcmV?d00001

diff --git a/Code-C/main.c b/Code-C/main.c
index 1bacd00..1a42c18 100644
--- a/Code-C/main.c
+++ b/Code-C/main.c
@@ -6,14 +6,15 @@
 #include "queue.h"
 #include "average.h"
 #include "growthRate.h"
+#include "getArray.h"
 
 bool rawDataWriteFlag;
-int nRow = 500;
+int nRow = 4;
 int nCol = 1;
 double freqEch = 250;
 
 Pqueue firstRawDataQueue;
-
+//Captor                     1    2    3     4     5     6    7     8
 bool selectionCaptors[] = {true,false,true,false,false,false,true,false};
 
 int cptData = 0; 
@@ -66,31 +67,29 @@ void *threadCalculGrowthRate(void * vargp){
     Pqueue rawDataQueue = firstRawDataQueue;
     char* fileName;
     //double pw[nCol-1];
-    double *dataLign[2][nCol-1];
+    double **dataLign = getDoubleArray(2,nCol);
+
     int i = 0;
     while(rawDataWriteFlag){
-        //add case first file traitement 
-        //(possibility 1: call twice powerfunction , following while strat after 2 file encountered)
-        //(posibilty 2 : initialiaze with a 0-fill array, move ligns for first encounter then let the following whil e do the job)
         while(queueGetNextE(rawDataQueue) != NULL){
             rawDataQueue = queueGetNextE(rawDataQueue);
             fileName = queueGetTabChar(rawDataQueue);
             if(i < 2){
                 if(i == 1){
-                    powerFunction(fileName, dataLign[1]);
-                    growthRateFunction(**dataLign);
+                    powerFunction(fileName, dataLign);
+                    growthRateFunction(dataLign);
                 }else{
-                    powerFunction(fileName, dataLign[0]);
+                    powerFunction(fileName, dataLign);
                 }
                 i++;
             }else{
                 for(int y = 0; y < (nCol-1); y++){
                     dataLign[0][y] = dataLign[1][y];
                 }
-                powerFunction(fileName, dataLign[1]);
-                growthRateFunction(**dataLign);
+                powerFunction(fileName, dataLign);
+                growthRateFunction(dataLign);
             }
-            remove(fileName);
+            //remove(fileName);
         }
     }
 }
diff --git a/Code-C/power.c b/Code-C/power.c
index 8f80d43..836bf0c 100644
--- a/Code-C/power.c
+++ b/Code-C/power.c
@@ -11,16 +11,21 @@
  * @param powerArray array where results are stocked
  */
 void powerCalculation(long **p, double powerArray[]){
+    printArrayData(p,nRow,nCol);
+    printf("data in pw : \n");
     for(int i = 1; i < nCol; i++){
         int j = 0;
         powerArray[i] = 0;
-        while(j < nRow-1){
-            double aire  = ( pow(p[j][i],2) + pow(p[j+1][i],2) ) / 2 * period;
-            powerArray[i] += aire;
+
+        while(j < nRow){
+
+            printf("%d , %d , %d , %d\n" , p[j][i], p[j+1][i] , i  ,j);
+            powerArray[i] + pow(p[j][i],2);
             j++;
         }
         powerArray[i] *= invTimeBandWidth;
     }
+
 }
 
 /**
@@ -30,7 +35,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){
@@ -38,7 +43,7 @@ void powerFunction(char* rawDataFileName, double **pw){
             powerCalculation(p,pww);
             appendDataInFile("powerData.csv",pww,nCol-1);
         }else{
-            powerCalculation(p,*pw);
+            powerCalculation(p,pw[1]);
             appendDataInFile("powerData.csv",*pw,nCol-1);
         }
         freeArray(p,nRow);
diff --git a/Code-C/power.h b/Code-C/power.h
index 043857c..e75d948 100644
--- a/Code-C/power.h
+++ b/Code-C/power.h
@@ -1,4 +1,4 @@
 #include <math.h>
 #include <stdbool.h>
 
-void powerFunction(char* rawDataFileName, double **pw);
\ No newline at end of file
+void powerFunction(char* rawDataFileName, double *pw[]);
\ No newline at end of file
diff --git a/Code-C/simulateFlux.c b/Code-C/simulateFlux.c
index 0016a69..9c0fbd9 100644
--- a/Code-C/simulateFlux.c
+++ b/Code-C/simulateFlux.c
@@ -119,7 +119,7 @@ bool writeOneRawData(FILE *rawDataFile){
         }
         int lastIndex = lastIndexCaptor();
         for (int i = 1; i < 9; i++){
-            if(selectionCaptors[i]){
+            if(selectionCaptors[i-1]){
                 value.octet1 = buff[3*i+1];
                 value.octet2 = buff[3*i+2];
                 value.octet3 = buff[3*i+3];
@@ -128,7 +128,7 @@ bool writeOneRawData(FILE *rawDataFile){
                 valbin[i] = buff[3*i+1]*256*256*256 + buff[3*i+2]*256*256 + buff[3*i+3]*256;
                 memcpy(&values[i], &valbin[i], sizeof(uint32_t));
 
-                if(i==lastIndex){
+                if(i-1==lastIndex){
                     fprintf(rawDataFile, "%d\n", values[i]/256);
                 }
                 else{