From 1c933464cc7344e73f804db56aa9d1126074fb5f Mon Sep 17 00:00:00 2001 From: Peter Schaefer Date: Tue, 3 Apr 2012 19:09:26 +0200 Subject: [PATCH] [src] test_sol angepasst [src] kappa & tmu2 eingebaut --- doc/doc.pdf | Bin 332349 -> 332621 bytes doc/doc.tex | 4 ++-- src/A_plots.m | 16 +++++++++----- src/A_step.m | 55 +++++++++++++++++++++++++++++++------------------ src/test_sol.m | 2 +- 5 files changed, 49 insertions(+), 28 deletions(-) diff --git a/doc/doc.pdf b/doc/doc.pdf index fd4522c2e6a8418399c8991bb1f80936ebcd07f2..5f9edd12531d941832268c01652762e707b63b29 100644 GIT binary patch delta 28069 zcmV(sK<&T1q!i7j6tIl80yQ$1@dE=Wf9+dKa~sDIzVlaX4ygi4!DQYKt8|bZ$yH8N zb|OYj>2UysCAlI1#&}4w4#{6n&+P1CX3+}}Bup`>Ttpb#=$Yxq*I)OnX>n1|;`}Z3Mon1xx`87-9W$L#z;>h8fi?q^N*@SX4$6;SnW`_#|7X+Dvq?OhY?6#P=k1^gr!k#AKi@VN#DcW4IVh$}5@%T>fADf!gCkP4 z&%b$sXz=l7TWfzt45WyLoFBpYT{P!8BoE_*YF;9H_C9t+BK>Bc2a%H-V%)xS^Bw1| z11&y(hdP0R-=&2Q;0zTLZZ=*FN8Apsh9kiSPp?Sx3GOi zVh@TGMM5ucebVyze=ADSxa7O)tBAqtt7s6KED@JDVu9&}x*+q#@`y9B9LR%o;ios? zi9{&E9T$1YYa{$`7&l;0n1F1ts6(Kemf97gL@J`=DYAC}e4g+91FDii$Rt90G)kAL zyh6IDe-##vj%b#^1*cf07MqBi zR3&KReDHcb^`c&bfV2h#kWDQ8J3N5Q*?cu#E~|^HA(Y9oS}o5)!7L}${9--7sP(Pq z6SUOoz}@2Nt6!mP>%l(}RswKlsNCA-w@`-O58?QibH6$)w;V54^|QE*<*J&^{KFMN z*%KvNwbdRee^ua12aB>elps%8P^BZ^JzG~( zI4vDakfoufYb}qmv<51vc((M@zau77L}-STciC6oPZ6_vKHmmrOtvq>k(7gm?C4z) zBNcI+Wt0G6ky1khYMryl-zo}#M1h4(p$yO~E$`wAe+D*oMh8{iyL-cc3E+T$PRCdq z+8@D9zV<uZ5mdh%W!t zk^jCLIytyrT$GpPeEG=};j&sHt5L5B7&#zY1^}PA6TCVC3<+A>}o)6g|ex7?MTn9-ynAHF-y%R)@;Bp(m z<5$-g)9SrP$;)P!alNN!)A9Pev-MX)3x;2um(W2j0CXGSFYD!aHme&0Fk-WGnp?D{O$`Ys_|@Wdj>u3juQBn9$1v$tl{>ap=+np)`ihdWJ;{ zfAGOgnRqZKQ`$Bs8|kKAr0M`wHXNS!3{apcbRy52a&bAHuIdy^Vg}#|;N1s}B#PQ3 zX!@fjk|)UddFhqQnpgQ%z8_D^#Z@)0NgxCEMWLtMhapQ(_Yf9p(+L$ct$N!WP|d`> zdC4!vZ|$WH0&&=P;wO$abeoRmCoDoufAFP( zWPI?ip6QrK=3ZPQX48OpA{*BRt39sJDHU`AiR3rfyepp&U+tU5{vOLtrdS%Ge+&&7 zx1m^wk0q3#T|SOIJ09b`jIke1%)^8ewhK=D*t-+w*mN>{N~TB(*_Hvs;GHnSY)O?( zBk1n%d}r9W8$54^6apIvYN{+jTAb+4VWfjnh=I(P9bOO;ihp%vuaf{ zr&%?xW}ZFO)#T6P#iW{#r_iH;M5rWiTt7+NTS@*$GcWja*hAca+Ff1#Y-mT`?Ck1x_PBvXdwsU@l3 zkTM#n;@zuhDC6K9-NF=?i2EQkrP%8}O)NKux(~gv5!!ViCv}U3-KhKVt|pAiz20Tv zKKLPcH=A5r*gY(mxmQu4e-H{Q{dC{I4kJBfV;kv@H*Ng`WVCJEi!n5DaiHWeEgjtS&pV}hGEMmMgXkJ#SUXvfi@rI*0l zlcVh!0&oSpgOWR9poAyq_=q$>U4|iu;@yLFxT&8?MFI0WEPgmKe_a)M7pFJT^{Bq3 zVX@3zK*%y^w%8!{Va-{4ngC6S-Lw=JoWVoKND`U)GNn6)3BpkJrnK^1L3%3O`@R>OXwL z{i`XP@7I@$Ro7`Uf1{(I%rMrS6*9Af!)DjZ+x>4o`msLx?El{6XD4yxKmjvI=^5&5 zX3O|g$L=5M^mn|y{mFfS*Il3sOh}E&xqqkHbwin19$yjhLDL(_{clQ?^^01IaI13d z(PV*ODwJ&1y~f+_MVNv=vHmWQ%QUFZo@8*~-i+>k4O@j5e-3xNINIn^5`M^<7ZCgf zg8x_dyo`_xY&Ln-c}v767VTA-&u08*?&t2ThaExD*6c+1+b1HOeTV>;U#;X0QqV^*nKII^YS?3WIyYGDAHZqH zlc4SgAO^tWfA8IdhR(;_0)h4qbmff)x-y3-**h<}+V_^=&Uvo>P3O7s%0eGwoclFy zhs+&shib(1b1;7MH~tijM&sKfTLD3I)-Tr?&aktAk>8Shiw}A1@=1ICm$qudd76)b)oRqrElTG2q`k zbn}wmvSdqi_G#o@D)t8AnDOrcQ*cf5vMkGK2s9S|0qeF~&68o=6NA>lx7NV{=a>;S zFbXeBWo~D5Xfhx%IW{&omtnL76$CLbIWdzVekXs7_XSW~+tz@K;=w%xrwJ0=-QC?C z8fe_z-Ccqc+}$lekl+M&CqRM(cfXx|_SxCzf3ND^>MDADbBr;E&e;nJQYBS-VN-i! zkc7RR3q3Of6AwUCUS67+3Bbg}%D}|Lf=EH3X6a%B`tLd-g$BsU+0x#Q=U*D4P9UHQ zcuju{=mJ)hx3>exxY_`iSpm!(Jj|RtOiTb4CMNFxD%v~o0K|Z9mZkuC27rvc9mpAx zLe$>D)5+4@!UgQ-zmEVa6KVi6H#aBUU*-T|Tac5b3D6E84|K5r*@8Wq0Br!O_9m7f z7tjCMf{M?=#l?Y#kh!moIy@*AXC8almQArThQOB zF(6U^)GREW|07VfH*;|ZI)MP-1sh8fkexHw#MRCeHP=+3s)a-y2!l{bdX^F|oIG0NQz4+L;5)ENwslMF}|u7Y`RY0MO3#w<3Sg z#@QZh4|D@s+5nBg3V($T1V{)g1AyQR|IW|Z#L3dZ#hJm`(&l%LjK8-5ms#A-RMg(q z7G&q*jQBe~F-s?q3ApZ_jDI)R+RonH&ilV)W@%??_B#zzR|iIQJ4;7bkhIv}I^ZS5 zKWpY77XUjG7bhzhI{@Se0C|{LF#dmbpyufS`b)_Cdl{UAueXD}1HcTN2FTaa3R zet-UX(g#nNslA^%V9^sFoZdKPvj05g9R8yA2R z{NVe)`}|XZ|7QQ!suIxh@96%uqqLowJ>a*C-&ugm_uow1{vHjgzt0Ud;D22x*nSQH!Cxl-|AznfF8{Cd|9>NraJ8}dOHTD4h5sKp(ALt%^KT9C zK)bqt4~M)x_=wp3UsX-ee@=h2Jjm42)%O2NrCosFLn3Tv4jy!RW;O;Uw*RbKI!jo3 zfJ~JvT}&+gGfe;O2H1yfK344 zoP7ai_D+bu8_mHEU=;qn^dI5`FpB;mE&!w0Kg0|Mkc|uQkC+?4DE<#|G65JR{t#GO z@(-~97^VIYD}Yh@53vClW&RKcfKl!bfy0siL*Q@}{t(!^;y=U%_OA4Yz}}Vr5ZJrQ z9|C(<{X<~yYJUh^4E29M1dc`H4}raF{X^W~$o^D<8!R*aL*U3v{vlRyEG8C~|1kS+ zDfpuMWBq#!7)}2mI3&~mDKsQ!s=fLj==gKWC4e74Rmn$cSPWV*!+J{0Ece_v^6&UR}5eq zo8NQ(#{^s#+dn4Y7;RnuNWpg#qun1VI48S*OTiB8|Ak;ckF?(E)rf+5ag$E7*&}zk0<2u7blKWc|Ib^S{dSyH$T`$O^8ClZ8F=zvN*1KPhs6 zqja_bI{&N8%;0~zy8Lkmj>rZ0ue^S<#s9SWckln5I9S8&UkI+D`yc$>7>|D;*uUq$ z5S)$I9|T9@^DSJh>LSp$VvxI+v)fWX3$Ge=3#Y^wqMv^8^ zwi|j9iG6bt%Op5uLOT_@d+N%v6;`qqEY=|D|g8mv1vgSFD>^7*xlH zx>Sr^aax>g>O)_%;DSMNn@?*uRH1X4EV^9j7M*GWC2QOv9RXSxtgI+*mgzO>sEo}& zhu#!x-*|hmy%bz_cZ*Q2IMznn#bU;IZ8fg>WESN=oH8?xG?`@7$B25l+iu_Z)Ou{| z%ME|o+86s)dv)5fA8##U#KqY207!Dm55Gi(#Nsk^M82-h*O6QF?T(BuVMq4btlW9ExF$56ol5|qF)-RrqPN19@Q z44bVm)A+LrYs*4A4fe2iH*En<$1Wx}#&&BCi zr1+Rl;jOFo+-U0AdV;xQ*I4SeYYpfj}h%^?; zW+(!;bBG0tx%N0X@xXjaJH3K23yJniv$=s~)WJSm>@mAC61v+%VAO0;s_{XWL|;w( zrq?P{WjM-@k@gX3cNgankeQCYZo{?frL~si&~TNT^+}3U<9033w1u=KfY*PNt)Nh? z7;TxIhi>P=4V%nFHlAi}ftS{6^6eX~m|9FP6l=h{BIVgl_W=Q*Lb8x3OB+c&8a$*!SbSYVSz7nIp=Go*u+etxw!Lqa|Ap zw>_ESX03xvpyeImkMTpZ8G<2F>^ki2l0 zA1hh~!5!~811=n0Lbl8ve0wl|On;=V3*&-ORhiXF@qCBLm+jodVwHcQXwbPG+<{Kf zgB-8^sUvjDTVO7hF|Vq&T-4TnbB}3)qRKblVfCHqPj3qxHi0P2s5b`gVkKunAmuAl z)UVKoKuc5YX@f1HFU^AR4#<|Qyo3#bA_MKhM4-n<;-RZ+Wa#Pl)v_iBKMWBOuLPb( zp49B8`HN@>vD!w)9z1`G0+z`fu;9ou)F%jmB>@cqpV*aIbXa+CKeZAJ^Jz;?do$m? znrKMl)j@Mm89K?ozl(GUH4O*|UawF%3S-eXu>@}4uwB0LpAvm@%J}T&xms&sK+`KF zyIk+1ceSJ)>5kQmFMm|XN=MWDHOh_Mm!YDVyX0M6?t2bkXA^%0DoFbv(>DQ{qnyQ8 zHkj$eP@@iK;&aTK)p5pBjvl|UNTv&Gk)(*h7n>6S_oNN2L%#1)+Oc(mFdlB0U0A*W zyAt%|S?{xqu}P_#y{yHQRg7JD8;4L2tcu#cezHc|RlvJptDc2pyU)1VY+QS_5~j&A zGM@=ouBdN7peld)$vYNUr;MUb5+dqgy#Gyz_pVjV*XoZi`CcvTbn24lud5ngj+3Tme;x6JjPY5g`a)|xS?hO(CZ8hK=D=vHOF-t37aCp@uRA|D)u*I7hgIVg z>*U%e?RtMtdylCWtD2KEkHQAku5@TU&t_`7cY+-`8v8*qEqKKDaI_@D1lGXNJY5{G zwQ;lb`T)$dMKfudbzL(&ekgfN9M)NCiq7-$8uVnxikcQhD%dZ2^g5(5c*HqJGX02B z)PzpTHwn!|%Ir6V{%W|!a=wfL6YSQQzi@B)mEC^}6kj-T*q3fDi9^J#G+@0yF-Kl9j|Z#Or)bDM&#-@cYRdX}1p2-V=P z@mgkG=}}vd@$3Ad+tA(oGS6!C?3VZ9VCNL9jx4XCS&?fqWL(Y}g3!fLRb5bu)lGlc ze|;0(7K#qCA7Dox@{kkz_C+{7s>T0#K%s(0;ZIfJRz#oJeppJj<_Y@^qOHVlh|Km1ti(dn4&!GxzHR{nvcZnqT=BNEg*2%X0GfY%rG5(=6J4g$TU z%NCIz+sP}5e4~eZN-^XP&{ov}J%b1C{Vt)jmSr23HAONNb2+U~bAT@O;(@U<=Y8#t zPgbG61&DfmE>fb0Z%#+=i8Fr~5yiIm0G~$@ye&arcaSrntN8M33Azs_pMClmywgnT z!fo!D;tC_EN^w^FMd8zFPrebh5iaYry*Q8qD9_P>rD#p$;uT2}Qn9uW)SK@%idbQ* zTsBqqEVFFLnJ7H0NuNbCcmlW&ei}~DfA+QZlbaU$d0f?pGA7CY)u4aUV(F0Dx;Crb z5`oKjHt|l#+2GxshvX|mVi-Ky5w7BGSf)<&iB;{61h(2*=M)MZIq)dhvCMj(R89Uj+gzEVb5AvN798^gLCOHPF;Gd zEqy^Fuiqb4K=wv%N!x!NPNwf11br3qie45@{b>1*~@oD_M&}941YR2E!zp(Gtwtji)V)@uD zQ29&N!;!VWcx%MarIe^T>mXsuXIgWSB7&)Kod1V-qEKzeM}|`tgDt1VVxC+oA}N6? z(Ef_j_0?B=Nq=G@Fy?4GTR#>KFKIxAMKV`kIb_36%~DD3;n zSGOA4SUW5PtUtv}cq%A9$+^@j*}rC5JzLc1CSZN4-$m|$b>zCb7y z-(Hb|;0sp|>0OxPmKgCFhZ2qX!cvmgo~<7l_q;dZ4yHfaXHDc`g~Wdwgx>j8%q2;I z5$Q0|IdU}EDu2#1oAO}y+(&Mr8gN9wPFj(a)4`My^?gMslr*))Z&wUqmW+!R*+2iH z>GjeXhA4m6;INM0T=c9Uu_?S=E<#-z;(U^c>$7LOx=Nt+7IXsm zjii5hG?Jk5Ti}&FfFkkG=4xl}b_&Ld{^6jW=E4ijRlQtw$P zVbkdyqL5@y{puo(cdl4T`CY{?l@C6q!Bc-7iSh9^TGf3WZGG5WmfpK(amFSpXxV%e zZZZ#(O#Y^zIuvCx$TIiTvFL=!5VpR#u_1_djF_va)|D z8)+{i3!U`HSrEu|7h@Pb793c+A&$OltLIla2E#VfKpEvZ0m^TMk*br`xdw*@>=S4l zoRJCdQN)XHY^+O-h3u4RFUNUR5@YOgk|J3=C&OxSLY7;6i9hel9xTafs0EW%a^Y{s zgRu#zCMU7l^HTR8~Lb3t>SUIk#vv%{nZ13eVSp# zqKpi_LreiUQ+B=SMS_=4`!#ne^vSQ^ z=c2N86JZVm>q?E$>C1BoN$sidSL4IZp~9xeU7IQj8?P&-OrW3h&tx$rQ3ZcCS-qC^ ziywmGWalb&Gq)-454Ww!)^o3tgjDk$d}NTQU$3O}U1ESKMi1Ah<} z3O3ON1#5!KdX(6*U>}ap_&$FnJGUJYaVHlyf|i(?%+8IUey)b`a-*Mc0J)~_< z12j5*55t(=N_E)J@jaXQ9E$MaRw64qT!Nv3AwT!TU_Zw~BogYzU3rRqT>kRLRPt-enpe$pb0n@|u5k9EXOp8WA3VUlAhWU>fBaiVUwy# zE;DsU(>%Hy^2KR_pRpe83njCJny!%B#B78@agJRg8%?KhJ~+0QJNZ#Ve0{Zu48ITQZnZ>@g36wVR$HM{ zIYrwa{dNlgl z9P+op(@>5Fs&Ri-a{2j(Cy=GjYz=;~N$KMvdD}wLFwoUzA@kT;(wj$EKxGyCYP?5- z!|&q~LqKOTlgxxUz+}^UTG#Z=nNOp`9inQ8dFJD7J_!v;yaZnGMH+os|7znshQF6R zd+>1aT4keY6*RiIfAZrfWj}V%65L>E(Vf3sJA)uo(>j0s{LYH1kti}m5T_$ZEwc!x zl+U;r$-%*;p><=X&P7dA(I|{?=8TO&(SzJ0$>S%+dEA{=8(v=3m^qiMB(~uRQzL2M zg8h|OiDvNG_?RWn5S(Tw?FCGlY-8t_rv)___n+m?j^S!m-4;O{PAi8Oa>?=k-oh6?nMv3voV}!yJ@s!mA5zQ18yE7*k9? zuGq)42Bhx`L~)?(d>!(mPDfBoLaZzw(sF12#`u10+s8Vij0%KD&WO7qQ`s$v)V7sO zJapZgL&zk5ZBJ(5Fs*h#@%3U4l3pc`N}!l#;sJkQrtH@+3wwk`=l*Gjk2C|jYp;>A zQsM1F5Vg^Hd785FP|E?;)0Ye+J=+`+nnZe9$oqvp;6;pb z9FTu%9&I5`Bs_?}Z5fnqKjpZidz4n|B@d5%nefb74nIkTluC85Z2$49VPcV^k84A@ zx|uK0zs}R%hpt$3Ku$Is1y3r*Qw9-@0TD_;4fO^ah`W;|?~(jGJoMgBF6}{(u+J2e z=7A@Cf$5Y*pVNctR|lH|QjHQy*e3joXW;!| z-aZg()yFP3nOcB3HNTFyyzq_Ty^n=kv>1EC`?pFiWFIcV_PX0686YJK#82Q=WhQ?$ z$5O^C@D)&mVRXJvJB8kV-yt2bKW@YMRrEeLSnFA0rl0=qf;ta|+&~+3yWYn#9-vY_ zC=A2?%8#LFe)vt)4#?$2Q^hg~df_Nqd!s)WM4at#X+oJ211T+I$KS zalDCj;X{znlefmG3Z=O{YO6G`BAtJC^?7%?e`)0BcRT5B!|gYj8NEoz-+YkCJRxYH~D|YnG$Z~L5MBD zqHRdZ)*4QtwtGKNA5*BR83|d+PZq)vl{E}OBcx34dyUv5d2&tv%O?*-OosP9T@W9~ zDh;l0isIP2PWd&2^H6Z)1JEJxDlIZ>}1YLxy5nCzF4dIs3AM3?831c z+)3o3A0Ulxs;)SW{6&8!nnckgX>o4)UMJ+(g()HVL;5FRg{h|Gx6r8mSW~yK8wnBV zQpk)?>cm)jx8mD1BCV|LWWMFaEZNhCm<`GW`$A92Yyiovta|v? zz}#qZAgimhM7atzXH{4OOD6%1)QY*fE|AQWqTv*NaXT);}DnP`bT zDf_}scFBLy7mq0zg&{fsziQsmS{OHwXpB59yw7Ps2h>Gqw3l*L+YVM|?Dn8=eaf*m zq`zA)Y@MH_luU29fh&mrES|s4{S{>3p-j)cgNp*3KR`|JzR3SL&w%UZObPb6Re)QuF+L2mB$fhtq z#o^h$kPDOEi1S}D{F>&y^s+XrF~SGA`}1@+S=VK(d!omDiHE`{ojj^J>{c?0;N)2k zueyJMwLHWtsK%|*R|IOA-@I=34_2aM2+kO4r3%_vSau za_xBvRX!?okM$fnz=0Y?BvFSpsTv#U%9WQ?MoHVFN%DB_a2LVVDc{j^O|P-@PHfjw z?)+iHH=OZ#)VpPL?Z#&Vaec2N6iF|WQ^J3t$(fUdz#A_*i`H8-vhT4TD(_5P)8Elf z)wDx56H4+uG#&)SY}?A$FmhBWnkj-I+pD3suJl7aw2h$L=_%mjXQ(k{wqizlk((f) z(4ns=d3fbtscy|q9Y2`-;5G{>h^WmQH%EEULm6-!BeMs6hpjoc;HoVVrC^j7DDHoh zmLkg?m6?P?#?0%kOjyk4bb3?FZPN3Qy{h*%ur7_Z0Wru7rYk2BCtOV6>mWTv{jtXF zK!^(xLCJi57`_uyI;(XtTv0s#s8UuhKL8E79u zMA56Ee#zNGkkWTbj~5;dnmh{sY_6!~3$Snc6(LKx25)lhlbAay`B{uEjs^+x4Mtq> zdGMWKyj;8PC>>-t%oDM@&1h=#WDFfw+tX{cHybM34QMISKerTB`FiPHlPG`nCF{;Y z)+wP|tX0PSRT1Y@AwK5N(%sRMhfasG7o4Yw%=-mCh>VhmDCwa7d6&$nL_@|bCNg#} z3h`4|kM5qsvh;jeX3RHBqNG#G4ENZ)lzqDaCFSA2l4^iZ-JhO6Aa5kx-Pp1suMWXX zHK|RA5R&4_J5?oQH0T@zrf+|_+bb=8r$%K-w2cg9j^pTC^qR31@BZm6G(n$DQ8PNC zs##Zzhwe@hN794ZqvpBsJowA)TA^@)wX|B{T%g0$QLFBeXm@V@j6^C_`wZ4rJB$6A ztX@d{S8~%;tj&czn_u@W1?|1iDr%pVqLNGqLv2>dlEgM`t5@6*}ekF;$MIr^9_DNXYEY44Eld8p^VKEg4IDNn{ zQ%sr6cDpphCB8IS`*z0*NXbBGQMYdu9xYGW=$k$kP?I_RjwF8&aDL4FwgyvzOuCs} zysUlhahFzmLqu}9?aD_JzB=By*QJri)(u;k09+;hG zVz`f336Z2R!|#7_JcNYGTx)oNmO4W&!q`^8s5#||M}#t@4v9MH$F=A~!x|jWMvU&H zmz_R>vhm28Y7q?+4!yw1&jD{R8ld+9p|5p7FV;i3Y&Tp!SPhT(HJFb&l0S2k#f%Ra zm}Wgx3q(ZiXdl4X%k3)sTo5*jeB^Z5I))nme5d$S<*68S~u?F?#xVn7)#vV zG(nX3b$Se(!r7N&5q#B!l@{7JlW9f~wA?wL#p$8v6S+9Cu$aw1u145@hJR41gN&D@ z-)bq~9nCo6?CyBMvvp2I+PKWPTqMg`)7-z0^D{jQCb^5yftIbu{_*KcppxQ(9qr&U z8&!ys1Ydt_BcQR{=LS3ITwC`fA}AkwrDF<9r7lYZCA{)TMX9x$%%KiHO>7x>VVR%j zv}m{W9z?ZrLGhfhEzFp}sUp}<;o>OqEh0PM?kAg0Nx(Q1)=lyoTKTNlr7t4ZjO=Hv z+B+h}oe=lXIpiOIwQ*SLn|SGgdOZt>c2grFi{;*fc^tCzah3gk z-#ZP%k?nyKFBZzm`0J56@&#K_?96vflvu7FmgiHRc2&(=-_(@j)`|=RvE>=MQJSh& zW4?9ax^YBGIZ$4FH6KNPYABuY$K1Gy{R*Px%249BAVH!6?6RZi-sU$uJGAAKVIr0t zZ&rWJyAQ2b3yu+Pa1g)kuiwaZoKF}bC0ZFYe}{2tVN?(6xs0Rp>~W;g6l5=b*6LPS zL&rmYI&rn5p93WM(j_VLe%ImgUX?n#Uo1y^x~I#iLVZ%qb&kM~!!3zkv%lL@La)X= zks!{b`6rmJce~i;7xAp%AZBj>oK|rS|acj-{vaSO|P%CN>7nXOqkG^+mGN2|q-I zU#+(O=?%jZ!)2hEjM)oHP%)ICLP-(wQZ)?PTM4uKQPXPLs#)p|u>?C@T5Aq>D)oPW zFn@#}){^3HG5detyF) z1`*ER!ZNLl{}H2*9B^E1V*{0SU)O6wb(hftNOP=$Wg)|nQsxm3pZCaB{Y4y#f$k+p zt?3zUBw^9UzLQ-*!SI>XTaWq4{=0voeWxg~O6v*!{9&m(Aq@dzp_nQEM_G2La>Oz# zqCQllyG%{q(_6W)0tI$6`Rep-))In`_QvXEj9by6)bdl6MFm+V4eGorB20OO>46-> z#T~eY0@8tywhftaZIF}UO%@Ow+guLbu|?0a@$u_Xb4D3T=%3^j?)eE&$T)vWh`jrv zBFQ6>JwZb319|2iz;u1ntjnh5pK)_oaxtKW$L`Q!1N61PK|cPc!q#&tF%%lIMonC= zNYgGAf|QMeZhpSHe#nXx%wlA(MQwak$hYr@3~r}dJoSZFw|?Ooa8)=zGpXtwlglM8 zh|qg~gVkV)bQjvSv`ZKm$~1pqj4)pJdm7O-5ti0kf<@GD7&<;;UUp@Ew;LY!?i;OK zAKp*eJS-7;0Q}R2cMfi=+A6yFqz~$pN`vAz*(Fpmnx{KUAPV$WSqe4DG|6xt8W0ER zrUMa9$f-8zo3^5>E7R`>ubY?V`R(nmmi_ zT(K08DR=p(1U{Y0GHq0&{1)+#(C`d|PO@C9MvI-(@ce*k3tjXCOMa@sd4)KnkO>=6 zTB7oDO4|HVQibp{fQ2RpsBpOK1rTazC|1R!yY);Bc{X>i<{rB0Po719$BLjVNq(#> z*LK`E7kd?7fhG{bfRcayd=3?%7>&b4%SGrPhq7>t#~P5Z&8he)2)AMkvbJ!>*hxBj zeLUe%m60sw{8xZla#E_bOAA(E=s_qx3 zn`lNeyqoy3i{w(DoXWR2b>451aqFtyd5!&>ud4%VaurJf#y@|}z1R~Jupq0lnD#f) zA$L4=EBM7GO5;*9SllR1l1HncYGA+Nvo=`mzua@d^uWRzpuGn?TWYaFS#%yF`1oQ$ zAY^@<%D8KbmDnOabF2ksa4KVd;kAK;W$XKn*$xW**lMY|!_?2gIItsfXI?qx&8sr~ zgGFVK0!PQ}aEX7(3?)F>-?OlY$((rv0$@u{mrQcWep`ODDRGZV&J~~L^ZqUBr?19# z4JoM*H6K?6_MY?34*8(iU@B-nRKP(ynX8|`*Vdw3qr2J_;+qb*l7o1yN^s-5&`({8 z>=|r=$el~nR3Ni(zd(Obv0T*}qIu9`H_oLncPV(&8;pM?nL6C^F3=*&2-1uu{mWM& zMfSPNtIM~taeU5hk~ur{L$@gN1>!RI5u909>!x*YSZ*_}7`aoUFZ?1hx009AbQml1 z?Er$(TX!c!B~3qMWYX;4zJ_Zey2pshey_IJd=0O-I}Fn ziLTv?U!|4m=rG~w>0j=sJsW=}cZntB3V!XPel8h2qLpy<_1)VKG4_BLV6XKmN>fv( zy@MpY3AHp)*L@X?_%OdS`zCLU7zp)ZIa@9!5%+&lxp1Q=)@j#x zlX5BK*+j}F4((DcB>QxPIxc@Z*YkPGQdk#=EA_yU|L3RuoLU;KnBw^`tX91=hzv8g z*ZNwY`l&2Eh~oeef~alJ`7NUQFqKzo_HUAntup4{PBbXJt2FecDzNrb<%7tF%48rTUJP4zI7rm!F0%+xpj{HvK8#fj z3(D|EV*7yH?(gPLUPYiFYc%o2$U;mClK~o5jbQgbg?dY&0vwjCNyf8D3in{e!aILp+S@l8e z?XZ6eiOll=ZX4{=wCVS99N1$#WJl(f9btV&dio)gkhr|=x%t?3%!e*XtuLvVZNg)r zSRWJ!3LL+UF%E~5EGrZ;lh+n*^ZqQh+}k?+p?Mu*_?6W3E#Ie@AZboZ-(i78X5oK( z5gSx0%x3r^l#qp=q`+KN9n6x`HZqRU(Y2u+Drok$PYBd$>rI z#Sr1H=mU}n)xbCqBDOVOah~qavbFQJm-s29V`;UYss74u0hiSVU>XL z8Xw*qE(OKV%ahX>4iq8W5sHU(-wUsCxHJs^`EU2P!}x_gABOj977xT>N*b!JT2*1{ zWRtQK1j7n$sEMU~I#9W{TQdxB82RO%82!$++cB)rSpy z6|cMr9?lA*{XpR_FRp(}0TssmMgy5h<*e+4$Xw%NwkRFA;v>0{Y)jE;@(d{n zs$0o!i9(ZC1mxoApKK0M(wSHirhVYe3UAD1i3f~f*_m2`a&;A z&WckHhIvulr8%++IDAHE9PZFX`^Ci1;f6Mel6yU%94*_Yq$JYFxHy27&D%3lDA<41F$bFRjz1Hsnd>EaZ?4VK0UH1 z6UC_ScsGTy1eSjU$*gcjC-|zq3&I-yK5S0elwDvQW^~)?9j%PAf;@8esH$%I)1V`N zSa46{+f0z4+o9p-l`jHACV`<;Qz&3={|et z`1-yHS_j{V-tb|5*QyZQ^F3A9rgAGD_2WN~$OAVDolBi9Oc1 zweD?R-qyl0mMtf@&-abl5uKjV_|E%*%GFhP?F)ZEuhM#@Qg6hP-%Kcq2>gV|dUNhl z^D4eY@OSYlEPr*&bzXCTOl&y**YuLaTDC=X;u=Ri+rDUSUfBo~Qr+BxW+pBl`ShzO zN5^yGcvH&8iQO5>f~cQRpUKJ_fo?39} zR-JzX|Mz33492mz+~g&`07P|}sZx!=iYEzyy@h-QRFvKJzH}qq5+W!)9i)VSbf`#4 zBM3-I3JgfMATS_O0s<0Bmr@c6NC{HX-6@U!U%z|5dw<_>*ZRX^&1S#n?0C*|_C7Of z-uJC5K01Y&JuNFydNzJ7e(tN4G?^Y&)b5il!)Jv1VgBmo$-_JHc)Bug)Vg0tR|=?~ zkG$6HHgg+%LQK!Qqq7;yH1ljSVj=(L_2U}+67!nv+><%~FUED+!5&bDpE4YA?4npw zj15FVY!6W%U8ZLl&dl4D)ICFae0DUB$iFrwnV+mHX+e)kaTfz!*q2zV~S?z?yiZ$39z4dj>OLzNW}f#lb)b9ZAH779 zecpw<_FPo%tcsajbRIpQTZQ7WZ7dDxa2zO#e9?EVys{c8%qu-V{DO&oUv>+ei7B`< zQD&q!_0(40T|y~tlD0qN2o!0#5nI}zVJ6%gk-NIerus{!e3Txg5S&$6J}NACIC1N~ zWXl$4;<(5~D28ETkz+;G-TZ#q=+96Se2M76HGLkc@$+&CbYNXUt*Rux+eAoxwTYz( z#4`X^_v=pUeHq+?I+HPNNdk zx;&#|38Lay$%XXEqIXNWSYFCXR9-h6SIa)!Wcwg~;67DhoX-{#-XTo{Wqg!m6MnD9 zn`r#@r3!Z2_y>7Eyak2kYS)b?b#4=KtnoiXX(Yfs9`u{=SP0Z~*snRNc=}>gYVJ=J z&g+zr6gfl+VRUaD`7EW2jv6I2&TSyFqm)0Y<(4XC{&aDmq?DBk_gQlf@9CIe^gSS! zX_OUotCJ-rGZ)87=LJ^2t7l{icB`f!iM?E9`zTB_R+^6gT~!rotJPEKwtrP7l25#; zluI*o^U>R(8F#7<{{2lh^*SyLb(Rw=&p37O_nzCDd*{kt=VEtDnqTE6-sM*ad@U8e zhGhCabG0Rc-w@`}k=T5?R`FQ)Ytoie{l`xV9tU}@yVO&bd-IApaZgc@PCoqG0suIMeYG;mQE$9}AtCOEDU z6I;l|(@M#YtqpebQE6kjAA~Vfq-CiqIeCBGk`LCNm>@wZ@X*mY!SD68r6`x*%5vtx zZx@uO<~&WU78=J)%6nZthgXEvT8e#HsX28q3_Z=hrORLNzJc0(y0k6L%L-|^UZw=&(l>^Vv?OyTJe=yIh}$WCo+sNMb1*^PGO!b zWFJqGf<5ct`A%W3hPJ#z_6zP^&SfgA`vS5a)voNX#N`4LFICS9VUZzg6O~?AT@J2$ z1yWLv)E5W?@9kmO?hc4_HAo8llGI5{*>1~tU2=M9RBNQ5Z*FcZgi<`lh;9`*_N9;i zc4hUi=Uazhsj4_JSt9>XOkXnnqN2miQsqQQWL>Ig2v|6c4q|?{STZar+HV^@SKcr+ zC}&U2AJ6v9I`q)fxcba=VAvMVYTBb@)g&VV%fg>!BJf6wnzf<7V#@8~gxj#qlOnF2 zp*&0Rb!~G6Q^ipu2imZW20f4OEg@=^Su&nWq9R+0J`EOPWz2&dm^+K5t->RFM31=k zn{qv>)q;QRxk|WV$htgJ<1~s5Nlt`yjY(W*Z%!%`-KMhS$5>vG$E5PHb*6N4N45+` zmF-xcQ2j`EIXd@FTE(OP=^nSb7=8H;#h#Y9d(`a!7GLD1*HN9((Dm=vS9}SU?p4i$ zr+HWJ61QE2{QMBY_W^^E;klJTRr(w#J3HP z&pIX0egR*C^QOMB3>oBA?xHWG87a{Y2MxdGb-0Z7#_l9@e;H^^sM|EvO5bXGK$+Hl zBE5@XcF(Jj)~T1saLgWN`O36G+8$SWp1Jjag(l}bt(Q$Vnr#&`UwpO@A9Qr`8vF2e z{dM#$c8!tm0fUP~r8D7o=yzQ;E2Gk5|TMYko)Oc|YWb6I5B z!D&_Z-~=k*^fr&{A8%;T*F@Q{n!uWt!}&gn1dw*{B~u0M4Bk;?aL0)4w6#97$Px)J zFc&W?$!sF~RN+U7Iog!f-y3ghwl#VxWi*uM`=+(qF6>ZjMdQw7>1x=U&+#wBZNK^Q zw-oeSHTWwk7l$Ba#1v`TZ_bNxK-*F|Io~W%+bsSUD(G$Xw8CmO1@ zY&xJt-;1mx@k`l${vgToM(qyabDB5i#ZmHkmv$GKJvX0_F=0fmA zT0c-E>{vnGywIjyO160RxZ9WRvaWjTw<_%I`zp#-)T^i3Zn@X9woE75$3D+y(tCP3 z?F8>#(_%?uJ(aw%_D$O{{2&nAtTBYiT^YQgtMb#b2EUubx9?`=R$)Hl0-MeY7Ty`5^MC*fkZ67yI~94E zMIc#%*71S#x0#%D)*&Ll-J8o+$AKGic72CjXjR`aVz)0JwgjP{q)Oj#leq`wr6!Y* z9g8FhO-mnZak42asA`um_Rdu^{%k4j&Ms-+CoMyCN;b}8JWW5om3&)eZF9THR+;S6 zfavyp^`^vc?$6tjf73z1L3?xw##b^aXAN(9ML6EOadR%xEU-5!E)e6T%H4B0b-Qpg z#Lu3{%?lhA`VMRw)a)L;D`N4(O+U~C)w6viol2AgG5&NOa;RB2BTf1#tUJ*>w$f^$ zq;Qfz?(k~v3rxl!E#gU-^6NX(UP3Oz(yRLqe)t3?D#uh6~4PAE}1uwG?N|;6??zS-q#x$FaWQjQ!g|_>~(32u{8gXv4dpVC@rQd)lWFT5zMOT0tjJ z3AD)DU~XeV+~Pdn_HLcId0o9giFuM)@80`5)wE5WtHgL)Hty?+x{PdD_68cKK7z6- z9nkCYBJ7wbe%^i$HWpeD4gr1qH}_U-pVsW+KUR{XlyC^{*2knyJjz3T_0~+;fJKf6f{*hbX_tvs${5*4RXrM*>esa z8lNCo;CQg&k?@$_gH+WZ$<{(p;6=Kv$$6d(#KJ^9OLD>e6S7T+wyPz5p*v{Y;=^xS z%s$+LNod;=`fyXx+P`->>emmQ`lJAjD|JDy?1qn`e~@Qo`fl4O1w}Sr^@?SN8%WRZ z-O62oVcIknV#>uk4Cw4=SCvlTS`J=rLM3&j+x8=uTHUkLiHc_Y6I5CY60c6cXNuo5 z5YTg34J*c_i|c+1OB{TgtMRL(vRP{QIaK1IO=3TVrT?c17bS`(^})Q1JE_C8@up#h zM%r?372MrF)YEv2mhA?|#?XIvq=Dp)P2??2PGA&AfzRqpJAJP^9#89I%~7)XnU9e_ zNtM7@=kTu~{p@rOYs)1N*_qp4v*H&B8TsND4xDd#=0xVO)Y-F^{~`_)mdS}C>6~7> zf6dtpMKZbfV!iC)D4tT>xJQzBWvEpZ{^Re9rGtcI>SEn54@vD z%y@!%N!0gH)YRkXb|?K_F}H5cqPRzgP+pxrC!I5T{JM4Tn0Gzh!fn{C_*1djZ6 z^!gb*K3Cs>;_ok1<$9cYSqg#F_V5^h8HD)aH~N=8^4sf$%a2FejmcRzq) z8sjAQxp*_fVOyzZRMRC>ME%lOSNKwtbF1B08%+Tz5qtb&kC_LZ`23dwnc79Nzvt}r zoe9yVQ)_-rEV{FQAV%%N65cacBKR*>3sLyl?q|W^G;eiYT zM#P|Y-6$qDf@t)d=`P}_?@^JZR$lEu|MY+bPmBv`Mk~ZOJZeS#ZHd&^{Y@ z`BPos!;hK4)z>D4?#bAAa0e=kQmf~_$hRNdw+*`GKB};zC0u-Sn%^*NF2EP`-Ra}G z29E;k_cx-7q25-Q_jo~V*(UP)Sr(&SA293~zje|?Ury4AEb1pb+TjcKa%{{PO@#uN812FYU9c6f=9AxFfz-VAl=+`(fE|yL@Jck(w z2{Ki2O(6AWVR9QJu^pKCDVfJR=&(LTBQhPkXLU}m?N5>^u^QNuyuWy(CHL?r1)h|N zzjJhVzbl(V)0dL)@BCXg-z}f5XMD)=v8dmz7fo%}tUEPi*>O1ao~GA+rcUZ!%+bQk zc&A3;HOEd4OM->aQ}`i&is)vXwbDxrOUAZjqFIqVlcuJnN~Eg>Y^sCpYmXCC_8xu4 z6h%W;-j#E73D@o8^mRv>_HY$3flWOs+4G5D^;7C%{)Ns0n(8Yc0^Ko%AVRzFw1gL@1V4zR@M^Vys&Kq}ctZEjO!Ih2HMx zsm{?Mj&rY%@4xhO|8%`U{`>RSo1t7E@A_*5le{>~kih3e(^$Lw>I%dW-`@}IfU ztdwju5U^&inw@eO)N3Y|RZty2?MmYe^u39yG!#04>B{wq<|kT=l(ksVtB+~1GKK3;rsiw2lt(=ck=8dcYa+E2 z)^+9&eNT49w774aZZfoEmI;#G;k-eim+PdejgU$q<#w&M?|6c;8--YPUWu@PkaAjq zcFi4eY~~522JMMdvu>@MVJA~+BGLVPmXmMKu={97_3ZlKsGfUQww%ynUI+>nEHSOHloZHv;9?Mvz1Hm zo{dRky%n}brCvV#@;q>Za9z!O!O&C_tkrqyEb*ZHCx;_mR53};4jaT&6qX#?uflja zjxH%O;7fgwc4&OH$j&Ezk21XvNA$OZf zHKd-5MVp-OIPaw&%=)G^jLLRf_QR~v59HOIZ0E0}Kx-47Ua*uiT02&V-HW7F{2a+5 zzY)h_G`qx4G8V1Hug}la)qRYaRCh&=-e*8A{ors@)*K7)~j5U+}qghb)DE0b$s4@Kr<9mWjVYdoBlsO?HjN`Qoww6(W; z`)AJX=lGLPO_f~Z7SKXFCWmoQzMQEthfTzThB4*5I={*3;R%i5(3+#Gbj)_VYN{_Z zgr>n>zaT0GCsLPs8^|%aZV3))t?cw(ZuptjH3wQQk8Y(2EIsG{YD)@gR+||feBC?W z7jwGQPAbO}CUrK*Rp%;M%5&0qeL2^cW=O8BGf%JmR_H6guNhm1?~U7**%Y`&py&5Y z*mSZ2Cgb@r8-45X`VuKMKOI}<*$ClQZ*SMFnT% zF{Os_trDFtEA8-!G$QG&W@?Wt^7i-L*?ICFT{=>55!a@`T0H512#--ptvz}z!&(W~ zYG6O^p+Ar2H<%1A_*Nq4owp)PW6%rfmary&dRow}D$<7!t}>7pO1DpTQEaHl60QsA z3*h2T?tcMdCt~}GvE1@bGlp_zDj7JlG`X=#5%}2t8HN8brmTST{P73v=w$n|JGestAKPjKc7_$W4JO=N?-z)a* zlS(}!xEn)!zJSq96x#105PSanZFUs}b~K~*Tgy9It2wT_oF7wxdjwFXtfRlYJE!Rg zM1MTVYk6oKs4O*q@71glH6?GMlNjrrSUqbwRBM5B*x(w$UH?91+cj3>^^$mwRTW7O z_#Lj4aWFa?oP7>D z?3KMYWDQz)83H`1z++9=`pqjTx7DmV{CP8z$IQ<^`|MA-+PGvrPhyT3ei9-w_)6k_ z@Z!8~lkTkYWHuepKGBYfTUmktR0@8JYUscy9ps z_)|1=N13jV!g~5FLM%S5a`XFV$1*>{O0Sq*#$?TJjxM{3F*mJCOSU_kevDgc4NA>f zAH5cpoBoh=MLZ$3Db;U2V|@v$A=GI&|FY50swlD4odBcLLpr^5UwTxNNIhv;Gj7Wx zF=*Ge6~u8ogZ*VZslP7?%ix_1ZPf%tEwlW_bZ~g=_*;=!CYHy0XKg#2AJU&2uPjbH z)}11HRw0l{y;)1jOmM!wvPsJKjjep9L7w+Lt){XVJfJ~(JQ@kbtO1(xqSsEHAlk`<<& zh4Gwe1MMCs@JVG{+DL1la;}!YzxBc??fC$+1{c5I_$PKV@X&&Xq4pglg`^$-<#XOh zs%xNt`(AeJKACl{H5MzYTrf(QtUYt8m@-rnnLhrLLbT$jT12!uRcL4V9pp<`wL@cC2iM7#jB zK5`3;VR*!x%yAl1sabR@h14;D{$^%%&2k}hl1(6=^_{JH{&k0l?+iMa()W2Hd+Ghb z#r18T4xAj~PxZUqQNe=I_Ai9WVaKzQk3YT~n-UGuYJqTlf0wqdqd0|B3zEg75g!$4 z+csElbRFCBn0KWR=Q6-JjDvwBm*qRpaf*txa2~Y7-wl*ApXQCw z^VQttyMsC%Mp9G`cB}-wdvkuT@uAP?1V2~@Zy6rL%LS_NefrEYa%IHn!=~;XO$FSe z_LJCzvV|KhiNzKWXs?R#RV`;W8gFk`M|RJqd?**YxH}ucRgI@~eJZ5^ zse*9B{0DSmgpPv0l*Gr{TRYrE>|j2a?Z7)&%Ee^*j+eXyaShPNmkRC^dmF8cWm3p8 zUf4IA7VLm53uBEi(pSac_BXw(jf53CKOc~B20Qmox=VwiIVw36gm?d;;*U;>>9(iy8!LNox=X1M_)2M9zWW@ZSS?i_(ZsZP zX<4UkQFG9oLcN<)-+04n8!`3e2}g}lukSC+3J7<)$1fqmMd5~ivLm@I(eK}c z@!=Ib^KKqdkiwt#x^!KHadcmmu_b%c{{?npn~3@Is?8N7`AC?5FnlV=s+~@}pw4I5 zHgk{9Cs<&-(7i2e3W}FUV)_uaZ{|q7n7Bb4q82~-6p#{esWw;727jTb{1ao!f_`DTGr}b5ynY9Rp)>%{X zZeLJTe2V^=Y{XT=aWum*K0XgJT!}c?Q*49oYt-R68eO~uuNa!|%JOW-NSJJ04pl?iCkeP1l9m(Yo4PvSLThVBbCSb@y1UR_Q}H5>eKiQt$2N~2o( zZ>;wBkekQfJRmbWYT2Z|%$4j)uZ_K=;>+K16w5{G(8heuI6Y>3!0G2WNfI4gFy|m( zpmX+1|NS=&&f2HS{RIp1S(qgEg8hEE(1S|yQ24`i^Lk3W&U0^OtJgLX&&RJ>JeG<* zE132Cu<`kr?B#&|8u{rfyUZI${$t`Zh08%jVdV6LGmtoXzUR=l1hqGsk0KOb{t}QA ze91HefBHl{`HbV$(`9QH;34PKJ=d>3Y*#jk7rpG`%zW&;K}d4%V21bsG1=F)rC|CO zJ1g|qs3Lq(Vs%Xobv<=#QxX0ik)4^7lUZn4F+T7|P+e9`n9KO@>($1e2;T^hKtX62 z0*c0d<-%hlgMt43S7d{5ng9ugVre!B)d>(VDAx3k0EeTo@BaueC?aVKc!eV2SS>)n z(rytlU`w|MB?ysVBnrE|MW_ZyAQ<-M_8$oZMq;0D0~63t5O(^HfC3?~pdCO!LLpdx zKmhgwAw70_hfoSTa{x%7T|y}W7~s%p_m3kK9NV-Dm_k4hU@XnvUlIs}HT_EhAraX3 z|B}G)r2RjBU|_7)9{~l1V>AB&gC{YwI};XeGM1A)OSlAgM9{~~#!y26tUL%BnPzY?$9~A=_^cp z)Xs!L|I<7Oj(~$9*lt21aeO2UOG-p!jt{JwJrR*91_FWLj)tL;C>$t22q+kiI}VA0 zA|SY)KyWb5l7V1o1Ppf^41omwv4>S9Cc2D)fi9QHe{>hZhl8L{ z)CC_ffOCHc2YmVS-_!#p11RDQL?OU1Ts+_i6zD=Q5I7WPu>pI){>8BeN8tFt05Jed zgunq93Wfq(>B3}S?_SiQE@X#5{*MiS09{A~iUwnFb_xg%gTO9i0)wD0Oom;cH!v5C z#0?q+(C0$XD4HY!1u!jKH zjiW;W9Rb%K3IfGka11D67nTG9xWP3D1tD=aBnkm!f;$cY2jMOW3W5GNwf|2gKmqOc zf@35Cw;=$v909*z5Q#?MqJTmHw80Gq1^OpH{@{;7A(0odL&Gj0g9d>wlpr(!0xr~O zFdT>PAmE|jaNG_ELPH>63~qz~1tGW)0=?}*2|~jl7oY}+j=%{943sq7<_GAIfA5I@ zR01#r4MpOVAV7!2%>)coNF+{xfIWa+|AZP0LjZ@03pz9ofdT%(VHem6=peY;35#($qVK|ohPbN^NU z{sSu%go5IBbU+8kEj~cWgF$hf0E|Q7upc61QXm1~1eF zKnJV|_7yFWtSJ}EQ@`_+MNEQs1=9c{b4l(AImL_+7 ZV&?M1)5Xe~915g`M3eLJ-O^Mf|9^JdjXD4T delta 27804 zcmV(sK<&THrWC!T6tIl80y8m}0R;mnf9+dKa~sDIzVlaX4yjs{g2}ueR_P!+lB=Ak z>_m*5(&GS$OL9d3jPa0U9g@GEp4r*O%%T?{M3`bFxri{f(KFMJufOhD)8e9_#rG#P z{{8m!*DmhaGE%33L_Lzf)}S}#jEl)f*#`+X{Tw2$%vEFv9xhhFBp$4Ku1;NKpa*u&9hC!edGd6UEpNr+B_*x+%tpq|PsG zOe-~Jlu~~R6@QS(3reJ91>-~u*ZklYFLmsfmZoXAU(V-(^02GxOQpd|L}Q ze&xg(+sN=n#wLn1(H6Xw2t7cS^Pj}N5IWcf@+;aP+x?gcPN+8n?9j=tkWKPlXA{?i zDtVl2VuEm{|DR1q&n9tyvPm-HoVSA}oW^we^mN-?5DU`E_Mn(9Nt|V^f57u?4US0F zJ^$trqQS?T?X3L?F_0n}a()cwchQ{ZkUWf+^}Is%?0xKtMEcb}4a z1X_Il26X}jze@`rz!@qg+-$NKjkztaMq|Os>TFxF(BRx|=cqv`QINp2uW?Himtgr0 zDt(9|+5{}uMZjQ;r~_hTf6U8=;6H83**bjC7pvASX{-5CHKX`Ybw}w?X$!>UU0vGxc8SRC@A3_?JHz-^l2mJnngsjoT1ty ztaH7RmN5nh6)r&0vDO1?sR+`TvmN-4d!VvAyu-of+_9T@hg+hh*-_j(U8#Ha7Pe1F z>_L&DNa*FQPg*{Ie?=)8mwZ=!6)|{y6%9gXOT;CPSYUdgF35bbJmO3&2l60Y_~{LJ zA`yyk$3JaFrrFO+Qk&5VeitHT#pXYmjkE&!4GKtV0jnZW* zuTXP9uw>Q%;I*yBG~|J}ORz8iD+lwK=+&P`R+QIo{jbWZf5yVm5zP{~;1sLWViS>P zbp_fuFJG>wUes$4kk)_zvW=yGhX;^3o3AFzWqolqf-+gvtL58JFsrkAezBfkH2T)_ zGqlv|z}@2Nt6!mP>++unD*-q&RBmnaTPQ>CMsWPgxnCWYTTT|M=2_gva#hb}{^5$C z?1>Vsx@wP)Q z?T_FlU;Cgq7^HdMjmCzTKL$kmG30M_{FBjGnes3GkNtvRa6n=W77RW9YoX^9qR)SI zTJ#W9Om#gMOwQT;b=Ob3e&vOrj>mW(XSp(qGJ3%xEF1Hap zesz5@t>1Z+yli(F*L(bSI$59hw*F#d!SIXo3OdLIfNm@NdA*#>W=&%NMr@Z(bL;aO zoFqhf4mG2iqmT@AP|EnLuK(v8=M$S0=V6ne^i7_#(q4`HD;olrs3s<+Jn)lA%( zm;7SNPdIOyYdO~)xK%$@3HJ;ilq_Cf6$O| z8;XVaSV9Tf<>T11<1s$S7>D7++)p@RyWqr+y?b$vO((;rWQwGaZ5coe-U}nlmQ>j^ zf^HAbw}y?|!Si-VA+Ujw=6(PH;K6_c5?ypB zeYR~0-=F%<{H|z@pJoRFmxople|Cg3fxvCE_~FiZfWoP}iat-_6mqr!ZW?zqt5*$k zn$`1q=GjwIP5wMtoYnKm6nZp}2$ckm>nDkOE6E=kh0E*6&7fgI@4*Zu+~0q4->7U` zI^8xWQYueJ+OGW{c>Y=8BV2@GMK+%i873IhS(8TzoecA#IvLJ=*^s91f46k#&xv-n z=Qf=`SBTqE6=OdYhCTaxd^hYi2mmO*Oou8Qq_8UOFO<{UGOjV?@kLsOWXjMywIuZ$ zQbr?Hyni(fWgMKNJDB1UaUX=H6novLiRI={_n|j7Lc1R1q;9dW8+A9{)r3*G*Sjp- z2j2(pW|M0RyN3lc_bMtBe?noUpYHqDVWg*QY$N^QrmcT~jIND)Hi0HC4wO8GJegL0 zI^&1o71~Q%FrnF3X4u_9FP)}D2s3DA_-O*?@sCc^bpFT#qk^OHImO%!wYC4@{^>!<1Us!Nvc zwY87oVry~MJP_slI!vb9acP)lUQaKo`B~ih^XAc_24&6j$$EKFoi_tn;pgjE{fBS3 ze>G+E{rYmT>N`zle{>X-8OHjvLS}Yw*z9_FyZ_BcKh{T|{qLLn>?E!nC}0LD14Et7 zY#E>G*!@GD;f{B(Ke;RLx(#%J38_&z_wQA^eke1`<0~RQXa*yB_)Uqbc~NT-ZdI;5 znk*1Zg_5oM*Ld5#2vhJU*53tknFjUQlMD{ro6+5`VXN@Mf8mZ7M_XM=!Vg*d0)oFl z@c-(8ml2YI%_gsUZ;ANCqJs+a*^Ga!NP)ZD?%{;+In+rdFR*NR#lLTkUJ zJwy9~I?Jm(@b>=zSo%2ogV@2h*ueq%m=QBF3NK7$ZfA68G9WQIG&3@nakT^$1TZ)> zHj^QKCx49h1yCGH)5i4h-Q9!R;;^{8ySux)LvRloJh*$XKyb+2a~{cg z-tVirx2OvHo9>?Oe^1XWp&(UKr582@83QFicFy$73`{%#QF(c3W+ng=6DtD~6AKar zg_@+-f&fqmMBY$VGqddqCAmd^KU}gm{bMP>8@-Q(0 zSeTf&|KkX9yTN(or1whTh(&=9YRgjspn~@_B0A8@MGy&Q< zfkRyEOo5I7aCU&Iv>ZUu9%%P3XSshl&;kD48~`%|^MAMdcl2KyS=#*_Y-C~rvb8s| z^RTot2bfvf00D{;atzMy&U64HJJY`$jel&MK;U>IS0hUsBV(|`-=!M?B!ra#Mqq{i z*5_p6Xld{4#NcFU^Or`(zqSFlS=`Q46l7})v~zYs`l~)MOGlsyxbGf}|2|x6JCK{5 z*MG>&($3WEFBzsT_KfOwmJTjJX|aF%fR~W|teFFy0qjg%oUB~z0H6Z^=x$=c_-$v%Ymcbf)yzD{t05h-*ppT^)5c~_t%gM+U2yk|E0s46T)9}B9#LNsZwKQ=C z7z52M?U4TL4qgVD{foip?`Y``&|w0PA2WdIub=;Z(gRPJDag*oX99i{Ee zK!Cqe{G|fizW*q3{r6~4{rlcf1O7LK0th_NKmgUB`PXG)XEFi*!~FmG7W^IZ|8MyJ z+U5Us{{L@95-v72f7_}43FKt? z*TVyjVq*F~Ht^juu?9abPT({9m&wQpd=Z`hKHk3=2)=6nTZg!v3CQ%XJITV%0WflO zH1a?Kj|7<50bb1D%WDdB|9^Yt0gMcGAZKt00Ib;uU+31g%8^9?3U*co}FiQL(u(#wNVgWEp{UKHWqx2tQ12D?`Ar1he z+#dp$BmalM?(f2~ zFarM|I5+Tr%wU0L|3I)p^FR0(8T~UK2iV>GuP-(5K_QM z!1^BumT&W4WC0gzV}E38Z2FIkV1YJ&P4%A;a2IUEX*4T1i~T=Nh6SwI{tvSLwXf4Z zn()`b{OK4gxCf3FAm)GA!SR12a)7IJvN3Y{M^l->U3C7F4u4#cv(Z0v{lymlbFhD% z^go4zJzW2R;2ygDK^Cwa_kSQbzsElitj6;Xf-CX-PpbcUPfT1K9l;OI-(Lja$L2r$ z_ZL18=ngbNT3G^_@CI8o26sKw3lq4}A58PXUZy{2vCvN!>mxW1i_aPCB&z0QugX?R ztPC+r=v(zFh<`Xdqew0WPX(8>v=>W>PLaSc4UA`f*&7v~)# z6k~(0(^xNj_mWrUdY1XY%UBXa=*)DuYKuDcp;6(d#s1GTJyNRVK28h6WV2}IAD)jy z==STFBI%C#MTe25Yq?=0&~2jvN%u|~gBg}m1DiFwMAkB4A+jE;3xl)UFD|la(mFFE z!BlkeiGTcW)%!nuS}3|>T+KqEIzQE?qVJ2-VrNsI`j`b3j*#2rtlv?EENHUmaHU&x zsR@*>a|d_&YoW8UqPkk9*Q&o`Y&jZ*FZq1u<;nI|c+=A(Lb>Wt7iAaoA=-1hY27=s zxZvrMnel70Nk&7osHdCl&YgFi`;MO6s4b|(r+>z?%a;9odl@}88pH!2$*nl;jSP;# zW$29fShKb_6#k328YghHBdujyK1;L8sS{@MnbmsYb29*g=K(cT4gsoJ< z;?0g3)Ioze=G8-6h~2sWfg61%#M0$p zMt?qeUs342g%3T(8dWT6Ag1zI&`DxkbeHh{-^RZ|23S`jgc`+LPa}a7@U1^JKa!l|LG8EZWDRT^{X79_OE@Ve zPzNU$Q(T&|{OoXxR`^mgFD39x%;(n?!hc#$!tm#mcU|Yr4a`>BwL#|P@8Eiwje9z> zWh#xl4SOtN%8dP}``tbVU7F}nGFc^GHEKa{dOVmTVWeB&skVg_Qo{QAe6W9`pF1^M z(gzqlFBu>W6WFkMYMk60+o{}a_YxVz0EOzBXrMg3Qxj4&Eg9OlC#hpcGqOZNRexi* zoLFb#;;i%O^rMXvGp{@+IuMMICqmc^MBw&LFjDZeGGJ=;+8qT!-Et?oCWN- zL}L7j7}0r5MA(mu?eTox{>~>jKg{`XRjV+l^CPGKPY37VZL=qzz7L19$<+0sTrjFC z^I9n$I3M`3otjy!QWW*Oc7K98(J1;*;?#3GL$=|2pyHQS84$W`BMZnW`lV!78oYIHzkiN!!)U>mKPzITqiGq6bY=Hps4U?w z#i`Hx$YIpgjQ$Sz`F|wSCmx!kg2hKRi0Q&WqaJ(uOEmo2Bx4yzpYKEj(@$%W#PE?f zn+pNA#7(VJzFn!$G4&%b?yeuYF?{^@CFsesK4uwXl2WyJT8k;G7(4SejlMgvDsCUj zu}0okz`J9snTKP0%(&fZTK{b&Op|43J{P83+0ckcReHob5r0sxjH*r&ENX9jJSD{Y zFenPa4LqTg(4mS#DtrlfjDIW@oW_lDgor#(@9H$Ttqp-kpEFZL(A@k6wO5{3?()dA zai0)l&6rqAS;zhh z`{d$#S|18V6Mxjx>hH4l0PUCe zK%HL3{vm;9N5F;wHImbMUAT6>Cn0!E!1{6 zf}Ob<$AL1fc*Kuzv?OB$)i@so@&S`7jDhvs-`ofqT!d z>{h7w#(#;;{t`|`z>~5}SEz1)EvF=y&@wuXqaa#)s49QAXnDDax+;+J4t-?iKKuaX z;F^_Y*@{eaFh}f^>o|)1=e+mNT+%j)$gT0o$o4n!oAUkH9zvPoZU#iD{xcG65*oRkTD>O%g{T#1`$sNBY-T_sjsc25W zn!3~<{k~r6#QqTXoskqA@kaz*&*IksG{(7k;$U$r4Op)n=7<~S$so0c6b-r8Id%_C zS$}U&|GAsJv@7q`R%iHPD-qL^5DorX&lT3yKD8wo->yTQ#-5h9MOMRC*ZenoJI5e( z6nPEJ%B(wvBaF}V&c_+)C-pUrES@=w*9pNT5MfAssvOSDEjB|NhF^>`4Y<5m#C0m@ zp{rlpe*AiO98@0Q*`0b2=U5B-{&aheyhNlDSq_RZGX99CT?G3VU}x;zq|%ONUC+o&h7ITS>jsTs12b_ zgJR)Nhs%B1o%4Mk5a`HNfQYsWVmi<@F?})iz7toqR}si?TI)%|uUxjmo-4nOvSCbn z6nH=GXZYg6%_f%j*l%^+rD3(>iGNC#zVlB~_?A~6JGAA^3y_7s;6UpW-nw5HH8kAs zvEpGwhQfi+Zi@5&{ia?@LG#2xp!@x%Rpf9dX*Gdw{PaL6n%o}Rs>Z)>Oz!^S~5a6i!u!z3L~5kWPCsMc7WbqTT*xPY$5G zMl&jV-%KuEnJ6I@V+%pOg|k`A3R~^GrE*}IWkb$H;ciX(Dw@IL&wX-aFiro($J$qJ zR_N%wdH{7ol7C3Q%HsPewSRS8R)-}bm+^eUgOHOx&V#$;Z-#_Wg!gB-iua+J+EEu) zb$b$+YU^F!Q0d5xE`#t@*6N>}%#xdniofgH`kPoYgBoZL{Sd+?54T|a`^8Vgj4z|; zL~9*x1}Q^dwJ^@4i?T)*(qSCCb=O)40>A#=eO3X23|$jIe!W*yblF)1l- zu|Rn&h!qq9D8ksTC4VqX$=MU)jv`G5=?cju8L6e|sW>D5c<U;*@1x5F!GAusn%O7x*Ddlp!qc8k zhEB#l@ww3}uY*qQLYM|gPF;*pX}Nn6;l07{6%)k=s?uyOqBbpbB@K-s?~I7KiZk2D zsWyL(>8rD1)H!4}FZ{aPyAg%`SoPbrH>&-Y)f+lQj9_}+)B1!k(R|w>J43k>PmkUx z0k0HGZAf>cpnpwobE667hy6DQrIPzwQXoRn+9|y=bL@9Ue8$lPW4_Ror1e+pXT}4s zEx41}WKh<0K1Oid)Clz6j}k6P3iJs3>8`J5BW?26Jo70}cCQ2ECaV5t1ni`hiMgFj zDUrLYLLsE7t-kwWi1TDzyeNJJKbwEwIKdF*>7Ul~n}3U56(%$XYzRIcpDlaB^TrtH zC3QLa0?hzYNp%s7&M(M=H^&^J6`Ydhu^s7;UkP~LKZ@{klmCLh9zDY)&!TDLNci?c;cP^oaQh?xvB5HJ#(L($32I9&T z3NIJXw&8lL&$fq5x83`tEyVdpgUZ57^b++P0cweTKB%1+(t*$CSu7iPwPdIzVK5nQ zr^7mk^=U|pi6$edP9S~yR?)R@-@1LsZKmMXHRXt zT)Sfj6`X1ECa*J1&>+Mz8Yb}*kqvmAtrs@wvvN+^h$za+T45#1A1^cTP#!kU8}T8b zg%$ItYD1mAnl(^Rtv}mOi>*n0WTAvjr?-zpmc8_?k1*c5WhLcz5x-GBNlrtcIuqmL zZGW<=9~y6a+FFr5cw}+HBr0s(dKPXr50y-sQcxX@G#O!8cpvB$Y68mliV zlD`vmT#w1W;TS(_nEDel+s2`WibZduH$wqPkRyAX%+y5mL$GYlo;mgY8d*Ok{pUMO zEObL9U)A3v{IVOr1;rl^WeKvfD4S?+B7ce;b;(%}$#s^a8QqucS$iPPc0a2ZR67L0 zw$MNs<~sr^?uC(SlGM3IMu$Q1H1}AWm;QAUQGX^g zYIyZES&Lf5{cg9G~%HN_JicB!ALg9&R`- zF)u%A7L4SYpPwfe8mFqGJ!NPkd4%<^{VAb6gqG6uC?k3Q@PSkuOtSGtp!ak7t#REh zaRxGqvs@V0Yy(r;c)0h>{ArIOG4+59j3=s@I zIuL{X5(AMys1tkZA(p(-3va7gYCx7~9A)%7)!U64UD|N_ea)byk&(C2h*LK`gG`Q^ z{YTa@VIj&x525aJ?A`f5MmV?-sJ>JX<6WHdEIzV3a>~sw+DU90&VSZRsTqR@WxiWn zW>+Ug|BfIv$fh;xMfNmPEvKhcsn~`fQ)^9o&I1Pzl{vmq;vIKn{#*KN*Kn8@$r((a z`p+E~o+cR!w{DTa2x}6u5sSn*_K9pX9mDwG*xnxG z#|`lH)WS1-Kc&0Y5j_hkdmveDhkVN|-bwcBM^5*quQ!oHg9*NER)ay$yYTcKLq~_o z{!JkFgdfVeP*NjFtKBz|ClcMT71?{r-;O{-IT@hFS;gh+7k?H{mOi&V(rc44z(s=9 zO42yo-EJZC+*aD3PgqE074zHpfCii2+c}zm&SWl$>BBISP1|LC^OO^xMyDG@^(gaP z@_hjb4N05?UeM1p`trfGrbl!?PY`?1Sjl=-lW8?Hnz&!m^El-oX5e?Yk+R|kKe-MD zL8j&n`o+CfRewWK6o^1h2cTMJF?JcBaS5`$y>nyR=3KqAnx>**DB;`{8-t=dxqG7f z5&CuPgH}6Ue)WVom#ieF!39$jX}}Wb*0WSI=xTDpl4lf7vy1j8Oqy&{SMSS`8jRaf zg_A>=T6K>_Acy1X=});NN*>M6kwaY>34lHG%HkF|Q-6Wzh%hKw_$*YRXVyjC3snUk zu7V=$Pv$TurCSK?c6v)Qu)FjdDk~M;F9lK?URR_kE04CGP`$w4vn8?PrZnl71#Jcg zV#;sTfx?Lp*d=KU1*cRbR*hiS(n^qThB%I3* zoqf=eL-J^MKx`K2Yb75P`h*ug&T&GjdA5x-9e@8Mj@3FM-EqnBK=&-I)=wT5^ET~~ zwGwub1SyqjZ`pD9Tf@X6R}a^Qa&4L7o_zn zF*iv6@RK?phFt&iyPXDa%Q%2a#fUHr`)_`9Mf1}sQ9C0pPnv3$8Q>d7@jATTLLhOr z{f!A_N;IUjj2(YBW1*}k*6iE@LN8t#<0_Qqpm%N3MwRKjYp?sWgMZ(@ z9_`vm_ZaNJXJ+&xqfB|Dkh_}&6p{=Gb+^p*6x2Pfm=fJeo-O_2bVs&QzvHSyQ4yW* z6)YA`@k7+&9*B9_=pjSs{t~1Np0dIDW$)I$rFz|HO{BPG6x}wx!D58CXMeT@Mr@S~ zDe`dZb7u@!0WQ5hlFD0UI_ABeX@9Ah>fipe+=hU}c$n&gYQ=OR3GK;M9I8?P z7oKk|h)L_s513^#)r`in1qEg8VXHf;JD=7u5opsrafZ>_ck0J>cbsvI#(#-5|TDzT5@e_Rvz>*?Zt6ex?4ELJnOAljcNqFHBhw)n$6_9)7E<&e{Z1TNN1D0!Pop zByy1$GPSO++b8@ZbuEvLFGoBZ!qu5%rixwkY0}%>t|)?o2rX>_s*>HbYtmJtq_=OX z_U?-xWX&Dv76ZI~hI;8Scts~zdE4(7#Mrq$ zqryi6+2sx$n!IewZP476Q52Hau;6m#x$Gt~8tnviN_Cvq^quXxU4P{o-O$x%aTo`j zu6Zq6)7PEy`)7YHrkUA)K$%p21jX9$Rd6u!xN{YCPz%cvP@f~ycbV6Wu5EjycD<2O zmXTDi7)zj(DH3l_>`p@=)Z?d*Q0-)^V{dCSc~ng~$%NztC7?LfX!}AjblK4r*~xtZ zL{LFV1y%#p1Qyd?t@Mn<>p$uD2Fs6c<@lT#IWI>+YhVY z3t28?y+YSW&CaX%#&1)vKW`NQp)cZDzS4KQBx8w?O75aslR-Cl_d}FCKc9RFuogL& zxM!lTmG^$LC_IBx*s8J1?Q%rn?)&N)r7W$AGfKCaXZ5IDL4W+p4kZ@!pm$`&>+3x? z9s;fJr}*J(S3-Xc8X7eq&?L(KGOU;ACv&by=~$u2qzO5f3gH|7S0N#w*$+%9`c-Jk zbqc0);XWjE7jzi{$PWwVCL%FqR~TR7(%QqY%2s(%;p;yUfMykxKS;mXAMV$D3X|F- zGNw>KDU*L)xPMjsG0&tCDS>|y?G%Jt?{7vFE3|km9$*^rNjYUq2Jw9@1~KBz;6^yT zW&;~vI(n#__qfA=Le*r6f-E`r2| zmXXLSZev_ggw_cQtuxNkWF*`l?iyBx6S9Jl_knRjOy1lI*cA7N;-gqU&m*!REd@SN z6l@_GK{Xul_b^kF6{K$SINy_?+qsD^h{g?^U2uPX(LUeW-;jKSjqF3NhO5ThYGzPZ zYJbx3SJ_}N)eQNbNU=o7m0i!qi_^|d$lyg?U5$U-v7E}r8qN3PEVi`aaa)pl>+L+K zo}+lw%nL5@hR>~395#9pQFLmtiN@WHN4uSjBlaOH;7Wz7#x!EMq&mlF>e{}HEq?sn zH|&BCCiu!lWrd_kc%(?d&Lm*2S_7i0b$@l(Z$Uk|N5+T_)gXNpaYCgn7U4DEICsDW ztFqOHnujA-R1f}C`?KHe?d!?IaD`+yXJ`4^lLHt(QhsLB>4r)I-hlZ$gus}bmDkPQh)ip zhW$@-Tc4W)m%L=UjV1F&K2)z|Cd!0RPqNM}N4by}M+#(kXDTTr0)l!H5KlJuR7FEx(%urD~P2 zrs>0q?;|pzgUB*7$8r15bod&6iS7}15y3rWp)cMiN5awSDzzmKyKnwyi&5n^S6>4T zCfE9|M%7y0-F)ESRY#GNr1{=XzOvTeIjQq0cPb%9gNed;bMy7|MQ(a=9Dnj)wX#rR zsDI8u*awqV%uicwH9|ak+-6=()GKa5upn72<;2UfWeoY1*Q44kn&F456C159RjbNH z=tcD-v={WYo;4P0^;eChTgY)rn}d+#l9&(TPZBT?^j2qM`|8Kh@vhPLYIN{nAcjT% ziv1#uutu0_NnCm4yNuN;kAET(RnQ#jsWVo=rDbZ!Is`uHRb;p97sIflriMH$IZZsc zCxn+o{SV^p-!(1R$W1ZIf_j zYRx$iSIEQB=~`?AkgkX4rl!Mt@490=;(<;&h|mB03aEu1#~2Dr&@vBF2YT zt6*=iB8o5&GUlXowtpBcc-knd>;*aLG=+@qn5w>o&SWFPa~d<-*U$zEX;b_jyW{Zk ze=Y^4o(#YH{cDDE<3i_U*LC1w@xZ-od#{<0?Q-HuwNoxL-)+dBTnU(TaC`9WH`UM1 zlK%aGyI(>bV%JW`6VceIAl=^UcGY|s|14>Ujy`y(jiLb=Nq?5!C(x1E9D9{lZ@&0x z&ipl|0dp3PzKA13uydAqQvAbXXR~VHtyk-7A-q{?DCP^r5_VhWgOo}H@K9f!vrf}M zi&~gIY7lxK35DGtV=2hqST9p@_lN4J`<+`rZJJD=yf+hmSJbXGQ(U;!W-m0cDfRix z{CFdMKISsfG=IXG?_LqtRrC9cI_aOE*&C~BA#dB}j`i3#ttc;Of{Z>FF>{VL6{ z>oS{YX8+{h&M`--V~sV<;**Dduj;Q+nG*cjtCMoMv=f%bk&fXT502G^usGqHom3VC zv41MT)1*x;u4+y_z&8$7W&-WpY{#ivb&f4xah}QhVK^K2s(m#fP@VO{l(k ziuRwBE`f&nAXjg%WyF<(iP{VViHOix=!yA1V>AK`UJO${8q^P;#L?qM2zhS zdmHqGe6F_k`D5|BuFxA-@5r75lR@D@O|yl&ZGWpTuB2r%%7P|dqpjmAPrdjnYx@yZ zYT26Cm$kPK%8Nr(qv5!68>3h0B<#dcGgv)Fb`mFk^lI#@(leYf0=t-A$|o;0M%e8c zLlaBg#9<-2&yllBL2ll{uW!ocs2PS`aFK+p(7mM}D4Uf#siA)djb3oWz8g}r|7wQ6 zNPm8?fEC;NBwxY_di9nQo*@z$aYIW2jazc?Zb$%oD3*`Zaid5L-DskdL+kX*imn#C5z|l{zY(fx2&EclZFf%BG}Td~pEvqPF?!KaESGY>^R)_@0|y)j(r%;(H{QsDD!o zk@kMPc)VrR`%u-7nqO^B;YZ5^IlH~Cy13xP1_wl}%AkO5uefud#6sh2;Yf(SJ^D>piL}c_SD%W6Vq)G(}j}x(4wK1Yr|hB{#%uI z>fqKM`CeXzlh;C~Y8`nW!nHfCriZmz%lD_M5Msr^RtP_mgagjtxUBK}wtsKZ&qqc& zN9jWwscaS2dFlbCcP!YrZR%e>e763Ml3vQ!gO$cAB)83qY~>cuxGlNbs?bS%85(+! zC~!&Vm$)1|^CHCqjbrHM>eQ=IbpR!H)HZP?eR}w_vf*S!Q~bV6Ugue0O_3C;==&Sl zJ^{iN*|E#)#OURX2ND*HbAMSg?%k(KJq@0+GuE{p{QlJ&nBP79rw^-qZ~N*(_e3gF ziYpw=He@u?Y%WM2nL;BT(W*jnD@fJ-aa{#ZZ~|e|B$kh(A>eio49EsjmS1ED$nGtq zI3dQUDA2@Hf6&ntR&{Uczgm)HnAk%!y<$GlG}j88Dit1$WNf;7M1PYypR-oAvS^3) z&P9JqgP&cSJiR4lyM3SHmmuIB+WTvhO-6vy6}R&0hQGAy`K?=zl~tX*w4|eS@3?W- z3pVjSUtu#Cv%-q4c<(07p59^AgZS%su2Y~3iEz;2aC6QA*||-fPXNBE8am25$kk}=)2O`cp}U9Plt!bVG&o7) zKF+y3Wt$xg&~*T}k`<=IjRx=HFs|yw+xk)O2o;bn6?x+gNgX&S8>oLpy+b+KOyJN< z^){D;kvibqp1|?3J;PFty|MQDA+oBmDd(>3JEako3Qy``JAdjF+^YG-7R0?vOQx^= zkx5B+G5o*nQplU-!bUy?&!0g&i9W^K_6Jd_$XW7qGAHWGRumvrUD~cy+bd_{%|1)s z`t5N4toQ*PE{~{y+Y66bsat>9F;aScrXY|{!Hu8!K0ujVHH-$v2PM|AqS}-FExJ#5 z*QU$=0whgMihpHY2P=@~9BQrMT7jf$AK+JaTg@Y0pDrXQBCvdjGD19*Eb|j`NzJom z)^{z7ati%0I&xt3%$l@HaKr~+?9pwBLGyA;VX{EXcOOyl`*2T{!-olfI+4MXDMVqh zh@A@ZJyV6URrF6I2P~&Szz=>p!7V4J;|KQoZiLx-%NDm6c&A!7$enK!|LzrYX zv#`0e&^Q%btKJk0nw2lm9;_z0QE|Kqzn+k{n))bJSqvLK-aP$aEA+*4alGfY<2_S8 z(L5YLS$_~#WWP;~tIR&;lqI~SRZNIPp|->;LLGu%lQf1Ehz`4~Er_rvIMwnp`?BPA zHrlVVK++NOLc-yMQ%CKCmcoT*FfRJI^`htPgb@Ep2Og^r;1~W}S?&EPUEJb{i2GBp zab3KN?o?Z7GTw^X@A6N8z!d2S{-Mi9&&|$$Z-3jtf+iLFyq`oyQ0vj(=^R3?^7aGY zxp~(McU{s#j&m`t(-KJtyM;15`#3~*Bf{S(dhw9oYP3~TrESuqnbZa(C%C6QmT5$i z$&lQ^@ngLEeHJy$#UAYV-tY2RP((^(JjghIbA;BtgnokhN@m>W^%^osynyrDtr27u z?SJIQ6j0E$@Nbep449`G)X?+NEsXe31KtT7lP^HUK$H%Na@tu+3 zuS&m!(`AC82tw~J@LhjKD%q@IbIpg?{(rm?hY1K8eAlv^!7A!mfY@mDZcyF${LyVPQ8sCSt0jj#ILX`HKvK(lM6+AI(!evcr?hk%Y?n%5Eg@z zS~xWQYpHE5-0fM3SQCs#n(sRM|ihDTe!BI5X7 z@&p(V6*8_9*u>?;L5uO8)|r=uOMgSP=7tF**E;2s$apuE6ui`C;;n*%m$<_a7w{0G ztra7IbBCLe*7|aFt*hqTa!fz{8-84hNynt8_TK(YMGpR>dGK=B6Dw0IY|iV|wA>)c zwKO_f0oG*{LCo+JVo6wNSluaX-%uTN^mnH+oqpg#D|S|TkjOxCjP@ug+J9oPXTdav z3{dC&tur=S)OY81wYZod2FQpBe^a2nGiJqgyi`hrmP=l7r(@7W2K@N*=z>B{4U&B1 z8zkvq(3cdioyRJd2fchr(&{m?E>Ikfe52k=>nzetz^&JG19ic3SnrKWY$vJ$0+V*J zRTx@pPE8#KRvN^?T$Z=KFn_{|_qMkdakpvyNB#oyph&l$@KeURYslq(8VRoW2et@f z^buo}P{Vbk<6P8zTDci6Ak2CQ9Gx~G$F}t5*tQ!LwkJXo4`l77)SWC~KxVVDxeCyZ z*_tc1k2d!6h)h;COL!CX{Sg2AC^;0~&y_#rZ<_FA8m3_RBMWAmokvYr~Udu%CV=~3@6%9c_1wh*jcd>>mr+-Jg78OL&^Zrq=e)nM@ zct7gVDv3}kYWfH%z`e5$`)w@a@`%aOp4UC;-l%5YU>xlywyAoge(7x_>lyw>6-grB zZcHO~rj%hEX{L(cw%8eT=h)ja?nx~oT%-n928n``dRkdlA~Ydg^}9}Rz$?BwggIJOR z`a8Myw~9ma@?yN!72F_Q;pZJ3HYpQF6t`FZAk+n`RDXZ9vV(O00B5$jUKQUTPw5_y z{a-fyBH}Y@xwaTlsLgE#XEYWSRPa+Pe-b|`msn?fkC58W#2vfHc=pd-ra5!QP< zIX=v3i5EC1trxRBC#pIv!mXhOc>p?7go!gEJimJxSPy-CcN}pS&W4|pgIqW;my_*I z4R4J;{%YgkdQ}oxgEq%}J9X2eV9yrKrFNtnylOc7uBYfi9KgwqCA zIwNLl6sODtM15d2_5)sgBk%Kz_IH0ZxReCbEQNdpR8`&9w&Ec~B&Cr~K{y>W3et_V zv>+)VDS9kI0cj2?-6Gwc3QB_@T>_FKo&Hom&zy5@#@>5J z+{QBCa=Hxqq8S?_lF(iFg(~Tk=K7uWU~+4Bae_M~6${UPhu>`+I(Rl-O)FWfvl#k} zCMQ152i>~=%iV37Yc7S{b_;I6sSdh}C9$1gPI*o@`9#`=o~2mF?ntTP^Nb&nfM&=T z73uN(E`??7nda2w1MZ=taTA3I9vTs-ujul%1aEh)@15Qjvdp=Tnk-XMeercqoLUEq zU)-K=*7(G$-zZ*5Qfg(JpPDclL8Z=-swn&_QD)X!!;-t+JHHoF#a_~jFx_`VVz=L@ zzRT@MW;|gcvhI2W9=JanR(!c;JB^90oye-rosjYRYrD0GaxcD3;P3VCCs_AMYsk#> z*jmJ7$7sZBe-Z`0Xg@I6UbYKW4;zR;l+M>Lb)pREoa|D@)91Gq#4(aTLHQcIb5jUH zh3)K6{qZVtQE$$Mk*C^Y3uQ6THEgZ}$HaiS*B!p27_B$q#qOnvRZ|-R>}9VM^Pk`0 z3r1WNu%0WX-hR1Dys(jRc|rb1fx53^QF|TXa7PqKXu9ID?W=Cno1sZlgQXqGAF;x| zk&mOx(xdiCi@INU3aYIvi(I5}j3?4~9-5^#$G`6bfr zIeqgZ^WT}GzNG_)v0Be_7Om6nVJFnPFm|HeGHqkl{I+jjP=*l}zo{19xcM`HSwDFiov?D#ttyDS1 zel95rJ8RxE-b}Rt}S*A;R?`YNX5rhLV>gXL;MZaQrS!bHB{K z?d`j0C%~zCjprk5ZeH?YJ9u!y{oO`LojNsh&2q(yUcBM{^l67>0sV}B;w`X*^W)&A zt?@Q!sxWQO{ zBToLSiRioO~}?zk%!q3ketn4q2j&TOi^8nq>ymSW}!>o?EfC^{RuXHg_yVmG5t_nPNiY+ zvUxski5e*HAAgb!k00aq2#r(1I-7#C^W@toHJkYWo$T zRtwG6qdc!0yrTEnRM%rgdfirJ5ae!WX_;yrf}g`4SzZw`S?n9_HjE7O)q14(xgzeN z6f3VpwxjsaahLOrDGk<`1jY8aE)NB%HThDeKw-Gf(pF$Yp>ejehduVdIR>Jaqa$w= zcpcSNaL-YrR`(U%x4pwH>H?~`8Co0Yrp?laQP9<{duXwDb&%JJn(p3)1(@)dg-w$%p-{h}3^r(`g$G=m|CqEP#^(~%C zi=LMx{>t4SJzbKWSB=FQ`iP1$S2b2N3*-~-)=hel*k8us zX;}V1CuyT`C_LrUchOq5+xd$!OBHehaVwm<)w!zGmAoVeJ(S?O^sGV)yP(0w2Lamb zcAp;a*znYUeEwN4YVK%n_f4*N{_f{_HTax2VKMin#P6jNOtLExcD|SURUky=dD?op z?<3|B!eai-lwij_NOo%33(slGnADRY5*v5NQ){-RwvTiZImG>u@6$7UxC(!e-hGxO zwSg8NxA~?o3+lDK*nPW3Wq8C>1IH9--s*U$p}P+rwczGWOHZNn{HlkPDC1NmO092} zqS(x^^%v!OToKOUjBTUuO?%!-OPXxAzQSf!aj{S~h&{S*c(qjWm4HCV{wmiM8&R_( zVtE&W!LlrK6M@IPg8q?jgY%nvRG@V)=OExp;*YPB*W9tLK~&c5qvX~@m_!2>Vi7j#E8XrET6`>)MB(M)J#YiXrOg5;qq$%iP)0o?Br-mn44eAShs74 zr?8lRB-;lfelbdSwFmoto>D=$3pLt5?*wyC>*Ad0y*7W1$>$^AW(6;7N8R-7+qgj# zaky!ctycu=Jy1)_e^)5+E5}L@%^j&u_Np|71e4RDGRzPkwQk;>rqDL;QmM0dl|jSo zc+k&wII(EtQ8iZeQ%%t)gRqyQA>F(lYcJw2zT29OqFoLo4cHh1oOB$>d+%O`npoo9%(gF|-6XenmLHI`SVw-6R2!{d8pl;=!1l;1OJPQ`xhvK#%uC<8?pm+Y756W>@=sBF_OM z_32F~(mcQ3_yBf==v--}%1=TsamfL1IEKGkw{&I&rwDK_!FAb%f^Sr(?l6f%n&b@V z2&fK*b;YsHv%`DKtX~Vvc$||>bdsR3k7d9Ldw{T^EflKe>!S&PB_os(eVhxU0jxmdQkRos6l!%O;m+4 zv0%x$;+{|SgrdgJB#pMmWSS5#wq+GAUsvFFYB3rgcw>2+ssA!H{o8bE zJ{k5`DcGkUGMOo?DSHPTyFYJVurqEa-9z;@sm6{Ct^P_eJY6epj{5l^JF*?|*jkam zF?`oJjpA+iz(z2$3g1U&Wk1nO3nwZB0inTJ$TL*}UVrA;d%m5oZA^#jOX8l*GLS&; zyygelTGepEKJzv>eSrj0dS|aghH~S%+m$6g3j1JPchZS6;n2w^he6*X)#n#anIREr zUTc;?S0k4TFXLJ&eoe-OI5m2cPpvu@hy{+z8fdwv(U?&@E)f%2{h+TMtr5|GLwHP9 zjW~voQi8$$^+I~4D@y>6y4AbOgn2(W0~zVO0*A#6C!kJ>Pugarsj3VX8^C;$5=OZRMQz^+fC}>KxGP1w(5@ag3&v#2nBylB8EnW^q^~NvQ_2ahIR%uayqX|cP3Y}QC zr>UWwU-Oc*-W5pvUL1Dsfzikm*HYcoStpGyf^=}Z)Wh~P@&5f>d-hkiL@4LkOLR!F zE+pRKUtJKL!zDbiAx4*9@gmXeyKyIHg_t;k$!V!fS8Ha|U(lI3lSK?0kE7PI%#wJ-IQkz7ceQeuyRCe$RXR2AzeFCVhb5WL zoDtFPSZyP66(aSW4jag}>p#7CMK|k|r0Vndb?#&cY~I42zPc|qPbgqI{I~jr3JS%a z$G$t2=--L4Y9|wk`eNeoshBG*JJw>41npD;@}W5occ{N_^U}pmVxx7mnwB5xeReYo zltSy25mfr$cI@8_VDmh^kRHOcu8TE7C|g^#@JgLnjndu`i(0sOzpUQnu3)5yYcp+H zsh-m0-F+rXJDsKnjFhxDY#1U!xqn3xHwy>&&6dIA1372RB#n+^zgWgR-C0^?4Dt@tqV@eics+G}jzN7hD>!wxlhs~m-$uO|u6+SpmS-*& zOLp7&G_Nt}Xy9a~!4#qS!=I?zR5C)cE3uQrl? z3t#cA@Z)~I>tD{+3x@F%+v1S!lO(rL8zZnhzKDCx_iaEoHR4q_J?4-mVC@;#IQtP5 z7YBFw1Rjl9bAK3TxqRu9oe53s*jG-Bh?dh_8|`~?+h}-oz!^g~&uAxM#(ld$YhROo^yh^X9#8iHw3O6XneQ*m_bhM=}TN zLg}2FN)a_0L$?~F-^mrbR-c`;>eRN-(q6>k-{l{xWdKB`jPnNqkJt}c`@YF*r2q8fdt9(lN1*?QIjd~fR|mQKM%fRFS?py9LWh_N6n8Fgo=$yyrM zKn~Hjtu^sz(}WIt2~BC5d+P3$AMfY2vTU!CM2kR%=U=t+zA;i9PN)@b#}J8@|C2**zqyCjM{*fMuJCO^r8BL*GX3HQODaeS^0oSPCgmD4ZUWu zFLGQR5Lfv4`u9U_vRvQO68l)Yw!>Ye38jQib=9(m;$XY!Al!RJZ7d|XO1;FuxZ``$ z@~>AjrH=ReX_aZ?grx$L6_|f96ux`VSG3DvFE*im!>M$#HrY?^_GK)Ua%W=TY`S&Q za;HT3TSh)Xn!d_Y2s#P&dfh5?-J#6vq4QjJpY+paMAV9%mWJ$@kFKk~bWBL!BH}K| z`lb>iDPmkzj*_(@h43Z04@Un&GSLx(-+b7os}u=&$3@d%nziXYgC;sRA0)>$r@--Ke_@@|VHJ|$4;%PZJF7{Af}gXxJ&TTqu%kLU=)Sdp{s zkw}IXS?%wqW2+$nwn5MHnb&?GhgxoKzii<*5`8RloBuQ9^DS(r&OXd`&i$no@Cfz5 zhya*eb(GHPcb&_PJ+`g|&phsX5_*YLuP<)bk{UX#GKJ~$Tc$bFv5d^*^j*|r5cA@X zoVBD!n3A7$3Is$`NLXf3K3XARx0}c@pnj^k6v9`Mo*;m^oDrFj{bY$3Qd_bhZ@OxD z)n)F8SL?l$haCN8r}bylwm1$8QkhMCbM`ne#nsfrtLBIuy#b?R15r+Q`6eEDKPKH}JBweEHk>k`Na;QB;4NrtJG$^EnxW7U)i+_N&Qy z3YXngwa4pEv@{EzS<|nh4_Sh~&5~6nm@3AcUi`}6lY^}gEb+f)EKW%*fhG`|Fz~`v zVM|83OvicoUg{Z-A(Gyo0ne}2jg^23o2n5@23)+!*w~ZXbBDI(k1ZllE4GsFMq*uV z))q9lrc5fk4!F*?eF2*U7<0lWLtU?g-mv~5+jH=<=?8Nw3w^kG3eWSq+9i;wXC#br zvBT7gN~w&?5irW&hceBK4$FPt84@-7zH4z#>NXldtpiEl^Sjy0d6#)bq)yWhDrPj^ zoZU7y6~G?mvi5zDT^(w6{5BP`EX>66o80GZ%+9becPAa+&byamLiyt3(rzIS+ za@zY&GNZeaS(bfZ<-<4IuSHtBp}uh3K{%7zCR)+RmZ{1JxTQ()(!Jj1N^a~dzHMyj z5c1=_dRRvG9lyl-U{WLi*&9bCZaxN2!<#J5eF@GX#Eh-GY z^&zR(lqKY$*Isd04QR}*#g37!(C@i|LAwwiLyHJHP5N$q8BN|mhs1E&_je!lgt+tD zHijxTcGo5+UE3azc&2@4)K4At{#wL+o6^247K_-ZGZ5$>Ro9`YW8`Vsa69vND)#Y0 zHrc$s@g??3BUi+Avm1s{13zejyM;N&K78zTx2^~(HG#FzB9&#*7HPRyOBu1USZJV!{~XZG)uZ=8JBFd91dbMdFAOzDDIno9}1 zzg-%eNJ~<=?_+jBVECvSOUMYYS$}wZG#Bw3(sZvzVwc(LULUptH)A!qUMWvgpZOQX8IHSol z?8=4YCr7CaxYb*_v!3k6BLRZHX7(pE1dm?l#L;|ugmWMzt>GJlReY?pM;*3q(z35| z_dlgmTrhvRAWra(Kzcu}Y9NFG^TyxB67cdeP{?RGC_@cY8wG>_f6B^C-^c6~@K+L2+x(oQLV$PHU&UdFEQv4!dKSkaGUcD<=`$o&!M z;yb}p!*3NJ_w4s$f;noubpD+cY9Ljmn^!-#V!Z-o z*n?J8b9?=h+|QuzWQkhMox(G0%<7&jIF;_u_rExHpyB`iJK1eih?0)(<(~4!iCvBm zsLz?x`h(@%58VT0)|*J*qX@*(Q_@m_1a?v%&*dQV7}r)~)_GrbF1HAqWO1}G>4vX+ zW=)Hl#P^1xRYT_s;TDz3Zwo~~EMl{FoWfPFg#Em z7e<~6*RZKoHtoB_9$7sCAZN2vT3U|vIL`vpQ&t91;{Dc00ph1)SDsECaC}ztU(4Fe zJfcPUF_0yY7k$VJ&0q)&T>E;s@KO+9wQHO2$?>(}jd)&}ahaUOOr;g<=Ho0; zeTHAWJ)~EabJ`>ndpZIiO*%3-%=HTISMaT$7&#NEWj_!pWYe79&2VySO4HTqw5dvo znO{AVor<{R>>S*iLkohPtt2f~FQ?4MBxP#7?>M;^Y(u@R%L12MxHBU@*cmZ3HFyWD zAro}fCF|HeI|=OcPLu1&Kv?8HR^tt7X)`!fk_rCkbXLyD)L1*M*0J|A5Blg}h~Dn; z7A46_7U$=+L(?H%(h&OzKdE_1hZi2AW~tNmHs6lVc0V}j1`T|ufLD_MAMdq+zCB~1 z&NFSU^yQbU$j8o-lEPPZ^)YTfKd`T_6M-1}RCSzmoK&o@VQ6aQ8Za5TYm_Ygd6RI&RpGW)pX^8EOb44rH^akF6c5EL4;oyJ$S7HM+a|sTI z;Mg}os+V9W2+sM>0E$N8n*IzRNH~u2=br%>hBNy!0K>qzqMslY5E27M#%%o=pn>mI zYWy9c;0RpG7DyF0`3po7^9v+#2?2xSq<;bP;7Bko?9Tv(M&f4v4j@n*c>B-VfVO9U z2M`#pZyQ(=fk9(%>^px)fYHu>M`$n{*YtOUg2Ql>yTD8&5`!}X1~{rMARbM-fI>JD z6SMbM4PeyY0T_i#0S37KO~9soKm?3I;iUHg888Hm3)=^@z%ft+Zszakk2BzdKNt#z z!96`v^Ye}2Ota3maw>-%#J20eETfkaiZAA%s4E+GKhbxuG$f-p23 zUkr*w!SFFS7=yrNoPf-wz!(hvMHm!z9>XBc>4#zvNZ?%n5D%apigP#xaf2Xe2pkvn z2ZEsyxVqE7k%QyN&OoYwY3El!L(X{&g~Fkgug*Yymq0Ky7+tAPK*F2DDaVI@#|IIvLG7v;OC1@lBiXRVv zgXnX{KoOAh#vso-g2dB-1vrmDBH(xd20{*G5}q*#Bod5oAb=sy%b}p>asv49?YXGJ z&=@S9AQ%CMLeF^uhoa6!0gi;7m%||N908+|K+@nlf`+31k;DEg@qs9yG3O>DP!N1d zp%G{({(=C8pK4$<0t3f)4vhqZ|C^3jfE^Gd>YNfJ1ce_~z?1VHW5E9h20c$W3>1lP z3}7Dujl~ZQ=3MOp1K*DW<3|Ps16;wI3lvq@|AC!z1OrVT~|bMY5M1ML)lG8zhlp4$>A3cpE!(J-9m z6+#*;5{|$pJQ@kdZx#SQG3U4qVCRw#4HQOvz5^Hxe^=2c4E&Gje-iuO+(x6p=SnDm zoudRAjlyp_KrsZG0$y-2K#jyNK^QO!@gJM;M?V;hfgtb*{})4?TM#(3;8P3(MdNo$ zFb1IAe<=AM9dPKmNMqnwfKdN!H3R~KA@QsRFcf~tfIu-Y1YW2C42ACq1ct=G@wOOP z4ublt1MvS};s0}_z~5s4L*sK6f`FgnCxAim69@uq5`Jj`Frcmdhr@s34?)9_=O&|J zn7{20xc+YgF(@b!F9#qPU_0@m00qNQP`qsduyaKh$Q}e1zhppxe8lf%P#}c()gB6j z5I;+xa0C*Jw{-x9#?Npl0tNf0vj1;m1wbVJF%JspMdB4*07Ky)NT4Vr8iB7Eg#zPe z1QZR&{Pz(Kh(8d8b9D&@M!=AGVFj>rieZ2Xd;Yy^|B=QpC=_7QKXU(StN;kyBH&L3 zk`X`Dz+J$(Er9|-z$X`g;nxK~2l794JaRR+#nn&~!pv@?We@- ErgebnisMatrix data = size(G_E,1); %save_* -> ZwischenSpeicherung - save_A = {}; - save_x = {}; - save_A_fine = {}; - save_x_fine = {}; +% save_A = {}; +% save_x = {}; +% save_A_fine = {}; +% save_x_fine = {}; %Alle MatrixBrechenungsArten mit dem selben Netz berechnen for i = 1:length(type) @@ -82,9 +82,9 @@ time = zeros(1,3); x_fine = D*y; con = cond(A); end + % \tilde \mu ( \Pi h -h + L_2 ) tmu = hmin.* b .* sum((x_fine(f2s)'-repmat(sum(x_fine(f2s)',1)/4,4,1)).^2)' /4; - %Fehlerschätzer 2 aufbauen V = mex_build_V(G_C,G_E,zeta,type(i)); @@ -116,19 +116,34 @@ time = zeros(1,3); % xe = x'*A*x; xe = b'*x; + %\tilde \mu 2 = ( ||\Pi h|| - ||h||) + tmu2 = hmin.* b.* (... + sum((x_fine(f2s)').^2)'-sum(repmat(sum(x_fine(f2s)',1)/4,4,1).^2)'... + ) /4; + eta = abs(xe_fine-xe); - save_A_fine{i} = V_fine; - save_x_fine{i} = x_fine; - - save_A{i} = V; - save_x{i} = x; +% save_A_fine{i} = V_fine; +% save_x_fine{i} = x_fine; +% +% save_A{i} = V; +% save_x{i} = x; - data = [data type(i) sqrt(sum(tmu)) sqrt(eta) xe sqrt(sum(mu))... - min(hmin)/max(hmax) min(hmax)/max(hmax) min(hmin./hmax) con]; + data = [data ... + type(i) ... + sqrt(sum(tmu))... + sqrt(eta) ... + xe ... + sqrt(sum(mu))... + min(hmin)/max(hmax)... + min(hmax)/max(hmax)... + min(hmin./hmax) con... + sqrt(sum(tmu2))... + xe_fine... + ]; end time(2) = toc; @@ -168,14 +183,14 @@ time = zeros(1,3); %Alle Relevanten zwischenInformationen Speichern - out = '_'; - if(length(varargin)~=0) - out = varargin{1}; - end - typeN = int2str(type); - save (['meshSave/fine_' out typeN(typeN~=' ') '_' int2str(size(G_T,1))],... - 'coo_fine', 'ele_fine','neigh_fine','f2s','data',... - 'save_A','save_x','save_A_fine','save_x_fine','b','b_fine') +% out = '_'; +% if(length(varargin)~=0) +% out = varargin{1}; +% end +% typeN = int2str(type); +% save (['meshSave/fine_' out typeN(typeN~=' ') '_' int2str(size(G_T,1))],... +% 'coo_fine', 'ele_fine','neigh_fine','f2s','data',... +% 'save_A','save_x','save_A_fine','save_x_fine','b','b_fine') % clear 'coo_fine' 'ele_fine' 'neigh_fine' 'f2s' % plotShape(G_C,G_E,''); diff --git a/src/test_sol.m b/src/test_sol.m index f2d7470..8d91108 100644 --- a/src/test_sol.m +++ b/src/test_sol.m @@ -7,7 +7,7 @@ mex mex_build_V.cpp slpRectangle.cpp % Test ausführen %Anzahl der Schritte -steps = 33; +steps = 30; %Art der Berechnungen type = [1]; -- 2.47.3