From 7b552ce0bad156fc4cd56bae44f8b76bc8830bab Mon Sep 17 00:00:00 2001 From: pasquadambra Date: Fri, 2 Apr 2021 10:57:02 +0200 Subject: [PATCH] update --- docs/amg4psblas_1.0-guide.pdf | Bin 1806938 -> 1807728 bytes docs/html/index.html | 4 +- docs/html/userhtml.css | 24 +- docs/html/userhtml.html | 4 +- docs/html/userhtmlli1.html | 4 +- docs/html/userhtmlli2.html | 4 +- docs/html/userhtmlli3.html | 4 +- docs/html/userhtmlli4.html | 37 +- docs/html/userhtmlli5.html | 4 +- docs/html/userhtmlse1.html | 4 +- docs/html/userhtmlse2.html | 4 +- docs/html/userhtmlse3.html | 4 +- docs/html/userhtmlse4.html | 4 +- docs/html/userhtmlse5.html | 4 +- docs/html/userhtmlse6.html | 38 +- docs/html/userhtmlse7.html | 4 +- docs/html/userhtmlse8.html | 140 +++--- docs/html/userhtmlse9.html | 4 +- docs/html/userhtmlsu1.html | 4 +- docs/html/userhtmlsu10.html | 66 +-- docs/html/userhtmlsu11.html | 68 +-- docs/html/userhtmlsu12.html | 66 +-- docs/html/userhtmlsu13.html | 26 +- docs/html/userhtmlsu14.html | 50 +-- docs/html/userhtmlsu15.html | 118 ++--- docs/html/userhtmlsu2.html | 4 +- docs/html/userhtmlsu3.html | 22 +- docs/html/userhtmlsu4.html | 4 +- docs/html/userhtmlsu5.html | 4 +- docs/html/userhtmlsu6.html | 181 ++++---- docs/html/userhtmlsu7.html | 4 +- docs/html/userhtmlsu8.html | 811 +++++++++++++++++----------------- docs/html/userhtmlsu9.html | 42 +- docs/src/userguide.pdf | 1 + docs/src/userinterface.tex | 13 +- 35 files changed, 901 insertions(+), 874 deletions(-) create mode 120000 docs/src/userguide.pdf diff --git a/docs/amg4psblas_1.0-guide.pdf b/docs/amg4psblas_1.0-guide.pdf index b34c817afcda937b8e54873d3978619de8aadc1d..ea0518e3a31f0c855fd118b141924c7b04a00b7e 100644 GIT binary patch delta 38569 zcmZ^K2Rzl^|G#S|o9nvfwIbuXdsrz^DZ3$&os>O_hHF-eGSXWrQbwsXD3wwv4Uy3h zQD{h#6;1qK?|b$6Hh%xlj>;Ng5cRJx zMEWdigF&`U%4#6c6hY3tmKCWg9!JN?)3O94{+p~^mCzSiCEjz@Ja}oQ2}9a5Vof_Br*cdqtIDOAcMeNQnj9gV1MqctFj!GaALVY%fnVE?h zFY^xxh+0ozlSoV?;hvl`qVie(MOAH;{5o}hG7XvACa)t;Cg~Z@HKNcMdNzRp(GhdW zEIohZ<|#QWlAkCqs=%TW$s{I$Ok)ulWHw=WFky|}ym^Rsusn`Kromn$0+~Q&(~0!I zXGN|L-5iOoOxKHA@6T&Q`bR~F@pcPc$LsRu{e*`28wL4q@(-Hn_45mi3fwgFnv3uH z&E64PcrS+WWb_S1h(6d^i0%;_79lK<6GrTku(PXlp3CpWV8{qrQ9-WiUc3B9Jce0S zIVL|>14FHPDUa1xLwc0t<*Jma*yvdpa#f}owoMmlxuGa0$YN2LRpGTb9Sl!^?ef+tOVADMj-nh${LG?h6KgBljvk**k4Y8pFu?o0m6h&5=%xb&&!KT(3nIv>%SMOHe#!U7+@$i5*dWO zEJUXh!OEmy>_uT3nMh|-kshL)0)h#_)(cZuL=ugLw1mJQolGNw&}pIAG7gzYA~R;p z%BB)wZ4wEDqVr^p-0bhqwI5yYm;L^;ZCGevNR*MSw{NIVAkT=Ap+TGcBWCRDv%%jt zia6JZN}}m43E_4A_im(jFwfKmdYiptd4t|yXQtka2ZZ`Wc?X90&n4^m{r4&&(~iX= zIG2pSQI3eTqMEe_obei|{CN#zC(i=qcU0Ybq$k%7TsL4m%3 zQL%H$B=)ai2!x0z-{>eb7+?~5E{HMS1M18qmR{(lKtF%LP z@`uz^8k0(7(U|{}8rAIYMu4!UlG%FeV|^k5VQm;jNJ2OT8+6_x@94-#BNu8XI&_gAj8U6JpayA^RBFrEWDQp(f5{|9?)dsgUMdF+wG`k|Ovwrn{ zjKJy%(-=f1ja8)-iCuybrb1?6v5>$hY`GAbO@s(ij>g^)p^+Fw5)0B=G`3EN4KYYV zPH)0C2s2qkHk*vZZN@$qX3!zZD2RIu_8K@Zk;Nh-YO&Z1VJeGArIL}#SnLNh20823`Gf5o0-Z#`c0_H{im#<;)zfgxZHwgT>-cE(Lyt!`t z-}@B3An(}FXl^G}FOWNpGIJn5e|Hjz62glFs@?%799{xl!hf#Ea$5qVCS*DetIV0% z3}S?VrJ>dh33O@X1&TqUv4utLwojb?~ZyA^hhF5G2c@^2Yz zX60YOLxDwqdA=YTG)T*Hc`V{#FE7EFF%L!0Cn(f+<8K*|qUY<)&2%ut)i6cRJ1h*0 z6BNijA`O(nwbQTj|8gSKb}$zWNVM^r@z4&kUa&V9)7L+eE2ILNt9Ym5Dj9qNeG3s! zz^V$f*hC6=ym12dA~NV8FOCc+V0RlpMrZ%3-P|NkW>AP^28)+Vc>B`1iVsAj+oAq( zlPsA=B0|9<%uGg;ELo94f(|l~%3wk)ts~Hg3@t z0Qh64UkMEahV8(Y*?Rtd(Y~mz(Cty((5XmyB39j+st3yT4nli=^EQ|TSO*FVYFrqA zf=r6^pLyTeKiE5BBX59$mTDrJ?iG;3PT<-$XwIfGNIbWu(%FAlpEtx+1UmG~@~Gi? zLp&Em4M1hF&^6EyLZ|*G4!X8)Xv9o_hJ|voE4m+?ib@w?>v@(%5?c6=+yL zJrQVh5{p0u*P=1#G-$vO3;c>orxIvvZV!z>CzBY^1C^WtZ&2xgOlY77DxFE7)46k4 z1Ui#M76U5)f}+8KkfG^pHk|-Pl0l-ep@U3dkkLJvBmxKEEG)>R5ExW4EXZI$1KyxA z*aQZR3@bCChfYRk(+CU(jUh&5vIq<&8!8i7e=fpla7lg6Q;LuiW)D?=|6_5|5ja6)KM2`pGpj0$niVxy-4J+Mh|ZWQRS zDHIw9T1)~PJua2aBCzQc&^qXb&46=*3fW957zX5LGg%-h8+ur92pSoB*fg3LjSPVg z^%32k1}zGRNF&jqMWdoCF`z|**3%ndkI~bM3Cg_o2CIi+30zsz0320CN z$>4}E4Aqtl*@h{mL!gonr|)uF2`p1O+Msr8>@H#?Jncq0teL z4#6B`Wd3n1QI5_VnI94o8Wl-c0|O+)CJ$?e@jx~n#X47g%EP7z$)G9l7gbtFSVLtX zUfEb}#H0`_1Dh->#A*YP2rtBXU`R7tGRgmL$>wgU&fC($0+kt@>d%##BxEENtA#vw zRhAG&w}K9jXE+Tctq3cLq&>mS5e7OyWdcHfhSL!M*ZZ|0^Pdfwbnb?l|C|xr6!byx z5StN#hGZVcX2^ll{v`|p@y^C-B2mY%azH9K6v3HYi{Q)*uaFJpj@{A~;<0;n-G{t*;B<#$0@NTC4U7MWXuwg1N?EX3&;=AtM*Cwut?RoC;DdjaNscG5FcQ%`FGzmLe}sU>%Um6EOa26eo&oJBghyi@L|} z9tAiSRka+1>^}qX(sEk-AQ?E=8UYAlgw%l3z(U&N%J`d8vw0qd`$u|+c9@HNFT;{W ze_0)w-zz7FG*@ET$j(z(dl0s#07N9524SC`!Y)Owp2pfE^9SW@5$7{lM?~}tb{X>N z4AvSsb_V8cFNb;5a+rrZi*=Cw%a2hh!H+eOyt7~r5^4`k#P2M)5ea#H7ON=$X#zbD zY8@iP&aE}d3J2P8;&R-Zfog3xY*~rcFSUXiF1He60 zfSI|#07#G-f17NKthxZ|Fs;DOQUJbzA|jM!ksvndY% z!vV6wAfV)c4Gu&DIuFb~{3T8sIai5wRYf~!6gmM)0N^>OofIMqG7k_#0162CFbKN0 z8z&*hb^G7E5>hzom0EIKQuEgUcqIYRx`dVEr?RV9m#`AExtAkVuD7s*QW#d1O$+v( z7O#sij+NGeD-^hDIjjUjk0=TqD!id;2L(4OI6%*CxN4c3q>xC4Ib= zpzzH&QwCclsf?>9!~IPAGTfLX{5Gz^*{Y%|$yXV~;&!Vc7#+NP)dK{VsgF(%t;G!p zBFBbsvdGJyiV7Sy{lBX&GKB(Bf&PPLLEt?8-oa3+&@iWw0U$uNf{G4R2*3vmcuF=3?Y&yUT6vUFa9cTlP8;1q}gG@#X7k~^h z;B6Gp0H&g#Z$l|V0R#<($x!N{fg)&B5b!E83oU$LL%yqyP~p?7)&u5 zi-UHfuR*1uap4?_9O&q4s4NWhH53xC06NeB6@eCD2nq)51w(+G;O-C|f{5UT92JJB zWR@5Wzy+PV5fz3QpdGlWK|>9F3mPmGe6Rq?&{@z=gT&7up%!ETC<7Iv!l2O&GJsDo z(Bpxefaxfx)PUFM=t2OH7)(@punmI+AwfkA&jRhx&`iRH1yP+aNpKdp3j<6CuhIZI zF+qFK0QVU^B2h;5)b{ReOk2z;RIm!Ps@R z8^{Z{iwwB(qyxuW`VgnC^UomUUXYj&5*-vYQ{UjST;>DL0@7HWy;HG zD9`|4{#7_3k`aw|oRQYQH3k%i)|g+5QmDxMNjcjCjkrZJD0=ukCqlJv#%aSgxA1rp za|saMaY79bHs8)PXB-CLLSXCe@Xy_jEmQF zn{aX)liy(^uk2$O#m0XYy= z1eyYQzdZbfCP4K0mtE+HQyIUNX6z7W6w&PYIo7-{bh@+C66yo23bNI(rKzg_ey7~qu z!w(Ort09+HD@!0A-{o*9@smUxU*J~$G8QDHzm>$~8Vhcncq9K)5|{gcJdmuoazPH~ z{c@4>e}0C9$;c}q1abUkM6m;B0mhU!FZ}OM`w-A27gGE{OO`xtd7f1P-hXCF1Xc0KkOQL7e-~L7*h^1fK~QC_U)2 z2ZDKtb3yF7L00~55bj+!&Jj7>jpI3Fd=Ji17IoJ@Bm(}3I^;`rytDumU7o_*UIL#I zPsQpWgi<)})?VD=UkZl|{EwVGh5suj)CiuOe+DwJ92t2$Ps6kx2wkH$IM;ua$mS^# zSq#QUSu#{~p7MZ3LE?#pctOKe9ZVGpOgth!oz40w9n(^q;?>_yBsi@PPL_qYl6Z z=(7hhyC3&Am4QS2U=#4W+}Y5Jg-vx&gx-!SPD;6`B#2`BB-_j z3n6|(U`h{!Fos)*#CPJB^8*m!lHa&th-ca`loH$s9He0w zkx=FLjpF7bIitUh6gl>98Urkb+`{0skS}9687MiZG0+nve8O!+lqR4!$34Xfb8DqM z9kAQq)BVvUyEN8DmP7a3ERlNF!?ZsC4~ zclm;wz@uN?u~lnj6oXWLQNfo+iu;9t?L|@*D^wBZP&v6Owv%Ey5I`PP?C?~4E`niJ zRRk!Oz^z}ETa4mjprE`(FYQzeQT+A5I;yCy2qlZlFgIRSTqgp{w%$;*RE5{%x)iUd z!(ii2MLiDXe@u~t!qdN)BG3*v8*p;qY7Dr4C7>6RbV$DN9RjEbmt4a$@Ejljbby>F zOM-i76yL)y4X`&{yrMUxOfjG&0GdgF&;bBL7NGuZ&c|3`p9(j z;u__POrQsJlqRx(^+B>Dv!Mf02>h7^TMS(aKDq$f2AoGB0StiaZ=j99#{ftJN#wo^ z+>pybfnQKiLW16qQ|NRKKzNib0tbQyFi9rRA3AVY3Y!Kb0KNVOB#%OOV2dCM9gq=~ zs{Be2L~ zD2an44QfESAsJ;ns4YMbK*d1>3c5174HeABfqNO$Zm={|UoaT_2el;1Yv`b9Is+w2 zC=q3m>0jO4g4A4Y31pZK24JFOl>u)s(a##NK715FCqgn{qFjyvypf4|7fSb; zKwIHt=w|}xW`H`GOfU|5>BMBAOdr(AWP!dZaC-QH0Ja5nGTE>fgA7!WMFQnhfbOwK zC}9S6$Rfki=qD%^8I%v#WAJ$ay>0-VuuyM8zd5jg*NZX0Em$;=fdLs7d=veeg0fZc zKa{b;g#!c(Y{+7xJ_;m}#fJMf^g@P>IxP5@3bf#?6tF!TwK9BHh6a-gm*$W^0PN3P zriNnhma>19Z{U<@hNVHh@>Y zdF{nQa0CjnT?jw(bpgFBg;EWFfWreEkTU)*4m7xkhEg||n|Q^H|A#jqe{j#h<&4hU z?7=+)azg~~jO^`%N?b1psori-&VqMMy+Z_Kr=g>GQP44(?tVuUAO+x9O(b0uFT)QH zM@8{e^!Mq^uh=XQn|o)?6`P+aq2UiVP@xt75gM`-DkPGDLJO`!L&ZRp zzi=vKHRP5Uo(`WKc+ddW@Xxq8{J^$&d%}&)-+MB7=l$0+xGWDmf=dJc=haN! zx0YXDGz_2?JAH(x)0X)^e3;48Mybqy+CXfi!JD=fVjTq_ZT}MMUsb?2D%AFBJXLVv zIGWx6JTug5I5YBE7M5qs7#Utfjf~uthkdx;bKsBR{y8;}DinZgLTctJMw%9y{VzvENkzACYD*2*i-XI6MoPwMJPMiC2*mK?ZU7y&&S> z2VrwB3^e7@>i4JcFjA)o!dEDAh5v!J`+03L9@^@7ZhrvUUcqZ0fi~}e#t1P7I~#F4 zhb0I=iR3vE1E9rU>W1>iRX19IxkLZ6GE@}wXy})2^cO9rykG9j|E$mZ7REiIgs8kG z@VN1LG7gv5nW}RG&p{nF5m5(~ zQjipN$cPlA8oK_jkLbQ5ykT7eca{eMUX0|un9b|T6p)ugIth50ztb)H_B1+Qg{myg z4-a$=co>A4u-!%t?s|p1pEelbeHiDOk}UvmVCKY7^SP$_AI$`8#62se#(#8`M zbh!B9Pu0Rg1d1|&M^P9QR+WGbew7S{R<(wKk2S_Hs-F7e`_$lT*7an(fjWj%)qfoC zs17v3#~*wjD51hg^1mk5yW2Xy4$BM{#^J-Jo^<tZg>LuvN9(7K6~4U(R7*^gYlS6kbcHca$>V!slYx9_CH)xD2Z zpFG^vugBTHnYm-I%jtcj!ISf!4{}1fQZ^rM+5KR^>Ul{-HtT$wF6;Hys_td34vi#X z8}XjB4tf8J%KJRKm!3MDaqRkVo;+jJKV0IorH9|)QzKQ!-d15x_I2uC=A%$d9%K z*UA`k4C@88wK*xCOI<`%$B%ZnkuR5h`{A(6UhTR|;F{9Hw~uUmqQ*a`zRvE@zpt#) zj>-9wCq^2-{rblvm$!rmL6QpBYnF`ZuPWMPFTOvX)-JX0@`0yV*0OE^;gxkxHh#91 zqN;5f>vmT48kiK{d~#5{PEPW{#7mZ0z%uQG;PUuS{v-u$FHUlT%e^H-67N1H?M&;x zec<3n1rftlExQcI6ULlRTVV>smaQ<`x-H}E`i{DbMitlhZ%%z9bttN#=jkKW2^w3v z{YRyyyYksDO%`|f$9MZl*8ecNXcpVDQPd-uZ~6hnQocrCbxp(R*O*UDI$PBa742G~ zd+?2X+)kyRhZ z1wJe#?qpO4k?hQTS|2`VY^rS0I?JC=ZTnH?<8C8R&&q3ig>5=jUs_+Q=s)hYuRU?T z>CrvIyS(t?-zM8Cv#i@!38>e-yX5ie{mVNtS({6)n4i_UKkoY@U80}k(X;Kj%ABi| z;Uf*xf-cIeQ>IZJCMB&m^sJ34cID&_)PAvJIqg=K-LCrZi;peSM0{3s*!A^>E{bQe z*{4z!9GmhsCDapocK-MzaO;>{XMhCp{I&oWlUddgh7LR875DC}#jt%eo@*~VG_b?{ zJ*8jOQcEp2&zVznXTp-Z((irANbZai1EXj;`T_f7JJ z2YU{m?_DCh^2HNGPxq{+9i@zz%qlY6AD``W#*Ode|OsIvh?a5W2N-L*UxV2 zcD|{5-4V7gscW|Ex{4HL884E04eNNf2ZLO(rVa|Q8u~(e zL9%s!MYvM1OZ!l7#v}i3OoQjOeVNy$I7-+L>rG@FT4@J^Z64d7!<1@$S{-StIKYu@ z+_wD1$**_9^Fup4+`b489~V0}WR>7_@9}}eN7}j7g8_R37^Jv3$+};#`iT19tat#8uxn(nr>c1vhpsNUWkMb|kkM zoG^~+u~@R8QvAkJ30&2~r{h}NRX*2D0%g(xQ{N?19 zldI?{%>msHzOCknzHk^w)K68ETg4JCJp90Vlb+<>r4`yEI1LT^$13S$B~|jlq0LRF zH&8jJUSw@Ew@ac=>bE{qL*G-IQ7x$7(Cy^MPgEw4?`*&?FuXZpH+OnOqcWlMB=ZnG? zWe#LCHydm*-K%e6eJF*Z!k2yQ_>#<*Z>&a^&v!nz^BLz3!b7WtZLcE@UYy_lw|^?b+M*ErtO zbIzPN8__Ju4&taERXML#rqILL8L{Wv6_tT`w5vu@S1O95xDY7*~Gyop^Q7BXqpcrR`5&ky#d zB$Z1l7bkPxM`|1ul^|xu-F~y%x`HLl7!EZbeHg4PRAf;YU(Yf5N+wzvda#ET5AK^Y z5NXyHE3{}@>%Ld*GI`Fy^Ys^FH4aNBRQgi(U2av~Z%#Gz=-P`^S1QbNvpTvXZ;9kY zMTchUVg9ty24=bMnw3#k5`OWn&%UtH{4cPv=Tk50Z z_}g!)<14qGZw$Ax!cUz^X*`FiUdcIhafS3qrrX+4^ACovYMxx%_He}shqz1AO15be z$x|Dey`^e|26_o{k8i*B*PfN2anD3i&~g>Tr_x99Gmne=sT;&zxr;AzEcJz zW1-9N#@%%7@Py*VPLs1VPTa8PmoHQKMfK<2x2&3|d3^31;ofLSf>-9nCmBC3Ja0Q5 zcyQO2CoM)DZ`{kKpRN&D&&N5r$E7mLy-R<~wl~Xdg7~Z7pYQwNdc|A*-qN6)w>pM3 zXHJ9dkPawk1_>-*o))Fx*1K*b$k6pqyxZE1Ph8cF+ zu<+Z4{BW72C)V+M-8y{Y>d-JndA$A54_B#$2@kERn0l8g8SBLd-(HwV-dQT!aPw1c z`OT-vhbD|q8xz&%4t0lrSWxru=?C?sgI{ZZ-qK1pkLw8j$!?mHby@Lv>e1t(?}heq zlm;F=uDD(6H@(3(dH&u86{l+g>s>xSK}3%4l?m+W8@Mo({=B*P#nHoz^Bz9}tQ=)z z*N8~J+@Jc=HSl=jNqhk2k(oioqeA|7=cf-SB8fMH&#v4#_ukMkg1Tsu`?Jif1w*;~ z=1aPaj0)Cnua+A;k11B^TwfAxw{!O=>Fzi69Gjxiv3%=8mtEF2dG1-lZ#@}%)HZjG z&@pS}8&`~%Y;C8i8dpXyQCj61XQDn5gGxm!Q^X?_J@)~AJPdD4ad`F?-rp_QU)u-+gg)$_@29I zp%w4SX>VQ-A24e;=1A96TgLY@*|+$7`n#vndjk^Ghua4gupYrmM=`TZcAdYrb06Qv zSC^CX#)DiabI#z>o}L&wainB}Zd_*S9`~R+)isM2DR-F2Zg}pHT-|8?tj%5eM^f9l zfg_I;v*mBxE!4}huC*PZ-n!=+sdRktDT5ZdC$-~uI3HgJ=~(>O?3is)(Hb~A@6fI| zVtS0xj#`z)GMBf8UuNxGV7B*W&QSh>eL23?REGuMqm}nZd^31u(3?PN4Cs5>IJkD_ zVc?kab1GB)dOsJ-Yt0>7Q@C4d|Dudc$yjcAA5PpqogTfOS35 zTdqH=wPs=0A=wVL!Cthn{u1`Cp10n*2JP z;1T0;cS1IXrX;PX7gC=4NjjL;C}EGiCWc8L%xduO;EQFhh^Uj}ykGFy;KySpb=#Nc z`!kxiJ=cjp+vqi#<205xHSc?i4!n_dYL1Wd&yH}mcEBWb4%sXW9E%3*= z*A-zUcs};^{iCXlLesk&+h)UoBQ^ADO@m5AV8Do05?^-&xm!@F{ck zia9AU?VByc4wMbcN;w&G=ES`1Ci_@gVM}NfLAfl=V|ga?aa|##^YfMYoJtpBS^H*1 zw3rn4_gKdyn*>IbCd`${a4lM35!LAGbw$WjJp7C4$nI6lmg>e>jnqQp2bYMM`!ckR z(v(%4Ln5gCCQoQ_vvUWPCf*!d(=eXokWgE2XX)Wpdv|A9I+4{jp5`zEt3q7xrI>;x zsw#WBS!av_oYg<#lv*X!3A^>Pp9gVvKOb>05mbD0a-)tz=!NcBCEFL$`zklN%pTv1apt?3 zOB1`jHg}eo?WZBllFXV9cTb5ev`kzrW!C6^MpI(vy;JXhYzt6zy>x#)w)RrU!JnUz z&pvOQGp+{qH|`QOCvv{(oUV)~ z-VbWdE)Ntryk^^jLftzzhVm?;8&?Kco!@-pNBfm!hVe@T_9~xOvR81FFE2fn-DNg{ z$VoA7&e>4kKXRn;r?J{{VtI!nS)mAH8J_y`jm4{2b-a{ZX~e14u`@3(9|;^Ud)r!i zygK%YL_*h=T>VCgjjNZYl^7T-Yxr>D9w+nc_;RLeWU$w^69X%WPQ!?4@istbu(W6j=_@GN~P8b;YOsa^>?TT>|wQIc&TsWCisf-n4C*$RI7Z4Tt2em(ls)Jvs3S<`-geY-a3jc@7RcC=z& z%@}{qDz(%f7gsAx2A<#4_ao`4=z==TKKqxw(;7~zx8+N1Yd`SbZmpD$w6K&PyT74C zwsXs)EiJ~2_UW`8V9Z4pEto%^jr(ltH5l1Z#nCrby49Iwude&yc#BfG%BYaTP3qV!-&6-lOjCcE4yyr{IM-vfvN93>osot z+`cijK0)A19d$)<-_3jt_0^*0OKkUC2zhitNm}IDu$N>&=2f+u+k)Ddrtb!uGOx#^ z@0)}7;Z(B1o@S*h@GG8bnKJkp&arwEpc8TB?0Uz3w?bXa-BfwFZ!~ozHf)Q(+8(Kk zmv_jhEs>k>&}%JtN;}$Qwy~r?dl_4X>_ZspK4|Ssv0goYrkn$XK1r>O$zTlP@u@%&Y*Vzq3QNa``jvDV!y^}SBz|M1(F)N8du zYUQOwxrJ0?WBd)YYkCD*YOx?>?R@zp*pfOT;1TOnKL4{JZuMPKwWGh5bJjwXS&UK0xwq zzn{5r)@kuc+?}K6S|^@G_+ol$sVkc#1}5F)8+`gM^maTd>8~I4T6xYqZ)xR?=}kwn z=uMGs!>*AkcS3>|Ka2M|7aT3k#;9Idw4mw4Fx6M-vvK*lS=ey3AYHzSeQlq+OEcCx4{z3b>Tj}kY| zY1vz9|M~HU8!H#?j%sp!&ABi2=&d1cR>Q6A&Hk(FJX*&hw?AHUf7`sATLxp~gG-I- zWIC^}egD()_S`M@qm@0D?h4n8-1#Ma)2bG4DBZeP_Iz*K2crYECF!0jdwSox%B!15 zC!K%1vPp3H{AYs`1#w@ihQv?B?2a(~IpFUxEw^TVhK6Cp29YtZJSmUF$Dy1>cb96u z3psrQ_qnGaX@i8$c($31NGf+c;56A#w;VzZ)}7T=t!=zRCAiOmy?UUY8i>zZde)!W|OetAIk zRYABzPu1>gG+kHquaZ`;`;U)l?sI#^p=m`hlp^12^+wPvwBKck*RN*YxHiWDr*%Bp zQayTKT6$~xlUAo|n3ewb$2^;DY<()@-`@>?wzryp_m2M4H@qXYu7odg$iGwnvowFm zX#Hk$CBx=0RwBRtI**S{;g^e+Jru4=*L_x5;p>|_RDDNhx@Mik(x7jb0>70;=|lSJyjUIO}VhGO1L0@{69` z$qs~`KII0l&hx%!OZ&BH(O+D(HmBWjaX8XY>>B2+GaYZ||8mDekNmKwlZkmL*M6>w z_qotzyW#B3+~`4%%fs}hF6yTVLjy}@-jR)^i__WLyy8xp7ke2r2YgvG^yX^RTA622 z%n(~nUP0VNH)p)X5Wl$jTq)`Nw4|=6ADcGcta<0!v~{n?=vCaYde+6>9hv?=+LeCGRK(VXV5jW!ef(b{xQC&1oQ88bM8EEtT;?f8{*$|pMRZa z-Ba_I&Yc(SPUPN|egAblcjf0~$DH-zORRHk2LyIi?6-bJ5;v9!_;(m>&~T{2 zb}je0O5(ia$QzAUmA_PZSdc!@y7G0OJLW`qea?vR{gF>q*al z=iwiR_d24w_6X+5Nbj}OC#PvHB=MuGbF570+J<7{Sn`iWq-P{#fKJRIcN%UnM!y{MU&71FL&0es8 zeOGa$T|XpQQ~DwPg;0}op2xYM*dA?OYZ}BD5?xVI~>h|S!_acJ%tjChZvr^Z+7uxWRs^*;8 zBx68HY5jU&+5O^J+7X5OV!_RxMbQ^@#x^(ct1NBIFW%eMc+hEuSx2*HRZq+*wY;)j zGMafV*AKou^HNnek%M`IttMKjhM6ApboE@nJ9dFKU!g5wm6rR^>^vEQ`~EM_qCNI; z+(MOdlR{7L#>TJf-=TD@KQ!(bZh9BJ$TeeJQ0Rl==NLYRRWemus>KO&(*|sVc4ysN z?fM{NSxaqIFPT2i-B2_7&STYReVNIk`-sJ96sWjINCF#WQ3m%_*Tv0?&HGP(Ld-erQL-h)hT;QuFk-B19Md+bx zaz|HJ&5Ih1GZ0_=;xeb_n2b}1PGRGg;&c6d>NSo^u~S#eMse$E+-kMEL{BfdOLk6d z-4S1OFm6ca?(Ppp$=WM)Z6b6=48xZtTWIf`CoL9bNIQ7Qs6Bq}WLt(qMMml^zI&_` zX*a_nHSgNWtEYwe96Vy&20Y!mG;flR;pA667umZi&Fbq?-I7s}51hL6#5}^owaH$U z4fbr$GWi>aqzrlp=Wd2v)_uAwXQM7w0pq@j95?QveqP?GMEl&+d@Y4!vq6!7TeGJU z9W{Ekm%Pnr?OvYyz1dr@Z)w)KDE3*KicNxR`-+#S!b`ibe43@>snc=xQ|R?C9S6Vb^tt4e6+4A`r`^DH_|oUz-&S#D*P)ghmJMTq8{Q3= znX46Gci-Rl=(b~7QmcyyKt4*i6BW`JruZ4Fa3N&Wk2EBtMW0%^?HUH=hf`P^cLUS3l^kLanyu& zUwypV@!EodlZi~udfVfHhKs#d2xs4lySA5&}C?ySC#lz2JBsk@yk6*QqBxolnwbp1R%_CA=Y4vv}+C8*Y9)yQ{?^}>=ig~lyalX0b z0}q*fX6t*bvI;7G@ZECJX*XH8>8IefuiuaP-mTIqw}0^#**?$JyP+dQ?^ey=LysW8 z*71_Bp+DxR%N)=XYTCDN4ZwQF>v^6#WxInXRy7VPSd z>#OyOkC~mmM@;F&#A3e&LtFT4=cV3~<%vsc<|I0P zcFtvexVAsHwWlE{)oPcYfym~jhv9Mf@8Xq3KdH@c+vxLJDFao<5}v%CEs;2SyWr66 z#-J@?FFA+4em%W>_vExJqOyr1*BGOeq567UeeD;4BQ41o<*gilp+YsAC10*>e&>4R zam9xjY*vnA#OjHO3sapGj!Lv&?x&TyK^N~8@1XHjx9L=iW^K*hfcrU^n3Oe!?EN7`D$`_VSS+2F zaNu@Xv7yjo!4G%Tbo&QlPIrB~G}WhCRh;tJQee z(|-A8Fykyj5RE*utz7hepbxx8>TIT$mqY`JGE-T2#@#)*9+^bMs0T-%#SIc=DRo)WoU5 zNBZiWi%e?v$-Br%r%ZZgJs0~>+vNH&G5@IbxmC@cR!SqBj*i3;Twi4S8nu-9zP&px z=xemi-4x(|c!#!n&#q&Xh1pJwVwJ=0u^w6Y7FI2F>PfTwY|hD7_XKTZ-w(K#tX&=L zaW>$(SLUoUXGXNO`Rr5Tv}`<-G^gTCerSYmXX&IoSb4zl#PRyOZa0S4sU@zs{+vHz zH4&Q-aQCX?p}|VmFUf1FmD0b3Pm=9|g{&xSs~^0SQ_Znjw8nRWYT zf%j+-*)nGBR@=>dmy9`2UFf+w(wkHd$E!;){bB?UKGCVWrju*DZ$#>o4f8u@UW)+L zNq>F~MLN55h20HteN8WU)08e*cZc0xqRJZmHr|6jI-4(l_C@j%`Eo?+@$;fHIuT2k zj#Deoi@5Fk{Nv$s6_FcqzIWc7={ms7JDale^u)78SN1)M67k;r6vw$CFeE0q0`!Kcd1x)?(wrR>ev@Ad(=zj&SiRwUaQW3@}=l`{5DQ*yD3d3aWm_Edra!fOcmR>~;n!JZEwBGtI=XjHnyTZ=&g*7bc zI7-~T3|fmtA#Lh0%jEc|Mtq$Aa?$whceBho-_~>?3uAUnY1_YH^qr~Q5$xVm9d+Zh zd77^637;#B{qqN}PDE^o6k=swtST%n-FBgg?GSga{k%-{fp^od>~+~shWgjq)^a45 zpD13E^Tg}Im696ARr0AyuX;@x?dSCcv=y&Vl!VrWeL8B`5-WP}yj}e=hf9}>%ZE(9 zKHlLydn$W_u(Xx>>&#mrUga(iu&Sk1`qwiLuj#qFc_3%`w3JSv`57xMib#}JqX6nzseiPHOb&&2U8(ym2VQD9( z9~*e*R1p=@PEzbL&3(di()B*s{Y=iCo|0C3QJwbgV8Y0mRWTPPrEa=c-qkA{ z&gN{syu8>bGow&G#3@OkYgJhMa^JebJK>iukk{C+59-(y`KobsQ07sM54%-1T|PSa zTD&Z<`04_6mS*(uoj0*czTuJ12fx0Yt0AH0a8c%XgzjZ6uaxCm_K4LfOy^i_e?{5X zb2rO1d9YP5uajuqWjZ|WZTIAL%hxWeH_t633qNs=;i|_h=eC@ic>5tZwL0^g`3IMh zO#|}AXI?!|>RfzP_3+OVZc3Hg6~$hkTHBeoZa~BN`2L$t8;!4|PM3}O{%A~)c09AB z=2mgTyV1v4NkgU!&t=>r4_x+BRSb`K%9N1x*q$ZzVeo?tHaSExgH(e!h74S7&?WtN zcJplX!AB;vZloL*FI0j7EYBeuW-g#Z%Wxc=hUo)tDVBf+p1NSUNX;6 zEOWNk6wEE3RZ#I}*@lhZK75m(y(2T6&Xlnd@9W=aRI=vhC#&vkRdeIZerml7C+6IO z|98e@ovpg4A2}f}@x$%UKM99@JEo^oLX0V0bxzGcmlq3gq)y$>T%&4ax5obY(nAN# z|6gHW9#3W0_3ygObI5R9W|ASpag0$iha^*E7DC2jDwI={sgjaTMH$j5%G5xEilR}< zP)b6mG>B+0*1Pu6^E~(SzQ51s_ecBMd+p&`d#$zCwNH2V_bWf8v@%I$vk!^W`dvot zy5!Az4@Pgv9y)#I*titq3EsnP&TA`uf0gLc)Y50}f>STjvEqfqnpI0p?k<#sr=6Y6 z5XmU@ef7+4XXM9-;mq;Iqvv9cmdf4!(lz|WYW20Icp>&@g&!*4&(-@pzg6!s@tn4L zA#S12!1`5st)cU{L`}?!%C7#Ni_WdO@DI?m-CRo7ZE~uUc|PtRYx~swUEur9{Bs^l z)L*0rSuUaIb*&`bv{gL&o6qE&vFG)-5#_kJ?mHS^_1>~^Q% zq%WhNXnM1ZNLRDv!re>#yu;6~Tv0V}_z%r2|IW+spWnKTtAvdyTHDqf!*EVS=&SVbr?by*L!U!`x9tT+qb9H*kkLr(Bcm}c6>j& zrj6FGAF9roYCPl~YUnhxu_~fMk!)GwJRPrj)@qg7staGVSF)aME@E@L5!z~SZ0iE! z!Q2YX-+?>wUXN}2DxM@$?)_l(dasGoqD$CPTW`gcO{f%Tbxm&1_+Dyf&HXb-uoV7V z^1cN<4(EtiSAI4^CA_8S_wY}ikapdV-#>g0jHONW@wN@VjK+@iJ#sIetDFd}FnwTt z?2O;cAhw6R5a=0fZZYYV#7I9drzbl9<#okiXu<3|kK|`vUFDfe!%l7t_|{>J={akv z4_22hm1>4)1@ zPRbu&wKl<0{=KDrg_~%GX1tl(Qz;GM7t1cLhz+_Ea3~|erDMgJO~WsgrY$xtX&s0= z^7Tk!Vu)XsugX}vHuXHCNZ1md-V^L$;jmCM%Ld9la=r_C4+k_-xQW=aG8c z=IhY4^^L3-X`Axrjc}IO`vI)n(c^kz#bfuJB%;HjyR$S+x95C4mh<6RE}vP-lOeN- zJI)UV-XBqzedl)E-*3Ace~`DGQ?7qc*V>3|`h7L4Z(sA4#FQw#D2$afFw!o$C zC9dBNoaj8nfAuMiCpJdXWOu;*WTEAqojrqE4h1j$?>fzI5WoNQ@?^^Dnv8SCL1Zbz zZI0Kh9Q6ifj!05vmyOI--yCU^jsM2NBRsS1u)rOxW;XI~Z0 z9XKH?Be>g5{aKQm=xo=FWaHO|uCbnntFO4%=ljuklPi#;W1abZmt-WK0@d2eIi z=04W^@9`IIoC;^Ni+l0zt4VT7<72@G@gMi;jBO8zNxsSw)!m*Em7k>-l6X}_Nl7$Aq2BhAmuWzpeAjR@(VTm)1Smz2ai1<0{T* zS(kQhl{D7U-5oKX_PWL^)8-C5`aaTn*eK0%Xv>-=-}x?2&#)vJ?Fgkdmf5ArT`D__RmaS<1zWvZo8YxgQ!pfNv<#|=FXE&S4 z*gnA?_caMq{ZAxE7JjOEbG!8ZzVba;ReR?*qw&hk7yQDVlDwRsG_8`N?a^bbJ7E>M zk@iN?r$)@{h)!P&Eo&rX?ndcoq@LEanzQt)vr&Cm$YpH9J-;`mVlx~%?;Q6(+9KDk zC7M%cpfb;Rx#9Kb%s%_Mf|X-QryJ(^<9~Y&x>TqaT;BDBcAcyBzM1dz_xWFrssa-^ zcUcmpt`(iI$Xc`S5{K3gJ0g`vPV$`(zWHWB_{akF#(6b?TWbpi3YAxcv8T|dMqj(_ zSifI4Z&isD>)|f@sco(f@w8|aj)j)G!%u^?ul^}{TB-i{VC7D3MSi1i%lEp*ZXpS{ zk8#{-^WY9?oJreqSa!PC<|)1VMhS;~oKfb^kC(eXoPJW(`&jOL`25XRK04~OJ=+`B zdVPEl+W1RitDE43?I&Bpt4p7$WE1V2>MfWz`_}$5j?#nz^}xdXrca;jD79aox(Q9a zFZ0B`WNqH9<`rHob;+H>X4HB*=tAbqU}6rJ(Npb7uLs+&Y*I9QU{@nt&Alh{tHWCj z#;2E^Yc7ijK1peOz_oXiOStUlm8S$v@7|=z_bI!gu@x1HCH8GLds{v}_@Jz*sQl5! zT{~O6^~xuK%duJ!CmkOhuHSI#&Zr32Iicsfe#U$mG;cp}j!@#oLq>atS84ts5ohk& zUmDtN>@9ZW>hP%9O{LE~fptcOYj$WA{VIwNPcq6qAE@c|nGohYn4xe}B zIHzs>nhOIDdYoK5HQ#Ve`8d^Xoy*@5zrtb3ttU;fy`CHL7?KAn*IM(%t*m+&8&Z0Q z_bG0FP3M}2xp%=DMNXxh9j&U0Jhehsr(3ucxKc?(;pP^e;>3z#T|a|mripO{ozGSo z#g@P8xpRG~P;;o?E4HiUpX?1aWqvv5AMvDdNk(3l$qLNXG;Mid^C@)IHy8b9$uzdj zy<_p0Wc}oK`SFB=Hudk{_hw1%@I~F4UpFGmxEzG=lB1(d(x>c*8>aDM_eEX(K8&nt z^OW?d_&oB+;bBCk3&XiQ=%mJYaKG8H6ISuo>9LCB$if2MHwrTTJ{R*gmU@&OKB9QB z{4=f4MtwzXUrt#(_5P*9Qn3#5E66STz7PBV77#pOVCEDU}h9K*&JKy>z%i{ORq~82joec1O&|tXI#+p zNixtwuCyJR zSC8tUWOs@^Hq6=cQ2MS=%+S4~a-3G7CrriL#WZe;Rn(8v%B$(>OFW||y4TE)U(%l6 z`s4ZMCliM^kGj>q`a>#<{Cr^4l4oP|8Vk14S>N#*)|-6-@0%M_!s;%23u=w8b}sW; z%HwEGvy`}1xQ+MJ;W1O&J(WD7X1+y-xXDM=9W)xXgr3{^cqIM!-CA|XL7<}IoOt?B zCpogDNnnq=a=c>v0TWy zFfc#w`t`fW3yrw*D=b>(!sfYZ4fef%$T|c6uw491lc!ysf8L)rYN}V?rM!#REOHCJ zoaip}Ox>WPIPJ*+ONYmM5W`1aYBEOyZb#8kNm(V@-cO_?yb}OsA-1GEL z*^Nl!U+ioewBecp;aG|bwQuIM#BIBBWm^L;Tj`{#AJYQ+q7+Z&V zp1B-V^c7I}Ar#dArGM{%Z8BuT-j|u0_Jh<-&osk^-pdx}W~DDTC{51s{@P;J8pXC& zE9T*eBTLx=lV#(s@h35sCx`^`1>HMvTU$G|seqSth?PQX$ATce>};+bbDD{5+g4h4 z#zw)}d#!d$&V;GCSpqxDcqQYjo1Dd#^_4gCtTm2F7>f$A_6Ww_>p!+HJ~UG~o@Kpv zmU316!+Iyv?3;{`$y;}0pJ=gI7tYBnGrioH&{$8VjAjyam27+6`oR^9(01JvvwU0m zeuZn!IW{Y34{-m^n1^4CUN3Pu+CG(8)8=$b=%)PB1j7pVspwVd?m4z8&T9`}@8aywj13u<5<{92qZ9HQc3!W4D;-j+I(mYtGWIWe-rW zvz-1Ry6<^*n28kl+;|f_Y^vE4TJO+mmgOGs`Dxsb7rd`p3XF!PzONq-Y1o!_blsVt z#`^o$r=qM!E{@W8V^v4{wDLnlOrDam#BS@4>&XdeQccJ9y;1tJPOzlBJLHJ>#H#!Y zr&%pSqZl5=TL(Yv5e!}7ea_TDAi({4Q|1eIo3~SM^UPB2ms%eAJhqSbT5L93bA4w} z>y0)2;f%t4?8aVmE@c1dZo}Az-_GvQJGjOCN=|0i((49uLxMDAkMQkvvMcRKoN>1l z-y}bBSh+^0KL73f6|&4^tB#R(M#*WuM^nojLyOYPsPV>9(qT#-9jp$q^WzZQpxe$%EzWHu7sl?hmu_TF2UB zBz}1nSGYS{^A_zG+0ymqVqV*UVa<~Bv9&t9*}X-E=UtBNoH)ghq|&xJvV2ZH9fET> zZpdn34=?O8>~o-M{t3>V<8pfOx|+>*8dKjjubr^Nh7mD++=r;jlYgF09y+hk@g&(- zjnLrvLalz4E-umA<8)PWPf*%T(}E9U+luGBr8L&9^IR&#mgs(JW7)?*-OAlZzfAb* z?Kb}sX!EDupsadphe1{YE!FRXvt6PQ-w=v%w4`v0q$#_g zeQpnTNAF2K#_R2yn5FKzOEG=&^UM9WI|PfoAM@q^QW3CQ=iu~C#PXa_#UuUS_7*&z zTQ#a5+8FC89FlZY5&S_NMO>;6E;%YX_{4_WZ4#^vj15mt@IE$}Uql(k7X6@keRSSMR`+ z+Ic0WN8&=FKY1VDb6WYpjfIeHA-&BW-t`%Cq!8~){rp~~vb&Xg{eE<<44=486Y3L8 zad% z&LJ)9WZ^`r;`uA{+RFC|hq=~&93J^i6B-QQu$L+?wmh1;Gq|-XSH6DM>@=&XKz^{u z{e|AW3-5mJZ8ylhUuzJW_Nph^zPx!hdhH`C&8t_uH5dh7gST~dvqoIQa_y=LZ#HK* ziVM%t$DZeG(^@#d@r-)ZHBYP0#LZZ}%I|uFT~bu>bL`z%PnneeLP=Yvn(UC?Q68F8 zaOXg4=MK@bC}$;UmHLQ~v&RBvqLhgh4}XoCl%Af6R~-{GX;78#*ZR&XvhD0W(PP&i zrXAn0w<~Ky zFaPa;kjRRRz;B|F#tyfxmL0qArm~6m>ruVolJi21s~e7oasFt37?P3mWlYymOq4s) zzUM@XlSj=&W#15~TSlwqamv7I+w}8idbw|PT+i>1dtDn7xn9n`!$V-3idnYH$G$hx zg5Dgot=Eh8+40#ez0Us+4t7=oSV-9^AIq|L~f0+eugw`%Xk#4JecZpZ^nIE;;GH0kfD{`Y7 zQ&(Nsbwb=ryX|oK+9}5Fj9$L?4)czb@;eHL+1AQ*$=_WgQta-2Fa4RihZ5g{^h%HK zv)#^*C*1{TWvm=Mdj&-{hn^9maIY)4k)YE5)1pJ`kmGT89glmeJ#0Fv<&X9pJ+-#^ zbVY8Ej%Rsy;+e#cr4kWwc+aYf5tF-59Pj^nV2j`0d%`YB7Ke%pHmP%wE>yDRQ{NeJ zZXGm>?$^r?EC^#y+c&&4J$qYz=c&gGWfX3c)|TEC`0%=7 z30B|!Nbw!{)AHc>M>4)mL6N`bU+I1;bxd`M`yTg1rabE@Yp%D^F^XfJRP*%DHwx#@ zEX4d2t!As|Ctakil8uN^RNqm>)BG^MzP%wp@TGvqpiuV3JuB3zUt3JFPNjNThwU`lvOx#AYdppSqte ze-t^DH$GD|d|u3zRTYCNdoTibRjRi%NWRAUoK{zrD1Lf;Z+)fb5`J}2!zMl9x9aP+ndZdg zC&%vS{$3=kx`gFBs|5SErI{RrxF6m~xgu|drieuj9$|r{{k1I-ztf$%9m^W#ZhrX>_ z<2n-ZCRl%&yRTZ(i~dn3JMT;Trt`0jH@pKYW{TG`xq@^bLyX`64l{qYhCnoIZf->bxU|72|>bZ{&0&aHS*rK zlMcfy$0D~Ld%t<-?f4qy!aFT1(`*t=yJVjCwHN$gJ%2o>QPcSv~bnL;^Z|342 z9T!N@ejP1!MfO%%25U_XYpfuT-7$NuS+2|zRzYItq*(R|aXt+4qkQEqK15@1;`K)w z)N9gao)_hJ8dS80$FAR+xY_p9>F5rZCP~K-olD|&)D=%3L<~CJCN><7!A{nHIB|?i zx`F?4wW8N;d%x-G>dS}svwWC1zg1^E$6J}>O2@e;IqV0j*p)LZzbReJ6OyreRKhZU zCM(ct-{`=h%EB?BPp^tbGxpBy*iK6r4*D1`dhP5w1F?BY(a;283*&EnD`oE%+CJX! zsH|y8`knVrEgnsM*4%gd`e8MnV;8?45~*9VB5O~joZ#_{x_)1Z`k{2^z$>wqtrg*a z^0$8w;`zDXcT%rI+ELv9Y<{k3a9BpYcO5_N`M_(Av+WOSKRo`ZYQpj%`?tN}aX%dz zY2{@7nd%FHY1Pm8m)w5PJNu5$*zDGitGjG_uH@ZPtvx#X*5Xs2KE-gXUv=|U@tF*Z z*d+Pu@-5q=BbRQ!$3>~SU$flc&*5||t@~lm!QhCS!3FresF*cZFK*(R>~id93~#yd z_UL#=WMuJn^j^ei15kSH{+zI|<2+V@Tp{ln~C^J)7X2h;H6=Q?lOYHl|_ zJ|a7^%%X|!_@#yGPq~N9j?RmG61n9w!Z25Vw*B;;656LZjWpNyKVu|k^IUS2a;j}Q zUJLh*pM2@6IoXC^MAu`L&OZ|HXDcfnktBSJo zeb)EWrEk?Mcuc?9Wv!U#>>~1Kzd>NM_2Ayo_H_59Z+2QeVdvSt>z*^FrqEk;;rkVh z4*Cb9=F%@QT{KCjHF4Rcx0Tipo#}8r${MMpa#E#%cr3PU?dtS{~Q!;nuqj(1Q=eX(Pveuvccjz^H$%)^%dE4QvT1@sGYV| zw&Cpg`QH6z(|voBNR2-6zqk`vT5CT9ta1+P&R>(MwryqiQQEJ3ozrVG<+YZa7|P%3 zZLcfzYj7`*?u)os+CtRNn3l<&>C;rsiG#_v*h{hnJPMXY@|YO)2gXG>@^zZ8xjvR8 zwDeHLIxF{hJ@1h%)8^ldC>IZ(Zn@u;y`i6&8Iw1byWOx${m0LkS5`P;``Po!Iq$MtELar4p+!o0#Tse(K}hn;Lcpa{Z4z8r!q%SLEA{TPt%8_|#@| zSqDn{j$}OY^F76OJne@0DxW{C?OlGQl3yNPvhi@~{Ss1ju`#ugtTXlOny^a84f3$k z3vEB*lH29xsy6%g`xOK<}>&g0;>!YjJ1}r%yM{}~5 zI67@5+Wi{4W|U37+V9MvZdjrFGEk@Q&I`pGn;ySlR|-5$$=zy`F4ot-MVyOo>fq{X zjU_hRDjdz9cN>pw-+#GsW5Bg>8ZCJ8D!*N7_NCbDHK+9KgGhTl=akbYE=-4NFy3EJ z%H!nPlO0!Jn3w6l!uLV*hp_Z{>d_M)$kl}wTfdziyzdv#)^cjnXsGxQ4t}Y-}l>K`p%QbvI+~LMA`XU z(OmB5JqosECdaTgFxnX}B5RA9uFU3JjUDs6!{GgBe(y;Whn|_#!7TjUlIHi-RQG<* zaEq3(i0um}ALMxU&$T7mIhqs)u$%*ME{`fJW=uSbWfJ6FwI^l zMN5yR>c`r}Bco3y+W%0q}ywcjU z0|u+rbtx>`2ZooqTS;HrBpH8rwWjLR->;t&t1oAay5BLD^yc1ty-q_Zv1j}1EsdQT zWaaE#`-E17Ud`WYfAVwj1NaLKkJ~>7F10^7xh^ zB?ABeNl|0uR^(HEVv~Ov~SLcWIK~Tot>dM0wqdX5A=P1xmo| zg-a_b76}siocrZV(;KkM9{yIJ6wK%ZY1lkhk{+-{e&%R)v$Dz)$5W{!nGd$o4Scok zV_7RYdFsZc-8s}|PQNQB&F4IR%)*lL-gtLqx=q;@-n_Fsdc`sZ+7Xd#!7Qe4&(wYP z+*Uzzl;>weT>5%Kmz^|R!Xdd^9cd-ELSk8iq0RfMJ&T2Y5* zz3#t>RQQ~F{a5_yGtrxB4piqgnrT5;g{8XLH#w%0NNS9ZUTnf-+x;rHuf&iK{*(h`xXc3ka}?heaQeju`6m=m5qbAlt%xfPgiEaTIMDK-I!Pz)v;=Xk1{c1%Oxt z9gFaf0cI93qm2N`7I^0X6d4d^5t$1RiNJIwQwZ?t=mUf=p!WhuE<#@h@Mx$2IxY8L^Gf^kD|@Q-Lx7h6;>42zwW_HU!9IeZZv!Iu*qD zg3y1_h7g=HQg#u$7L?IP%uon2Ss&r_8p2e4z&KXDNqXFtIU?bB&8AE-96b&ykGyo$a%+9_$nLVT4eK$)G-P zJHf7CRbbu&v}FL?MT~nq6u`+f0svjWO@<|aH3qo=akC)YT!hLE#4WG^1d~kGrs&dO zLjbV{ge?e-7xVz^X}~De2k!!LB0{PLv{bM-310%92v;KVb7o^2=&00U=x_7tHXnK3Esvq$}(bya0oB5U(MG73gRH@XTm1 zf<=bFBhz%y2B87a$Y3C(G*BUd&(LUupN)8EfT|dr2crSW8qp|0g^ZAnVK4>7GTJ7E z0vtYQe?asEE`=vT^+!Vy(6v55EtA1D$g;3=AbCVG7-az@iV%Aczu;v6kOod6s396e z+D|mkTB_wfyMI_9> zvjEB%z2}Y{arVCy&|6gCUvQD<>Vo*ZaOLpI>iat{lN4~>I-z$ACz(XxsPj7EgdC_U3(Y$U=!9T==ezq^q@ zL_0g(2?KyGz2Xj`qeEZj!F&{W5S9qAT8v)oPKa9nMJu%r3j&f*3ZzPeM~cW{7oW_C z03`drawgI*kb6Knfd0P_!;qa3_NnOqV?hM$j}}5qbIj)=JURL}Ew=vxVIxYne=G=) z&i`J}^1n3Hg_OfwXb}h+^4_AR@R3`TGf{>S{RlJ2{{@0fZ?0QJ5@th6DgYT_cO=B#{@kQa6a?-S3yXPs7?Q)2ci1J3*zaR7r`gG zn8UOdy@#$SFvf>2EZ99f5|VHA6cXj`oc051L*l>-oy$~JzyUE`~JWChp!Ml!<&$8c<)Wb z@<60-3W2lcpF}LB7l#nqOb~2mIThRZ;7Cm}SrzuWXoJ5J9VA(_++PU=NhtJ^G(t*v zaTZ1OuQlNlL=Q+JMEDlvx@e>>-6fQu(J7$>NyR`Ff_2fii;Ee80KRc(3yT6HBru?h zgc0i?3Hyf;E9o7fVC0%G!Vb>jY#5PBFAO6{0yeR0gYCkA#TzC-(kk9c#6VVProz}| z;gHm~Z-ukHIh;@?Em}w=Ff2Y;bAyGoC2%9aB>4Q%OT&pBbgOVen!cP0nOd8^feN1h z9eM_p*tv+u3p{)cIuS%4uBL;oFFHo*+B(a%wU?{%XgNd&Sq8)h%c)qd1g=52a-aeY zp&LdM9^%?^Dz3pEa_hnpgAqH1k(Oof=CFV$dPy`P4vFDvG@(v{Bw`&F8y62Be%*gx zmVv-imp&g&NTCZHqT7Up;R&1&=)Mc?LhPPE`iv+HAP>PuT^k}Chyvj@fvzX$#)NKi z`f!aw7atzDI|73T#F##uO0ppg)CZ~wxZR*j4P3JklZZaLhQQ?jD(FIsZfr(Ta3Mjg zgZjWw0rwr42CFg&8uj7Q1D8#G9(2`)y9X={*FxZ|M1$dCgP1n-(e)JFZu9{w9LOW! zDgt=kkeq-F0$pE#djTpyqCo-jL-?oyGad~>68y|gmy3l(jr-6278;<0m2vK2uR4nh>!-hrKv$M;thp8S7PX1i$EFIsvx( z#gVWG&}7q}#u4%K<8kozk_3W1fD8eSc(<8B850l^tVw1L{Cfx%C3_USRBz_MCrE^p&%y`lIuZx#DDXTaiBfAW&cA0+9Q6Ne@H-kbnE|z1hfY}8M?_f z!kFHj2-zWT8^PiAPaog!qWCC}538zk=ayJX-4mBB2TSl1@n{wu5Oeq!S7>_!Xk% z5|$hc9q`$3m-B-B3`i9~mK7Zx&rGV!$+l5jqM^l9oJY$R=<7_Z&8)1gwapBz&9x0F zmS)zLW|ju}mSl=G+1kR~f?~B=?*Bi;hyUN7DvBOp>3i;hzw4oT<18#Vh6)@N4mb-= z^+GjPRIowOq=aOisL75B4kXh@vQ4PTi3%zExYKov@DUyjJnKNpNqCyPG zbdk&nHKCUUmq0R0B-? zNCqausTQb7!D&z;OOe1JIMo~p^-zJ_OM-pkR5R2xK!qWafm%3q4Qj&whr)udL^4nT zr<$UsF)B=u3|7IZt5I_mDxijwn6brdG6$no=nENyq8JtJi3*!w!C_B0mARl+gMKDl z4cAhUVW?mtoXXr1_!*}{jNnvc1d<;rz|S}pVgx6-q5@(NCqZ!HBwr|~5Xm?dr+>{N zd~m87-65NZAV`~0cOdG4c*dy^**Fyf5vOAG?rdVc8VUJ~1aX0r0+5&)%^+iG^2J0h z9AZH=^lqvueO(SA#X%)wqM0uAs2rl2lL3Pw)E8p<~XoBPK9W}C}3HPvIZ3^v87Xsy zSTWItaptPvGmHV52V+1aU<^)7^v(jagBvv=L@)*)Cfd0u<3~+!Fvfs5!5EMfFrWey z{n^T#A&Q!iH*hMbij&~{V+@oD7!sHmX%r2iF(6%F3<8Z|swsn-aO^RLG$y9Im${5A zYRaJ@aHuc_1TDr;z{H&VnZ3)<(vUJR1_}oT3JeAuRg9s6iS2p990HM#G2kFz3^;)p zLmLyT4`g=e(oh>k3j?{GVStH!DQEV;5yTj4P_YUVm#SmRP%JS}05IU}V+<4-3=|^_ zNG%w{3=?0U#~jBzZA@iQBrxnzfkK1j=F?=xb<5A2GbcRUG8pIJ!Vnw~1 zG4U1@I5-KzB!2&4)@a~uoP^FQi4zrks9?jyXJVP%5C%92T#l2LpaPv_5;zqn@naGq z=+NQ7LI!>Tv*ILhK28F6;v_LlLJP$QPLeZOq&q;L{bf94qQEKWiWB*DqXNe~7& z34Dx`j4+8Zcjkf+DL4s|3r>Q7z)6rQaMBu7n4tn(j*~1fiR(Gc!4Mia2{Hvvf^fk} z;4qwIk4fBNWA?6t_`^vM5jY9li<7|lI0^iYlfcn9mFW^MOk%)_Iob!wSTKT|jG9z9 z`8bseBPqXRHsN^UROCM@97>#u97aW(r9yJRsgN9SDsmlF3>C;%RB#Rrr=s&ig|vZF zk?*JkMl#f9E+~bOLMoXxX^fO{m|2s-NM~*_YqA*7-!N-(n51kkvnG#8s_`;w3Yg?B zDzm1DNsioR){u*-kcx3CvMLpMgQ|gyNo8&f!WgHbxTX4|0samqJkoNA?hA}s-HI^`g@=cMcz?!_>dN0>=sY|m z#OpJLT!$OSB%>cNLa4j=4wJ{jL*{;jPhghuAp758lo8@2MyhV;14f$%v7Nvu*TJ_1 z4^s63qk&WmVub3BPhoU<5Q$0HkADWEz-cPinS8vfvX?yZsu2L41E;cJNYI2v3!aT%M=&T*T`+MhIc}$|>bhq@{u3 zp9552_qRED4bm=oIno*ki$_`+MD9?Kl|!~V%3#vB^GMfi=aKnbgLEX~r1_aF8U;zd zDZQ%h)p;3H0wO;xBZ;tjWMtkl^h8Bzl9(F)zA-`JVcb6N818ROpe9M< zpLyuY=$|9MgoO>(5oyGL@F?zz(AbceAaqa=RRxWT$}7zijTE1f#>&#^R6Qnz!DT5t zEF=Mzk_?fF^Kx5SM~TY`{AGY!*jj4Uw&`IAP@M2CUG}&WjtcX+N80 z!C1?m&8%Q-odA^z`y(quuonaIwhs61F^i0+W_t^jkrl+Y(A<;A} z!{b5YQJM@o$FgC!pXWhgYeWZyMuY_U2E`<3l1P+a<6zxUF}|@esD;CHZ?FIs_n?6= z#~i}rg8cmbP|qNM3@plTGSDe{)PGDX9NmMc0nWi~{FiYI>S;8RMqq+ZRFEHM9NGFw z7RQDK%e`Zxqjeqqfldf#8d1YLj63H)^GO;$-oC+LhTuY|nQ|tOelCV?M@?PSL`fRl zZBTwmdPE}mH-4!sJ=QO~3)1t~(?N!!uyummZUhWv<*)W+G*&~9LI=-b)TzW^SMUgu z$$CUG3yF-ymh%&tdQ?O$4%;L|VNqdA5;7Wxy)MWk!F)1O6OX+sNN0eFQ;@Wc*j_;z zRgVm|h(T zSQ-`N6AKK1<qq~fN_q|`=izk12KV|cm~}+*ON6uyc5DvzEP8f89|(xq@NP``L6>5 z3IpqhL6U~IuP;P55YOe6tFB()srIR&QIN_I&q@Cj+xV`(P)WU3~+?VH(0?BCUehY{Q-t zpp*4T^!yZT79waOErm#K$BGiTr2S&KnF`L46KVblkz~YSN=8-=4l~@>I|TNHIP;GU z(8x$zDi*(#tg*^J)H^DeJ4{9s4(OLcrfOJlQ`W!R{ifw#fse^TE-#f<0^j*V$^Q^T zZIw*pY=fFG+T;f3XykjXjGV|XzQK@?$SJH0VqYi6kGN0C%@v}6FVi7&uT8^hA?v2( z<_ST<)+0gB_wCsEa1Np6uu;OvGS^9rXCZ*os%CV`aijZ zxQyl)q)Zk!pUmNuh(-b5K$KHBa;}cKA7HAeaYA=QI9EADgpiKiz(=JK>oU@@!8~jV zfl7hcN}-boR2qp+pwOT}r$Y~oOrSD}ENH+S7N>_wpphVIP-rv)jZ9}kgHE7PNi=9M z2s9dJ4wFD*5Lv<$8jC<S1SYx}D1%AmOkqI-7Ep;K0+R-NQ6R}O8FZ>J6TlvGd>i%B68K$k2gpbjdLN`TDG zfCdp7ERrx4j)O!b!cr=U3@tK;J`Gw_3Ki+tg_Q<=sh|i@2FOf>^CZC-4W zNzmdPH=I9-L?hCK7ZNB$gr@^XLOdVgM3JyPSpB*kyRkl!JQT#?C{`H@M?@s764qbk z9~~YW1@QrNPezmrv3k;U?ucnvSa?h{VJ!@h5T8QqN(^cVu!!I;Y-EH|fPNlD=@rP~ zPMiYb!6QScTXG0%C4gKg#9Aqmm}JOM1QK948W;=-&={B{fdmdlqa#wekZ7M4W9KSS zH3_iW&kbp?Ke(#!;>E)HW`wmAI&ga$tAaF?U@^e?y%MY%GI|gzinxD~RzvLDaXNfp zdZh1XN209AB^O9JfkcM`BX>lIEF={pW2$oaiEWgxUw8V}na!^-f` zkjQ;lbtr4HWTcg`6b9=rDhwoY8mFuJi;VuidIJmiueumW%QVjFUu_|vEm3Lkfj<~x za%7f4o#L;NH^)H(y2r4hf~a}`9IKqbs$u@e7OSBjxUPJLBO)Hhu!27Y;RW!+Ie!iW z{#{*wm(cT9`Rn|_NB%1ZlZ2RlleR+6oCKTr^Tc3(P^Su_dx~>nCZ~X3hcD8^{AnB> zIeQAbgolOL!*OB3HvT*>i0mi>8h_1`Q^^Nds{U@vj7T5j=wxNBXC5?tBnlIC5Yi#{Z#iDv5$rVx)DjR5Hh{h5xl& z2t`Ot6_$pfl72gSwEtg^|D|pb~aMG-s=q9ua@MF8Y8 zK=)(@6_GlRB_XA?SSuwekY+%LrTv--frt#e{AZ>bXJ5qRJa*-LC;`dnY$$9Yg;Iz_ zJvu~J$c}mt4<*rEAD@Ra+jAahQ^#-nfP^o~`q2d7AM5M9`YpqD7(e$Wf z$Tno)4pIuBMm;9vMhXRHBHllsVg)@&Aa`HjM5Q5E|9*UMR$z<-q)lE1!%L+i&*f#L z>-aBVzY4=GLmlQS_N@f>W&+cRtyM+)9!+Bz+J`bHca&tcJ4vBMHQ+M1J_gtKN`*FY> z$52Gc4JG1sFM#mKVJxKpm!tg~5Br*w#C!1S1 zB{pO($YWr31QH5XC@caBL8x$uRPZt?6*4guz9|3`INyLIAjwjxXqqMpQvoNiAfiIi z1p7cvr9$F|4FTK0DiR42Fn}VU2WcIE5P?L21ynTipdLx3QQ44t$S6pH^>m2i=s7`> zq6t%>o`O7zvce#cVLHs=R9HX*OcbyJa0nQRjAnW&unu?(j+V>>{?W;F7-FIvLrx{5 zh=~e$lFR}OMMWVNCs{)dVpB*cRRA^=3XrCvAOg)NR6t4;Dh1Xc zQ6K>cfW~1B5)h4pKL8=oNN{ujJP0)42-brmP*I^NFa&50RslL>qSAx90UZKG*q+7$ z=t6}}=|uDZnJj>b0Kll=R&=;5f+2_mR5ptVsz$*e3QS>~PDeKaKt*Su>;oM-6N~_* zOozKd=s_(2cBBD}0#pb*05D^aA-qrl_ApR50sI80agGxc@MV?gXX0FByp4giaB^;KhhqoD)LMe&*Gy(07c2l{GIq< z9-8>n{s&z!1rA*h1eDX&!2RHRIB6EG^s%E5Mvk^3)k63Dqz02f}K#wzo| z^)3oqu#lbpBmtN2sML%9ht%MZsMKgg< zau=Z5hYxWp5o{0ig>>O;Fn{U?_~PhiAzeWVidR%as*o;zDM%dAki%>9L2%`;F9Y@7 zpX}2)sOx`$6*3bMdklPd_Tbh?p;f?tZG+CiO)_X~{$qf|(IXu)H;R$rl*Ey_pWt5+ zy*LytE6d@@yj1$nBZFZ4H$xoUPlfc(8TtqJqh$`N#$Rfq|AYINBcJ+k7@nU1AI$)N z@WkX;4xMwh?ugWJtTeK(9VdweAODj)gL|CX|ARh+3zR;ZJ^srF87ao%7sK0a4XP9EwDXl zD4hRf7O}_T^$h+8tB|v}ta5QCnz6XwoNDeL6(UEl_;pA{CytI_p5m6|e*~VdJ;B)_ zXP)9#A^uNs%aP;H!S+0!;)Iaz@X4tL|7!)#!miCV~7j-J6e0dxPfmX2YZ>t7mW{nY60ejJ9EMdDa4lu^JKCn_xgA42sh0~qda z#v(=Vy2`)Q$*IiI{QtYKEE=i?8V`*lGmH9c0~7TG(_Wz!uYDdOdkh_hHS5`GRVh&&y~nIKhF zaJG_ycv&R!9nSJE-hnqx760}fgeyw^d)z-oydh*UF5BE>@!uOHKr9@#|GGf}e1u3$ z;B+x4;Q24>(9g*znG-lmG)h_{Djz@!p^v{jr5^qCF2*kAxcsV*a0txFU-}>Y2!ppL zA&7CjDrFKjF`qz3(GVv)g|quJ4*_E1hJk`95R40h_)xY1sce9o2~;vHjT<+AmyCl{ zR8cg>Ijeuo>h~N8DUwsFLOCypNX_6({{$<5<+!R%n*k13pjxQke!`X+oF6#NPagjA z8pz+^HG0Z_vjA_Os8VS^gD29Wj>o}9Bqm8_+(949JUy){hh*PeZFk^gb1>I!QG2x*G_QpwNH@tYQKxLQ54A z3nB)TBV;0tEDS7=0h9v`DEdg~1qN(MCWA-9?KxD4fcD|ufdG&mT?`c>4QK;mC$P!r zy)+GQ8JPi5i-s1B4A_{4){`tMQ5a=_LPTpYs9E6t0Ig}E7(#WxxtHf?fm14?wE>hq z6dK4z17Jd7z+@VL4yYz!G8Bez>kUHyFQBky3Bz?hRDbY>0bOWyh?WUZ(?G3=&W2iy zQ{KTnJ5(LOBBvyUG6mEI)g#IxdXvpS1%d@oHv)@LGeQXlOrupK30#1N-lLPyf1wp3 z8SonoJr(q-p9FObRNv^TH1vz33Sb7+7i>m_t)amJ?ZDdu)RVBB z4iTJA0*dJM03Fp39ZrFcY6uEXsP}*~w4|h?m+SBx01Y?_Iw+qG^nnmshQjd*(?PLx z^fsRc1s{V1`haUFs3*~Gum=WOG18$9WT0nC2Y+LrUugZsprdswD38Geim2f5$N+pm z*#{*iSs1-%VzA&S;BN}hQWNe0Kzb%rY-}(WC^1oDbfCiG^A$SR1@tfUION5O@Gb9}6J1w20uX5QxLek&aIghdEcK7(To)=8v=t zwIw%%Zs)^eAj)vAP2su#y*6FT2mbjpJ%A9LNcnHM2cE}C51N0)P;h2W00U>{{{GAo zh@k*}74n52Z;5mY;Flm?=ddzJjsV^iK#%e#czIYLh`0YMc>xQYCvtfdpQlrKmd^v z#hW08;&{&E0@yxRjoU@>(mZrDQSs8CEaIFp`Yyu(V4Sk*cfxROiYiR&FUsJ{9EE|6 z|CfVbuSlq-z@7ksS|QuT@wR`R>#xTkAqimU7oY+S`~*}eewUL`#?m|2hZfRcjRi@s`EI3(ZOH3yl;qi^VfBr~&>LD5C^UV^sl5BufKchG#cX#8^&- z021z#kygRdxYf#U*3CrkOY{V}gpn*9o_k+{c;i3-CMp1O4+rkgL{8)IYM4JNSV&UH zbLdA1?KnLov>Qj@r2#BMU-;l}+WwLf{?rJFY5grFx3qzL#SHT34fceJI zj|xOY+Eno(NVtqNffsH!IEDewDt}4^L*Ntokjwr!I!Ixtd!k7a{Y7&m_m?x1;EsG# z#G4@p=HvIHRRl*a6(m#59q`EVRMj|IUXm2=GNZ~ z0v)-k1VN5tfx4XY;e+gpv?$<}v42-fkPy+lqWjAl#FX)Ph!31{e}KSomZ{ohW8rhW8e zkY>69$H~M{m8&AR)$nP5IvHREZhT8!h?ho;Pr)r|J+y%8y@kLO8R1vw{1|Y5BGlm} z^@ciLhKu#zs^iV2(13FsDhdlGK0qHL4t!wiuB^vDSOD3G^HG+DmsNsV13v7PKyRGs z3$Nqr+_-CH>L!%r=L>V!%4oyilKgrm0xU!%Qvdk`9sC*r*hHr>5uRarQr#9Uc_RVA zRraP9=1c3Y=*z$5XDf_8$G%!N5ovO}ItYJxZv`uGPTo?t-kN^jK%;^4*kg=KCm(mL zyd})49&B#7E@G1+B~58Qvzg>NhsV={RGn78d6nUV)t%9rk@m%n$7TJF?;A$^QGyCz{yxIxrj>5B~E#D1=ruS>nAKe5omOuEuMTV>-3n3oc!GL z%{@+MSAYGiUV3O%OwHFRnH!axw?B!IKJ0Pm&}QE3ETakK6-G4U&-zuZ6YRNpLq^HZ zvo#5alU?F(3LLn0XirS#nTzsgYUZCyW3G_E62NUEbgwkI-tLLzcSldJKs|y zP)2G_P}qkEgP1sygj_8r$Wv4(p|aI`#0IEU9>gHc8n$deT}C?at4l1&0~jsS(BMP z{g|@Aa9M=)xAP(Qq*Uz&9vKxMqt2KOp2tp4W`F7|!_>XV-kWZ=B>Y)W><;EGzKiw! zbd1YU%o*y(V;{om_1uy_R_|ta@}@mIRXwDibopw)OHnz~F!kNUr|Ppz#f)@Cv3!?s z@fDICk&DX26FTDEmPS5SS#bBuLf^t~XHQp!$w-EI&hw(_#7eqhRALuAZ7|v$g{(I` zZfSGl3eVH@_Y(p)_iyWlKk3{$v*6BfUH-s(p`d`l{(Ubkn}cq>FCm0hrHrswoLv5) zV(gCW1vP&ckq6FG_l0k@>iJq%8sF@SR}e1>TvPGVHT%fX<&`l7x)-{ZUcNzU6+9&N z@%+c5cXh7b$8CGH(tN6+wwyRn`+U_~e*{I2T z+l1@#m4)LoHHGUzKH7oZ#505J~^4#WP zFvAo{IB9R)wDv@$yDNR$lD%*4AMQJ+O1s-{uWg&u?l3!G_{7;xYgwO~XJO=*ZyN+l z*-aj46K^G-Sh#P=xW_-o&!?wRk@ya~e7*hX3yoy~H_AR#WhZ^nq+WF2;_ozf>*gEm zY-W!2(v&Nd))>34H>on#ngKKVk5oqu679cUzjo?Jxma)>OI1DIBxjS`$2!g0EjTIf zm^zC*vL^W|ktZmy%(y+YeN2As`w!&2_6o0c=k4#fBwT+j;+=WI;d|3#(YMEQ+`nO` zT9RE_y5a`clr5~^`sk~S>Q?rRQd^4L`lN%$GNVPc_NX_Hv>3+ZgqwY> zK6vP0+>fh@yf=iD+2^WFHT4T*Ex(;}pjA;cOo(BRKAJA1wG0Z6iI`^dIuLHqL%cSh z9XHBvl;5(4gbNbbc>dxm%HXb?fcYJH4(IuLjLRh=Y|qj3JA`N7w@VbL-3c(?E%m0V zl6@{9Wq6*;1A&7wR?E({3_1?x>~|LoPSc9;-hQ(6R_L}qi9&@{En6eP$hTvfZe(~` zOYQ8f&nwq-N(c|zbDFX=M(kxx=xe*(i9@%B0(@2{4lMPLEOw0jRGL}2SJBW>JTr3l z?UP~0>t|~w8ZL?$Rt$wxb{;9A1dnUrHrsuj!@js{=BCo2;2mNcmY3-$thG}}!8M;% zWz5+v`DJw(jxj#(L-v#P-n7}Z8LG~jx>x^%$XjnjGmkRa2T%E#FE~Wl-AUs2b@c^3;F89TaL+j^3o)6LxSS8fqOU^Yc!&DRMO)1=d}TGjOZMJomcqX zmH3ipFcp=%*}?{{I}^{gCw67jS0uH(8n^vSIym?|>*~7{jJ4nsDUVLGnTxRrtCA~Y zjNjkbf92Fy;RJSE+_$f$FRyA|)(WQ9nr{iSlu=y$T32lsA#Cg0>wDr4No;f;Z&~qq z-F744`W3PKMX4Y6l_O>i>1H|P1F0IWmivz_lk{r#IUSm{@Zhn~vnCF!YZNy=v1+sY zBqC;ciM-eHl=PN;_C~%kn}(Jp;ya?hw^@h1A9*t8d~AoA!G&OUhO@X^?SrQkuTI?C zIhvIH?8a1fo%LjT|La6${`O5p3wh-Q^55!C=e};3@2{VnJ$4^6W#?2TUNURohVYiO zRF{(HZQp(b*HC?B#!lK@59{`v_4Altnz~*zLfHT zYKIf&H=WtzX4YjDIL&@Fljt6`>3gX-kEZ62&a0{_4jb0%j09b`Qy37uBz)(iiCAUp z)o(Z4lrx|CH?HyFS$oO#$;P)igFEMsDy){qxy2s~n`75H9q&9=G9j#Mv?BcICBp#Q z&Z>oOl&o8a9z0CHPYL~`b*=2Ct?53*O@`OBd7aszjKP=NS(frenK#&KKg=G4NHbJf zCFgHQmAo^adJ|?7U9;W)N15=3;rQKJ`|Wq<&i~x?d{s|qQKIA<2R5JQEl*R4g12?r z#*)4w$4VAnjHbFuHExh$?m?aX?<_hudceOU!m=uf9PbD^rX=imvzyx!U}vblz}6C3aB3CvuLd6YdksM&sv{i6ECi2B+3 zyJ0V|e20Y#Ud}Gcb8XK&^|@Hr;==6ljAPBW#@2MZJMKB;_E00GFeIfleC~(1B!SQ5 z>}gyQul3TT`#K_gCX?D{A{;%_h%+AU`}iE+T|Mbs7wBCv2UO>d)tK*PM=Vi|zrHu=40G*Rx0tfJd8B&+`%U34JMmEyYQI;^s@td8dErNI z`h-MOZV;g6eQ0Ih>>jtpI5%lB5S2}Z({6lg%dbgIYdt$R_wkKv!k({+c3u~3UgMM$ z>x=tF+}$FE%giiQ1iSgE+uuKJxw1YXV`8q!D9KE;M7L?1J@%LZU+uQS4M)gN9d$1& zhfh7-G+Vg;&G_=mX3C`=!BP8iS_Ih25^Yz^9ea|i&F1}zFR(L-Az`Nr59+v+Dl&uJ#@XW?<^vsCTGMp zzs3xR^AsK%HS#y*DI$L-w13o7V{U7GmFc+hlC#{C+7j=&r%MC9mq=WApc9n8-}vg$ zK@6s|FFz|G=Ig3OCf7f?BpwX*(do#0d)lI)Vf5hIT_F)up-oo`bQ(UC>>DmSwrPZ~ zBr_3*iP^be=}AM25$TJE)PzaHY1&H<%70~3x_2{QB%cbg%r*~f5m8#cxz}Nbi|4#O z>v(tZnBNe5_O#Q9d<;h_2!6Oa_tZQ8Bc&@EH_p29GcF#PHXNMJHlLH_^y0Km_{yD$ zyI;31d#hqOu}5DmU}yJ*#Kr61M69^2=3bN`efS~eQkH36y3O~VW`{_+d3?+@zmL63 z`M%F#Un%R%Pl{-JUa(Zs%}g)j8h@9j$TOR2F(;eetd|{G@;jW;>+?@Vv1ABMZ|9w*Wl(CIJ%+Fo-XRXkfvQo=)&ZF!)6R<_zW}Wgzz1k;U7d|0k>z+Tj zrKWW(VZFCyip=-J#DZ-&DJrr#B;TlHrty%w_pK84%X+(b5o%b&g@TFcJyi$ojyD$Z zZ1$Sv$rjJQlM%v`;Ij~S99geHx{`Kl&S{g`5XvIUXD7r{&b?vF8=a0CID8379t)Mb z{#~oovOIeC#PhW(t^!}=`6u_^$oOD3_f1N~!kM>w1mr5;eUwjKnw+}-Qs_nUk)F+_ zTXa4gzG7v|PMnxA_fox++xM;mw|POf$&T)97w_8IYqPq&zI^g4MM8dPHqb(S`VQte zP~IQm+iG-0v|-;JC#!@~>lOo1@31w_W#ZV2-BX|x!efnZY z?sesuukYrhR(T3*s9zqQkNvicnp7?&n}0lrZ6mOjwdN}uUp5fz_Q`Z&gB0f8=e-s| zIQlZ}zHL0k<-tu1iSj`=xrdo!nj?WtyHA#UeLx$E%^M3!v){d;U=e=IFF8U^Xp`f5 zxuPeciwit9+2x3?eY&ObvgYi|sqmb;$tT}kJMhwN+2EQf!tvd9XV2P}%83=V8BS!2 zw-+!<$CPGrR(Ge!S^NvKR?XwX3X^O?&=+p5MF6L zWj5!_>r*dw2z7Nk+i53fYh29OKD}edsWj(#xl|u>>&V<&ckEBg9?*TY@AIe>r^Y+Gt{d}!L4rW@xChw9pDeGRO6$L6R zuR8wp>Qr0b=Aht#mz#^)ZjYr1^eT$Z{d{)&qbaW#`S;nSd}q3!7JX0?{IZC#GQ~35 zZMaoT7#li6KFYg^Robe+=9zJRUt0O@>mD1+g^R>yMvrYPw_873i=z;Ex72n`x2$nJ zFPquA&O_F&*sJ-Hzi9MBuhm5assZ0VlgTuElV>j@=bD=b`$p)fWpK1g;p#dMek)w8 zqqNx_W~|SOWgDjamfGLGn3*tqQ0r@|!HZLsx{C)dKW~rCzt+Yxw4;cvaJnGF=tXn= z$1&Lu!=1V-q8ghqn0N5hZQB4^ znSbAvrRBi~8?(O<;$0i?*3!!*>7r1j+R-8Lbv9?O1)*Se>9 zR`U5MmDP&j5ouYgFRw}Y#{S`1=%rV~#wr&MSnf+1nd^^VV3BS#(A_nj;COKd>p;kU z*9&s%cTTO_v6J9;oGpd$NU!+k~ zeD#`h@2j5aAbaQ48y?a(Y!P}x!VP9to$%=sDf6K|Qor9a$L;)GyN&Z^+4lHkQ7cIw zWWr-h*ruYTaUWmaJd^lX`$V0m9be{jd771Sgr>h|uF9)uQT3Ce^S)~QaPZ%>E4%Q% z&a?N&U8+ysPisqSyCq?ok@mWdzjF(|^EM^xbcIbxr_O-m`-aG|cZ$ZR3&*}Jh;5k2 zZ8;@QV7gnUAP&Ks_e9~6S&3q7Yh#MbX|p2>+S`x!2UXSWle!S{g5I;H>{!pz*&1o< zMoF73Zp}_%GpY9{9*h;1*DCVPzeUQbJR^7IU}l+N?n9EJl)#xBnr7GiI~5PSirbf; zs~LJm8ylY|iYwErQoYi+K$9A{&M2w8H)MOQkHEcBhizx}NOUw63|}>;C@f>Yl9#hH z-JyR9N~D!vC+4iWBhhH3Rab7gb)!t)aC!J0hUt{-u$Pc}x6?k~1y|;38BPV>dHUta zOg(*ZVXdsu>EUaQeX<)DOXuWlz1w2$cv78jfczxpo0+py&HNwhBenJRlxn)pOb5NE zQK#1A%JX^hD~_BXTzlq5>u}V*`I_Bx$wuilHG1;mgJoe}d$Sn~ZH-6v*GlixSx@7~ z48D->8kscKr`4a$t!#SHwDmb{+e_E5a+mRoH-E&h$<0hxKGMn;F-?CG)VCUUZ@*vf z5zndd<3>7TamDX=9$IEOImGnNHF`G?x;6d=h9sH$;*e)kS7J;|r`M*MXT>K=%-Ccr z#^4Xy7C)O9gLkV#PmJ%*Y;bzN$3bh{mQ0RZa^LG{dh5`onG^d$Pw5e2rZo9Xj?Ez| zJh<`HY&P8E=(|0Q%>ntHhR&5QURmfrf9^_Zu4t)^k0Z5fU*2UE|EM6n9#htFSx~fk znNL#9bIu}?f&J9IbZ-tNf2je;TL#|5ls1r5@rnp37Hh_ZjpFO{g*Nw&WRZeAV}m zR9t(+D*BvN%GGmx0;N-oM_c0M#Es44JswA~eKU-cC+pQ$<-Bbriw!XPZuFY2uoNGR!yq?PH&qsh6W`=*M$+F?_>W0d^XMwj^>ymq{jdOzdb*3Tx3 z%k{pNe=8iQzF&zb;7im;cjozI#=lg|DBbklUF`9*FPpEvJ2x&dr{R_Urdy9c(dIdf zUAoTLT6{8N%{S@9vb@>(2cMOm{Jc)luiaTT&9bMMy+fvG^>nMC#o^ix1CIi^*g31& z0goS%yIzk@J3ZHd|2jcxa^Ta>Jr?)aC-<$W3AQcH7u<7Cv+M9GCGsyqP(( z42ROOj=m%Xt0g#%%{!aVB=1?7mh~+4g6YdsF30-ssUW#-F*A~J2+3(JaY$8_a=+K1 z!TD6;>f$kWtrIEe;ts3%r&~7Ggq*+kK=EV(TV*!Xr1)$U(NnNA9d}uFPJO~4e?*T* zwPJ8=&#G5miwi!!ZoAhPKH(ZLXRB94C+3rAv|}`qWnqzd;;yiuWbZ9nMRzmaYf;Ol zy|2;*8~3iwF0i#qj%brfs=vL+c-gVGxGSyY#L+p+!vxtAO&&Y8R%u&kEhuCiowKy* zu-YlprunKZNs>LH8J}5eZjcUODqe3I2tE}mo5-_x|CX*#0rSah*B$)UPhKB+u_8R( zyV-Xa(ORi_eHF5zrYPokZd1AF%&k6VreGl!LS zx@_5d#_4^A)zD3YHOp88q-pJ1(KWzo>ZR!+y((p~;ncn?hb# z9bEOX{?4H|S;0?g^Vs2;a?=*4FVB%(hMS}e++ySbF^;#fElD zUhVSU;pEPugK-N)zS=5B@*O>DHF3g!ed+@Fn0=|uUi#8ilQP1>o1bkl(Om0t2Bjjs;6Zc@hjF`|&LY-8YZrF*s`{L9h`uReE(@l~7spk9z)cbI|A z)mf0f@1d5hgZV9^2H|(vr`Ubeo&F}Ww{_Q6{LplJ$JL9Pa|uJ z7E>^@vD|rf&T-}1Tm`F``$Y?`eIU#){Z`5wJ)6X@KY21VNd7Fcb^SL6+sxSMYQMLt zFooB4B~O|9%)!rr<0NH1{cX9CGdM3vi~S>)R2v&jd@nn%IJs_+d3p%1!sGlY%CmGd zDqN^hKZS2#k=fwmj%O z!M}^TJ$r-ad6RoL)O5lxF4R5FmT_HyQ;Ev*y<`z~UBG40R+X6C)`IJ%C6{`|%&c%evH}&R@PV7(ce=Kc1vi8}k;XJeL#q|e_9avd&*;k%USNSe^ z_4dj2%(BH#au$CpJ!~#0@^*#x=8x%A%l71flTBi04&2^QeeuKowF+mq1oowP*U5O_ z-SE}krD*8=u`AbZ5v>;^8*8a0{QPTl+(;)cX)nDjw>jRqHuBPu3U@or=Ip8CaV|13 zHdW>8s`50JFUep>CwlP5JoB3y(a6?sapI-qyKkkx?DBnEKB;D#cR5(&w%f!~y(w$e zCw59%Cv3L#Zx@vheY0fw<7UARFAGmm{V67KPbYL==>->@2`I&rm1bErs|?N(Zj79) zeOuz4=DJPwx_HgTGp3|Bd6em_`&V?k%AS-a71bJO`$SisXH@qXE3H=Ut9dS;8<&9s?`z)fKzoaU0ymL}cuz67!v*NMAH@Pg|%|l1#>HPiz-B+EA*kNPQ^}K|gYCO#`WvkRGc=Nr_k9X|6PVMV* zwwU96K7N;Sw(q29RL}mKuQTLRE3dEbteZEs-%8Sp5w!k=^QXCIHPmvF-kIHaTT=3* zJ9rA4xJGlleKs$#Z@{cA?8o|>M}PSG7F|4iF@D$QCmxsPJKSL8RXy2m+gA6;cFFSv z*VskhRu4D$sh@NnKpO4d?M>X)aikpMJ}oMGcfa4__XT^y6_V(@W)h0NM|}ot)t~VZ z?DpiH651QNq5Zl~%$!OaX(S;Vv*q%-<$+mhEhj3>S2Jt+CwNzqbHq9h_J}z#zZ{RF z@{B1w4(xu=7M~&T6+hds%i!RbWEO=#{JpCydr0N6Yt%sx(*vG4(sA2)F~kIM5p(rW zH~vEM$!In6&bO;SjT_Hwd%c!7;=r>UV;u%dX`RZY&@Vn;?Up2b3So8=7g(0m6?m2A zvkPYSdr4HjIkj_9zk8cL4<%cYA%3eww_ryUjxq8eT(0p@Qt4lOL(*7`75pTN4vH<8w=Gq;@Y@xd$I`lfT}Tq2!Yw+#o-mtrGm{v+FXR3nhvriLw*-%`MLH zyLEYpwCbuFT3+nvtokHZE&O3q3R`^Fxf>39J2dByi5x9CJKMuPX{*=o<+jti^2~$Y z2stlIP3cH%lqCtr+_mXC?u1))zeehr*kxo@m$eh^?0_y6fi1E(iP^Y^yWD~bYM&~0 z1}v>TRlmb{ga7FJ3|d^%Rq-eTz9DwN@#rhY9Y@>VTAewui1zJcotge+*;_VSav zU1Fd@`RNPGc6NE!1YKSyxUzPpx9#Hx4KdQ{630ZJxgp+>fuu`w%w@hLFFJa{aXpD~ z+vCsz@s}FYi0kcF&y_yswWTT8;un8^CirG}QC9&owT;}tHY{!4{rKtlL&hxw!~TRR zevebb`p@p?x-Z`!IQF2hsI%etLEG;u%KWLMoX<-l-8^$|jmWy(e0wsMcjN9Q`Nj79 z%Mb2M<<&aWd%1E{AcW9BTM~SGqp;?O=Hs2&X=CzlaUZXQJ4*3h37odLb;dl=z98ZI z{X3tnVhlCAS1zTGu?uS(mpuCZb%($8r>}dHd~h?(F)3$8`!s~=)h6aI87nUc4pGXO zX^RqFrq8F7!#?0F5Y!iy8gVn+zDBsdBf~REEbTic`0(@Gq47l@UyLVtM=uTDw|x$O z<|AyYkbvAd&l@UpMK(~snY#L03AmcC^V@TUhZN!SI!xc@<}HJ5>`+Ik-n&mVhn8zB zA|ve!uX=v5kWmlbk<@WiEwIMSEA;YgxkU9leD+1PHTN(V-B@8g(NZ0*?yVlIV}ca6 z%jm6;GjXjJ?hh(?Ha{10yewGNeM7LlqxIdI74LKmCqjm{1;}}?t~#N5QukHd(CtIZ zZmW+>zpp<{QKfI}W^R*br|w*~;VGYyas$0f^;-xSES=MlhQ0f%ng} z_g3iqU@|>MEWY(O9ep3acI9A-l=b$xI;ndEQ<-yDGN(>_iHU7kmtn$w<8s<+;p*Ay z6|J!+^D|!Vr3^_gCWH^$9<=G#kMn3B-zQ+!ZoX>K78?;t_ZC;t0~>7>#zXcFE$c3| zy=rgsY~GErrovohJF8=@B2`iKJjIE^Iu~cR?B4ds<n>J74G;^K_&U5{| zv;FjF`j4>W?R!YO58S&o$F+>cHp;>@wU;#o4M|NrYUvNe&`$-5YOS}Ru~u&~>2iEE z%%AY)?bg?yHSkM2d!~D~pVMdFUk(2^G+!C4>bKkYHQdz*$YZ}O4Qq>QJbAO37F*?a zTaReA*UFFku6A$;5bM_36}2$q!#SI^DGBGc49UJe`r1NiEG}YxvQ72&=7DTB=N0Bf zby$R*n_j`mJNLZv1`mpVs{C$f>zl78&!g#*(Jgy*_i^V>wzHAEQ}a`Ks|pmFOfR?< znLg4?U-yCa{c_g^SCyA-+pn(U$rQ+l_6dE`X}R~}>BslfgB;~2M|7w5#-{PzBVTiw zEWdQxu()8I#w12muzsLTw>?d2dP)Y`Kq#*h8hdw_L)X`xlO1z^EHIIZAwTw zQeGJy9O?Z2nDS-~N7L4+gjT7Ed4-+Ci{pxcIyZ3<>N1J*<%eX&Lllyk9Ha!-0(ah zyuW+obX)hm+fy&n$D%$idS)b?tMKLqAG>@LGjjw=y5m36QiF>ttC|Kkd>)kQoL|$S zcKz(d!l%@AGDa6Zl6`8zRWkFkbrVUG*N*q|le{UsaDwMUIW-+*X;qtDBL| zemJ8A4y2x$o zy6V1J(&V9Aq(^z3ujh?wZ%|M_G*M--|IiM0Sfzfs`*G=_jHb)R#8p<4hac!H&v2kz zeC54jp^dtNve^bXTN7tVkN&7wHXoLweeukO3lYRRdkYYaoe2(O89CEMvvvk z3wO*VC@YIC)m4f-drc|gfqmn-o7fdorz5^~$CW61PRN}RiD;@2AXQ$uQ@z&gut>(l znxf!^A9m{8k@8RUlc71Um+f6QpmI)s=M|QAj~m-bNUmv(SHO6e%y+r)bx!2RrG3wc zkE-#`87?*p&fW=}%dhc9$8KWjn`ZL?5dx@o<^^{R}aWhv8zjWKIC20xCHdgrP7 z^$L49+Pqc#&0X7f39GJ_w?!8gww!)QI9C*Q^q#jB<7~!{{4m?G;w05PXUi)@|7(5{ zt??f%$D$b2ab;Yet{ymix_< z9!F?dx$a7V(T?BA8-)q8?&GHiq;?D^EbPTA@!6;A?pvOTJ-NI{d9Q#`C;O6k6`Q~8 zf#VNK=9@J&0m`N&@3Mk7=$$+7aDQmu$-)bkRztP(%+fs!*nIu=lk+WHb4(qG211mJ zqC3*Pg^ggKRZ!A6q^;;j`$iU%Kc>PbFghlf7Oyef8r zS3qh*U5atb8V#aTnfF6svEWO4HbuHTU=R1|#Jp=2Csh|d8+}P$*m20`+Q9WQBW~$w z`&aF;5huU6x%k32W4_ejp(gPZ_XC;&yR~AjW@b`%EEFER-f(Y@^Enxd&RfO4s-o@& zR=MVGgLi6kLX697s(i#^)aNLBe@-dS8(#W)Q`th@Uc%f3Oyv*iy${zG2_xn)TST0E z*j1MmBu%g1*kfKNAGs+x#Y%O`Hd%|OYViGnwO91E>bG4xvUJDp*_&PbcQ34L^~RL+ zy%yZplQ3L*tU_TiL)>yk>|6Lglfexq6-}k?m|9I0ia%BY%;vD43a6K{Mo zG#WhG?5Zz){P-ib<0qu^{7Rt%hNnuiH0vo>rIsI2)zwGGBD_$qx&Qdh%8C^(%-)fukw+ z8+V>6T+i5drKm*t{^s2=+Uf1`SxQgVCg_~r*{`{Ga%al}tszpClPr7sQ@@Jw(@KA9 z<-wh6cd_l`L*Lq;-Zd_hsraJo!i;TKEFyCEn6sR^&VxwYz}fyMO-m9j%Kda>oynr$0C61pPsM>>_VT=WYlyUH35dXumxmhDw+S3bO*t^~Qa9M*^T({-Lw!kmO$~^j%X8lY%m7U^= z-hNTN5}7L*VG7HG?9F@%N|!V}Ji^l$oG25zxVc~a1X~c%itIb*Y+k(oY=!eWyTI?; zm8XOrz9A+(nYV4B+FnBgp^*<$xbR4Q&+@N!&+ON#P82ak*(VoF>X!%1!$vOnzrx-; zo~rNrA9s!+$ymmF$Sm`8%~UFxGtXrnN}0N3jBBn)BDyLKBDXXc5*5uEDkYjIluSuQ zR6_Dy`|9<2_5S?+`2Ep7d#|(iK5MVF_S&0=p3m`-7C|TT{NzIHjt|$JsZ} z$Z2&v3EJkpy;S{L!kp#y@=-IVimN9-e*VB)p;7zf#qz8dyF@xOd5Rw}jnO|8om}0p zxmt!&DK(SxAb6`m_{uHO&vklrq*^tOB<{^qXmjH}+q3J|%*Mja0_~TAU)x7fHG)<9 zXZdO_1RbX3|13Lbwx5+FTHtZ5G|4YEAc*Nvg;RijdaLxi?zrSh?s1>? zH-?=%|MT5G*}!ef%-8ry`HW1*9cob-kj|reNcc2opU`yYJ1G5_$`iQdR5kTCMf*9~ z&y_`bJmdPBe6@g5BJ+br$s1t{OFQbV;JPBKo5qH0b%Ry+*s^Pi-yJDt=hOUY{w00o z8P(f&@0$tIZ|FHaJ7&8*+c?aT?X8lJRE-wp);>X|Lx;SoW4-qmmL7TgtXY{w;ZDv( zl9t+6p>UnKPXU9;Ic{f66AhSG<+_~OttDHgtaD`Y=*|iQ0k#I3Y|5+5JIfzS*Xd9P zP3mmd&3Oup+c+Fs=o#IAaEIzfV}s$l5ij34ny-2_T4R?qHui&-IIFEQ z**#zNRB4nav0-qAbitrKZ`T3*x2(dwJ-kDO)2$sT-;ww3> zMrjW23BTN&F+fhOGMp)`Td{q^HC^ew15_iyO9~?qJ8TxN5CQQd-Sc`0y+$H_WrM3c z=|$RePG8p5oPL~@Gv{!kUsmafh1O}gx~)Fnd}18lX2#sskDqyZeI~{wm&-MAQYNzX z#G8{QHs@ZiisxwTPI&kEfpf5XL`%u_AHJq+GDdzb>O9fcHnnF2pS#`hE&5Bc+^zeVBMi?jA^AqMHsp1;1fZROe6q&T1Uhu`t~QZ|3yoAG?! zVV-4TzVmkz9jll2?aB*jY}{UJCBy#Xgsjk#!m*0Oxt)zvHJN>n{OkODY7l}2kg<9U|e^@hx>EIxei z8=JGqwT+e2&RG*3A0Ov?MJlsWf>vnv2D{ZJ@@+NA2-|k_{sUW$@yiM-?20N*JA*m1 z{4X9LR`(?u7!GS)k(x8vG)Qd8NZq_3obvwQREqj4_`iq}a}UmUJ)78{zvD*h@e`8T z*Ay$IsYjl^-lQ_SC2Fmwr_RfCKF2R#H{CI_9(}dtidp^6iHVUFL8SHGiBwZYUla0-d2v+hgtG^n$_O+oLnvaP&nrT&(TMhV%9d@doM1P zarqZ(@f>?+qJm73$L{0r4Bpy=viF_5RNivey1)5$&Jt>(7S}T(l{-Um2dwQ}_i|iT{_4K|Sr9Y7g?q=osJnZ=9r{ph z`lRd9Q@19Z_o}mH)Y10K)0;kiwTOFoBp|DWCd@m*(NU?v);j zjIH&x%3~Jw1ypD4$(g%*+^>8QJM7=Vo_RR;{+^B9@}|}D4dH%4fyi!V$hM!#^;OeD-LyY2>q=%NyGEst2@= z&~{t1bxOp{%^mu^!$dD}=H%4R#Pl-S1DVzUVaqkiN0vrbp0koJ{T|lXrJ9~wz|_6$ z-ZDWwtt;e#<+J+i3h{>})fbkBekJ2_VTYb@KR6kQuX2%*zgSjW`BH2w?~vVR;X$v^ z!%0N1F?Y^6!}gGd*yGyelU-ip{kcBY+Xt?)oJwKpE#lFta+}Dd2F)xSE#5(VI<+9O zU9R@X>)5Y4?e!u}q-5o^9Wz=6Zr5KY5mS_eD6yaEjB8q zs`gHQXi(1aF4$M^tIBg-{u_I4da6B}>l^VZ;DD^&Ld#R{cW6=ycLJ25NQr#Lq3PH*aIYY_SJ zlvT`Mv}I0Ke8qw5&sE~8-%{f>R%`7081-%Tj=B>0?KRS&4PzF?sXMzgEl3B)I<6`| zw)^B5EFURwxw#?b^U`kxA~$_lLXYY{`hs(;d!tPb8W-eyA34^!#U?jgaZ-1S>kHB8 zIK%G&M}uQNVB&KjHT)4WdY75SBbjnb);_)eVCAf%ATM*jp7-~Km$l|YCe*_-CnWQ- zXS#fP9?p@@`Fi-Pbn4^BYu65z*9xpB>n-W^@)kJ}BGnM2P@JV8sWjC*zwxvrORAGe zDHlgQY5xi@%XS5wAx}zke9Ns;E|#{g&<&lG zLxa?{9ObbG2gw|KLr3TLmKsS=+fH*)rfs5~i*a>~LN$n7neD!;z`|zu?o;qjl37*y5 z$F}!T*v7M|H;f2zqt#NzmsC7My$vzVH&DJ`hZf ze5{-;p#Q9_x%gDYbgQZLs!E?7*{7B+ohfjeS}7o|5pNvc_1Igr*D+e~gnaB5V;yc! zn|;fLjUVq$c~gZqn+N0TpeP6VB&!x9cjtBT| z8e&R$(x9DN7_nSthn~E&&L_8z22>*cqb{nYL>BO|JBY1u|#-tF}J#I|sI2 ze_12N(|@Zlo_E8sy506$-JIl`{H32-_e57lx!;xGInH#&M_YgQ^o8s!qqGglob?-5 z`aX7{p7{p<`z?ZUnoyzo#v?KMI=f zsVuxzo5po|p*MljIAMdoIcRG$u22=Ln#%0x)1pP)(5IUdqCeWbLPXS9`yqMNl*|tW z!=T@DcS4NY4&RXKeRiQ?y9)0!?fT%H5YebW(M<;q=NNtBaHZ0oH}xL9$oo*~yYa=G z-j+PykFKxc*jC&_=HO~Jn)Q% z^YoF~qS0UerZKYywCEz+t9RFp1hH0xe+qsR$s1fBFDA<4DL^)R&$G5!dBYi&`xjk) zY8xB>;Ou!mYGGF>I_wozqYz4}G0v7AFT6Hq+G8-=7st6bkUc*pZ}6MJ$VGBt!SFjz zU*(+Q=q(0Y=cs{)pA=CChUG4uO}*{&-TY^YL`}&@QOnCkZ&}1o=$z{9_^x|U&H6pn zV5D*@7Whk9+vZF`C4RB5Q^3x`?80;5+};+t(SXkPZ0v5P(gPj3U`HOWrTxTwBIk?z zeaxr3`Q@wHQod~ovCE%KJa>Z9E;%u-ee`VR`4qe7Ih*s}P0d1q(S~Dzn~D3Z>R}GCYfllf(QQNPj;4G$T+5qxX74F^ zmixw~ne4UE)Jj%B`^QGDFYqkS7n2Ko!tBY*ouSA>litMUIJ~8w(@wXV_mrjG;|n`V zt?M(Ga(ErW1F1*aU0h?YsjN7f?r1z#WiI&5xm-V-XGaRf+q`J>{8y(7YXx8C8}2FF zam2^%>eQ!~UB5rbx=&m_x{0JvoMJoJApF7i+sdZ-SFe{*{XV)|RJM0igzcGTdA0Pi z=)1zcGkXWcxu>Lu9Tnf)^=5C3JaN;1mTUTmsonQ0tM_)eKT-db^^R#_V*EFA+voBw z;WQV!rV4IJN7w0|=aDU!`li0W-1w8JqkB#1nhWEhdE=jB()xYyvX3T~chzIoPA+77 z@#@`)!K7&O*9z8B0~=wckMy46M*`-eczY;ZmU!MRYIM;_QO+u8rO@L zFO(Fo1$Ho%R(u|$$6WEveQJ|1@2kA(5X(c-bpMN=nu9wZa%b;59O1GhiX1Vuym9Ka z({{gjaiO49MQ-``>uT!Ot~tcyO1s*ma&764Z99as-+CK|u&EDG1%HKR9rI6qI5grp zCDJSx3T*rOYQA|Ky5qp_cW`D+`=`3*{s)p3!9(2X)AU{YOs96z zE5v71&CVAmopdlWTFyRg*K1(rs`|`Z%XBk$;d7M$_=6(WS<88UEkzMXB+&73R264_>7=P6mt3zc~Wblf;l>i*&E?A;|{!q295&q^BE=A57olsQ!Id|2?& z((3oi-FmT>dJ^=xt)_A{1#buY`5z5v`09(Wt}6LFmFz0FHS=*~p+N023zpI`;YzK9rC(Zo8GR%dTcKe^Vxh&!S#dr6q*mx`)v-X4MB2Q&lq z!@=tFDc_Ilrbl&U9%reTDC{2e-m-2_MM~54y#{yV22MOX(=kZs8)QPMOX$6nlMG3D2gErU%qF9 zxu)AniHDjA!~^q*@}-8eIRj2dP2Tj-D?D!wP~2`@ulicbt^NF4>q1Px+Tm=QkJO?A zZ#Dah`A02lm(BLHIj`$FU_Xw1d}Npy;Wu?rE&gyq%Ce>(*PP!9N;pL9bI+3aCXu3B z#?$+uE@)rqPNvT4vul-;C)`rG)k{5w4vamL8=%c>>Imp=xpymhd3nOSwyzbACzL&; zt7g{T?9Y$qk6gjw;i)a|z4gU{caM?#=68Y>O*7MlD$(v`#-wdYF9|!sG%smnkx>axsrmKz+G$H={;(hF#5sk|;qoDnnS~Y6 zQFD44?nYiY#q^dnc8Nao}r+e zJAx{n6b%RHZB5E}mRYq@tv$l}RFU64q3XS2w!v-p{cBc^Fbx_$D5 z%U51`HgxAPXRz@Xz0hr+=%h8V2b-6*NwSXol)NnwuO*WoikFj9udk2bmreZ zrtEX%^v(AN#Xl5B-z!&&`#i`Kefr~Ej$qNZLNA$3Go_y`j9&)_E`;9NRT{#7bxcxX z|F+L2)Ie!@Ve97jpe0>458Px2TJDv}4t>kGy?rb+L{wnQW);q50x_W5Mh-BJ?J@;egt(I>~TgAB=>zBW*`ds>@ zktB1j7Axaw+4jQG%i+}xjiXd4x8J$jC$>hXB^feP#l(HJOczqwl}x{>*a&u22;aXn zvMWgZVvJ>ww#%1kZ*BpeV!zj>^$H_OBXP1DmE)Lu8x7KZ!+6K-ey(%;F|Ay-piz1t z$*SU7pL6SZ(w5j86Y8&TsQSbSRNPEw#90=dyMks*%-nz9`EbrPYEHO) zO;{*eZ8Aw+Szy?4M`xvP=;f>G(cL@Nkgu^>`FwHPKhSQ?bJ^eFQM+C7;F%{1-qB^{ zd+$HG`Ez9N-t8F+CQ__Bw>}tOh^(6Wd0=JICWpqPkL+>yvPSjvE=v1YquiH)>g^>BXPc#qX6kkv?>~E@&RH#wWXH=RIkL5d{PUMAgWd?RVdz!HA z^K?Evc~$qp!Eeu}rixGLJJ{7YYl&^U{$}dyCk4Vpcf(9a=p`K)s=ny$sh{?vdAk!V z-+kY2p5FZCMU!817zlRZ1vSg?;S9Y+_r5Y0FNY;{A<1z22l31nAES=_o2qiDI&+~(-4BeJZ#+=2bzlwd-#qxlsL5T3_BKmm;A2nG zT=7b=c=fam#`n~lmu1Dt{E*Yd*ZgK1&gJ~UnPd@EJ3Cggx5(Z`^Wqh|PmyQqZ^xXZ zu(ZY6FU8({kY7ui(=`0lxY0k$;Y#w_T}5xV|QzX50&tXBnA`J{QMC*o|9pD(#u2&B$k{hXRlc(YYvmC zbiOxIDc}{{pk)$pVj_yjAsxsTyhMH2Ze4OSXY|%JVjYXTmY#8vsA$$`L0$Z$rAZY1w{XNm61V zAtIQ%G>hFWa3J$Z$0ytJ3K4$)y?(@M4yKSNHdf)kjd&=H6$rolSkUHny zO-lOtuW$t(7CGi~)tfyhRKf~570qY(sUj)O8KJ!O_OcwR`C?wGC z_lUG~@8BMNUUiAF(J7hnU6yMnw}zg2ZV|pPV%SymI$OT~)nWQl@o1U-d-jmZdOUTE z`3fJeG3L0jXUoN)-ZcHVW%pK=y5HiG+8O9zNqyI#9k^{YZ%bLW)73L{@e9vm$l9Z% zp?oDHjN+tnOv z{Q;+eh<7I0Kb6v!TyBp2UY8blsiHyfoM!1*;+h88``j9dbLDE^PklJ0JuGnGuwcSt zj+Wps>alNht2%msW0Fuh+dcb8!Gp4m>~`XlJGez0wz|wmT{iADE{^NHkbhuK<3_5I zq{Zsplg+wDXIPS*heOiaTSM|Zm+n;W`JDEU^=yTx!gIs11HV>y)i~V>oU_bbD7bYC;WdWNuGz zhkJBW6&}s)#nr2VHiwGE7>Y0bp4u|eN$ux#kv#fpob868q>*y>2br-(*CyX zeV)akp832~l^1tvjijO_BkfL#3S8_D4Yt*LIrpoLYwgpd)GEI77sqA7c~1W{O}*YT zFF$f{;52XSsWgEr2TnJ|iCrDk?OPjn#|od288i2({@6Fnedo(T?U<^c-%3CErYh~f z6xX41H;MTRS4pDD`K-zHK2Ig{t(QdSP1Z_%`x$KYa)%}#wJV;(ZuRJ%0H)>R61`7q zbVVm4+28rBo%D^oE4EQ7QYKBYH(-0>PN6={dnbujR#W1hMw*t@;Md9qewIg~60v4V ze$M2cPD`m3J~9bErW(QAA~$nP2d3(6-%E3oZpFIBv>%V@sS-*IcemnB^$?{UsB5G8n=ZwR6e!FVr$vdNm zl{W~5T~?h|$evhbZPcQHf0uuFC}e}OU1o5j@{W`3{(PYqBno7y>%B8|PBVL~6*&J) zUNIo*^@eBr@7$nU?7Jgr&9u8%QIs{adLw;SLq_O{hwAib;i2T;zXg5{iF*B>*=E%{ zLN79`^KX6kIQCPbpjN81$EU_3mF+_#hc7;P@Pcjdr(x4{=hL{(HE9M~oQMb@Iult+JdT(Cp^WWU0vv2#`e{aW5bJH&RhUxJg z8at-4S@i8mwoX2-*ZWEzZcsVY_(EuJ+qN@H9P(fKtIV$|_K_KQyKZF)Df3j<%^u~P z=mcQQ~G0?)%8F?+%*0?6COzHHNhsvS*Ehz%r zMdh-m-o^;uOKC3`TIPCR2z z0Rk21)DTuKni8S*0yPV8$^ptM!gK|cR1Jisssp&d8Vqn(gpjL;xJCdyR1>k)F+Q#X zP?$OZ@e5jL0XZEY00T}g=m2n}TIvX@6+mJU&MgPx_yEolm=WRfY5*Ru4uZI4tO$kz ziX?m$;lBd63shhvLK_Fp5=2e{DuRV{5f6k8AP?*6!57J}ETEIo;?P9+u81*ClM3e3 zLx{#eD}j178SUG|))&!7D#6Y8u=w&ouR-oblI9CRV zh&mA5!AN~zEac4pFw9b*aSR(FxLL<9kduNLh$p|CJU}xzlM;Z? z&ctdsY;HIcO7IwUCcMDa<}QRHZO#Q8m*_$`Ak5+p7eWEixczd04w|&hu7m=hKxVi? zON%z>N+@wH;*&Dg4Wujo%q>c@af7Zpv{*Mnne#8707)5;7twsI2o|J+=pPf%M%@Tm z(16FCP+YMX$Z8wI{KFC3FUr(3Yjk6~4_FY8 z7Pu2eaIhZrAanuGveN@%i%eU;nenLEOqkGqmI7F9qc^dF7U)Tc=mTLT@-dtRi##Ui zoMf~Yv2y{!ei5=0F^q!zAC?eJpv~|mRxZwg7(vk(M3BaKE+T0&XwDdIIGq?z$Sv^v zXAVFi{&x;FL_Egmvj|rVDUi_@o{LDe=onzUhrTk^Dj|L;MjMVDc>aTXs||!Ou#CkyMAYEmV03~+$$0*;jW#4nM*E+S08~1f zgcc?OJQe>bA)>oTBm*bc;$jd+JQ@i2u#3+>+CgSyOt8q=1Sru}-UPsqqJ;?aq1y?X z{m)By%dpZsTB|l~nK#hY{V@q7Xf!(n^x12PUZtKa}&LjA-3z&KP< zS9GMJ=bxPU?>As;dJ~G=i?%~2D|GUto!Ub~^}uTH|&?oq@;f zAR%uN6&MPv-wLqA#rXLn(L|EP>HbK_NTNk+$RI>%Yj-dbto;r!y*4dk2SKG>4I%hd z7_;drz-hPWy~Uw=AYf2}{qsG2B+#R&hY;4VClf;m16n~SL85hp5F24D3560#w4g9T zgmx(uY^FzR3?(7}@%3vsWWM!bgbZvG`1A@|d>El5vS<)_Ah!=xwcB870Ak!O?+m0K@*c z5lINamtI7|NAx+=OoJ(rQNXV7AEyNHxMjz(j{RcS#riL34 zBn8wNco^Wm0TsBPB4$d)%>{8d14A+-0k|^iLiPY^10(@L2DsRxvyicy&}m9j+(q!w zT*KkC)1W1X!+A=k)rAv@wBcQZ7|n1Ocub$>vx{&7b26Oxa2Fv#tBrreF8jLP~f^t!!4Sm2TfAxW|L1c!1Ng$LFXY*eY6w5%fBN+#w_&ElS zBYG?$vP6Rs(K{-ZKRc13xh80iI353*0W?Q^jekkdCj4^({Luty{wJ;dkszh@8Ea7j zzy0$?QHwu*V0?@W3qbsa|4mJX1t5aLza+2#_;KuS37A`(42%iLOs`RqA5Vzbl0kDI z==i%MXpRU}{+bLlM;s`BNkDT%Z1R@`GzZQSnrl4VE_D;Y;p#~6Cn^|vYazYi;{4wr zEu{B9X>YL)(i^tiUwx3?e-6w=BWZ!kKnJq863GN`cy|IMCc6ZJ?SJP2wGm6lUvq)l zh#liE38)SCTUt&6VNFZe4bw~}5c_D_iLehCt}fV3h=M1ocM}3ae=HCg=pPjv9;T>@ zLJThawBFssDqQ1_LPR2DO~$scNF>T>6-i)emL#a9PziC`gCwG!=9&z&i>IkX9nB(z z5T->W6WuhURQRH23ZaSkh*AgvHWa165rFEn+bIMq>{EsuMW!C34++w+Z3eN5CY=g; z%cMa^`&0*i$+CB5qSl2B3hF+Z8(GAX8ZHCQyGLf8mo_#QlTx) zBsc|40wbd$8OacaEgtKT2WcFcgp?G+d?38g=1_B$*bEt4LIZ>j;o(t??8+Ey5f&5oe5V$g))*ctYC(kDyUFJGO!0uF-1)^RH!2v zmU<;hG0SUoyOgM5c36_Xc)}y8_D)f*HX+&9v zn);|PKr-}oiZN;$qQVHt81bUqCp>YA5^dXlB7`9MBOwziAb@cS1T{{9$ipc(ZG1mr!$(3sBSBo?BtIlx zMzhHwgk&gS0H)g@!N@R}>9#MnqGcPS#+XnO6Ufk|?aLu9Eun+oFuE~q%>lxRN=5@< zAs879gOTMi{<_^T3L_)qQy6c-o){UdgOR~H7}*#VU|x)@i}80?Gdd#=P#{<^GT0a+ zgBviiA;v$o_!bO}k->!+8H|XL*JA?wZy3W&a4Bd*RAFRGR6s0YWQbafY=#Mt*%*DS zXyjZ%3YiZhgZVJB4JKf8fbq%&HQi|8@Se(8IP71H&WxlUm;n5@D^!EkV{}xbW2nH9 zrm!}Qjx?pO!h~Af7;hniV04I7j1HlR(IFNwI>a7EXY3Q$9T**=9HT>&V{`}$j1Do2 z(UFDdS`e-n9od8qKEvoNm_VI2V^;7UMqi2whzX3&jtSh^z<9xdnh-1)oeLA_Sd?+2 zCb$}-L&RWolo{xJn85dCj1hvU$w(-mDo%n-fYBidU~~~oP?)O9=mUuYql=;?GSn1D zO(|5Y!~``TF}g^grX=bD`w62%_+oTvOwhKR@m3B^4aozeqi~=@Oks4`uNYkc6HHsd z=mL8Hqr;ZK=&<=Px&|g#dya{Us*DLyw9t{q>3W#p$B&F5u=_AN3JUrfOo%kikfB(j z!)C|mkOeS0iVQl65jqM7`g%;r`X!?uy5$Ler}mZy^kD61W^E z@u33kX%aXUCvjuKd}fR;C{s|ttT+jrkCVWiI7tW-R%vCt6+ul%`8Y`$72tfFggi+S z#i+vM2u2@pIZlFT#Ys|_u)PdJ1`&*tP@W>e7R5;r#yCk16ZXE(c&mV#;7^>Sj0tb8 zV#pxAaS|h`L44yRbxb(cmGJ^%9;cw)PGTfea5hfT!i4uwPcdFY1mGlw{)|50S)7C% zNP?Y>lOPOm68IP=8DPSvq!=S0Qg9L@2}3~OBx6)Sq~N6Wr~sGaBoj>dsxzZAga%H6 zOo5XiTyPRN3@2G*!nfK{{NW@B3!DTIfs??!I0>AOlfdsd2^@`67%tg@3HS0c=JG-^ z3?ozxP?G{XAEzwEL}cL*fEsL1oPzvEflY~1ki#fwu@p!SI7I*z$aNGURDiFjI0bo! zg7yyu(gsdJzN3g@B07vHp@?B3ft`%0S7IXR*BLc&Or&&}QIo(#?&&gWk{H9oFjfj9 z$)i+`Q=~Ccq6VWTgOU0fp+-S2ra;EQDafi6732;IGA4zwFbHFu;)4obR6s!E6v%}5 qX(o&azfY?zAoik#ofgH25H*mcouh3mBqTXBbagpoWjETIa{ND - - + + diff --git a/docs/html/userhtml.css b/docs/html/userhtml.css index 57ef3dd4..c0dacb61 100644 --- a/docs/html/userhtml.css +++ b/docs/html/userhtml.css @@ -20,6 +20,7 @@ .cmr-8{font-size:72%;} .cmmi-12{font-size:109%;font-style: italic;} .cmmi-8{font-size:72%;font-style: italic;} +.cmsy-10x-x-120{font-size:109%;} .cmsy-8{font-size:72%;} .tctt-1200{font-size:109%;font-family: monospace,monospace;} .cmmi-10x-x-109{font-style: italic;} @@ -63,7 +64,7 @@ div.obeylines-v p { margin-top:0; margin-bottom:0; } td.displaylines {text-align:center; white-space:nowrap;} .centerline {text-align:center;} .rightline {text-align:right;} -pre.verbatim {font-family: monospace,monospace; text-align:left; clear:both; } +div.verbatim {font-family: monospace,monospace; white-space: nowrap; text-align:left; clear:both; } .fbox {padding-left:3.0pt; padding-right:3.0pt; text-indent:0pt; border:solid black 0.4pt; } div.fbox {display:table} div.center div.fbox {text-align:center; clear:both; padding-left:3.0pt; padding-right:3.0pt; text-indent:0pt; border:solid black 0.4pt; } @@ -149,6 +150,7 @@ div.abstract {width:100%;} .shadowbox { padding-left:3pt; padding-right:3pt; border:solid thin; border-right:solid thick; border-bottom:solid thick; } .doublebox { padding-left:3pt; padding-right:3pt; border-style:double; border:solid thick; } .rotatebox{display: inline-block;} +.figure img.graphics {margin-left:10%;} .lstlisting .label{margin-right:0.5em; } div.lstlisting{font-family: monospace,monospace; white-space: nowrap; margin-top:0.5em; margin-bottom:0.5em; } div.lstinputlisting{ font-family: monospace,monospace; white-space: nowrap; } @@ -185,9 +187,9 @@ div.lstinputlisting{ font-family: monospace,monospace; white-space: nowrap; } #TBL-5{border-collapse:collapse;} #TBL-5 colgroup{border-left: 1px solid black;border-right:1px solid black;} #TBL-5{border-collapse:collapse;} -td#TBL-5-10-5{border-left:solid black 0.4pt;border-right:solid black 0.4pt;} -td#TBL-5-11-5{border-left:solid black 0.4pt;border-right:solid black 0.4pt;} -td#TBL-5-12-5{border-left:solid black 0.4pt;border-right:solid black 0.4pt;} +td#TBL-5-10-1{border-left:solid black 0.4pt;border-right:solid black 0.4pt;} +td#TBL-5-11-1{border-left:solid black 0.4pt;border-right:solid black 0.4pt;} +td#TBL-5-12-1{border-left:solid black 0.4pt;border-right:solid black 0.4pt;} #TBL-6 colgroup{border-left: 1px solid black;border-right:1px solid black;} #TBL-6{border-collapse:collapse;} #TBL-6 colgroup{border-left: 1px solid black;border-right:1px solid black;} @@ -200,8 +202,8 @@ td#TBL-5-12-5{border-left:solid black 0.4pt;border-right:solid black 0.4pt;} #TBL-6{border-collapse:collapse;} #TBL-6 colgroup{border-left: 1px solid black;border-right:1px solid black;} #TBL-6{border-collapse:collapse;} -td#TBL-6-5-5{border-left:solid black 0.4pt;border-right:solid black 0.4pt;} -td#TBL-6-6-5{border-left:solid black 0.4pt;border-right:solid black 0.4pt;} +td#TBL-6-5-1{border-left:solid black 0.4pt;border-right:solid black 0.4pt;} +td#TBL-6-6-1{border-left:solid black 0.4pt;border-right:solid black 0.4pt;} #TBL-7 colgroup{border-left: 1px solid black;border-right:1px solid black;} #TBL-7{border-collapse:collapse;} #TBL-7 colgroup{border-left: 1px solid black;border-right:1px solid black;} @@ -214,11 +216,11 @@ td#TBL-6-6-5{border-left:solid black 0.4pt;border-right:solid black 0.4pt;} #TBL-7{border-collapse:collapse;} #TBL-7 colgroup{border-left: 1px solid black;border-right:1px solid black;} #TBL-7{border-collapse:collapse;} -td#TBL-7-5-5{border-left:solid black 0.4pt;border-right:solid black 0.4pt;} -td#TBL-7-6-5{border-left:solid black 0.4pt;border-right:solid black 0.4pt;} -td#TBL-7-7-5{border-left:solid black 0.4pt;border-right:solid black 0.4pt;} -td#TBL-7-12-5{border-left:solid black 0.4pt;border-right:solid black 0.4pt;} -td#TBL-7-13-5{border-left:solid black 0.4pt;border-right:solid black 0.4pt;} +td#TBL-7-5-1{border-left:solid black 0.4pt;border-right:solid black 0.4pt;} +td#TBL-7-6-1{border-left:solid black 0.4pt;border-right:solid black 0.4pt;} +td#TBL-7-7-1{border-left:solid black 0.4pt;border-right:solid black 0.4pt;} +td#TBL-7-12-1{border-left:solid black 0.4pt;border-right:solid black 0.4pt;} +td#TBL-7-13-1{border-left:solid black 0.4pt;border-right:solid black 0.4pt;} #TBL-8 colgroup{border-left: 1px solid black;border-right:1px solid black;} #TBL-8{border-collapse:collapse;} #TBL-8 colgroup{border-left: 1px solid black;border-right:1px solid black;} diff --git a/docs/html/userhtml.html b/docs/html/userhtml.html index 2d1ac7b7..43741db1 100644 --- a/docs/html/userhtml.html +++ b/docs/html/userhtml.html @@ -3,8 +3,8 @@ - - + + diff --git a/docs/html/userhtmlli1.html b/docs/html/userhtmlli1.html index fc40c5d3..c8ce8f48 100644 --- a/docs/html/userhtmlli1.html +++ b/docs/html/userhtmlli1.html @@ -3,8 +3,8 @@ Abstract - - + + diff --git a/docs/html/userhtmlli2.html b/docs/html/userhtmlli2.html index 07c2f5d9..b52fe6ee 100644 --- a/docs/html/userhtmlli2.html +++ b/docs/html/userhtmlli2.html @@ -3,8 +3,8 @@ Contents - - + + diff --git a/docs/html/userhtmlli3.html b/docs/html/userhtmlli3.html index e358a8db..a102ce70 100644 --- a/docs/html/userhtmlli3.html +++ b/docs/html/userhtmlli3.html @@ -3,8 +3,8 @@ Contributors - - + + diff --git a/docs/html/userhtmlli4.html b/docs/html/userhtmlli4.html index 2a339025..dba92584 100644 --- a/docs/html/userhtmlli4.html +++ b/docs/html/userhtmlli4.html @@ -3,8 +3,8 @@ Citing AMG4PSBLAS - - + + @@ -32,23 +32,24 @@ class="cmr-12">When use the library, please cite the following: - 
+   

 @article{DDF2021,
-       author = {D’Ambra, Pasqua and Durastante, Fabio and Filippone, Salvatore},
-       title = {{{AMG Preconditioners for Linear Solvers towards Extreme Scale}},
-       journal = {arXiv e-preprints},
-       eprint = {2006.16147v3},
-       archivePrefix = {arXiv},
-       year={2021}
-     }
-@Misc{psctoolkit-web-page,
-       author = {D’Ambra, Pasqua and Durastante, Fabio and Filippone, Salvatore},
-       title =  {{PSCToolkit} {W}eb page},
-       url =    {https://psctoolkit.github.io/},
-       howpublished = {\url{https://psctoolkit.github.io/}},
-       year = {2021}
-     }
-
+ 
       author = {D’Ambra, Pasqua and Durastante, Fabio and Filippone, Salvatore}, + 
       title = {{{AMG Preconditioners for Linear Solvers towards Extreme Scale}}, + 
       journal = {arXiv e-preprints}, + 
       eprint = {2006.16147v3}, + 
       archivePrefix = {arXiv}, + 
       year={2021} + 
     } + 
+ 
@Misc{psctoolkit-web-page, + 
       author = {D’Ambra, Pasqua and Durastante, Fabio and Filippone, Salvatore}, + 
       title =  {{PSCToolkit} {W}eb page}, + 
       url =    {https://psctoolkit.github.io/}, + 
       howpublished = {\url{https://psctoolkit.github.io/}}, + 
       year = {2021} + 
     } +

diff --git a/docs/html/userhtmlli5.html b/docs/html/userhtmlli5.html index b9362d17..a6336dcb 100644 --- a/docs/html/userhtmlli5.html +++ b/docs/html/userhtmlli5.html @@ -3,8 +3,8 @@ References - - + + diff --git a/docs/html/userhtmlse1.html b/docs/html/userhtmlse1.html index f7af8224..2d0fbf00 100644 --- a/docs/html/userhtmlse1.html +++ b/docs/html/userhtmlse1.html @@ -3,8 +3,8 @@ General Overview - - + + diff --git a/docs/html/userhtmlse2.html b/docs/html/userhtmlse2.html index 8c0e20d1..474befde 100644 --- a/docs/html/userhtmlse2.html +++ b/docs/html/userhtmlse2.html @@ -3,8 +3,8 @@ Code Distribution - - + + diff --git a/docs/html/userhtmlse3.html b/docs/html/userhtmlse3.html index 3184abfd..c806825f 100644 --- a/docs/html/userhtmlse3.html +++ b/docs/html/userhtmlse3.html @@ -3,8 +3,8 @@ Configuring and Building AMG4PSBLAS - - + + diff --git a/docs/html/userhtmlse4.html b/docs/html/userhtmlse4.html index 6e7558c6..b932dc14 100644 --- a/docs/html/userhtmlse4.html +++ b/docs/html/userhtmlse4.html @@ -3,8 +3,8 @@ Getting Started - - + + diff --git a/docs/html/userhtmlse5.html b/docs/html/userhtmlse5.html index 8f6ba4e7..b4a1075b 100644 --- a/docs/html/userhtmlse5.html +++ b/docs/html/userhtmlse5.html @@ -3,8 +3,8 @@ User Interface - - + + diff --git a/docs/html/userhtmlse6.html b/docs/html/userhtmlse6.html index 5cbf8ab8..15f0ee45 100644 --- a/docs/html/userhtmlse6.html +++ b/docs/html/userhtmlse6.html @@ -3,8 +3,8 @@ Adding new smoother and solver objects to AMG4PSBLAS - - + + @@ -215,22 +215,24 @@ class="cmr-12">pass it as follows: - 

-  ! sparse matrix and preconditioner
-  type(psb_dspmat_type) :: a
-  type(amg_dprec_type)  :: prec
-  type(amg_d_tlu_solver_type) :: tlusv
-......
-  !
-  !  prepare the preconditioner: an ML with defaults, but with TLU solver at
-  !  intermediate levels. All other parameters are at default values.
-  !
-  call prec%init(’ML’,       info)
-  call prec%hierarchy_build(a,desc_a,info)
-  nlv = prec%get_nlevs()
-  call prec%set(tlusv,   info,ilev=1,ilmax=max(1,nlv-1))
-  call prec%smoothers_build(a,desc_a,info)
-
+
+ 
  ! sparse matrix and preconditioner + 
  type(psb_dspmat_type) :: a + 
  type(amg_dprec_type)  :: prec + 
  type(amg_d_tlu_solver_type) :: tlusv + 
+ 
...... + 
  ! + 
  !  prepare the preconditioner: an ML with defaults, but with TLU solver at + 
  !  intermediate levels. All other parameters are at default values. + 
  ! + 
  call prec%init(’ML’,       info) + 
  call prec%hierarchy_build(a,desc_a,info) + 
  nlv = prec%get_nlevs() + 
  call prec%set(tlusv,   info,ilev=1,ilmax=max(1,nlv-1)) + 
  call prec%smoothers_build(a,desc_a,info) + 
+

diff --git a/docs/html/userhtmlse7.html b/docs/html/userhtmlse7.html index 20a74c45..be3f3a40 100644 --- a/docs/html/userhtmlse7.html +++ b/docs/html/userhtmlse7.html @@ -3,8 +3,8 @@ Error Handling - - + + diff --git a/docs/html/userhtmlse8.html b/docs/html/userhtmlse8.html index 01bbaa41..bab484f0 100644 --- a/docs/html/userhtmlse8.html +++ b/docs/html/userhtmlse8.html @@ -3,8 +3,8 @@ License - - + + @@ -36,37 +36,42 @@ class="cmr-12">AMG4PSBLAS is freely distributable under the following copyright - 

-                           AMG4PSBLAS  version 1.0
-              Algebraic MultiGrid Preconditioners Package
-             based on PSBLAS (Parallel Sparse BLAS version 3.7)
-  (C) Copyright 2021
-  Pasqua D’Ambra         IAC-CNR, IT
-  Fabio Durastante       University of Pisa and IAC-CNR, IT
-  Salvatore Filippone    University of Rome Tor-Vergata and IAC-CNR, IT
-  Redistribution and use in source and binary forms, with or without
-  modification, are permitted provided that the following conditions
-  are met:
-    1. Redistributions of source code must retain the above copyright
-       notice, this list of conditions and the following disclaimer.
-    2. Redistributions in binary form must reproduce the above copyright
-       notice, this list of conditions, and the following disclaimer in the
-       documentation and/or other materials provided with the distribution.
-    3. The name of the MLD2P4 group or the names of its contributors may
-       not be used to endorse or promote products derived from this
-       software without specific written permission.
-  THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
-  ‘‘AS IS’’ AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
-  TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
-  PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE MLD2P4 GROUP OR ITS CONTRIBUTORS
-  BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
-  CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
-  SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
-  INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
-  CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
-  ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
-  POSSIBILITY OF SUCH DAMAGE.
-
+
+ 
                           AMG4PSBLAS  version 1.0 + 
              Algebraic MultiGrid Preconditioners Package + 
             based on PSBLAS (Parallel Sparse BLAS version 3.7) + 
+ 
  (C) Copyright 2021 + 
+ 
  Pasqua D’Ambra         IAC-CNR, IT + 
  Fabio Durastante       University of Pisa and IAC-CNR, IT + 
  Salvatore Filippone    University of Rome Tor-Vergata and IAC-CNR, IT + 
+ 
  Redistribution and use in source and binary forms, with or without + 
  modification, are permitted provided that the following conditions + 
  are met: + 
    1. Redistributions of source code must retain the above copyright + 
       notice, this list of conditions and the following disclaimer. + 
    2. Redistributions in binary form must reproduce the above copyright + 
       notice, this list of conditions, and the following disclaimer in the + 
       documentation and/or other materials provided with the distribution. + 
    3. The name of the MLD2P4 group or the names of its contributors may + 
       not be used to endorse or promote products derived from this + 
       software without specific written permission. + 
+ 
  THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS + 
  ‘‘AS IS’’ AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED + 
  TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR + 
  PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE MLD2P4 GROUP OR ITS CONTRIBUTORS + 
  BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR + 
  CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF + 
  SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS + 
  INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN + 
  CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) + 
  ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE + 
  POSSIBILITY OF SUCH DAMAGE. + 
+

@@ -78,37 +83,44 @@ class="cmr-12">abide by its terms: - 

-                           MLD2P4  version 2.2
-  MultiLevel Domain Decomposition Parallel Preconditioners Package
-             based on PSBLAS (Parallel Sparse BLAS version 3.5)
-  (C) Copyright 2008-2018
-      Salvatore Filippone
-      Pasqua D’Ambra
-      Daniela di Serafino
-  Redistribution and use in source and binary forms, with or without
-  modification, are permitted provided that the following conditions
-  are met:
-    1. Redistributions of source code must retain the above copyright
-       notice, this list of conditions and the following disclaimer.
-    2. Redistributions in binary form must reproduce the above copyright
-       notice, this list of conditions, and the following disclaimer in the
-       documentation and/or other materials provided with the distribution.
-    3. The name of the MLD2P4 group or the names of its contributors may
-       not be used to endorse or promote products derived from this
-       software without specific written permission.
-  THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
-  ‘‘AS IS’’ AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
-  TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
-  PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE MLD2P4 GROUP OR ITS CONTRIBUTORS
-  BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
-  CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
-  SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
-  INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
-  CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
-  ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
-  POSSIBILITY OF SUCH DAMAGE.
-
+
+ 
+ 
                           MLD2P4  version 2.2 + 
  MultiLevel Domain Decomposition Parallel Preconditioners Package + 
             based on PSBLAS (Parallel Sparse BLAS version 3.5) + 
+ 
  (C) Copyright 2008-2018 + 
+ 
      Salvatore Filippone + 
      Pasqua D’Ambra + 
      Daniela di Serafino + 
+ 
+ 
  Redistribution and use in source and binary forms, with or without + 
  modification, are permitted provided that the following conditions + 
  are met: + 
    1. Redistributions of source code must retain the above copyright + 
       notice, this list of conditions and the following disclaimer. + 
    2. Redistributions in binary form must reproduce the above copyright + 
       notice, this list of conditions, and the following disclaimer in the + 
       documentation and/or other materials provided with the distribution. + 
    3. The name of the MLD2P4 group or the names of its contributors may + 
       not be used to endorse or promote products derived from this + 
       software without specific written permission. + 
+ 
  THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS + 
  ‘‘AS IS’’ AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED + 
  TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR + 
  PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE MLD2P4 GROUP OR ITS CONTRIBUTORS + 
  BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR + 
  CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF + 
  SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS + 
  INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN + 
  CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) + 
  ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE + 
  POSSIBILITY OF SUCH DAMAGE. + 
+

diff --git a/docs/html/userhtmlse9.html b/docs/html/userhtmlse9.html index 5a1864d1..21676acb 100644 --- a/docs/html/userhtmlse9.html +++ b/docs/html/userhtmlse9.html @@ -3,8 +3,8 @@ Contributor Covenant Code of Conduct - - + + diff --git a/docs/html/userhtmlsu1.html b/docs/html/userhtmlsu1.html index d3d37856..8498220a 100644 --- a/docs/html/userhtmlsu1.html +++ b/docs/html/userhtmlsu1.html @@ -3,8 +3,8 @@ Prerequisites - - + + diff --git a/docs/html/userhtmlsu10.html b/docs/html/userhtmlsu10.html index c8ffd9e1..0340d20e 100644 --- a/docs/html/userhtmlsu10.html +++ b/docs/html/userhtmlsu10.html @@ -3,14 +3,14 @@ Method smoothers_build - - + + -

[nextMethod smoothers_build

-

-

+

call p,imold])

-

This method builds the smoothers and the coarsest-level solvers for the multilevel preconditioner  5.3). -

Arguments

@@ -262,7 +262,7 @@ class="td11">

a

type(psb_).

The sparse matrix structure containing the local part of the matrix to be preconditioned. Note that ].

desc_a

type(psb_desc_type).

The communication descriptor of a. See the PSBLAS User’s Guide @@ -150,9 +150,9 @@ class="cmr-12">].

info

integer, ).

Error code. If no error, 0 is returned. See Section 7 for details.

amold

class(psb_optional.

The desired dynamic type for internal matrix components; this allows e.g. running on GPUs; it needs not be the same on all @@ -205,9 +205,9 @@ class="cmr-12">].

vmold

class(psb_optional.

The desired dynamic type for internal vector components; this allows e.g. running on GPUs.

imold

class(psb_i_base_vect_typeoptional.

The desired dynamic type for internal integer vector components; this allows e.g. running on GPUs.

-

[next] [up]

-

diff --git a/docs/html/userhtmlsu11.html b/docs/html/userhtmlsu11.html index e00f2580..d4117714 100644 --- a/docs/html/userhtmlsu11.html +++ b/docs/html/userhtmlsu11.html @@ -3,14 +3,14 @@ Method build - - + + -

[next5.5 Method build

-

-

+

call p,imold])

-

This method builds the preconditioner p according to the requirements made by the @@ -83,7 +83,7 @@ class="cmtt-12">smoothers_build, whose nomenclature would however be somewhat unnatural when dealing with simple one-level preconditioners. -

Arguments

a

type(psb_).

The sparse matrix structure containing the local part of the matrix to be preconditioned. Note that ].

desc_a

type(psb_desc_type).

The communication descriptor of a. See the PSBLAS User’s Guide @@ -158,9 +158,9 @@ class="cmr-12">].

info

integer, ).

Error code. If no error, 0 is returned. See Section 7 for details.

amold

class(psb_optional.

The desired dynamic type for internal matrix components; this allows e.g. running on GPUs; it needs not be the same on all @@ -213,9 +213,9 @@ class="cmr-12">].

vmold

class(psb_optional.

The desired dynamic type for internal vector components; this allows e.g. running on GPUs.

imold

class(psb_i_base_vect_typeoptional.

The desired dynamic type for internal integer vector components; this allows e.g. running on GPUs.

-

The method can be used to build multilevel preconditioners too. @@ -272,7 +272,7 @@ class="cmr-12">The method can be used to build multilevel preconditioners too.

[next] [up]

-

diff --git a/docs/html/userhtmlsu12.html b/docs/html/userhtmlsu12.html index 9bc621c1..7c8cb956 100644 --- a/docs/html/userhtmlsu12.html +++ b/docs/html/userhtmlsu12.html @@ -3,14 +3,14 @@ Method apply - - + + -

[next5.6 Method apply

-

-

+

call p,work])

-

This method computes y = psb_krylov and hence it is completely transparent to the user. -

Arguments

x

type(kind_parameter), dimension(:), intent(in)—.

The local part of the vector x. Note that , single/double class="cmr-12">precision version of AMG4PSBLAS under use.

y

type(kind_parameter), dimension(:), intent(out)—.

The local part of the vector y. Note that , single/double class="cmr-12">precision version of AMG4PSBLAS under use.

desc_a

type(psb_desc_type).

The communication descriptor associated to the matrix to be preconditioned.

info

integer, ).

Error code. If no error, 0 is returned. See Section 7 for details.

trans

character(lenin).

If trans = 1).

work

type(kind_parameter), dimension(:), optional, target—.

Workspace. Its size should be at least 4 

-

[next] [up]

-

diff --git a/docs/html/userhtmlsu13.html b/docs/html/userhtmlsu13.html index b5d808cd..f80e1a19 100644 --- a/docs/html/userhtmlsu13.html +++ b/docs/html/userhtmlsu13.html @@ -3,14 +3,14 @@ Method free - - + + -

[next5.7 Method free

-

-

+

call p,info)

-

This method deallocates the preconditioner data structure p. -

Arguments

info

integer, ).

Error code. If no error, 0 is returned. See Section 

-

[next] [up]

-

diff --git a/docs/html/userhtmlsu14.html b/docs/html/userhtmlsu14.html index 4e4894be..705fc60f 100644 --- a/docs/html/userhtmlsu14.html +++ b/docs/html/userhtmlsu14.html @@ -3,14 +3,14 @@ Method descr - - + + -

[next5.8 Method descr

-

-

+

call p verbosity])

-

This method prints a description of the preconditioner p to the standard output or to a @@ -68,7 +68,7 @@ class="cmtt-12">build, have been called. -

Arguments

info

integer, ).

Error code. If no error, 0 is returned. See Section 7 for details.

iout

integer, optional.

The id of the file where the preconditioner description will be printed; the default is the standard output.

root

integer, optional.

The id of the process where the preconditioner description will be printed; the default is psb_root_.

verbosity

integer, optional.

The verbosity level of the description. Default value is 0. For values higher than 0, it prints out further information, e.g., for @@ -174,7 +174,7 @@ class="td11">

-

[next] [up]

-

diff --git a/docs/html/userhtmlsu15.html b/docs/html/userhtmlsu15.html index 427962c7..ddbe5fb8 100644 --- a/docs/html/userhtmlsu15.html +++ b/docs/html/userhtmlsu15.html @@ -3,14 +3,14 @@ Auxiliary Methods - - + + -

[prev]

class="cmr-12">5.9 Auxiliary Methods -

Various functionalities are implemented as additional methods of the preconditioner object. -

+

5.9.1 Method: dump
-

-

+

call p,global_num])

-

Dump on file. -

Arguments

info

integer, ).

Error code. If no error, 0 is returned. See Section 7 for details.

amold

class(psb_optional.

The desired dynamic type for internal matrix components; this allows e.g. running on GPUs; it needs not be the same on all @@ -134,15 +134,15 @@ class="cmr-12">.

-

+

5.9.2 Method: clone
-

-

+

call p)
-

Create a (deep) copy of the preconditioner object. -

Arguments

pout

type(amg_).

The copy of the preconditioner data structure. Note that x must @@ -199,9 +199,9 @@ class="cmr-12">, single/double precision class="cmr-12">version of AMG4PSBLAS under use.

info

integer, ).

Error code. If no error, 0 is returned. See Section for details.

-

+

5.9.3 Method: sizeof
-

-

+

sz =

global

logical, optional.

Whether the global or local preconditioner memory occupation is desired. Default: .

Return memory footprint in bytes. -

+

5.9.4 Method: allocate_wrk
-

-

+

call p vmold])

-

Allocate internal work vectors. Each application of the preconditioner uses a number of work vectors which are allocated internally as necessary; therefore allocation and @@ -326,7 +326,7 @@ class="cmr-12">argument to the apply method. -

Arguments

info

integer, ).

Error code. If no error, 0 is returned. See Section 7 for details.

vmold

class(psb_optional.

The desired dynamic type for internal vector components; this allows e.g. running on GPUs.

-

+

5.9.5 Method: free_wrk
-

-

+

call p(info)

-

Deallocate internal work vectors. -

Arguments

single precision version –
double precision version –
For a first use it is suggested to use the default options obtained by simply selecting the solver type.
-

info

integer, ).

Error code. If no error, 0 is returned. See Section  Optional third party libraries - - + + diff --git a/docs/html/userhtmlsu3.html b/docs/html/userhtmlsu3.html index 22f31663..d6f3c0dc 100644 --- a/docs/html/userhtmlsu3.html +++ b/docs/html/userhtmlsu3.html @@ -3,8 +3,8 @@ Configuration options - - + + @@ -44,9 +44,9 @@ class="cmr-12">As a minimal example consider the following: - 

+   

 ./configure --with-psblas=PSB-INSTALL-DIR
-
+

which assumes that the various MPI compilers and support libraries are available in might be configured with: - 

+   

 ./configure --with-psblas=/opt/psblas-3.7/ \
---with-umfpackincdir=/usr/include/suitesparse/
-
+ 
--with-umfpackincdir=/usr/include/suitesparse/ +

Once the configure script has completed execution, it will have generated the file To build the library the user will now enter - 

+   

 make
-
+

followed (optionally) by - 

+   

 make install
-
+

diff --git a/docs/html/userhtmlsu4.html b/docs/html/userhtmlsu4.html index 0ae348d6..f6e35e2e 100644 --- a/docs/html/userhtmlsu4.html +++ b/docs/html/userhtmlsu4.html @@ -3,8 +3,8 @@ Bug reporting - - + + diff --git a/docs/html/userhtmlsu5.html b/docs/html/userhtmlsu5.html index 128fca5e..dac9747d 100644 --- a/docs/html/userhtmlsu5.html +++ b/docs/html/userhtmlsu5.html @@ -3,8 +3,8 @@ Example and test programs - - + + diff --git a/docs/html/userhtmlsu6.html b/docs/html/userhtmlsu6.html index 921424cf..6adab472 100644 --- a/docs/html/userhtmlsu6.html +++ b/docs/html/userhtmlsu6.html @@ -3,8 +3,8 @@ Examples - - + + @@ -128,58 +128,59 @@ class="cmr-12">. -

+

   use psb_base_mod
-  use amg_prec_mod
-  use psb_krylov_mod
-... ...
-!
-! sparse matrix
-  type(psb_dspmat_type) :: A
-! sparse matrix descriptor
-  type(psb_desc_type)   :: desc_A
-! preconditioner
-  type(amg_dprec_type)  :: P
-! right-hand side and solution vectors
-  type(psb_d_vect_type) :: b, x
-... ...
-!
-! initialize the parallel environment
-  call psb_init(ctxt)
-  call psb_info(ctxt,iam,np)
-... ...
-!
-! read and assemble the spd matrix A and the right-hand side b
-! using PSBLAS routines for sparse matrix / vector management
-... ...
-!
-! initialize the default multilevel preconditioner, i.e. V-cycle
-! with basic smoothed aggregation, 1 hybrid forward/backward
-! GS sweep as pre/post-smoother and UMFPACK as coarsest-level
-! solver
-  call P%init(’ML’,info)
-!
-! build the preconditioner
-  call P%hierarchy_build(A,desc_A,info)
-  call P%smoothers_build(A,desc_A,info)
-!
-! set the solver parameters and the initial guess
-  ... ...
-!
-! solve Ax=b with preconditioned CG
-  call psb_krylov(’CG’,A,P,b,x,tol,desc_A,info)
-  ... ...
-!
-! deallocate the preconditioner
-  call P%free(info)
-!
-! deallocate other data structures
-  ... ...
-!
-! exit the parallel environment
-  call psb_exit(ctxt)
-  stop
-
+ 
  use amg_prec_mod + 
  use psb_krylov_mod + 
... ... + 
! + 
! sparse matrix + 
  type(psb_dspmat_type) :: A + 
! sparse matrix descriptor + 
  type(psb_desc_type)   :: desc_A + 
! preconditioner + 
  type(amg_dprec_type)  :: P + 
! right-hand side and solution vectors + 
  type(psb_d_vect_type) :: b, x + 
... ... + 
! + 
! initialize the parallel environment + 
  call psb_init(ctxt) + 
  call psb_info(ctxt,iam,np) + 
... ... + 
! + 
! read and assemble the spd matrix A and the right-hand side b + 
! using PSBLAS routines for sparse matrix / vector management + 
... ... + 
! + 
! initialize the default multilevel preconditioner, i.e. V-cycle + 
! with basic smoothed aggregation, 1 hybrid forward/backward + 
! GS sweep as pre/post-smoother and UMFPACK as coarsest-level + 
! solver + 
  call P%init(’ML’,info) + 
! + 
! build the preconditioner + 
  call P%hierarchy_build(A,desc_A,info) + 
  call P%smoothers_build(A,desc_A,info) + 
+ 
! + 
! set the solver parameters and the initial guess + 
  ... ... + 
! + 
! solve Ax=b with preconditioned CG + 
  call psb_krylov(’CG’,A,P,b,x,tol,desc_A,info) + 
  ... ... + 
! + 
! deallocate the preconditioner + 
  call P%free(info) + 
! + 
! deallocate other data structures + 
  ... ... + 
! + 
! exit the parallel environment + 
  call psb_exit(ctxt) + 
  stop +

@@ -288,19 +289,19 @@ class="cmr-12">.

-

-... ...
-! build a V-cycle preconditioner with 1 block-Jacobi sweep (with
-! ILU(0) on the blocks) as pre- and post-smoother, and 8  block-Jacobi
-! sweeps (with ILU(0) on the blocks) as coarsest-level solver
-  call P%init(’ML’,info)
-  call P%set(’SMOOTHER_TYPE’,’BJAC’,info)
-  call P%set(’COARSE_SOLVE’,’BJAC’,info)
-  call P%set(’COARSE_SWEEPS’,8,info)
-  call P%hierarchy_build(A,desc_A,info)
-  call P%smoothers_build(A,desc_A,info)
+

 ... ...
-
+ 
! build a V-cycle preconditioner with 1 block-Jacobi sweep (with + 
! ILU(0) on the blocks) as pre- and post-smoother, and 8  block-Jacobi + 
! sweeps (with ILU(0) on the blocks) as coarsest-level solver + 
  call P%init(’ML’,info) + 
  call P%set(’SMOOTHER_TYPE’,’BJAC’,info) + 
  call P%set(’COARSE_SOLVE’,’BJAC’,info) + 
  call P%set(’COARSE_SWEEPS’,8,info) + 
  call P%hierarchy_build(A,desc_A,info) + 
  call P%smoothers_build(A,desc_A,info) + 
... ... +


Listing 2: setup of a multilevel preconditioner based on the default decoup

-

-... ...
-! build a W-cycle preconditioner with 2 hybrid Gauss-Seidel sweeps
-! as pre- and post-smoother, a distributed coarsest
-! matrix, and MUMPS as coarsest-level solver
-  call P%init(’ML’,info)
-  call P%set(’PAR_AGGR_ALG’,’COUPLED’,info)
-call P%set(’AGGR_TYPE’,’MATCHBOXP’,info)
-call P%set(’AGGR_SIZE’,8,info)
-  call P%set(’ML_CYCLE’,’WCYCLE’,info)
-  call P%set(’SMOOTHER_TYPE’,’FBGS’,info)
-  call P%set(’SMOOTHER_SWEEPS’,2,info)
-  call P%set(’COARSE_SOLVE’,’KRM’,info)
-  call P%hierarchy_build(A,desc_A,info)
-  call P%smoothers_build(A,desc_A,info)
+

 ... ...
-
+ 
! build a W-cycle preconditioner with 2 hybrid Gauss-Seidel sweeps + 
! as pre- and post-smoother, a distributed coarsest + 
! matrix, and MUMPS as coarsest-level solver + 
  call P%init(’ML’,info) + 
  call P%set(’PAR_AGGR_ALG’,’COUPLED’,info) + 
call P%set(’AGGR_TYPE’,’MATCHBOXP’,info) + 
call P%set(’AGGR_SIZE’,8,info) + 
  call P%set(’ML_CYCLE’,’WCYCLE’,info) + 
  call P%set(’SMOOTHER_TYPE’,’FBGS’,info) + 
  call P%set(’SMOOTHER_SWEEPS’,2,info) + 
  call P%set(’COARSE_SOLVE’,’KRM’,info) + 
  call P%hierarchy_build(A,desc_A,info) + 
  call P%smoothers_build(A,desc_A,info) + 
... ... +


Listing 3:

-

+

 ... ...
-! set RAS with overlap 2 and ILU(0) on the local blocks
-  call P%init(’AS’,info)
-  call P%set(’SUB_OVR’,2,info)
-  call P%bld(A,desc_A,info)
-... ...
-! solve Ax=b with preconditioned BiCGSTAB
-  call psb_krylov(’BICGSTAB’,A,P,b,x,tol,desc_A,info)
-
+ 
! set RAS with overlap 2 and ILU(0) on the local blocks + 
  call P%init(’AS’,info) + 
  call P%set(’SUB_OVR’,2,info) + 
  call P%bld(A,desc_A,info) + 
... ... + 
! solve Ax=b with preconditioned BiCGSTAB + 
  call psb_krylov(’BICGSTAB’,A,P,b,x,tol,desc_A,info) +


Listing 4: Method init - - + + diff --git a/docs/html/userhtmlsu8.html b/docs/html/userhtmlsu8.html index 0bf399ce..b6230e14 100644 --- a/docs/html/userhtmlsu8.html +++ b/docs/html/userhtmlsu8.html @@ -3,8 +3,8 @@ Method set - - + + @@ -344,16 +344,14 @@ class="cmr-12">the coarsening algorithm;
  • the coarse-space correction at the coarsest level (for multilevel - preconditioners only); +class="cmr-12">the solver at the coarsest level (for multilevel preconditioners only);
    • - - -
  • the smoother of the multilevel preconditioners, or the one-level + + + preconditioner.

    • the Gauss-Seidel solver object, while the point-Jacobi smoother i of combining the block-Jacobi smoother object with a single sweep of the point-Jacobi solver object. However, for simplicity, shortcuts are provided to set +class="cmr-12">point-Jacobi solver object. In the same way are obtained the ℓ1-versions of the +smoothers. However, for simplicity, shortcuts are provided to set all versions of point-Jacobi, hybrid (forward) Gauss-Seidel, and hybrid backward Gauss-Seidel, i.e., ’ as well. -

    The smoother and solver objects are arranged in a hierarchical manner. When specifying a smoother object, its parameters, including the local solver, are set to @@ -442,7 +445,7 @@ class="cmti-12">after the call to set the smoother. -

    Similar considerations apply to the point-Jacobi, Gauss-Seidel and block-Jacobi coarsest-level solvers, and shortcuts are available in this case too (see Table5).
    -

    Remark 3. Many of the coarsest-level solvers cannot be used with both replicated and distributed coarsest-matrix layouts; therefore, setting the s class="cmr-12">may change the layout. Similarly, setting the layout after the solver may change the solver. -

    More precisely, UMFPACK and SuperLU require the coarsest-level matrix to be replicated, while SuperLU_Dist requires it to be distributed. In these cases, setting the +class="cmr-12">_Dist and KRM requires it to be distributed. In these cases, coarsest-level solver implies that the layout is redefined according to the solver, +class="cmr-12">setting the coarsest-level solver implies that the layout is redefined according to the ovverriding any previous settings. MUMPS, point-Jacobi, hybrid Gauss-Seidel and +class="cmr-12">solver, ovverriding any previous settings. MUMPS, point-Jacobi, hybrid Gauss-Seidel block-Jacobi can be applied to replicated and distributed matrices, thus their choice +class="cmr-12">and block-Jacobi can be applied to replicated and distributed matrices, thus their does not modify any previously specified layout. It is worth noting that, when the +class="cmr-12">choice does not modify any previously specified layout. It is worth noting that, when matrix is replicated, the point-Jacobi, hybrid Gauss-Seidel and block-Jacobi +class="cmr-12">the matrix is replicated, the point-Jacobi, hybrid Gauss-Seidel and block-Jacobi solvers solvers reduce to the corresponding local solver objects (see Remark 2). For the +class="cmr-12">and their ℓ1- versions reduce to the corresponding local solver objects (see Remark 2). point-Jacobi and Gauss-Seidel solvers, these objects correspond to a For the point-Jacobi and Gauss-Seidel solvers, these objects correspond to a single @@ -491,7 +498,7 @@ class="cmti-12">single Gauss-Seidel sweep, respectively, which are very poor solvers. -

    On the other hand, the distributed layout can be used with any solver but UMFPACK and SuperLU; therefore, if any of these two solvers has already been @@ -507,7 +514,7 @@ class="cmr-12">_Dist has been class="cmr-12">previously set, the coarsest-level solver is changed to the default sequential solver. -

    Remark 4. The argument idx . -

    @@ -540,7 +547,7 @@ class="cmr-12">.
    -

    +

    +class="td11">

    1

    the minimum coarsening ratio, or the maximum number of levels is reached.
    -
    Therefore, the actual number of levels may be smaller than the specified maximum number of levels.
    -

    what

    		type

    val

    deault

    cots

    ’ML_CYCLE’

    character(len=*)

    ’VCYCLE’ -

    ’WCYCLE’ -

    ’KCYCLE’ -

    ’ADD’

    ’VCYCLE’

    Multilevel cycle: V-cycle, W-cycle, K-cycle, +class="td11">

    Multilevel cycle: V-cycle, W-cycle, K-cycle, and additive composition.

    ’CYCLE_SWEEPS’

    integer

    Any integer -

    number

    Any integer +

    number ≥ 1

    1

    Number of multilevel cycles.

    Number of multilevel cycles.

    -

    @@ -656,7 +663,7 @@ applied.
    -

    +

    @@ -670,7 +677,7 @@ id="TBL-5-4">

    what

    		type

    val

    deault

    cots

    ’MIN_COARSE_SIZE_PER_PROCESS’

    integer

    Any number -

    Any number +

    > 0

    200

    Coarse size threshold per process. The +class="td11">

    200

    Coarse size threshold per process. The aggregation stops if the global number of variables of the computed coarsest matrix is lower than or equal to this threshold @@ -724,17 +731,17 @@ multiplied by the number of processes

    ’MIN_COARSE_SIZE’

    integer

    Any number -

    Any number +

    > 0

    -1

    Coarse size threshold. The aggregation +class="td11">

    -1

    Coarse size threshold. The aggregation stops if the global number of variables of the computed coarsest matrix is lower than or equal @@ -747,17 +754,17 @@ class="cmtt-10x-x-109">’.

    ’MIN_CR_RATIO’

    real

    Any number -

    Any number +

    > 1

    1.5

    Minimum coarsening +class="td11">

    1.5

    Minimum coarsening ratio. The aggregation stops if the ratio between the global matrix dimensions at two consecutive levels is lower than or @@ -765,41 +772,41 @@ equal to this threshold (see Note).

    ’MAX_LEVS’

    integer

    Any integer -

    number

    Any integer +

    number > 1

    20

    Maximum number of levels. The +class="td11">

    20

    Maximum number of levels. The aggregation stops if the number of levels reaches this value (see Note).

    ’PAR_AGGR_ALG’

    character(len=*)

    ’DEC’, ’SYMDEC’, ’COUPLED’

    ’DEC’

    Parallel aggregation algorithm. -

    the

    Parallel aggregation algorithm. +

    the SYMDEC option applies decoupled aggregation to the sparsity pattern of T .

    ’AGGR_TYPE’

    character(len=*)

    ’SOC1’, @@ -830,11 +837,11 @@ class="cmtt-10x-x-109">’, class="cmtt-10x-x-109">’MATCHBOXP’

    ’SOC1’

    Type of aggregation algorithm: currently, +class="td11">

    Type of aggregation algorithm: currently, for the decoupled aggregation we implement two measures of strength of connection, the @@ -850,22 +857,22 @@ href="userhtmlli5.html#XMatchBoxP">9].

    ’AGGR_SIZE’

    integer

    Any integer -

    power of 2, +class="td11">

    Any integer +

    power of 2, with aggr_size ≥ 2

    4

    Maximum size of aggregates when the +class="td11">

    4

    Maximum size of aggregates when the coupled aggregation based on matching is applied. For aggressive coarsening with size of aggregate larger than 8 @@ -877,16 +884,16 @@ class="cmtt-10x-x-109">’MATCHBOXP’

    ’AGGR_PROL’

    character(len=*)

    ’SMOOTHED’, @@ -894,25 +901,25 @@ class="cmtt-10x-x-109">’, class="cmtt-10x-x-109">’UNSMOOTHED’

    ’SMOOTHED’

    Prolongator used by the aggregation +class="td11">

    Prolongator used by the aggregation algorithm: smoothed or unsmoothed (i.e., tentative prolongator).

    Note. The aggregation algorithm stops when at least one of the following criteria is met: the coarse size threshold,
    -
    the minimum coarsening ratio, or the maximum number of levels is reached.
    +
    Therefore, the actual number of levels may be smaller than the specified maximum number of levels.
    +
    @@ -933,7 +940,7 @@ class="content">Parameters defining the aggregation algorithm.

    @@ -941,7 +948,7 @@ class="content">Parameters defining the aggregation algorithm.

    +

    tine

    what

    		type

    val

    deault

    cots

    ’AGGR_ORD’

    character(len=*)

    ’NATURAL’ -

    ’DEGREE’

    ’NATURAL’

    Initial ordering of indices for +class="td11">

    Initial ordering of indices for the decoupled aggregation algorithm: either natural ordering or sorted by descending degrees of the nodes in the @@ -1010,57 +1017,57 @@ matrix graph.

    ’AGGR_THRESH’

    real(kind_parameter)

    Any real -

    number 

    Any real +

    number ∈ [0,1]

    0.01

    The threshold

    0.01

    The threshold θ in the strength of connection algorithm. See also the note at the bottom of this table.

    ’AGGR_FILTER’

    character(len=*)

    ’FILTER’ -

    ’NOFILTER’

    ’NOFILTER’

    Matrix used in computing the smoothed +class="td11">

    Matrix used in computing the smoothed prolongator: filtered or unfiltered.

    Note. Different thresholds at different levels, such as those used in [33, Section 5.1], can be easily set by invoking the rou-
    -
    tine set with the parameter ilev.
    -
    @@ -1076,7 +1083,7 @@ class="content">Parameters defining the aggregation algorithm (continued).

    @@ -1084,7 +1091,7 @@ class="content">Parameters defining the aggregation algorithm (continued). -

    +

    @@ -1098,7 +1105,7 @@ id="TBL-7-4">

    what

    		type

    val

    deault

    cots

    ’COARSE_MAT’

    character(len=*)

    ’DIST’ -

    ’REPL’

    ’REPL’

    Coarsest matrix layout: distributed among the +class="td11">

    Coarsest matrix layout: distributed among the processes or replicated on each of them.

    ’COARSE_SOLVE’

    character(len=*)

    ’MUMPS’ -

    ’UMF’ -

    ’SLU’ -

    ’SLUDIST’ -

    ’JACOBI’ -

    ’GS’ -

    ’BJAC’ -

    ’KRM’

    See Note.

    Solver used at the coarsest level: sequential +class="td11">

    See Note.

    Solver used at the coarsest level: sequential LU from MUMPS, UMFPACK, or SuperLU (plus triangular solve); distributed LU from MUMPS or SuperLU_Dist (plus triangular solve); @@ -1231,53 +1238,53 @@ double precision.

    ’COARSE_SUBSOLVE’

    character(len=*)

    ’ILU’ -

    ’ILUT’ -

    ’MILU’ -

    ’MUMPS’ -

    ’SLU’ -

    ’UMF’ -

    ’INVT’ -

    ’INVK’ -

    ’AINV’

    See Note.

    Solver for the diagonal blocks of the coarsest +class="td11">

    See Note.

    Solver for the diagonal blocks of the coarsest matrix, in case the block Jacobi solver is chosen as coarsest-level solver: ILU(p), ILU(

    Note. Defaults for COARSE_SOLVE and COARSE_SUBSOLVE are chosen in the following order:
    -
    single precision version – MUMPS if installed, then SLU if installed, ILU otherwise;
    -
    double precision version – UMF if installed, then MUMPS if installed, then SLU if installed, ILU otherwise.
    -

    what

    		type

    val

    deault

    cos

    ’COARSE_SWEEPS’

    integer

    Any +class="td11">

    Any integer -

    number

    number > 0

    10

    Number of sweeps when

    10

    Number of sweeps when JACOBI, GS or BJAC is @@ -1382,52 +1389,52 @@ chosen as coarsest-level solver.

    ’COARSE_FILLIN’

    integer

    Any +class="td11">

    Any integer -

    number

    number ≥ 0

    0

    Fill-in level

    0

    Fill-in level p of the ILU factorizations and first fill-in for the approximate inverses.

    ’COARSE_ILUTHRS’

    real(kind_parameter)

    Any real -

    number

    Any real +

    number ≥ 0

    0

    Drop tolerance

    0

    Drop tolerance t in the ILU(p,t) factorization and first drop-tolerance for the approximate inverses.

    Note. Further options for coarse solvers are contained in Table 6.
    -
    For a first use it is suggested to use the default options obtained by simply selecting the solver type.
    +
    @@ -1447,7 +1454,7 @@ class="content">Parameters defining the solver at the coarsest level (continued) -

    @@ -1455,7 +1462,7 @@ class="content">Parameters defining the solver at the coarsest level (continued)
    -

    +

    @@ -1469,7 +1476,7 @@ id="TBL-8-4">

    what type

    val

    deault

    cots

    ’BJAC_STOP’

    character(len=*)

    ’FALSE’ -

    ’TRUE’

    ’FALSE’

    Select whether to use a stopping criterion for +class="td11">

    Select whether to use a stopping criterion for the Block-Jacobi method used as a coarse solver.

    ’BJAC_TRACE’

    character(len=*)

    ’FALSE’ -

    ’TRUE’

    ’FALSE’

    Select whether to print a trace for the +class="td11">

    Select whether to print a trace for the calculated residual for the Block-Jacobi method used as a coarse solver.

    ’BJAC_ITRACE’

    integer

    Any integer -

    Any integer +

    > 0

    -1

    Number of iterations after which a trace is to +class="td11">

    -1

    Number of iterations after which a trace is to be printed.

    ’BJAC_RESCHECK’

    integer

    Any integer -

    Any integer +

    > 0

    -1

    Number of iterations after which a residual is +class="td11">

    -1

    Number of iterations after which a residual is to be calculated.

    ’BJAC_STOPTOL’

    real(kind_parameter)

    Any real -

    Any real +

    < 1

    0

    Tolerance for the stopping criterion on the +class="td11">

    0

    Tolerance for the stopping criterion on the residual.

    ’KRM_METHOD’

    character(len=*)

    ’CG’ -

    ’FCG’ -

    ’CGS’ -

    ’CGR’ -

    ’BICG’ -

    ’BICGSTAB’ -

    ’BICGSTABL’ -

    ’RGMRES’

    ’FCG’

    A string that defines the iterative method to +class="td11">

    A string that defines the iterative method to be used when employing a Krylov method ’

    ’KRM_KPREC’

    character(len=*)

    Table 

    Table 1

    ’BJAC’

    The one-level +class="td11">

    The one-level preconditioners from the Table 1 can be used for the coarse Krylov solver.

    ’KRM_SUB_SOLVE’

    character(len=*)

    Table 

    Table 5

    ’ILU’

    Solver for the diagonal blocks of the coarsest +class="td11">

    Solver for the diagonal blocks of the coarsest matrix preconditioner, in case the block Jacobi solver is chosen as 5 applies here.

    ’KRM_GLOBAL’

    character(len=*)

    ’TRUE’, @@ -1762,46 +1769,46 @@ class="cmtt-10x-x-109">’, class="cmtt-10x-x-109">’FALSE’

    ’FALSE’

    Choose between a global Krylov solver, all +class="td11">

    Choose between a global Krylov solver, all unknowns on a single node, or a distributed one. The default choice is the distributed solver.

    ’KRM_EPS’

    real(kind_parameter)

    Real

    Real < 1

    10

    10-6

    The stopping tolerance. +class="td11">

    The stopping tolerance.

    ’KRM_IRST’

    integer

    Integer -

    Integer +

    ≥ 1

    30

    An integer specifying the restart parameter. +class="td11">

    30

    An integer specifying the restart parameter. This is employed for the BiCGSTABL or RGMRES @@ -1809,16 +1816,16 @@ methods, otherwise it is ignored.

    ’KRM_ISTOPC’

    integer

    Integers +class="td11">

    Integers 1,2,3

    2

    If

    2

    If 1 then the method uses the normwise backward error in the infinity norm; if 2, the @@ -1831,32 +1838,32 @@ for the details.

    ’KRM_ITMAX’

    integer

    Integer -

    Integer +

    ≥ 1

    40

    The maximum number of iterations to +class="td11">

    40

    The maximum number of iterations to perform.

    ’KRM_ITRACE’

    integer

    Integer -

    Integer +

    ≥ 0

    -1

    If

    -1

    If > 0 print out an informational message about convergence every

    ’KRM_FILLIN’

    integer

    Integer -

    Integer +

    ≥ 0

    0

    Fill-in level

    0

    Fill-in level p of the ILU factorizations and first fill-in for the approximate inverses. Additional parameters defining the solver at the coarsest level. -

    @@ -1907,7 +1914,7 @@ class="content">Additional parameters defining the solver at the coarsest level.
    -

    +

    @@ -2241,7 +2248,7 @@ preconditioner. -

    @@ -2249,7 +2256,7 @@ preconditioner.
    -

    +

    what

    		type

    val

    deault

    cots

    ’SMOOTHER_TYPE’

    character(len=*)

    ’JACOBI’ -

    ’GS’ -

    ’BGS’ -

    ’BJAC’ -

    ’AS’ -

    ’L1-JACOBI’ -

    ’L1-BJAC’ -

    ’L1-FBGS’

    ’FBGS’

    Type of smoother used in the multilevel preconditioner: point-Jacobi, hybrid @@ -2023,21 +2030,21 @@ class="cmr-10">–hybrid (forward) Gauss-Seidel, class="cmmi-10">ℓ1-point-Jacobi and Additive Schwarz. -

    It is ignored by one-level preconditioners.

    ’SUB_SOLVE’

    character(len=*)

    ’JACOBI’ @@ -2045,47 +2052,47 @@ class="cmtt-10">’ class="cmtt-10">’GS’ -

    ’BGS’ -

    ’ILU’ -

    ’ILUT’ -

    ’MILU’ -

    ’MUMPS’ -

    ’SLU’ -

    ’UMF’ -

    ’INVT’ -

    ’INVK’ -

    ’AINV’

    GS and BGS multilevel class="cmr-10">preconditioners, respectively -

    ILU for block-Jacobi and Additive Schwarz class="cmr-10">one-level preconditioners

    The local solver to be used with the smoother or one-level preconditioner (see @@ -2163,22 +2170,22 @@ class="cmr-10">Note for details on hybrid Gauss-Seidel.

    ’SMOOTHER_SWEEPS’

    integer

    Any integer -

    number ≥ 0

    1

    1

    Number of sweeps of the smoother or one-level preconditioner. In the multilevel @@ -2205,22 +2212,22 @@ class="cmr-10">respectively.

    ’SUB_OVR’

    integer

    Any integer -

    number ≥ 0

    1

    1

    Number of overlap layers, for Additive Schwarz only.

    what

    		type

    val

    deault

    cots

    ’SUB_RESTR’

    character(len=*)

    ’HALO’ -

    ’NONE’

    ’HALO’

    Type of restriction operator, for Additive Schwarz only: ’NONE’ for neglecting it. -

    Note that HALO must be chosen for the @@ -2339,28 +2346,28 @@ class="cmr-10">its RAS variant.

    ’SUB_PROL’

    character(len=*)

    ’SUM’ -

    ’NONE’

    ’NONE’

    Type of prolongation operator, for Additive Schwarz only: ’ for neglecting them. -

    Note that ’SUMfor its RAS variant.

    ’SUB_FILLIN’

    integer

    Any integer -

    number ≥ 0

    0

    0

    Fill-in level p of the incomplete LU @@ -2416,25 +2423,25 @@ class="cmr-10">factorizations.

    ’SUB_ILUTHRS’

    real(kind_parameter)

    Any real number ≥ 0

    0

    0

    Drop tolerance t in the ILU(factorization.

    ’MUMPS_LOC_GLOB’

    character(len=*)

    ’LOCAL_SOLVER’ -

    ’GLOBAL_SOLVER’

    ’GLOBAL_SOLVER’

    Whether MUMPS should be used as a distributed solver, or as a serial solver acting @@ -2477,19 +2484,19 @@ class="cmr-10">process.

    ’MUMPS_IPAR_ENTRY’

    integer

    Any integer number

    0

    0

    Set an entry in the MUMPS integer control array, as chosen via the argument.

    ’MUMPS_RPAR_ENTRY’

    real

    Any real number

    0

    		Any real number

    0

    Set an entry in the MUMPS real control array, as chosen via the Parameters defining the smoother or the details of the one-level -

    [next] [up]

    -

    diff --git a/docs/html/userhtmlsu9.html b/docs/html/userhtmlsu9.html index 03afcbcc..83e77773 100644 --- a/docs/html/userhtmlsu9.html +++ b/docs/html/userhtmlsu9.html @@ -3,14 +3,14 @@ Method hierarchy_build - - + + -

    [nextMethod hierarchy_build

    -

    -

    +

    call p,info)

    -

    This method builds the hierarchy of matrices and restriction/prolongation operators for the multilevel preconditioner init and set. -

    Arguments

    a

    type(psb_).

    The sparse matrix structure containing the local part of the matrix to be preconditioned. Note that ].

    desc_a

    type(psb_desc_type).

    The communication descriptor of a. See the PSBLAS User’s Guide @@ -136,9 +136,9 @@ class="cmr-12">].

    info

    integer, ).

    Error code. If no error, 0 is returned. See Section 

    -

    [next] [up]

    -

    diff --git a/docs/src/userguide.pdf b/docs/src/userguide.pdf new file mode 120000 index 00000000..7b032aa3 --- /dev/null +++ b/docs/src/userguide.pdf @@ -0,0 +1 @@ +tmp/userguide.pdf \ No newline at end of file diff --git a/docs/src/userinterface.tex b/docs/src/userinterface.tex index 833d1959..683f8056 100644 --- a/docs/src/userinterface.tex +++ b/docs/src/userinterface.tex @@ -137,7 +137,7 @@ logically divided into four groups, i.e., parameters defining \begin{enumerate} \item the type of multilevel cycle and how many cycles must be applied; \item the coarsening algorithm; - \item the coarse-space correction at the coarsest level (for multilevel + \item the solver at the coarsest level (for multilevel preconditioners only); \item the smoother of the multilevel preconditioners, or the one-level preconditioner. @@ -155,8 +155,9 @@ the hybrid Gauss-Seidel smoother (see Note in Table~\ref{tab:p_smoother}) is obtained by combining the block-Jacobi smoother object with a single sweep of the Gauss-Seidel solver object, while the point-Jacobi smoother is the result of combining the block-Jacobi smoother object with a single sweep -of the point-Jacobi solver object. However, for simplicity, shortcuts are -provided to set point-Jacobi, hybrid (forward) Gauss-Seidel, and +of the point-Jacobi solver object. In the same way are obtained the $\ell_1$-versions of +the smoothers. However, for simplicity, shortcuts are +provided to set all versions of point-Jacobi, hybrid (forward) Gauss-Seidel, and hybrid backward Gauss-Seidel, i.e., the previous smoothers can be defined just by setting \fortinline|'SMOOTHER_TYPE'| to certain specific values (see Tables~\ref{tab:p_smoother}), without the need to set @@ -176,12 +177,12 @@ coarsest-level solvers, and shortcuts are available in this case too (see Table~\ref{tab:p_coarse_1}). \\ \textbf{Remark 3.} Many of the coarsest-level solvers cannot be used -with both replicated and distributed coarsest-matrix layouts; +with both replicated and distributed coarsest-matrix layouts; therefore, setting the solver after the layout may change the layout. Similarly, setting the layout after the solver may change the solver. More precisely, UMFPACK and SuperLU require the coarsest-level -matrix to be replicated, while SuperLU\_Dist requires it to be distributed. +matrix to be replicated, while SuperLU\_Dist and KRM requires it to be distributed. In these cases, setting the coarsest-level solver implies that the layout is redefined according to the solver, ovverriding any previous settings. MUMPS, point-Jacobi, @@ -189,7 +190,7 @@ hybrid Gauss-Seidel and block-Jacobi can be applied to replicated and distributed matrices, thus their choice does not modify any previously specified layout. It is worth noting that, when the matrix is replicated, -the point-Jacobi, hybrid Gauss-Seidel and block-Jacobi solvers +the point-Jacobi, hybrid Gauss-Seidel and block-Jacobi solvers and their $\ell_1-$ versions reduce to the corresponding local solver objects (see Remark~2). For the point-Jacobi and Gauss-Seidel solvers, these objects correspond to a \emph{single} point-Jacobi sweep and a \emph{single}