From f67d0375f9aaaf9c9c7cbab614dceb6937cf5641 Mon Sep 17 00:00:00 2001
From: KirillTepechin <kir.tep18@mail.ru>
Date: Wed, 13 Dec 2023 19:48:21 +0400
Subject: [PATCH] tepechin_kirill_lab_6

---
 tepechin_kirill_lab_6/README.md               | 106 ++++++++++++++
 tepechin_kirill_lab_6/result.png              | Bin 0 -> 19320 bytes
 .../src/MatrixDeterminant.java                | 138 ++++++++++++++++++
 3 files changed, 244 insertions(+)
 create mode 100644 tepechin_kirill_lab_6/README.md
 create mode 100644 tepechin_kirill_lab_6/result.png
 create mode 100644 tepechin_kirill_lab_6/src/MatrixDeterminant.java

diff --git a/tepechin_kirill_lab_6/README.md b/tepechin_kirill_lab_6/README.md
new file mode 100644
index 0000000..1825437
--- /dev/null
+++ b/tepechin_kirill_lab_6/README.md
@@ -0,0 +1,106 @@
+## Лабораторная работа №6, Тепечин Кирилл
+
+### Код
+
+#### Обычный код
+````java
+private static BigDecimal findDeterminantGauss(double[][] matrix) {
+        int n = matrix.length;
+        BigDecimal det = BigDecimal.ONE;
+
+        for (int i = 0; i < n; i++) {
+            int maxRow = i;
+            for (int j = i + 1; j < n; j++) {
+                if (Math.abs(matrix[j][i]) > Math.abs(matrix[maxRow][i])) {
+                    maxRow = j;
+                }
+            }
+
+            if (maxRow != i) {
+                double[] temp = matrix[i];
+                matrix[i] = matrix[maxRow];
+                matrix[maxRow] = temp;
+
+                det = det.multiply(BigDecimal.valueOf(-1));
+            }
+
+            for (int j = i + 1; j < n; j++) {
+                double factor = matrix[j][i] / matrix[i][i];
+                for (int k = i; k < n; k++) {
+                    matrix[j][k] -= factor * matrix[i][k];
+                }
+            }
+        }
+
+        for (int i = 0; i < n; i++) {
+            det = det.multiply(BigDecimal.valueOf(matrix[i][i]));
+        }
+
+        return det;
+
+    }
+````
+
+#### Параллельный код
+````java
+private static BigDecimal findDeterminantGaussParallel(double[][] matrix, int threadsCount) {
+        int n = matrix.length;
+        final BigDecimal[] det = {BigDecimal.ONE};
+
+        ExecutorService executor = Executors.newFixedThreadPool(threadsCount);
+
+        for (int i = 0; i < n; i++) {
+            final int rowIdx = i;
+
+            int maxRow = rowIdx;
+            for (int j = rowIdx + 1; j < n; j++) {
+                if (Math.abs(matrix[j][rowIdx]) > Math.abs(matrix[maxRow][rowIdx])) {
+                    maxRow = j;
+                }
+            }
+
+            if (maxRow != rowIdx) {
+                double[] temp = matrix[rowIdx];
+                matrix[rowIdx] = matrix[maxRow];
+                matrix[maxRow] = temp;
+                det[0] = det[0].multiply(BigDecimal.valueOf(-1));
+            }
+            executor.execute(() -> {
+                for (int j = rowIdx + 1; j < n; j++) {
+                    double factor = matrix[j][rowIdx] / matrix[rowIdx][rowIdx];
+                    for (int k = rowIdx; k < n; k++) {
+                        matrix[j][k] -= factor * matrix[rowIdx][k];
+                    }
+                }
+            });
+            det[0] = det[0].multiply(BigDecimal.valueOf(matrix[rowIdx][rowIdx]));
+        }
+
+        executor.shutdown();
+
+        try {
+            executor.awaitTermination(1, TimeUnit.DAYS);
+        } catch (InterruptedException e) {
+            e.printStackTrace();
+        }
+
+        return det[0];
+
+    }
+````
+
+### Тесты
+
+Тесты проводились на 16 потоках
+
+![result](result.png)
+
+### Выводы
+
+* На матрицах от 100 до 500 последовательное вычисление быстрее.
+
+* Параллельное вычисление работает быстрее на матрицах от 1000 n.
+
+### Ссылка на видео
+
+https://youtu.be/ZK1qSvbUNTk
\ No newline at end of file
diff --git a/tepechin_kirill_lab_6/result.png b/tepechin_kirill_lab_6/result.png
new file mode 100644
index 0000000000000000000000000000000000000000..9f2497996be1c9437939ee1a303d53f6093c3ffa
GIT binary patch
literal 19320
zcmd74cU)6h-an2yDuN<XL<B@dKmyW>bd{zwrAv(g0U`9zAz+~h2n<3Xp$bUvy$9(E
zq4zE=bOMAHT7H)qcXxJnc6a7^J<s>``#1NVbMC$8oO?g-mVjrf3YW;PlaY{+Tzc|Y
zR-J_8lq?C!X+_f0z%QegS)L>$O43hcA8NW9t&SgP^bRpgA9qTB<UK3T&2pFJ6YnD(
zUB79>Bj2lxl#Qkvr|YN5D(_x3&B-pjsSa|WHE6uqMKi({Ot*hb5T8tt^6q3$EdE@q
zA`_+*+KFnp{IM5CD)jVYcm#f5+VRnC63XvHADh#!U43}ZP@aaIyZKZ7#RsH5^sgyN
z1xRjDnvp!Zi6QZMfILO|1~}oiG14jYd1n;yXi#_h;3G)L`>NEqWX=Wn++_ue!nR7Y
zmA6rJR!$rbC&QAc1sv~nHeWN-%ID!{!A^QVnRI{Sz)3P3QIGK$5KW!ux;F=ztXDRT
zwJj05Z)4Hb>f<Mx_z50+(D9sEM0*e!?J-%nyz*FMp!tTp_1T*-#>mM|l})QC)qRJ_
zdV`*JF_egk9KFS$bs}WIspxQ8!>|Y2+?RQ(D!(hbmXLBX6itJ~aXx&O6Te5L0(GXr
z`g|^C;aIFwr>552ZWj=wVL}+uUpEU>&G1;AJI^3^9DdK~uvBGlGuwD(<Y58rYu7+X
z{NjCmTPN;agr)OUvv+e@_a`G~eEa=VGWXB<ki*;?-bXQ~kB)romH>Hr8!eI-sdJKS
z0K<An^Q?~4+L(8y%rg0kN)oGcw~ToWZK5;Lr^q0wJ;89L_|VYpfT?=M(qD&J#-#Ns
z=lP^ay)$En9N&tlkJYm{p7kVtFJ!YqHk6#E<gr8wDSUXB6>?li{`9JrzncgQOXACG
zYkh$Zo1a6{7<9u7H2g|;=HS?sKTizfA#6_J5wslRI!;*vi~y=vX4iQ>w$e8iyTps@
zjzv)r*8?niZM-tZ$R3m~K&}$w(#wV<mrst$R}Q>L7Ywj;ms+p)U%S3IIVU%a9`Jd)
zb79fsR_JvVp;WV?U_+PcnsCm8uJ@0C=d9_AQtaC;H1Y!1uN+)k6jlxTHVA^hC@M~|
z{h(`bP}u5qV-?ELF|?F*GTM8{xz~Gd|IBjxm-Xn(U=}F#5m^jpAz7(ihom~Jsgg}}
z)pBALqC>tO4<sbvqbJ;5cUuO)JBPT~sRz$EOtEsfo}eYc><zmGIa1hWhH*`za*tH;
zlIgy!_2K9<?T8pk8f^Hb1}HQHrT;wIUzR4$;PB@uZR--<J!0^gsRqUe-BI!R6BP%w
ze!{0VsK~@qEAmQX2#G3oq)pe%GB**U-b@`u-U;Cexf(XrSBK)o<8UkSzICaO4E?Q}
ziX^s`y{+&<J=0S2dI!x>+GAiQPRz2xC4xfV@pRfkM0NB;7(vUwO;TiY-8i|N0T)wL
zR>|C7A^~sXdcxxE1-0OjMEC3Hk$wZ8e|lZf*8c0&M^x-Yw|?v(j3QEher(NmajLYC
zQEc6~N!yeRp3yfFrHXsB>ziX6@yX+hYl?5HNNMb<!L&|%^hL?!gmSG7s#oRu3E^A}
zq8qQi<rpzE-2bU%GpW5tJG$hkDs9QboT+|u>zy%C%pat-m4PjL3oR{rR#X&X_ssYb
z3_GmEKsO<>KWqot!ARH&4eY8`8sZMbv!H9WkeZ_14jVLfYutWU%mnT@UbYfHUbWTP
zpS<)o%q|r=z_PPy!`t6C=Ef&zewhsRaP&@D#|J4JK`ou}g;H)0off^*<oj|LU`7WK
z^UGF5!|*9ZT~m_--nOb79g$2UIcu%dE}iO-d&{uR>^b9`7Aa=(c8&>8H{l$<DAQ)Y
zgd-!C8Xc4DLD7lI5?P2*P<sjb{s=6&vNE#(`ksIm?zb0ng>PaiuV}A3?09e0Eyv0`
zWlyxaM^-v=0P}gj=^B<~70zaWtlgsjOurq+fwMD7b^SU(;Wo%6KS|!<=!Ip<+d3|?
zcdzMRM7_djI;F?i5WiOW1s~7w*RDz%UBwpT%a|oziN&li6DotL*{%x!dn5O{-Gw5N
zBN;UgQNl>QGVz(39)e@Hf7_EvA~I7jBS5wCh0z?U#r*J1MEwe?a;i36bV8zuochvn
zLi(o}YTUK-_1f@&Ii8i%y++rHb*7OVy6gg>J}ZdZpz&7)9;jMQ4dIF0pc-zkW5}CI
zaJk*Hox7NLqZ+>?n5~NI`Mi`MNG|=>IqjkBe9}YBmW&4l<m1V1Z&pESPM-{T#ha2d
zc_|;ykkEOe?pZ58R^16U80HDhazJK`Me=QA+pJntx~F{pE~~yyP^r*DIe-HPO55EB
z{c_WKu7fhFqqfnP4eIymt#F)(ReWalV`~EJO~N8+Y<g4Hnvd*81o|DEJSa(nN1^~9
zBEm+reeNU5&G-Q6qJc9WRV(3hit-grE5u(as?TYY*Jd&kBgpEu;Q0!w9C^N2B*2|(
z?pkoNE9SH&?X;ssn>NA+LoKK{_jRXg#nVcoJz-{sKJcrPPrUX+8J(&UZ^$X(cjlxt
z<eTjiT5HEi@l!4&kjM^95jGIHWEfjfd$)s^XYZ+$LQM0Ss9_p-HAMD1Ql#v|zWQXo
zOE$Dnrgo1<yms&NUCIxebkKUwATH6E(eMy6i|yd6&2y+KCkOMK*4x`D{MY@<w+qBa
zIC)%cs#%qG%s~FbE@l(SytYK&AYNwo^OR;59T*9Y&GR2yyW(p-xfRAQt=_3IeX>5Z
z2L2`;k1+IXkeDNjMRDv;C)45%&lrAh${5V%9!~Vmc6&SDvUEh<F*YC?&n&!;lSr|+
zq^0ACdNVZZq)a{*-j=C{kXz+68qnT>jm?(#6Rh{-BdhZqkmNYN>KU-idaRqi=#9x(
zQ4@9;@Nw(BSK+u;9~mPfQSZ42tniE06*!YSJ>x4{%hr<$iXC){)U~f5sZS^4?7UG&
zX~`Y01p7)QP>!$<6?fG7-#qFbo(r8H3QjY&ah}H3Onr$(B?&<NB_zaq$35Ha<}Be^
zbK8AXoV7noP4LdHbBCfBRQgl(BQvycY3Ao9FWpefwCK83E>E-+*1q=Hqa$xMwoQ9J
zK@e01(&*KJ72t}#&xt}++`F9TczfOMnb;7=SC)Imt>!T`Z<!rM^eJH2a!M1V%zRA`
zVyO9ppd~V89y320AIhWDi5i|Oh%($C%P~(=3m6aQCoeWHZV7pHcjnWfjW8m`)<dZ3
zJLP~jlg63d4*CS{^HF_Icr1lSkxfE%=so*&$STM1>Z=WxH*4Cp)^&^DYi!D(5cl#L
zE}flN-3%>MG+OhKTVBJmt8NTqwV^1cN4#RgAWoN7I6lk8;C=x?ne)YrC_V}iCv6@n
zgPeI3!OQFrP-5!}*{meCiiRAGNb-~?XfY0>nI4qvk;9BbZ1C@<b-b42x9kepwd@ro
z>f)CnZU=b$681U?qzQ}+6?(dtyChoX;dLliqrG1)XOR5yV&~AO%TA{oZywiQB!S%e
zO)>T_U&r59YsO^|`5cbdGPgaEYEZHWwdj??T>7f-b7%Dpzg}l+&B<4|?L%LZdXpbC
z<C)c8XRq{5n@unm_}!)6U-V{A;lp(S&PB60OiXIhhvMbZu5wUW{X_sCCS526JwTIF
zu0zJcJff)S=I56byeS5m$s8l|U?1R39*9`d(`fZb;oXGf%FZP->da~HO&MQ=J@9cV
zQiECsZ5qkZ6P{@&kv$8^a^O6)^Pn+HNpbcr)zgS!ksack3%d+q1&Pj`40jpWz7&SA
zn3a<bX`N;pv;!1XMS*KvZ(w?W^v=}NP}>D9ef7_wLc}^4AoQpF#b5IhHA_;H*0lP<
zF6;ljA377@8Gwq7G@q;ODz=|A(E5Wi$5Ic1YrKZ-a?2B>mCebpi}-OaV(kjUel7DJ
zQOm}+hME=Q&nQ}w8j(mvn@;HsVapqc;MYt#jhN=a%=SmC`zdOEQa|*iC#iq*qip^!
zGsVoU>0nc!TD6*&{NNJkh;DqCIzW0uU%?DHQ}BvHE{)c}L%1UsZ;N|**r@OllYSUt
zF*Q0p8qb;L)g#RlbPZb+rxfFzvT`L-x?Q3cdOvq~_)4U0^`P^N@d%nZ>b#bqxrj86
zTw1I}+)K%oXC8Ck`c`v3os(lvp)*?boS%)%*>)4(BX`?QkKo5Axh%?+9D0k+l(~!^
zyXvF2Mw+9;A0p!SDmo5>wF6acuc?7fvmVRpnleQOq|ok`X*t_(;j|4?y`R=^Ymo`y
z@+bNjS}&Fi;`F5*mYHuDs#p1p`@x?-<Ia8uiY|&t(1GpLd4#r%7%b?I-!PVEC5}iP
zS+9-{w9ZPC$+6Q(zbx0%x7P0Z8f`_KoL#!KpBjaa-wv(u%lT-P`p^hTY%3<54p=Le
z&QiZguqe9K_1X3)YP--|XP@Go?Dz~(ViJ44>bi1&bLyH%LTz5V<K95vg;5<o%&`{N
zT`@He9Cc?f(z&o4V?(j#N`vp@yuo^oyt9+2wI?%JcVxk8wC|zzf;qx!60LN>m#KG5
zXSDsWi;WqSTrWX%8g>_dT56><-48Kx)nTrl=*2F>-$c-wDBquKBkhP+is|(!lnMba
z95;Du-ZAn*a}D8$Aaa-;SPK#<7NyI0x>x%<k>4Qzl}~%eZR!Bux&TiZAb%R=n{v<$
zt)b3tXo5%6qdho0iVT$c>P5nWuwC7O6<vXHFKPDWkKB@z7J@=oCwyx=e4&B04iaSW
z-Nf%5%}Cm(OcN++uM5_JIgFR@!Iw4KZO;x`nL~WOu^|k`gr{_42wC@q7S29#VnJ(N
zYJ!@n*n#ca_*PS3BUWHO`3*Ry_gZA0lUv4Mdu^=?|9a_iJa2tZjR2a?e=$o8!$7}0
zbnP)!MdL-MS;xg@cN=f?Q7cci)zDs9Vvj4Eojp}W9+5%xn8BmB`?<buZ;!wf){r{5
zf`LVW(=+**pgEph{Vq1?xvYGm?N)*aP&BL|TA<Jl;CfsBeiJSnRtei93G}7bX(onW
zo;m)bE|cnc)bVycj^OMoEQ6LznE=J3n2RrZ+j-85?pEiNZNA(iNFcndzs$ug-lrZ)
zBBYqKu6n}L1SKCPe*<~*&K2m#G#Kp@_5my^#A<fi&>6Hlz9$bp&Ui9;$X_(Pv=Jh#
z1QnOR=SfXjNy58&5OA9grpA3cd2Eoqr<0+-xpR(*u<H6sLn3k}pvNVYe9h{jN2aM>
z>2<4T8eDU(>5j=wYo^D!TRC%;IaKdDQIp_fffprkN7z_vZzlDd)x}_2x2e(Ay<X@X
z?LO47L1{TRpfn#(Fj1N5`M4<hP=&s7M`Xw8d;{;=w(lEk{rBGEQZ%2Ru1Wyv@%b``
z$za=JPh2)psq(XRt4)>uog?v)%7jx_8W<~(ZP;c~xtD!g!CQFw6iTcsS4;Ak3p3f;
zT_4eNG;!`z0dX$Xe1uLBNJBNjSY@YQN7XB%tZ_10v#)EQ+nvJ$f7W7KTjwkBQ=eM(
z64J;ecp^d?@%ba!GEvuot9`3<@+g_uP&m^Yp1I0d)_Twkyc$!SWt1?oV&#Yfo7(s&
z!YN2dHE)i?71Uk0i0nSFLe;PP-}Vu`(Y7JGS~YmqrSV;r6B}k*&+4vTbU)87il0B^
zK=XJtgvKna`ThcxS@csM+mt|NI2EH|<!VkJ&t+*+uu;F6jdfyp$C!Ef(?|eruua?M
zW=AKot{vrcUFSMRXErc;o>4w|Jhn9slV>i#`fl@+!Oz@2s#a6G(9a_=@vZ7?qKdsv
zh=R&w4R^UrcT?PlU07#ISL(IB`s$;tnacul=2!a#2jkkt-Qm7_Qr=eW<h9|N-%R8y
zdvElCu-baB`}v;TDE`^Qex4oYf6!ETuP3<Y+dyyFekNJbo_^wb=O+aXOGS*V`!eF-
z#xrCxO_<sO1Lfj;a9|SoT8DcaallqGb}OjLLBeu>&7zJeL&prC&VGncy*TH3S4^^O
zwx@zz--cqNke8X*UfK~{!iZ;Qt%V!KM`#4+$)DFNegkf~BwWT@vARZuFd&Rg8PbGW
z4WD|@E$7YcI;=COeg@ZEkb0#<eBvWrMxHUU%=^{s=0MZOKq>)ImkHRI0zNp4!NsLb
z_3RomXlGmnNgIgm3Ouj)<4GM?gJ#E%zS!ez<44evYpomUD^)JQWR+Xu?%9o1og*3w
zRv>X^ZPqadj6xLL7Y<)ejz?!kS5}=V%H8&z2+V1zH9wA>LRaH(=~?`Q`0?yReP&nZ
znQ-j`GnIUotF}aWLYcYtS({B;K92t8Wp6zxS6BZWZddd@Tow(sos%pyzrMPx82s&{
zsg_CrJ*XnRz}u#q03S6_t-#sTGm0@pp-a3c(mg-jr3+JAtE%9u@60-fa9G?(!$EN#
zJm$XD!6y?eHz^4!XSs1%G?fm2Shg54ou7~rl*8C5JZk&0E92}SrU8~kvBp`oH`RZ%
zWH9`_$6&W`D-6S$wpaV*12U$#uD15>eq(UA?VC>5b0VhkHQnM{@un4dP>O@u?<N{`
zSdI96IEx?s4vH<Wu5cZpZ|XO~?Q3Lv@Uun;0Abgjq5KX8SpC1FgkuUH+1L{=QetI%
z6ghZ{^(MB2T6r42a<&vHoOMapXpfviW82a+C!E+X*Lx0y>4cm=>#{@opUAd<ts1Zm
zXG>CdzIHwf2*dzo7L65NKr<=uni*mGC!i*4w(&CRiynbW%dmn{c>TE;eDpLFo$>gj
zfxhFLRa6LLq$&A+;((?mAF(=Js3d_jS2hKN>FG)rqD{2zOArCoTdH&i0%aRNNbPJW
zeQI=s11qQ>C?PXQ;mO46Zc1&aqFWw~ygX~Bp>|RrAEtFX$@7=jsay_6g=@A@O0|yK
zFGqok`f-$LtQD>XRph<_Rk04Cgycmi0IkP=IiAR=nLMER`PXqda6q}H768I3YGyWk
z(9^iGf9~XMayw%{v#%yl47j2ug8)S>poUXSw0&^W=%b~I+;~<xOEQ%f29ZwFP<<rB
z0^CWe?SsYyowLR#M>&OnqX^Fj+lcgt`f9iK*LoV=`!+To*72DSHFG$cYSaaQj=ve@
zzf;t_(s1!$Wj<3V<su^8_$@E<ZRsaD9H!~Av;sbae3oA;yg>TDHeJ6XPycm3|4&Yk
zZRB0PQ4r5VJjg3*+-23WJm&JO^w6W;Z)N5Lin{4t#7r9TYobUkj~u(%f?j+uA13n!
z_@z9imiQ4}qSTs0@=^mR<8<ed@kI;|N6L7ig-)jK??TawbJgZYCV5lur7;)IuFj0P
zH1C%eR5tgG9Yh+MD?}Iem7M-rg%5+9XLBo6Vux4QLMT4*gD`k-R8ATGV<|c#H)f75
zVd+p`J?q>`{UhixHtLtbF5>qZcEa*7GfILv!)i$Ck3AuuKc81pQGV7x=&#6a6SJP(
zey9}-s{K1r`zC8`rp&rN%AA&AkRlKuKQ%7mm0O10vtSi@nngvt3|!!{PLE${+I<?L
z>HOJFar`QowfqYfxipavu7o$0Irr$>lk0cq5olgcfL(pchbh{{_v9WeJMS1MPN*S!
z@tw<xJm-4nQZts8Vc6X$;>G9I%;St@o5f2$wk9`oH7g5jzu`9ur1%Vw4W+;I0E;6y
zMlT}r`^>vm?qH<JQof=0TfgyaNAB6wpPi2WZ-nkHK<EbBcd=pw6W)8>$^od->IF#Z
z$KXgIG%l1K4}A6TC)hj4GJ3gV&c-^!wSdjjsvX}Xv~n-kcG^j%Pn%>qgH#~S0AgRH
z0fYTGKa^;uKi{kuVr`_ujxFoNkG<zGEwU~5f2tg#tp#wmN&dQrE=d!{yV*dCY6gYi
zcS@(Fhh552GgH^rzGOq;iDi$D6T|+6`_0OweQlLlugtpSq03Xt14&)zdke)*RcM{x
zMHn8gFL^oTuE?<ToEKKo-)}}XmGES(%y*+Br!TsJQ?yA6<t_@4G}8S$o#8)jm4Z=T
zjjgc8&#MflUJ!jR3m{-F>gzW8`wDtT=ASw8RbOxMi{&lm_`7&p9hr6#JLq?)<=D-2
z4hhD_rcFov5c6#KJ=VRpbs3&V$ndAMwgc}inDw^#kIICVnrY_Pe)1)}Y*-(OHo;zv
z#NJY>oY?$^pKHiEe<}64MsM<95ZJJ2T<+T*SgH>|$*+@mBq24hUGGcQr`6RQRq=b3
zu<W_U0e59@RdO)6-<tN$Yni*vyc&qSvEyaC_@jf^L&rFWnUXI%V~d{VyIdEAX_v&T
zIOFFs%=b-V^1`m}iVlv5pzj#g=kVq+3(YdC9dITjX>3}$D`0h)xOd0CHUtGv+f>_j
zOe}e64sLGUE4mYWNXSX4+1k?vTim+i%_L!Qx2EZqKIxO#!+EV3aHr^Vn>bh%U8?_R
ztHZ*VAX<F{$5x6u6f9OlK2BGPR(rd4Y6fu?G&$r@3JU<cS@Lgv@S4_cQu;Jy16M7y
znr%#GEpxZNI77+U76i<u#lQ|OJ`)$`MB3JSInXqVha>q_&G1y*QPJblxgTWKLeSLu
zdm#+UoX6V0b#lSH!FN~yU)}$FQhd?528YYetswescT9*^5u$rtTE>?@Y<D|N7J4}2
zdZcwiGs)%aQMtOILuLBahglVP<5wLwXs~ZfpV-9F!^d7_=2muF8}8(Wl>BscD~JkP
zy~EKtIDA~i)_^E%ZQg-Iy6ETK(CW`rPRzS6^uI+FAJeJ=;OirI^DktOYaYLs$IM1Y
zF;oV8(O=;MgraN0yxy0fl<+WGZ4$LfAz~KjZ?*cfOS|z_&3EN_^N9^boHRq)pnhC@
zK(n%z|H$G_Xf<?bL7E#sdOTwaQ5E+-sC$=Mbev(h@5I2aENEusmBDCMY59c}GqH79
zE@gHWN-l7wNOn%!v=#0onebC+3P{B=DV!Ur<A&X!_r96J>yc`+$72~H@^qtK_8pNO
zw6NDGehzeqdudD{uliC<cWxDN-%N4f8;0eU2Q6wyZqRyOawTij?=U-?%UWx6;X2Qd
zcu-j_6$L_)tSC(<FL7Z^pL(JWx@e4qo~r>Qo!z%&L0_Ia#gfrfbxB6k4dqg4)s??+
z-o#%fKl4VA!mciAxTo1?_T}VP#L;Fs$RE{eO;6pqG!UPfrK3oudgpAyWCXDGc2izi
zk7m1AFSvtp5<j+j=;98(uzujd?4Z?$rSe(>`y!?3_D~^T0`#hvR3FnkT@)WGdz?-)
zz6z;=nXr|J9oS&joPupt5~u9MUh}pX;2t>^46^!F+_c+M)Ns*7!a(Aayqsf|RxkS@
zso_r}oFtg3tCSS-d!-Im+|_!tRU=O5(-=4`aE3jPx%ZCtR#CmTi>kE{0gI_5J>Wuz
zA8s*=E;*Z|HIA!xa)@~dOXEh~%RCo~dNN31hXiP+6ajp8h4K70x2*XZPoMyX{5ZU`
zJ9XoWsjwWg?X@fSvF-Pg(S0H-Ov}+LT~V!#+)^B5KP4YFMe|PHR_eg+d`{NYgD3mi
zyJ2nO#&oj^AHy>~G~cp;uY~NhlE*XHe*0L1YK82?p?2?Lq%?bcZS0RFdS<iuPWum^
z^3jbJ5O^AR_~tWz@jOpLo0^M^_!C=QDv8nqoy|~lQ@TmsrxmBov0dW1e~nL8>drtM
zF^Y)VIC(Lk!&q@yjM#dQi@3c(N(u7`?T{7`^Nyi!vAW#51I2ADAr)>i6ck0*8e$Xg
z?ZQ9ER>ss=)~p^p{ziCYF-><5v$a@7y~a~SbnK%K_8W<054Q>7i@K9Zc23SbN;=*|
z1nclWAd~%Qyj9gz-ci&E?dud3)ZZM3d0EdJl!_6pMS*Q$xo^q61HW1=VV?8k@<=YD
znOL4FX0<j3Zde(!8e5^0zg#tZO0)7?y*pikPOJ3Jd4FwS-^T4KQuhVZTlDb)Y;tFr
zBr`ZNUZAMWgr{n&LcdcrnW~E2^y;h5SD{SZ&zc_+#;rg~lRK8qb+@u}#8S_##R*6m
zQgCO48}(V0@&!TP!Z&8VH;iZ$j&@JuPLJn!jlS)z^BA`otlzTp(9K;DEmF_o7lLQ_
zY#J$a_1<Kj<5j}VHp86Ew1!-BC+w#hcQh87VOSncakq{YmiV4dNB&rkq+5H(w;A=+
z>OH@*4xp!ezbzkg&fnjG?pA;I(zbtj;7rmF5RU<(mF<@U!Gs)-o`!GnnoS%VKY8T(
z&#iWLX2n2Bd}zHPQUelGRhfWZNZ6p;9x3;J`G}Q1zThGA(ogRn{x`OAh@Ch|Bznsw
zj!r2&L_!j7VmvNE{NQhSaY8=yI~VG!2Pw^l*h%q5qt+)>v*sh@uHN1WLb>~OZ*_d<
z49i&MTTMljTX)q~l9+Y!1;_;SN<C2KQ=-^LiR^Rk@BUU?LPmQY{X>9_8NGt7xhjP@
z6&88vG^H9nEtZhiHZbaT3rfXG<E#t=TZ`b&=u-mcoUpr!@Ei_3djOzF#Uj#Ar_WB8
zUg#XbfSWQ#57hsmWA^*-;!2*MV2n(t3I`5R|JA>UO6zqoMa6eFbF^XNxDnIl^BTl9
zm(!y9EW^q8Cdf0lvtxqwH1s2OL9Ji?^Qg4G{0qsWS=nv8%MY-Sk$e<1w;m#QM{fA)
zT4anT<P7v-Nw4fEy3ID|1}KK~yadDNW~1Ddgwa>z(u`b<=b3Dkk0(1w@o3T4y#Bcr
zD$Y3^a{%PooXVDEP2;1e^g&^%dqk36S!maenYVtuzotPbV}MLBpQpqoWc~_|HX&KE
za@{C7p?0Tvhq6-5_qf@Xf$AJUj4}=8#RZ)q#Xr@8f|iWrm<ah=NmSmFN59mBo+-19
zS-XcjZt)C!)?;5kEo_Z_RGR0ust|2t!5pwO0~4tEswO=@X$=v5!lh>>H`lSgNm!}U
z*U-{%ad|ONg6`74l5v9;i#cFx=JoIJNXg>pizZWW%j^5#IRsJ*pd<es1Rwi4egzBv
zJA-~jEr1MI3mbH1o0cjZEU8qqz*71_qsuF^f_&Wa09oXg>aCOY_Fudy|FLE~F|Yg|
zpKyf_Q>l<W-$p!M=GuTiK}1WFC#=N}_@Fmil0KY?{Z3o!L?S(idTQ$h4phZ-l}Ip{
zAhqhyUq2wdiL7vFRD}Miuc<U<6Ga_-lI9|Bi)|sPD(|X!iBaF4rSA5|NljmdJ$iU>
ze>bi`$|gIPPiPlyIyfAb%H?zd*e>D^DxSoz_B;Wv=_9H0DDd@F`=sWYh9Ye|QvIyf
zrK#GNhJ1Yo4%>MUC&OIBZuqaVfKe?`uk?zNdD-uzFg<y4IWijgDbrScoRnv1NnMa8
z?sNA~PKgT<m4ilJ`46efm!!N#**x8}3+z`v4Us*NXyWa3jLj(%i#kvg5+^UJ+tiX;
zL||*`mm1Wy;2KLg(<s8jw)OKZF7dF=QFqR#YzNVI#cCjpZI}Db8Y!C!hLTMt=<L^L
z_x;+}M5j{KPTuG+jKfU<qS`5)$60etIVvVK*LFtVXs?w2cuUXC_gHUmy{8hIbZ=0^
zcDDtgD6(_UO-gEj-H#kL1i&8-Ki$gUSTlE4u5Efj@e6pHJG4T@PQ^y86Nmgc6H8PR
zBIf5?gi3l2zY&h>OjU+>{h?=X8ujr)L$00IVLCFUi|5Luz}sLkqtBW`r7nbU8*Y(X
z%mE$C0Wia=IC}%`ErHTu!j5B5Wr>P61Wq4(YNQ*?FDpj}3|>;xKlMKCHpbDxx2avV
zy*=ZXj?%7!gDHw%9i?p~pYVCZ$O?k#zXe|xenYCEqT7Y5T;?YClM6$Al^-*cyY;=w
z4AUESY8e$rkD363eN6>nG?UO#+nWbefX<7l9D|^q>pAR8eH{KgzVoKNClbwPX2047
zI&dq_OS#XIo_FHax$Q?!GEM$JA~J@oksSMvjeB;2N%D#HKz&ysPH(_s2HQ>uOMPH}
zVbRpls8|0o?s#awEuI)Bup*!sZa5%hLRXD8#k`Atc5dNWASHHIqba5zGUBgYqfE?c
zSLBncxbP4nHOyOGibYoVrnr55^jk+5z3FR;v~_?;SWtTU95;3ZmWg@&f^|GQ+^ttL
zBFuUG>UFFgG$-+5S6vVvMstp{j3^dsPLyJoHcT6Hy_YZPA~rcJN&Zg7w~}Hs(A;ea
zbv%THehw;z<-li@IM?z*XtCSLeupfUp=SbBzhx~lp}<b%V=$J>Ji>uh1r?d4ooa?{
zy=9n>f}8VYkvlyz`tTIrpAVhZKPPTo^|!ip$9^+G$J$GB{zZjT$Ep$Xh7D(1;u#)o
zE5NA|rI-U>+~Rfw5ADles678cM=#yeq;0r^)oU)}TbICofIcnH8^5Y;(U@&XVR=;v
zJ%_BpW?qWAg6daA#RxNN;#$W3gsy+S#L=gr1WBRIi<jgSQU?9~MC(E6ch5Mk|4-dG
zUp~THsX83n-5b&4mr|pn#d-X&D-R*uuVQ8gP~9!TsY{s*;{&F+rRPd79@gH4bU+g}
zJ>ynuq5K9MSyooeW^QHMOU<X4O0-aqfL9s!|CSraTnAC_OTQe>tj^~r`oJRrFt+i)
z&soUB&0EFCkQ_G(OJpX5g1UF#3dAM2wQ$xDtzwOAo!r?@2Q7lc$B)BYE-@UnXnT!%
zS94)3=53{2yBvjwr8oaQo;s3|v(QQ%(;37+S%03{^O;R-@vR%e{6@mYd3CZP^CdrG
z64n(@!=o&?=biL6<Y~3Zr+T6)TV!oMxv~eO%P7Dk^<^AJf~rtw>8{VoN{#um$SK%#
zExGb<&?y2`bl%WRiGliKIGjair)?&a(F{~(sw|PE?-8AA#KCbtjK<ST!v?t2sH76>
zqKT5k<&0uKMG(%k26rZk7tK3*plEtf4T_2T1Ly=$pC9-i_+w7_xQT3YbMI6(^3wf8
z0k&A;CWon-&NG<#>PH1<i(wo$0}4atO^jU@jb8s%sKg=+e8ooq)|3M9y)0rRpt<ie
zCYn$OLDLuUk3LPv^}2rNk4Xv@`9R8u$407wml$?JoK&0_=m2VtE`jFkTeV|j6{Kw#
zt2xZ0r<*GW-CbWU0V`8o%huBUjMkK&p)8YfR`8fR+ZkN;HDTlKFJ7BnUjLd?e*6#4
z<4Eec+EFNWx#I`w^dtDF&bjoe#JJC8To-PP9&r*{iVA~lk!c2R1O(*3p2gE~-k->;
z?{7!<vx;v4hKu#X=QZB-s8O*Ne1kkOE=@2eaVd~n!0t8nH)fkZ0gm?UVP;~7$r-DY
zbUWY=kdCg%11b62<LyQz6FczFv!0298$nmi_P53_>_!cHi!FVSm$7aOe<IR?RHD*i
zBBO~jJa!EkcLib1*li{lC@?4vu_EWM4MhUu?SmI}atF%vjGo>r*d$+^xg<k@+eB7$
zYg$*DqV%@C+|3?u%-`u3T%a^a_7-B`yx1c-Kj_?HnuVW?JXtlVd068~^SU<JK=k)e
zbq;@22b3OmfHMO}t(l3NuhAa9uNq#|dzaXACv%bTg`Bw(`q?{*&h$V}?c9PLK;$Rl
zjV#_)^05S8H(M0N6H4c2+InDS&TcEFQBp3;!xoyxyPq+qSLs&^SgyC;zgrulV!AkL
z)Xsh1r$noz?@JA=IjyaZ)u|__Q6^N^fO;i}`8@ML{KxT-tKM`AHru0DV<{$<ovW?_
z1XY`ZX(tMyhsan|VvX1ps}6}O`BpcIXjnOTKG$2;(*hFto%M!E%6vf5Scd5R?)W@#
zx$6N`SQ+kuol88*M!R6qmI83q{Mp^2K$g+6GwD~oA*z^eYp-XAt+wOgH&k~%gH!|9
zD`Apcw%%Sw;|6DLX77>X3$Pa^wV3CfJ40*J2Njvi-a|D%PSG=LC91x@)(A$Py2<sI
zNXr#rz%Fyv?y=8h0ZWxP{3(i>Js&e3p1%xKf}zr+qaHtC5TLY-2!18Ud!syCI?s(;
z5JX9K9`Es<m(>#{^)hQkq*MJ;)4BKgj=Avp0j8b4!RRE~T}8;OD*zJ|GZ{9lKQG(~
zHap)?&)sz0tS?&K@K>=Wq(xCPDZ)L=)y%7l|FCW&Zhj13I?);~{yhJR@$aF{2H>uV
znAnS?L+5;uuF*D1tQ`=lE?$d%)$pkL;6-SCc5&ZYr#93Q=0oa{BH=IwU<uYA1=db&
z>ckQmLHcct-9I~a_A+H@13_L}PwX5QBaF^SH=gn#q-6@xI^`vWFoF^3kJMO0>EFOV
zt+kRn&yds~YV=13R_drJgj6YN{wSW#)Y>sv*+A6iB(DqwUps+riQ7r_z*&~^Aw`Os
zzbTps@!EhAnYU?ic18YIbg%L7Z-9i+1(9CJ?ed&0GBS4h0CS@L?U8&#zXw)KNGxD4
zdKN&Y`M#DnTnDO14^n%h_yd#xV%jr~zwW*IKL#6rh-Ys5>F&PY_9RSt$FZV$GSu<r
zQ@!K8(Oo*nb`fk<+0bTgSGPLtD5@<c4lzEN!A+|V5bd)Ew*`!kkD|76=vGo`Z%zPC
zT0jX<?$p&2o8G;L(QGDYXT6(=JeQ;V3Q0{i(o%u@ux)SR6r549d#|b^M9|cFWr|z>
z)O6fmv**v*=>5{Bg}F4;)AR;7$)Y`^KKADb=R2&D&n=DLW};7AWyw99G7ZYz98P&P
zhcDhWdyA!nOqN~%$7?+XucgwyCi$!O{2=H`4#`i?F5*{?s*9n-$~eiWk1jOWtwVfe
zPf(p>g$3a1zw3kg2dW5yUi`Izfw@%NKVZe>L>g?LBA$Kz;Cz%2WTq0YO>$cHQW5sO
zy<LAd*m>$(5@4e@>tOsuiOsn;HGN%%ljKyAKT1(pZ(t8C9goSPz{j_v8YUX|*^ES1
z3V}QM1Ao5jFaG?B812`tbds3qTI%hqRpxOT0Qh)3c>N-RO--fqywP5lA<I!!_1)3I
z=x1G0x8UnDYRMVR02@`X98%`uV9HD$Pyx(p2AKdN9VizL(LTewd_C5&SHttP_8*&B
zC{11~D&H|h97&7}s>c@eyNKkFY;YR`uQ~K@b02?tHNPMCLnVdYVY8kjB{|^nL}tWR
zg})Wq!UYn-QdtUXv1tp{sbx1NjtaJ;B#GkYLr-bqVhuOug8~q=yE^`H1v@1XKOLF6
zIt<KFv^7qfr<tP&xXIiR4Fx+A#%fZCs+3l?NW_yL=popfp*7Hf`fa^uA;i_XoexGH
zg!`G*(3UVu*qrKrcTOWXW3hX=DKB?;Eoru$IQRMEnr_DSv5-c4{{USt)Tr55xgq-A
z`FozoGk5oUs0YMdh=W~hH)-YO1w-@1Tl7?|a>}>ni`Tcc44dVl4yHsp3l%4#Lzb>c
z+I>j&4Xr&2y$_I75V+!Ka7yA}SCp_QHrEV&vs%$h=vdReW(DeV8q$_%)9Bo45eBo|
zn}p?-*m^#x^wEbHNn3u^frrU=bGQlBu;8;OqC9M7r<&yLMF50q>m;2Cb=!~B{V0EY
zIM8F{?wOTWO!={vB3>|JV*2NMK6g})@tN&;^Nd_&9Tnr`y$#l;qK@C20&`kv_SG{=
zc;Y5}*U0rJ?>^jpy`e~su(c8$t!T`X5TbXNm$uKEX6>>olVaBuc9`z|s9z)OXHOb3
zCl`~FBM3ckdD_3*A?7o;c)sZvaWOgR_rlJ2bFF}rh?qVnd6^!Ml>U{g#4r8B4}<=y
z(o^sLSCyW~m4Bo3*l5-Sj+|&cu92L728fXUsaOE$(;sR*Tw=pwZB_oknJxa3Q~)hl
zyo<7f9QT-09>iA*%rV(uXf6+r+`>gSEqK(SIy$Ro5EI_Qs(EE%lGMSdObXxLRc5!r
z0sHXD9cCVK*a%_t@(sbb?We0fd879sNqi=FCBOr%b%q<Vlj?mJ#L5oVM=mmh6tno}
z(3SB|-?peD*o3;D&|$J9@OvMme^y6IYKR=t3AMzexk8fah}5)>p3_WcXpdd^3l4n^
zUyP8y+fif`U|oKR4g0w!P<bV3xfFk6r{7w}OItg#G}GEKm~JJmX%x&MY7g%!I0o-F
z$vV}BKyWK+NXZ)aOLFrQ(d!gON^}hd0t=V2^K~It_A){=-WbzWD&S9#Y9%c_``_1p
zTD$%L#QdTDGsm-xyc%0H-gN&S^J#z$9q+3uA6lkrFrd8enUYb(g%L%>*&=Ce4227q
zczeIGw~EDw28^e$4E@zqL4n4!P1-z7oMcRbi*@6z?C7@7Q2$SKcdiPe?H!bCEgi~s
zeb3nczs(sOJbOY{*+jrU?k)I^d&C-BzUhh5kE`5PQhdd8pYl(~59VA9;Ocq0hJ`Z3
z=lGvGL;DZdV=HoC64F*%lWvN=*E4e$_|Fu4qt0WlV2{!&0R@djC24RTh}t-eEgcGb
zvM5dHh+=$quT5;^-E!@y&iNi<+Ywt3E?Q);mSYj6zzD<p#a@<6TdtJKwoj<!e(C4@
z<{hwe1d8X-p<$0*iga?y;LIGN;JYe^Nxmp4xJ&(!=}tn3SaSEX(ao1V?R?mevRZYY
z7FIS}U@OfA9qp=P)gp{Y0S!m9)t~*%%`0Q_=yMHnY)h)+^^E7}dRiF#7$@fVV96OH
zj&1xUOGjInpQk&Y{80p<_Fmdprm_mSy$Rj<HflawsO~8g(V#qE@}c6$r6^Cyr6^8d
zm&_S1BK}%YME5yv8m~xbIhJi%me=;l!zBn?dqupG30pzYGrqspz$k`8u%$By^rS-o
z0V%=OLki6$X>JzPV&uFCbFpK5|36J9EN(3RWf@7fiz8yfv*Xazl)V!LZ`9aFsrtd5
z3jQ<+{w%Ts+Vb`;Cde+_%uh)2A^fbA4R0V1A3W0OuwNNmbE`Ybp8mSCN|FiC?$j{H
zm?-D}npar9|7%{s`G6(7T;|VfO8=ZJJ*hAnZBA-?5`WoWq@(nQDG?~T8YQ00GbI)$
z#POZ|urnUF9-RnGNr?-y)BYFBEVJ4_GE=D|0V6mBDwjCmt;=)Zk^&egSK<KD@O`lr
z<%36%7Qx+qX8|zUUI0^wx2q=@$=mH1)ZYeqaGxj~5&L=?#nVO{A3X1Ksu7<3T_~dN
z3rh+QvEYg8H9N@t-k6YkSfP);#<>#rd+i9(S>VpOth>}Brz;qvSS$Q}v%@j=g7sTw
zp!@T(&qWwYQcLgRv?9VTgIPftI)BrNu7W1C&1Bb}_m&q=d#fVZ_NIW)mZH+?u~kW`
zu)_O%omLj+dPU8DxOlyZ@?-+6Mhwb%Rr0oGK7w^d_y1Wi`b7llsfJf{MugdqrES`}
z?C}ye*WazrY2C{-S_%K147BO%c91voGj9;xde44c$QuWz@GC<#zdI@EH=2<D1(?Ge
z%-=;#g2HRdmiqmj-w8zKteLO#Vc6mk5op*_QNzUhU#NpcJ@S&&dbJ6_jCXQB4vBzM
zh+0Ixb^~xrElu3syLGY(S<-%uP|*H!gaTFl1D6HX@}7C{{ErgJf4VaLx0&Yuy&m4b
zhX7=RB=~q2Y#xA9%IrJT-E-V6X4^tYYd72C8gZO0Pl;hN|4smM?j@td{Vg56BLDNw
zJxsKfA!_1#vaP!)PM0_Yz5j1Cpls%Dpeaqmp_k1^^D4)lsDrCmlD{Vwcr2D<mMV3<
zHZnSDSxYRxaB~%7ZM}X07HR)M2+F(skL-X`7M)yLU(L?uyGzkuw01W}g{Rtz$QIki
zi%(NN;KK+onJ?#b%%RoWJL;|SJ;$YRwA##(nM=zWhf#5K5;#m5XEFVoV4AoI-s!pX
zOQoeRe~3TR*F6Ae^dIxmQs!RXgwEIfJKgQ}oKLe<_k_=#ev-m)5GldTyiDEDB;*Hx
zue*M3xdI8Q0yR{T`xm3*d<Y7Og2OssQ9{ZD)uc-KgPR(h?Fq3sW`hBr;N=eeQ>2ta
zm}r*fIxPk7VPEzj!1#6`3VZ#-kS<%RIid0|hz*7T!gL*6nHIU?QB@@=(@oLq1wx&4
zU-Lx@PQ;%otFt#>2LfdV?Nt$RA(1pO3|}V@?%T61b+XBx8^Ss{8G9W_0MQe5;R8R9
z>c(OV?+2U@s)5w@#$<qd(0EbDiC*+aW@|EnzdYOL%U%+eg8fK96mlNl*il^zeRcM|
z%YV1uDe(V`{cdsyfa`{hz=raC_~-ITzmra@E(ypnKT#e*Um@86@9}Dgy0D+cm8N)q
z?)R3DKE3?YW1i%H8^b0LtL$HsOG%U-Z^2D+gnU^PDh?@E<b8PA6F?K$S9*U0vm-xn
zw5pJj;lKG!l?b(9?s1`IZUah#*H^o&(vj@%L~!p-9t_cRB>VQA?zkOOC)iZuKShDG
z!5ZAsKUG};tOGRg#S6=ERYd43YD%8Oe3R${%PLf@wa1SPPwvWu8PZh8?!vqWecdo@
z@#~I`ALmI^6J$oKjt$ZGYUeRr)vHpqVngqR$C@L$$m8m{7&0fOQrgbVK7pfx_QA((
zyRpz!2C$yJw^R1tOnZ}0{(p?^#4#7;5-hzq%gw~eg+|5M`h(Zs=5Wi&bU%syuKl-Q
z&?If_UQv{^3my}m^xZeIA(Hv6Vc-v>$1(mVJ$?U*^s4GKHDL2iU9pmPF=wsA$HWBo
z%eyrVosvf@R}7h;7r2J$STSy^8u58WM=mOUoTiB2)dntHsgY!6=8536Ay%5*`K>52
z=*?$DG<{t+T<Y6y3XNGLJ%iVdNAQ+WaDmah>+-;DAGdqLgZtC3I)j&%GMjU6VGl+$
zTNb|$y&MS6e;ZXqOy142%k25_=}z$xPU+Kc62QNUoVpx>lbDY$+t!XB+oxqwNC7k0
zxtWV~SdHE&1(vG%Z_R!y9-AEzF6De05-*|t)zGWSxU=SB()bV^cjEB<nu(<?z{r;=
z)~8zLJ!=Sv3-+jg0)8qxp`PQ-7pQ~jW59r0?*C~@oS36z>)A#=sok7l#F^IxUcPRm
z=b}>VYdF6JR3vQzE}Z5b8bg~h`AR@I^J$(@^|`hhp(&>7NN&Fhv##Z0q$f}M@51E6
z$3&)qvqc&bnPg|)MA&4V`!AYS@oM3-g8v>TzxX3g{`qnB(p8=}i*YOZ&wo=CHQ12D
z)YWA46_!(0+6UA6^$l+k$r)y5OMa9|Z53N@zn|O(01%POhF!N5`kVU<ov+PlI-?=l
z<xaz2GaHoI6Yeg%O0hXdt)0EUTGW|K$7($%o1vxKhaLu~9<R-PJ{4m*Y0g8J`9AG%
z+m7s@{LX3d3c5}wAk<X$k5E%&?dm8Pe7F$)eWAzNHpR<r@8RBNrm#S<((7tq>GS~R
zJu!jKXijG%z_W-wPP*Sv{au>CXnp;%ftkkEhToMpCKF@)OG@yDiGtj36mX09(mzh3
z%U4QCR1|5OhV84<2g{7@4>7GG=p;lYKC*fU0Ah)_08R~AMy%H&g#nSK7tB9_^#h`I
z=iQ&%Uwy_M^Wkn3u#*Anq!P?*cvS@?<8i!M!RRDX9KIFFT4z$_%}R9DcHX#zkKnX3
zO9cXl+ST~Wj&JDFjHr2tr4&ad?bH&|+PDNeY-z{FVScywJckN2y)vooVk29giQe~k
z1MU7yE=%x6VQFE)3|mW1hi!vM;Vc}{VBi`$#YDXRb{;cIwzh}G(2Ju0Am?uLEd@qb
zI>tCX2a)REvvm9<$FPSOc>yC2A6;{Q_6^a=G$FTf|7T`~?$*n7Hb6Z0B*9)$>nP&U
z?>@TA_Zsip)5xQTT<mpW9|M~a7t;l<Un;rxlx>{Rg!_r+Ht;6?eS%&0Fb6C?J2i}2
z!$SQyb(Zz-<)uYFC-?$$aLeDt4Fi~rg~tn(bKIpR1@>&#-~XZO@KomS;NxXmL<3jH
zla{gLs^P(us^Jcl$64tGo6@n;VU>E6w0SuX@IM(#nab$-E1W<+%;QGb-_ladiobjX
ze^o4WjQ<UcY>*xzRa%LU(ETKtEa@yOLw|nGr6_|7KAJUr)<3^=33<KTrDM-gwEVNG
zdzkA2U(mE#=`ShwjPINe{-YVue{iEvSk$%30h!|BFjFB`1hSIQ!Ij$?<5AjQel;X|
zcW{**X`hWlI#o3A7ZlIZOg=Gd9Q%Cfm`S#q-R6TAxo5LL`IC&3kNF?R;H`*An%C|d
zbl65lqQ6Jp;eh1g)XS%t(-E#o;cK+R2?;tPgCwM!$Mv4sZuKr&ofYHw8J?w*o|b#l
z*AXvWo6o`1sc)VGL>UrN^1oyS1IiUkl~asVQzzi^(}$jyQBLhcXY*`@Jl(<_ggAV~
zwA2*Nnf|K7(If29{Y&4=Qyy|sU$3~^|M~V{29Q3wOU{ga{kZL1U1#>^7$<()=s^CS
z+)<j_bfwV|kK#aj+g~XGGbxa@-!4gg1;mQ35xTUKT7f?pvW^&3dvEaBQk~2xQ#0gZ
zSmSS{IVm8w<Oy_=5gpMgrbWgga|iG}+^nx02Eti#Zh1)`JnwR<oxLl|tuDlYaLS~S
z_j~+@5Ya)6g-%aJ(KIeIokb@73@L0Odi3RBwE6^)j$$qaAeZIsGy@fodK^fIr9t?2
zwn0|9>lJa@;bi)0ffOSAvQ8wB<BT84s$beBh6MqAtHfpIv&d-vqC7rYmh%PFyI&oC
zDbP!Pis9+0^JjQG7jnan4UHmJp<20YDF$hZZZDq!0^`XaH)vA+9RdQ^Q>kFJKt!()
zb$}aphy%@qJ=8r}GP8hac@mv}bm;qA-1l2O?oUpABk>Gb*DW26bBiy5f7Bv*BBv^w
J`{>2n{|~Qrp{W1>

literal 0
HcmV?d00001

diff --git a/tepechin_kirill_lab_6/src/MatrixDeterminant.java b/tepechin_kirill_lab_6/src/MatrixDeterminant.java
new file mode 100644
index 0000000..4a5b4ad
--- /dev/null
+++ b/tepechin_kirill_lab_6/src/MatrixDeterminant.java
@@ -0,0 +1,138 @@
+import java.math.BigDecimal;
+import java.text.DecimalFormat;
+import java.util.Arrays;
+import java.util.concurrent.ExecutorService;
+import java.util.concurrent.Executors;
+import java.util.concurrent.TimeUnit;
+import java.util.concurrent.atomic.AtomicLong;
+
+public class MatrixDeterminant {
+
+    public static double[][] generateMatrix(int n) {
+        double[][] matrix = new double[n][n];
+
+        for (int i = 0; i < n; i++) {
+            for (int j = 0; j < n; j++) {
+                matrix[i][j] = Math.round((Math.random() * 5));
+            }
+        }
+
+        return matrix;
+    }
+
+
+    private static BigDecimal findDeterminantGauss(double[][] matrix) {
+        int n = matrix.length;
+        BigDecimal det = BigDecimal.ONE;
+
+        for (int i = 0; i < n; i++) {
+            int maxRow = i;
+            for (int j = i + 1; j < n; j++) {
+                if (Math.abs(matrix[j][i]) > Math.abs(matrix[maxRow][i])) {
+                    maxRow = j;
+                }
+            }
+
+            if (maxRow != i) {
+                double[] temp = matrix[i];
+                matrix[i] = matrix[maxRow];
+                matrix[maxRow] = temp;
+
+                det = det.multiply(BigDecimal.valueOf(-1));
+            }
+
+            for (int j = i + 1; j < n; j++) {
+                double factor = matrix[j][i] / matrix[i][i];
+                for (int k = i; k < n; k++) {
+                    matrix[j][k] -= factor * matrix[i][k];
+                }
+            }
+        }
+
+        for (int i = 0; i < n; i++) {
+            det = det.multiply(BigDecimal.valueOf(matrix[i][i]));
+        }
+
+        return det;
+
+    }
+
+    private static BigDecimal findDeterminantGaussParallel(double[][] matrix, int threadsCount) {
+        int n = matrix.length;
+        final BigDecimal[] det = {BigDecimal.ONE};
+
+        ExecutorService executor = Executors.newFixedThreadPool(threadsCount);
+
+        for (int i = 0; i < n; i++) {
+            final int rowIdx = i;
+
+            int maxRow = rowIdx;
+            for (int j = rowIdx + 1; j < n; j++) {
+                if (Math.abs(matrix[j][rowIdx]) > Math.abs(matrix[maxRow][rowIdx])) {
+                    maxRow = j;
+                }
+            }
+
+            if (maxRow != rowIdx) {
+                double[] temp = matrix[rowIdx];
+                matrix[rowIdx] = matrix[maxRow];
+                matrix[maxRow] = temp;
+                det[0] = det[0].multiply(BigDecimal.valueOf(-1));
+            }
+            executor.execute(() -> {
+                for (int j = rowIdx + 1; j < n; j++) {
+                    double factor = matrix[j][rowIdx] / matrix[rowIdx][rowIdx];
+                    for (int k = rowIdx; k < n; k++) {
+                        matrix[j][k] -= factor * matrix[rowIdx][k];
+                    }
+                }
+            });
+            det[0] = det[0].multiply(BigDecimal.valueOf(matrix[rowIdx][rowIdx]));
+        }
+
+        executor.shutdown();
+
+        try {
+            executor.awaitTermination(1, TimeUnit.DAYS);
+        } catch (InterruptedException e) {
+            e.printStackTrace();
+        }
+
+        return det[0];
+
+    }
+    public static void main(String[] args) {
+        run(100);
+        run(300);
+        run(500);
+        run(1000);
+        run(3000);
+    }
+
+    public static void run(int n) {
+        System.out.println("N = " + n);
+        double[][] matrix = generateMatrix(n);
+        double[][] matrixClone = Arrays.copyOf(matrix, n);
+
+        long time = System.currentTimeMillis();
+        BigDecimal determinantGauss = findDeterminantGauss(matrix);
+        System.out.println("Time gauss sync: " + (System.currentTimeMillis() - time));
+
+        time = System.currentTimeMillis();
+        BigDecimal determinantGaussAsync = findDeterminantGaussParallel(matrixClone, Runtime.getRuntime().availableProcessors());
+        System.out.println("Time gauss async: " + (System.currentTimeMillis() - time));
+
+
+        if (n < 101) {
+            System.out.println("Determinant (gauss sync): " + new DecimalFormat("#0.##").format(determinantGauss));
+            System.out.println("Determinant (gauss async): " + new DecimalFormat("#0.##").format(determinantGaussAsync));
+        }
+
+        System.out.println();
+
+    }
+
+
+}
+
+