From 7f19cbe523b6408c38fe372b2cdfb11d62d92dd5 Mon Sep 17 00:00:00 2001 From: Peter Schaefer Date: Wed, 28 Mar 2012 22:24:37 +0200 Subject: [PATCH] [src] nutzlose Funktionen entfernt [src] area und mex_build_AU umbenannt [src] Variablen hoffentlich komplett einheitlich --- doc/doc.pdf | Bin 330346 -> 330452 bytes doc/doc.tex | 5 +- src/A_plot.m | 35 - src/A_plots.m | 4 +- src/A_run.m | 8 +- src/A_step.m | 16 +- src/Aelements.m | 19 - src/{area.m => areaQuad.m} | 2 +- src/computeEstSlpMuTilde.m | 34 - src/expmlAA_2DQuad_error.eps | 824 --------------------- src/expmlAA_2DQuad_hminmax.eps | 865 ---------------------- src/expmlAA_2DQuad_norm.eps | 552 -------------- src/{mex_build_AU.cpp => mex_build_V.cpp} | 23 +- src/test_sol.m | 2 +- 14 files changed, 30 insertions(+), 2359 deletions(-) delete mode 100644 src/A_plot.m delete mode 100644 src/Aelements.m rename src/{area.m => areaQuad.m} (92%) delete mode 100644 src/computeEstSlpMuTilde.m delete mode 100644 src/expmlAA_2DQuad_error.eps delete mode 100644 src/expmlAA_2DQuad_hminmax.eps delete mode 100644 src/expmlAA_2DQuad_norm.eps rename src/{mex_build_AU.cpp => mex_build_V.cpp} (88%) diff --git a/doc/doc.pdf b/doc/doc.pdf index 8597ce2bb69b53a29b1b282f29e53f680935b08f..add9a02eabd437711c882856f64871fb358e9ade 100644 GIT binary patch delta 17963 zcmajG1yqz>*9J@}jWp7TfHafL4Bg!lA`J>6(w!qEAV^3IB@!wvB^?S#mw-x#(o)hO z0{_kP{@?$-&+x7Fv(_x`b?tNRbN03OzRo_w;#U&w*JRpqb7KN92nkKRrNoE6r|v$g zLfLX^L`GP%5&Fzmy5J`cX>jg%>DvvHw$|jxd*L>Oeb0(m=F6PpY#$7%*mNycbXbW* zE!rO}te-c0<6Nbt+V_$XyRgu_KBQHkL~qhTn;xK)eoN>_iao4eCck~hBM{18_e1_U zpSO{jmPoe@-^S(z;cJ%O9|4HAO!(5dPEzYhP9zpN!@jZ9DL&$b~?Sw**#eEqR%=<{7+u%C2PDgU6LalFlxQ{2-!*cqf z0TR;33J)@Q$8&3hy`K8cy}R{j{LPDzb2GK6jj|`4P&RwPY|<6}MvKXWFYR3n#x$C_ zx7Ib*=D6pG6APHNONO{6Acg7-vSbpt_q;K&tc14xsW|pvNT4IMO|ex%N8c8!KBRlJ z-2;l{6e4^^C>`Xpp;-D<)tJE{HoHiqRWL2pu=3Q44`P_=W^)jXaL= zd|o3bEJP9U%E6VO|8XiK>zAjr{A?mj(>DZlBDzSs6Ks6lB6fD~Nl;x7Ykw}4Na&(p zI=e$9?L@spn4z$tdqNe9NEwZq+S3Z>+ zv?1<-$wZCnH$Rr;N0;xJ8t!?H2Z_F(31>4raBF5#6*Dz&{J810?@iqI?WeC@m74F$ zEM+ULr$F|kmTBN4W_5o105(G@Bf5r(<2n_Dlv?_Wva|Y$N@mqs-(f~vW#C- z9_GT|11U*e^kL1EpoHCuKZYTiWv;cR|L0aKc*2mi;rnOpVhc-LHIt8 zO5qmam()l(H4sj4?p|^0+MN*L`7cYs@+X1L3T1fAMCo}u^LZD7l*-`W73t*o3;Dyd z?(SL%I@7PGz!uj0M&kn;-ls$}{b}>?GN@t2T^%=_dD4-iB*Bsk%k$LNRf=0+ z^UHWHt~+kP#}Wd1G}ZPg`=g7-b+h=YXdzn8zapgecJU~(*q+q(`DaBsL}-~FI$r zSxu4IgWs(MH*ma7QpsgV7twfbUiij%LDghBP;Z>G2_`}we3_}v?@Kp2)C=hPXQmkm zp5j}Ya1i~(^L7ot7m|tld*dnQIX{l4GUeh8n`CLSaa=S7;irMtGw(X-X7shn@V{Q7 zgm^g=vCEX&CXw(@z4T;w1$~7spCl|daNHx8MC9GeNolv7w``qAMQhPdng0sUuyEWd zBnCxzAKRgYW$a9pQ%NJDg^rElGC#I`>JlWaNO~rQe98f0=(jlN85awE`n%rW?cBHS zkg`Dh^QlNrFTJ(>K{2_yk+K$RWND&+ite&DGcxSthAnvr?r~k%@ySUYF2^05Zv02> zCj!AW`L|bVcM)XUtz%PYs$WFqA+@->)zx^d(x7Iu{uvpZ)GC4QvD3Nd7iy`~(>UrT z9AuWbl{Mj(=QOMop(W?>KUTJ%tmA?fLq_A6KOL+e93gCw?SeU!JNO5yAN{^f{u)%IVb3^Wu2{ zQAG$|*!d$ospgMNuBVBN{pk1Pr1Kvczh=8gh2r#Ie3^@kjw2;PFPx*zL*|4p(9;)m zOBZSv7o{qLbY$nHcIWN$7kLQp@@xL*_P8rm0%R@BMLNgNWX{L*Z|rx;6^95qnz5c6 z6%#=twj;SdGY#(9p{IvXby*j|Qw`k{eoqDAPpo7zr0e#_tUkFvb-cwvUJ$)m=1;iQ z(QzOxcKi!>^M`jd>EG_4Vwv9qf1J&Yi6|eRI9&5RKE_=h99+&yb-DhDTk0tS5nv=F z%l46vVt(qL#o3?bAW;q8EG~Ja5-d?0xQOa2DE16D;xMM~9qc&BfywLMg>{iyM^j zduVa`5F=a*h2xM`=Xwi|<0*VLR`@IqWkso#Sie8Zxf7J^GD(McTOyJw^vAr7IE>c zg$HGWId6Gw%gxSxRnu#bo;r(cXqjN7vI9{L653=A?bYbY_Qk*UdI_m2;bhkQc8`}d zBYiLD!I>X)%L@H9k*+Q!H=a~Jo^?D$EHFj-?^$jalY$NL+wv?*l^5fk(oth?I^j?S z(FbZ@Ga5x-|4jFa=VMgQO@1`GUD8(L>BYgnkbWe z)U=w;lV^)+ZoJb}y*ZBaVyoWqW};MO!=GCYE*!~6O)~Z9vS%OpW8yQr)uTl1U{;Zt zNGbKgiEF{S&$rbY1W~cEP2IO*T*-ShEH}F}S}pzmQkFU7Ri}H*+jxGg;2HDdoy&Zr zwLOVSUekJa=&kjnG9`?N`EXL-kM9|idRTG7heRfG7bkWIF*#0ki~EVJ=k@z~eva88 zbzf=2ml^sU(3J>PFB|$-R6_$cnSvzrER#M09c8%hTv@6|Od9+X$T;3ua(vZM(95fy z@?Z_@9o*U6eCil}o~Y2hB-rvfIo3C~`%ecbmdDg5weO>pechD=<28UE5UqNJH^4$2gzg>tn zKc?s!i5Tq0Ek#|gj^rV|?s6ieoh9{h_%C&|QSI5isw45|s+r+?e0Od%M7cSw2rySs z!+E&J_mETX^6DD#`9jw|-3_j(`OeAb>>QGuBz>_~F!zZRenVtBOxric(=4XKl4FB> z^nUle^4d(jcYo2^t~h<|kLvEdS2D_pp-R;>UF=I*dJr5u(|6%Iw`Mo_zu_p|a(+E8$K_LBL$h)6ZI7w@aB@s`5OrG*<^C9K+p3VSuel5Zzu~=KJE0=N`lB1l@|QCJ~XB-+;;yrDVeoh z7s)>SfwsYxZ;ht!eRiYj6j{2L&ER@f$#*?b&u*kO%r?F}{UBfA{NCx{YnG0~!)~5e z4D9)}or4@`^K(m$NiBTo1b+_ju2HM<*UBI7-zD@_AO+s1|9IY1kmZS#6Nx)rgpq>*z}@b6&Z@k)(=NuCBBM(8D1{PR4<@ z&GI`OQ8nmtylr*<*+;1yDP#GKBDR$DCj1lvJ`+&&l_&M<-pG7Sn##cgXfk(1Kq+6o3j=)LGsN5JYfOH{Bbxe6H&!-mtMv zFWo1lGACn~cQ=^Ui#g* za`-Fv4Z2HOfNR+Uhxo04<)Hz2!A1srMAHAW_F4%NZ&gd>nfiwUi}%hMT?`bqG1;54 zJbU{CF}jTmE1kz7^NZuJY%zwN@$qbfb)`ZyDP;FIXU=ilH%LDZ^p-^@9%OI@vx=wl;j=EXLR=&3g(HlJS>Ej z@RdcF?qn-OMv`=?AV zx>jnbD)7vWeX`*@cyDAQjCcqSpN%1RapsI6BHk;kz;$}al0Cw5wLe1S;ZP6Udn&S6 zZdbC3IFQ-~Go3}5S$@(%w>+|cr$-deMylCFK~0BMB7LQ#A(sagA{53s!}k=@Mh}0( zEY(;TIFp~aa@4_p9&$z}gy^bE5Vv@VhP@Z(j~cN1nBtpn;C=W+Qm^F$wIJo2CyIs2 zi`C1S4t|1v~R0>wn3z}LsckTgMEDf{9u_TKXQ5e_8tT+IWR<Mz{Q~9Sp3< zIbv)pVIR{kD5aEEqB$bw2BqH_k1E&or>S*s?zBuZeW_8tRcv4@N%qwUy`^aS-O+mY zet;*KN45EL2=0j=e5WqR@HCJ*$ii>BviDdlc4WBc#Eyn>YQ(?&=Upod!-Rm4MDqtb z1qZP*=RgRRURG|<$|yOiGJPaU!efMJzquvjkSZ{;_vuB}+S1*l;Z=D_mmW$nft;wF znF-AY9dB%_ypqT)%N}l1p${sG=YOK|~Um4m-yZb)rrCXh6 zQZdK1cgfWg^ikci%1ab&l{p;SZNo&}_?!JGWER@MkK8f7pma)&H~=ImnUs0BA~Z3= zbw16lnCHKyjOOPu`JzN^cjx&@=RFo3*Z9rNSUsEPeO>-$Tu0yh-l2=neQWgpv}!Fg z^o)>%3{d2I)8yyXX@74VdzXdCRUBy&j>-7jYgTHdcu1VL8xpv+BqMyWbm&i8B(c&) z-G2Oy4N*m2a`eJeI#*Ss@vr3jr8RapKKbVqcaQ;lCSB}t+;Ggx_Zl9(fg0L)-@iZn zOc?dj;NafyODVrLu}NL@U~q5!cI+g34)=JcU`5?APdh2|_{go{{LX^a`1Gnjh8nLA zL**f>PZX?KG(z%@WDunSr7kXsLdQK>9-?eZ;VYUc=>j#r)RR`a_FXv2JM+Yx8l^b} z(k%UZ6lERedc((X(!yWl@6n@uAT`5~Hj-c(~rt z8?hC><0~D1j%BQVdEObI?Qj)pZU147GNGo-=g(eKY6^#(fM&~iMYx=5Yzt&Z4a*SW{dU}3UPzp`U)OI7yo zrD$@n@%{2&V5?JUNsAql)h=nQcb3QPYuK&KudU5T1O@5K#{`gkc--@`wABHl^j$px z$=-b2nh}`%>zJU#I1y7&k6{abQemvtzvi=3u)`UrCO%wv7wx1*IkCS|!WAUPuo$Vc zM>8oO9O@Q>$is9;PrPUODig8iC~p#6LHzPrUP-N;MhVxE)*sx-IF6NahNLRp!Sm`Q`fG|@uF#R7K~9a`Rge2c!Lrcglzd(9p6h_99OK`V4wC|{)&mH9SZ(bCf3 zz##iqp2~%3(rr%16E7(F!|{oV(594(?ky=Jir08$0xUOw$&sr|rPwvhz7$oncYh%J z-nw4yEbZYRF8q&A-#Vq#ABa{J>{JvzX128nB&r$BX12MD9_X&#|07(+B$wDLfbsKo zBsbdjp~f6J(o`duP|o974~>$S)|*sI6Hni{=gE$}mtdpsG~5ZA%%4 za!^%G%NgOYc8ufsD*x$8_OYNa$rHn~&Em)j!JYH0MXrVr7kSUwcfoko*SxJ43>(|x zy`%kq$EQTKqV=z%D8y~qnoZa8X|JVqxSNrj=hBr@dirDJC%-zCcPm_$0PJURwl9W0- z=lp&o3+MzB#JVz}{oZ$EJDV1)n_RYSWib7t_EcIfun-udt z(;ITT*Vfd!eDcMRvL0&EyDJN)1TD|Dx^Y%;3{q3`A8FoP_sR~3>E*vB_9A@pRBxnS zn7#k22UI5$9pJswy}Ijn%R5z#n{tO2X>mJ27w3ArVv}3)5w&hjP>TncPV@H2Mz7pa zz4@Zo(nB}3GLx>=wHbu{N5386xAR_u1Sqi}*4X>P6@mJn)obV;4x~{EyEP=FJ?9(F z`vk=dIdrwM>2NO$oVZPbkHmux3@GLhDBSWqU9{9+^eEm;>2RfAjoRPb<=&z#@yEg4 zeyo`6m|yjkvO{Mx3Fiv1&1qiy7I9y@bs^TTK@u(CKjUMCx-nAK6)(&ZR)czqcbCLT zLr{eg<0tTfBB#dney=y>FXuNsaSkUGa2Gd_qdwCx zkzX-lzZO2pvWTqNQ&af2yYMeW)c>U}nlRGizRf#m?r zbD!U6GLig1sb_WFFWX-HeSvFzvq-l0+t8gY9p$#Zp|uZdm1j&7cy2#+Dx-qL$xU-g z_xKz7+13>byEh7Uf)!tsfiawVJn0p`9!53IHQsy6w)dT|QEoC}HDip9(3OJSP4kkH; zwl8?^ql9Ym{|p*#Y|wf8?uUIY(F^EH$Ms1>oM~d3%EvV6ZU73lWYrNU z<$>lV_s3W)EQ8!5o{>k{>2o@|QQ43YMXnpuX9wUNeqNWRc?~(vNWTyF26W7G2i`&Q&XQ2nY|W`7$s=`{q=ie{#GKNC3mf%;O?c zwo2-XH*Log%;z?b@|uQMev66t8ST%as+X(} z=w08?oT03*Wv+C8z<8swlG!9lreiel9el^Y=P;`5q$ z#->b|VFhG&mnN`!`9?aQeR5w&^f5hk9LziB#h~>so%)n&4O*oq(hc#f6ABe^R#KXr zXxr{O{CiH)H=VKszdd4R!(Y*Pp&cRM-Z|;8i=21ZeqsUVS3DMN!nW5T~gE`d&&mP{M0y zwW}X5^j7auGUEB=XSV!Ks-)z_ABFgRu6dQKlGJ?e@iOMTAHwci-|o5cg(&Cf_D?pn zS^5P$aea)6-THpT!&!Dy^5Hju+jHsthIbMw`rk?KIT^+|mKwZnTxuQV-*o(RaO~m| zcFtw$j&n9**6RleT0cAMJ3Y+QXgG0mdDpR)j8SngOV3>?N!XB38(Osr%X2Sm<~OV& zC{hrKtPqbwpvb2Y)VTh>&0|9{tUmo{w@+&YaF42*aWgID&e7H9X*{HA^Nm*P7A2m1 zzMg{$!CLxb_jiRHZRR*l_gFiY>^IdxvhQa;(Q(9Pwvi6JV8L}3vh0fs4%M0HQ+dHF z4S8_kbegVkNckd=`Ar}6g-DSa3HO6re)I(dEF&|#v{~=I7%!$hrea>EczqGCf{vJa zTR4`&P!e?E$bXJ1ZWxObrLgP^X9X3uJ>78%d}w7A0y^x=Sxvx^eTi%A*ZyrFVpQHs zRjNDkXM^bjh+WF@%)JMHcxo!@NtHvHRiMBTj!v>wsC z`#kbqFhk1HjPH%^_gRG+r81>w9TtDz9BMIg81J4Hep&iWa$R#W(nOs0X5p%7+|ajv z#>ZQxAIj-Wdi>s+-;+yZpiX-H1?7??8xf&@L;^Bhf-qwZ1NGv>mfr;6~$q* z=tY0xw|jcBspU^tqB8O16zg;bBgTH+N2wV#J^YF?C^@(fKF`rNIUUFSF#7o(OV^Mo zBJ)MqF}F|c=FPa;cyq0kBJDD8G`rp{de!V2dkKF`SXn z^*{AON0##DF1(!>PkPy~FjN_?kw$ik!!9;A^Hk z+xTouU{b(7(_TnPvJOx2I`nt!% z*4MN7sH=K7dfGb~KUY@yl_xqv_quJVall!x$bO>I&k1TZmEAy+z={k^p)ikN*x$Ad_K_4aN_4Kyb8K?nj762)l57Cc7&znw*=BYi8%ysv&?TlwjSNx2=*O zj}ZCgU~p~G-*G}{DY97k4w)}=6x2g?GRIu}-fI!pzOxS}gYkDtW24BW4;LmaSqggt z-0tGO+ox^Co)8wKfUNO-CAxM!v9Pf2tZ zNKm|eA)fFOnWc4K^GuKtedMn|`uH>V1R`^Z$}h`{8$5I0$#(VV*qxf488b@l_+fvS zf>k@5%5QWZ?@aN@8dr1^{`8~IJQR+u#mC=`$IS0@WpwDcB(eN4w2ThF_kG(qfv$dZ zu&PVHge&dU&z44>rAMFT=K{%6W|=6Ztb(BB-}{)1MO9!*&y=N^Xwh0nqe@$M{^+|z zGgc^$GAOn#;M+6T=O{hw8(z-$yj{apoOI*CW)YJEu~X?p4ep#nkvL0i^pO4=`QOR2 zj(8+Q#VcHdg`F&-aVe_HF4ofOi9F2YY61Jz#o9!R?T?0{W44nwlKeZ%GZsOdnUZ&B zOp7?A6Rfl$-6JL)gmCnUN8jM?0^hIn+>LUU_3@I;TO-*{@wJiG!*5ez#KZ-C_jUIa z6Sa8^lq7oCo;>1pYntB>F#jl}xy=(Bl`8TbH8jvlX6plE+%`rZWv2dMmMUK8>zi-v zEGbhx&wQ-f`tpOc_)sO)Uw3XP9gz*~$v1+hZ?|V)WnstmOQ%^LydRj2+eM{$Srs`hkP5r*wU}Owj{YPG7Cxjn*k8bPi0RfD=La|5 zQY&6{`B#6pNl)4J)|o61d9~ervcS*T-q02vvs%ZAk9!Kxm~jQ}L@0vXx>u2;UCFl8vq zI<(wTzWE(l$?Q?32XgxItxgM?*yb4^kFFH=5UG~A5z?*ap7ULDKozvX)7R9=^u9oo~D1(|)} z_lQ<}w>(suF$pc_8Fx6WDsH!Kj7rtAHzyYaz1B5@Ozwpa5qv9Iv+Demo$iOnJQ40m zY|5+VlYJNv)pV2bI6uc6)%^7`;S5bF2B&*;BDyj zL`hjU5~b6iq{B8s*V9)5s25K@cb<8?^Em1@<>cA)q~cu}(1c`y>UL5&YHu-p&CS|NEuI3}}lLW&szKO7;sXA8lbL-F5qskP|WCg>>7GhOG zL1C!C$=hAbzH;`C;F)$)&|NQ0f zC8*=-e(}wJBjDAG(wrXT9$ywTwzf|GZT^yO)`-Kyhd(pxFmfMi;VoGdw#3?~p+(0I zr?Fgo@{TqP+Vkjz;uEKu?1P|KozdF1q;mF5V=0YYuAp9Hk4atiz^6RSB24Gv6c&qR zS`jk09HhLAnvD8bjXF0f>Jo*CBxx7tcsrckkOR9U`6mHJo}CSWB;=D(XXEO;Djzj+ z!CO+vF*y^F2XCz-dXP}x=}quzdf44KM4++Fgy)eaTGOAFM__mK+Do+boQZ5hLpS-P z*KZ%bVM5{+4H@yEAhuH39#}a`jCq;bs0zzp{HhfuVtEfySA1Do#Uu3=e~6xcpLpxw z_t$-2ozwgmI7r_P9Xba@J*rx;P-L6no*qd++R-(LKWC7cTdDpY4w)`QP>= zd(#gROM%Z6TPm2+6ZUwg_tm;+WPLfyi1ZiiV`tFcV7G`CgiDk}ZLTLoX_#tDMRBuB zyUA#lX}_lT&1>5oACLTgv;R+l(u5ber&~OsNcgh`v9O>NHbPeW)CE3A598`IY+CLq)m{tV}`ro%M#E9~LmLM1EWAlM^iL%Gss33Blj-G6P>G@fQSW%;} zh{?PVeRp+VwvD;>$&a8o0k25s?+cr^=*xFuv5-GoeTgNKDgHK2#;Q)cvtaFpHVt`$ zr$q^6bi5naPTiN%+hZHgFx@$CUo31o8^4TlX5T(uLhs&f@H~098RCF%wi4LBBD?w| zTgsQ6eoj9`-js-w+4VUQ?vy62swRKksMv6?82Xh|Al+j72TrO{iwzbp6a*K$DTw(}9H z;NH6E1Hl_Y>YwTJ7;hd2P7mu%Q+U6ORx{>r?_9w+4H4qcFf+>aKgIt#A~5oL@j7Xr zi@XqeK)I^iyHov@JLmCoQ9!w}A8&=)@1ZQzT{EGP)K2#w(LTieCA$;;uz;bPH)+Q*Cy79-ET+eNDM! z&PLCjxhgzBc|yDU>8E3BvR<`}{LlLe0p^g|BnoopKBpq~J2TLOH1-HE`?%}#CiRp_ zDcM&P`?I7RpIV=wj;MrR|M}A68b95pB-h_0JQi5_X8n%zdXm87h>#%3e(m3RiwlwK z{AK7dq8_y;_+w6Sd$A#y?0bh3MZfpd$%`^QIbxr%QY{LmhYb9c_v+4#09_<|llLSM z>+*|0G~Dj(&$x#_f@o>^7d6cPnmXfgE0>?UUA8fpdn4}A2W5$b_pYEs%6*p4A(|2% z#Dl5Y+j=Cy=DEiYoH7af9Hvdi-myKaZ0oQ?r#yOnFMjxr$tOcx-l)Ei7S5vpU2<7kJXS4zV>2_-!}#V{8i zqUpayDm!iULWt8LhT4;H4jBwZqm0jvK*>VFdpm?Y5-ggV(qN(4kS-PTaqVDI^YMF} zYkqgfZqrX9`2O3EWMso{eM}vZT{`%5@qbYnN zydN`xJt#k!F1b(8`M6MuvDeVElpO8Lbsn41nQmb4D=)xmIb5`A_|$(L<+A>Hsb|H~ zVjGkgpS{)2kk&y$Y;)b&{Iz>mv5HtOn*-&P?91`aO_ikrSuJG8_3dIl>eJWTS6X*m^9~|@W~e*QSj7w(ATdijEb>9E~CMxO3+u1 z?^St1mGIugR{cFB8My%&xh|nsE%c`T`5p4e-@-(tD*@eojZI1)$R&dweXEL2NM|fx zXc(j3OVSn6vxGVxD_K0%ScJJ%qobnihRDbjL|@rc;YxiGdfr>XZMRHrW7|e1O4#kX ztTj-%Scb=x%H~h(IDPCpu={E#2k!2=(#L=a@OE1f;qTj7Wctv04fcDQf;096Y_V(0 z|K2pwv^VcDkT{W6+U8HZ3_7j0n_AU9(r>sRdd4Y-TVjf>b)s%gnXfa}pnswXgc&NH z1ibL(uqQpCXf#x+m-)abG*#*^SCuW{sxe)Q5jpxRLQ-mMiY`E;-A_^?CF|guvW3sZ z9e}8Vh%)TA3d#MQvcOmX1EKLzIH%AxjCqcfSG0yY^vLyHK>ZmCS$ZOs;#-Q`@Zq|;_6WHL=SUs?}v3ern z_R#jpLmzifYxjp9wJJ9W5JlJi8IK8TBS7IoYFpX}yl!7dfe?!4=Q8)-Qyh?y#07D^HjBtb~=GNM!5sAEFgdtHF*GnS=2D@qmK`^bCMkuh_ zB>CIo@RQ2n|z!bPtdK`@Nfr4a?aY6L?;n2Jjy9CUR}5L^_4yMAeeqp)8KLqWk9 z-Akh=7=Fcwgd;FH{}>VAD@FtqirKj|A|O}jgNY&_7}Q7o zRust1if!PvP{@_2K;b~rdgSB+2<9Fy~p5q`x_C=@JOyYriX z6(1k?=LZbcJ^{A?7M{Ie-I=L$jyBMD`_O`%6q}6t2F{aT-_WP0fJ+24+yw{ z<-r)f%Lh0>uI}Xn7z)Gq9so6a8UG&O0od6U?|_=ULL?vr;9A*30!4g$IFQRxCj_2A ziC{$v1c8AdU<}hQjDsPeSTzlT0?b^&z+e>C`+#vUDeU|KL6IN`HvK@}BeBW?1P8<6NX)@0 z5JM<79UwRYdUdBrU|;{H4Fp9&ukId!#1=mg6b`z|5E9T*>=;5(;Hv{sU{S0p0TzZ} z7y2c3C3yfvL18Bk91a8jcdUso;{XFAFVCnmfDRY}n|v4qh$fcnz&oJ-{g!`%iGp9j zfVtQQA_(vtTN*(yI1KvF8vmWrFd$0UQ3vz{d_mVeMUmK;CJKOpJ<(^LXwf`DVc5I8VF*xmsc z0^2*F#IAN(AfT(k+6)p72K_r{0R90GyYmAi5p@+q0?B~&LZI2WBGhms(Ar!X2S;Hi zH5>_OEE;Pd0s<8G6%2{p@q!Q#6qY9-5CX^r?D9lFVgCpAzgmI*e|vghGSKK_FASW! z*c2iFR!HHrWFYmz0yGdOkWjzATzN`5zweBiiw0H zuXHI$I0^#g_(%W#R|Alu@GIg2VA$<6=u%m+JOMD&zohl@i~z7J`2id~|68)?|6Hkn z691!FV}}mF5ZD_8PTMP`hz6AP%GrtnB7z-H;3jm1YT&Mg zefR@QK(PBL029To2r!^gmnZK({Qr*z0H7$gCt#qVM_|Pm42Ugu^9P2(f$PS=y)Xi}V`0@Xz#j~|WC83-0t0ME@uL5Cl+R*b5_o*kMr(K|)2bTQdMdVc!^_K!b-} zIRN&rg)gfUC^<9|y8{J86fpdoh0DtU_Q?huyOXB?o_ zS3J3TQMgo!E35(y9O4RFNF*A&D*<|v|KrLJ9Pe042#fZ~u3x0N7#yVE=EQ{$F8&01X@VW`J`IyMqSiBCypFz>xo|{g=|u0?-cu0e2oO zPk^o;_|~9@_?oOG7$GVPe0u}5yTHW`rl=@)dFum<$|Drz6~RyiuoQ>H|L+hp4k;-T acMmH!4gNf8D3Dh4A9ucC^FxGfBxx1K$-3|jF>iUf$ z5XpV7u&tfSSQMyHf|N6BuN5UjHEfw3_>Rbm@W+DURNM*&U~-iGi)qb{l+Mb>bB2D5WRDt@Yzfa<7~?e@D%H zAR@n1Upu_BMOGObC#BIswEmPw+Z8Uq)e+tC$mQrunoJ(6(0#vaTD+`8BrB|KLZ*!P zV2;;99Qk)%4k`x9-@iL{O;z~AL!v+dXU?~tw@|XMhWGd?e2L?GnCp^mZt+qX62ga< zwQN{Fr!tY)ONu*2rU!VXrPK7@AWwB7E4e;i{bkqbHL;h9jCpQsRSY?$tr@49rX5%K zrxb82A1$?qJHvR&@d~8rwBK3rBb>%K{HWustP3i&G->%x zUnx7uU3qVlbE>w6;lh-!?RQ!5WX5J(&vO(gMQRGOS9u+?Ta{da$)(!t6nznvyneX< z1sR+5jPuGg`>U2{Bi+YnLRHG-6v}I!_nux;_iB*K{urD5i7dLVjXQ*@eD$*Dr<-@9 zg!jRUcv6a;(~p0WKfC?PQMscqcvD5M%8@RXdZUZ4D?ipt%X+{ZlYT56M0VY2tWV_2 zS+}W|glLloV=s5&9D?QXTwdXPIKk9Yt*Xgwf&Kn0ZSHP#>|q0^arv;PoDDk*?MA50 zkO{>n=qQv}Rr>zYr&@X~@pOp|dqIV-`z`HwYx8xE?C0dFbI(!g$G7Lx+}Vo_?_75u z4nM4Qd+Xn+Vm&#XjtP=U=1(8w=RQ)sKYlgATk!rZmEe8?272VkT>D2&1-AG7kS>9Y zQ!&9s;{|lllg>MVR(D(C>i4FdHce9v+|^5Mnv_arcTMSMy|Yw=TgYCiIyH>nEj<6O z>XiCSMf7}5$3p8lO;laM%Eo;DYr#^@*Ke1Yyx(RhO?t$48=hm#b!|f2vo2HLXQW$& zI38eEMIYdg>Rxita|pdY5Kuit5zwAo<7i1a`Z6ULK8U4H;Ur{yk*uO`H+K`#tK-^M zTo53P2+d^Y@Pr%gTVZqsbSTm`=S%8RkYU5lzquH>bx7a4xt9Mn;J4rpLB%JO`U!w# zuv~qMB~~jt8Fm(1YAFFO^5{jYc2uG7p&{hKWIkr8O@G+hKYO?^QbJ!P@GNIH$Su!( zexz9HW#r@W=}GBNXP-P@+_m!yeqMCLdV9t0>x27KM?320VuMiifeZ}jMFgAKJ?=&x z^rngPE(3N(l@Z&>k9Q+*Sif>uACHVUbk!6O_{FbWH6<7hrh5u zP;B2a0X-fVJ6%n{jrl~RHL$^exT?Ux#U=V;9r%%Z9cxN}xRoVcm8~o#L71hdc=xVk zm9)eyymt5BYd_90J;!fSOV{ew>(W+z#e_3+`ufeW6~PY)XjqNtg^?s!QVDJSk)&jb zhL)Dj#+H^qMYWcD#1=g&SU#@c2W2ycKN_Njil{=5rd-bc$jp!Y1#t8v%-4c*HI=({sDq+cO`X z{hcMguW==fUBVbdgGl_YQjtGBARNc{Ex+9R(HZY->=wywQHgJj*M5ge(V39!lsJ>% z|Gd=JM@$;hq%YD?A>n>WL|8)S(UtUNa>8@t{QXcohM(W2FMUFPI=_05b5JA^{qskU zHUdFKr^EX(hGUf3Xps~jA4W5XGMWM*{Afn^u^BX^0`1aNUNq?ObcDT z9OB`BdKUMKv9ZFfe*EK}CY?QIuxP+QEhT<>6XPeHWt_VR_zEx_G8~-mTt9I;I-fl& z;y%~;eeNJj&dl2DAUeQ8e9lsv7|GGLs78TG6E)hLydTW!jyB+Sjc9fyAN3DY8gpl-T6hhRE}f=#uvOy0%!7BXI+luOAE2TgT3S zs{ZcPpL3U_jU4GN)qK6=$Ml71XQYvw6c07;D_PHlk+i>*Gb0{%$xn05L6Nj@;7opo zz;Wf2->JEJz6^gtuUc;mMKFylVESb!bcR_yOK&k$dSfW0qvIZi?CmD#@vnC@FhH%Z zU)}m$bR2_*hPSN-^tix;f)_ zu8G*dq-|o9e}?GTUb#cF$wKOjNO4>=)EN?|7rE9uam|2_e&MjF)#ImT+4yd$xmX?-alExPK#QOTqn*_X$?{%6M3nLj~D`~t5E@> z>F0z&1k9{SnY<}dnB&G~fB*e@7Qb_-^rakLj5i^b-r?za-V(;ysX2G{Y8e4uZ3gpj zFQz+?g!%Fc>5n{&3bU%JD&@!7B8s_h0%y{u43zHf9w$EzuB6Oj+Ur?(CDM_-#4h&4 zCLFgGcE$2fKghlQ9_eF4?P^}IWB%63kXk%y+QvUE@Rsh&nXKffY`dX>X&v3t`D-_z zXP7y@P@S~eG)xznfvQeNC|lcsT^nO5-PJRZV@t}clHeC(F=TbmKKo-5vQ)xttRieX zl1xx83<9Iw)jRgphlMw)OxbHF9q|$je)eU>H_v442pzZsHSTP84D|dPqpt4ngaOQtgNTDNh4LKTbe;gegkGKC+YTsCzXa~P z^4Vk&4Xn)=HoZ456_j;nt~Y)>Y142#viR`6p5MxEB(u>#8bPnT)1r_m&lH4_6DjgK zN@y~H>TTS)vp=vfhrMbGFIqBN%CoyB5CVz*;3;_|+V(~3(qL6jFv;Fsg|TF1Q+O8U zsP$nPT}I8q#H?Bo-)u?TC(o|_FOpIOc}3=gWfEG2&2Aw(l_pk~@FCCe`RC!w2f~}DqCwT)E@>oG_VU^%6o81 zzFcCY+NPXDKG^7U*Y zRgdh?3K|srJrqhkUf-u0i45K{l|ip?-#9m2teNKip-?5zLAVPkx4}Rh zn`ib&%|jJ2BC@TaTN|QGQC7lJN>@hc+P)GFOGmI1lo@lSceE0oX^&<5zQz82cdBi@ zw%EP2KFQ(sd{LamCD|Zi+m5*UHO8hFGT;y#}QT)xNYZqd66K z=&AG)S;dftH0ig?Nd%3Jh1XneaOPr!)E<{PDA5ftJ~wnPSbFB&CtWPltE=3Y3;Ai| za`c{4MslfDGt%8DaX>-yXM7ASD{p_ZMMbkE@g2dr6htc&@Z(7~N+t z#^ON1PkozTRK?=G)%IZZWquzLlQr9L8m+|k?RQq0HJ=~+dJw;tIWaf?qa_O|%FNBO zOaKO3n@MkDr$^=1 zYx0WKa5vHp*Q;xf9i>!kKQYy%Z4c_fi^z$#9k(TOj;=@{xP7-LA%RG5>e7}Y=Bj<@ z&W#zAjc{*t*>>Vfw1%wyyOGDVlDP%j>2XES_--zT_XjHM4|a_<-6z?6_~(U$#fQGl zXWlQ_{2H3sIB zIKIR}rv9y)<8Rv;wPQ|SDmz+x#WCMuf(8ySy*linl`_8E7IRxb0e$E_I;3sb)#8;eN;x)MA!$6dZl;I z>lsdxgG3GJZ@)QYBx?m1f0%a{q`d;#kG)6BDEYKf{@o7j{SU#SWU+6)$tjGO& zlYecvnNjk9>mBUe>YZN{>P#vYbH{^C(}AFt2b$rQpNjCFDM@LFy5X823ThCvev7%! zgW*Yqhp=~EBiSna-nvrC>)~je>1*veUWBigDMC2i=Fj}Q)Z_p~>DRmMs*;(6GILoi zXgpb1^zRkO+vY3XaflnRjfuUgT3{TUo^Y8OFqoBNx6Bt$E9GFE6BBctzk#8E(S z%WgzDDkgjm%fF6~U1O0|Ub)v_@4DyktVe@PSbV=R=8pY+rF|zFu6Tx%Dm(=zcF~_b zH|fyBWKgf_1iQ67ZeB%YmKn;&JX>#%w3EensF07}X6`+=ug;a7%b2CuF_N!2DIIg* zs)%y%s;Fuk)4RojnUVa-iP<3NP5VM3)8iWZwMQUM#O>zQBE)B|Rep`;YqyRa#ps?u zy@x@>Y$?*u7*Y$~X>DgeaNRn|2gjgiNnYzQuZ|Egn5LE0^;N{3*vj0)r`U@f6;^i3 z6@nGsA%5>Z7qy5sEuW>E2TdXaTRv69DVNL1v=QFRXw8_e#rE6^xq@kBxk_!(QrAv6 zn7tu1!oQndu*3-JR*QMH)73xg#HRRT%Ez6|Mmp>?s+wZkoh2N?J79 zOYzD-)8arIil=j7d$J}@VhYcUq?_b|&wrFOV1?aMWd0XX@x@+37!jt4|q^6#jTZ2Sud0HZ7om}PGZa{5D+E65S zbGbpuAl)VZd9rM`>{L>No}E~=+`@!UuZa47aH9HKpC2#=_d2QfOp%P}?9(dQ>uHIN!?Jp{;{zv_Nze0!kWQ+|HD)?); z$L1w^3`Qe|v+ANqF$oY-uFgi@=~3Tx>-vB1O(ibt+=- zw7fWcDh?eZC#yP^C#SF5hs~<<6|ohNYI3Xg8IHp311RN`F>-mVIZ|Xk9}LF#m7NOF zeBT!cQEZIRKK0LO)dw~#c{gjkoz9l0_dAG^sm%6nPDH)g6eh9DH8;`w+#{h#OICd zZEt*B$4L{h!HpStKWTvA`lXl1V7<|Q9HJun!1D*hM(dvTU^j*4{`AV4celUD-Y@es ztBz8Wj{BdVY#U0&>FejwJmmRt@?O7ZQ^a{BrR4N7EuksRteRNKjjQIlx6kc4w=-Sj zP0~e`gUp{kdftE43Ksod9`&k`XG}36u|MYhbKdKLp56vRn6bJXxkPipQH62_tH$)$ zKKGtETbb9+Z*CHJdw@e8APYv`MEr<<^(lSdlStnDeCg@Pt|;p#y9Bb(=PGwIoz_+oCX9BWN|=7!~8);Y9-s+dSoDwg)jOb>kffu z=WuzrTeaEy%RPp@Ocb4T<=5{Bg`2h}DdO|g&;Rz!(9-*f+1@dD$;~*f<~p=B-(luE z(RC(p`fj^jH+15-c-HvwN+M5 z^!#_+6gSm3uNsRdR_3J&7T5Q9zcp9m!I~LpcRV`$Ia~D@ z?h=0A%c;B!yq7WIO*@V-)OXWRdBHR&;yId}XQWczm`iR=&FekgytK{4^D?)df$x`+ zXeP0ml=rrBHQNCiOI`FCxWibYAQ&@dn9C^o}WKKezuyh%VZD$y+*;{*`Gi zT;{lf$alM4Mw&|$Y|!YEW%KHg@u>{4jqp{ok0}&}PlTR45rWp8TC>lYGE^@#tyi%^jF7WtjRNsZ%(}y0#KPy|&hcZ5y_62t(9d>5j*9rS{KJtuoB zzPjO?yssY|Fdq?fZfHs>pjxMDnd+s|PLyX`m9SAD?d$Ewcrv8Jo-v z6!+w$EXB=>2xN76Hk-OF@(+<%*{sUyhaZLy&viz9rcdXc;$L(V zh9u@Y(4y`Rt*|D4lTBr_Y}=;s|RM5`OH-b`eQRg5|&y} z|6$K%E6+o&Ytj>XU#$zvwoM+6x#V?F?3o%2@tpM5_rJB%JY~Grs5>TR^_`5G^Tpf; zGo@+?_cE?K-!$^*%H8E|oE-+WKYS$$rZqcS|C;fB9@f`Wv5TsTQ#h>-fSYnD*{mmo z>>;x8F!dHiy?9J6v*COijEF9LpXkN5eYHOC+Xm%S*8Bs!VU@v8DR(<% z-Wqi2Q}cZ2u|^Np_P$5j$!Rl}j7E~j$%+qi?CZZNYaVT+N>);PUVBBZY%WkdasKlY zE&7R#;g?H%1W!Ttouqpg2-51wnhZyNI#uwN2w?&}-ji-G5YMa)r9CTLb4?}rAVB|n z{=FLij*_oR7_W{v{P?bk0*_MIhGHsxG}hmwOd#U=ZRymgt&*-VR)1+5l}B2TG`S)a zd4vBDA=la2lGG2|9evA^!yxrt(P4cE<}`GP^?=niDJaInm9SIpWBISl9=e55rjKII zRtK1IDJs+12ZmPTYd2?hb~?vSB+1K-vQB%Y3aE)C=`YhQH%a&DXbrt_IJ0GsRw`=` zeC1ZDy8BEW$zMBBYI=CMANng&s9hkj7tv!zq>JYheM=ymNh`}`SoXnME=WLQJR`S* zZ0q#rjivQ%$@fT^7TYx&qVBMf0>6dp@)Gwj22{2f)C-pLPi;m%pK7A__@? ze^=YKhkU+vyjdV$?!m!#U1Mi<|HNSX{7WCzS|b|nkc?TRb^k7lmAjl4V$*XwR`&Bx z$1j+==pv!kUho>{&2N5|VQA3_x%wxU$Cm9D3)z!ej#zI=JojQedn^hXu?SH#qfNTc z>v@g2Skfr2T|WDkIZOECse^(~)!rmn2_A&}h8>U5>%NtMmvx3iKQ#2D2~o#?jTUN-CIaAZj-mh9!^SLt}>V@mqWF zX~}^hoHl5XwL+Dmw+L0^5qmubnidYjuWP-LI2z{}*GkVJn|>?D@Vlm$L)-mBc1h=!9*?d!AD zGRNQ%wLHsG@$sBCnZ;J$_2dDeC~>+(@Q7+JpW5{4SPOe#iyqwyp1qgMNZH%(F{KAD z95TYmYbs9y6jLP+j9Va_e)>rY$)%NM80Ffl^a*#+$4*9$4khpKV#znEn-@4rmo!#q z1){1aD{7S@MVX7k^ZIsuwascaKoRQ^w+4z+#hth7nyv&yZ`lOSLcitrbm)nMG+lRF#+z5URI`;g)#_sK(9CpR+6I%M zDJx?YBX0NTbpObB=CZ|EL;Bsm1l=-eF^!t78*hac9+B3kL;E^>I-tKLVA*qvq_z&V z;>T&yc^OkATuPoM{3-G~e>XG!!U1t>51Kx5_D1N*|ah6%;b|`PcD3R z=QwO2@AgPARe`gU!U4sk@EXGhvPb1RZ%eZ*Og1hHMh*37bT|et%-d2E|Dr21PR4BC zy`FjP?D^Z5ZJ&H%L>-oj-Y@3u-J`U%xZfQu<)_eB7E!d@o-(7&qu%WOHpckP)Brnk z#kk@Z=G$VlK^uFvR+>BpdZp{$gHfPlq6c0@5OR-Hds6yZY|Cn#W_G^}$0@xPw%bf_ z-rPdR`UU#CT@~G-iTF#s=|zD=bRi5=w5S*5@7itqt=&38j#=j0*QB};y=UE>Bh8Op zl!J!936xw@DfRJbjmuen;s4!%E+yy4aBM{ZQ)#wBf!$W7mfr#=ED1%ACwx-%)*pQ z(7Jk--FQ@kvRq6v8;1Q>IMEbC;cryjoBff6feno`h}x5qW4P7rvO&w9$2&ilYnBVY z?w#%0v7`zqRa2gRS8nwcv8p{)Jt*WCOy0L4;-2_@*rZdIjgTID^853NPC`txm!SUS z4^*DGeVkO{3Ngf+`=?6hrM=6wL6~I&`-6T>Vipi)_UX{gPr7f8?@Jyh@46- zpv2I@XJTXHr=vY0YV=W;52E-^{3{gG#|Lg> z+S}74=eNgJ7~4p`Uzm&u=MtN)(JOG&hI!7!3r8c)NOz^_h?B=tcoIrB z#=AOIS(n9KrrvH4yV!kpFVy~kNl?E_#Fasq#$lx;x)BFkG!^MOvd|YBMaZ4{tIN3d z_je9i+>2oq^|x?poSEJV)*$xme(+TJyT((Q{xQ*P!)Gj8-VL=H*hsR=aaFHB)!ad8 zjeU+wC`v{44x*o|aTh0y_P9%$0sG*D4q@0I^N=H5?j$}-TxPm&kr)$(=UkiEi(ACT*cs}5B725Kjrd|~um4wC1 z#GtQ^ifT`$QBWvs*`vUW`55tjCQK&!a~J3G`xFZXDJ`4FF6kiFshzH4C;k=60_7)C95n(^dJM_*q)!9-Ev@|BaF8*$ZX8UdgT)HJoFd%EkG1@$3 z=fiQXU{Ad-or+(TA8#oW8aPG~!uHEO>vN#_z_hEGTUE<0f?kKOZVJP(hG2R+kdm!3 zam8sRkNZ}U;l|~UFKd4FmNUM(-}kguFs=J#r}NeAoW0S5F4IS2QrQDVk$MQo=X~D@ z*Xz^==IZ1(Y29BFfX9z_%^I%O3J(z7bHJDIrGZ8bQr#^dop_yrnCZM3a-P}Af5u}O z!{O98A75+O=PYYg^i>CVY0G8(jyWW!Qbq&et(rYf%~ILax>=bhTmEAw@O4{nB|!+n z&V}+h1A}-LcfHs2R5$yMP9LY}rn73Zq2a+~Z)Lg6IjS6Uw^n{2i^2xo5;VWbDi4g6efk1499 z8L69mv?ZH)NjCKj`YWH>{|;Ddrt&c!?iI;7<7qVE=DEw&8Sp-squXL;F)n{`^+|%5 zsQL$b3-5t4n3Vgz<%ipK3FDQsxmw@pT-@2kC`i7%NqbNsWZH-BZeRIAKNAvf4>R0Z zFd~NN-s9&hc>3e97#*e)b6%k)Vpd=3COAA=!WCMt^Pr2{l{*~%IqFkyJ?c_*dq$&k z=2MFlW0HvWjQ(w*vN3;c{<9|^CKZ^}!r>`e5_ zJzSlkNbJ2U_w7^ z<1$oU66p&3fMfF_6e{&xl6hae+>cqJ(AIosskztWgJjuwGpk+GJ1Z{9Tq(>7VOnc( zZ5Y+jIm|MmJNWhGK3S?O-P$Hu>6dIlY>UMdadpsg=+|322uTYvLwiv|1<=-LgSoee z=Bp8+%f&>DJb9IGG3MbSsjzW5w6lS`!f=AZ&NeuOM(~PHO2>@&l85fwAJg9yo+J3o zvQSI(w7)(xN3uSBLQhN*A=IygJtXqfP>fnb&uDKIedJr!&3V)s@074cq-CNUG&XEE zA<7TD4D6sa{G4R%8FhRX7nF%d=|e%mSv_^f`3v}Yu=6zD_7`yCvEq%hVGumbih>^Uxs_y-xb7LKd-fk_Ak>8 zq)QDxgt?Owu5gW*mDWDsGhYg4`+cc zjrQ7 zKE^(>UG=kM7}&>}F!uVX=Goi;4VihEW1g$Mk-hxlVYk*j*sv8e3tUT{6Qmx$oTu>{D zn3XH3T1I-^@mm>VaYxFcJ{vLdzC_niYBhjrS+YOLVAPCWbgdSMRDJb(c21@Mo@-cJN+F!!95ENa(|6h|arIacOPBO{G|?2R^Uha6(@TJOAh z8{?JRTSd79j&Anwq1nX6!@77Cr2_AyECuuXw$&(nv1%0CnMz2#8A`wc8F;d0ISzwK z6-hlAh8YlI!M3!5ya7$$#jejhkLu)HX-&xV?M+Dbe8VR?HXlUowf4hDDL8qy zz3=9S*C(IKBZ@2Stoa&jpZgBPO%*~CUgCFQTms`TpY9)AS*F-^ITl@xO*D>@eETj_ zUaod!*Yx&tjc^JESCM0V)vkE2dw~&Zo+8H|d1YEKwnWn|36}`Uf4!rqE_Tlo;^Klm zPs_rD1eQ%$vg;n#CRB`j#Y}xX^|%Wa_MlGl@{fK8g5O;nl9P!F6_g;kLFWC9@`O?F zu}CcMG~%kd+>kpJG-XjjR5sUl{scj+Q%fIg_&J4Fxy)vRbh7vz9aIaCOSq=R#c7`b z8_kNq+{n(T@pijcBO&?vmuMq*a$=(SdQNIX;qmjKPA$Ivy_Zi@*WV54<8?MzypD|i zBAWYZz)R`2-QbqaaL;@IZ^ZGVv!=bIqW$cddp+RF*(0;eUB>y z1XfIQ>56me=*=qT{k*Q{L(Rqg3SKnG)yY$6q+!zJ7}vcg&E7xQ`$ZNi6|1J2@aE}x z_Y2RSZ!>Cr0=#9PlXWOG)Ltikb0^xqQ*W8inVYhi|J6d5jtx>Df1h7L<>6384ciOa zhH%p>d3FMW7~dyd&|QrNp{51}`EXYxUp859Zgf(}xa&}TV<5ll&s9k#S7wR#XH zpR`?8$K9-*c_W9D3+C@(hE-7PF5D+w%aj=*&cY9iYfI~!Q$IONS}<-Kl<_9)wP%7W zkeYH$9kU_qKY5pY?`Jy0I!j|5M$w-#8y?!rEVf->`7PXQsxIY|KKCD-ALfIyS5-tGWg3$y2y0nw z%XY_5RUE2S&`ly6OCQBMCykNk!NmmGbkwc=>EF=48kFPAyxo5(z(%Xq$X8@JdR*B+ z{0PIe_o1ISlfjWTO!)T`{LjHkzcfxGTB|eqrMq<$ci5lFY2Q_9o7M{QA}w;P%!m$9 z$kbB25B)x}2D?X7;+m$lSv0}i@-Ub|-#OGeW}Y}WcrSQj>TvT12ObT#W|RucueNIo zQWPV(RpClH^2#sPWWnSgZjs9#2fs7NYP$1a0-PGO(qsI*cxYR_RJHHznovKl6fxig zUfMJUZ`9FI9@aa!J9HmkhrcyUquFf}AbQBKlw82*pXv~_ejas=LVo`(II8mXz_UdL z#o~p7OuPntx6t`@!pzEN8Y*nx0#Eb9lHG(Tuj7AH;a2q#J8pXLWiQS*k`&LA{~*sG z1CyiqRgf!7kM{KjU!tOv?$(=M6T1ttAA}4x7bBTkYWsQa`~vsyv5SenKWCpr2plQU zORqJu$NQ6*LehN3q=oYs*%<@F}}IJ1(#`dR7LY7rhe4x z|1`>EG3L>4t}6SO6r&kd(FVnMh&70^7r4B+{;nr!>M4_v^d36UcW_w9a%C&D&&?Qf z`AprgWh_FoKCDjN3^Oxw!|jHIc5??Dx0qiCi!sinaVwa- z8_i2n3L_XYk8e0KiQ_o&FYCR3zos@Gtx)xu?K&k-ow-7!skkB?eXxf%;yQ()&*7=l->sZ~-b;V}-hr#jmd0L}gp{D02{RXkW4<9#& zRu5gyhY($@lOc&Z-3`;tL8ljO+TVTLX}|qqsZI zd+g?nbYCev7*}A4Vm{#Ve)hIbEyv}v?O_1{(W*B)NBbR#%aXKK$?uz(l0Wj!EhT0| ze4o{O1fNWGCqG#lcCls8A!rKd<^GQ0x_8=r+NqBK=P}K;lis)O zTATeB?4iy#>Oj)I?Uj#j9x85tYtEW(geSeEpFjlIP-xuxq_fdK`L)|L_`#)XMeaMX zS2=ZA?&=%)<>h+Py%45W$$HrQsARv$Cw(TE<#cZQyUGJaMlCsb-yz%tbjVHOgDHW> zd%;B~isiV9b$rMDG^KBUg;HJ1jLtp8niw^iJrtWmF1yn_Y>>2aAymHWE7LATlguse z&79~i!q2V!7T;vIL;OtDW&XuDu+7xf%(noCL zk||3M79I-O+HkJl%>PhKc+)a7$+_cGcgxxv8XHK*_o8nux1N-V9XnaY$;^J)dRCor z^6G3cszk@<>4>+-ulVVxRt>+eSNIu6PXpOAHgZREwT@y>`MKmmd2J2lB`*EOTuE|S z;FsO;=<>gfzdWu~1qTy7%aFXhBdNyDnG0REc#OYT(Z zx>seG<@b2a=olNuHtdFyyR1rZ6f0*qj5eDw>stLw{cF7}9TSyH%=9m+lk*e3_R3O& z_hmsh*H_aD{BzwB3OUbSo3utqV%#~^1!2Y#eVgFYZmssTUyS~9Up+%Ntw8l9Z>)cRv#Ks4MM$f<=~Yd?&>7WtqK)2Yz@%FV$CdGx zE%1F0?#ftu!ZJ|IVerd5NbQQb2~2G-Qt4e;k-^wr*et^}Kobj6Rd|w~@s;eS-FX^Z0YHxO1rqO!V zEa~Kwo`t}(Rk!JD4(nC18-p)KdTUnOJV;UK8(SfO0(_e4mL_v+Ei%7cZ;Q_hftZ%DNZx^==(BO@n4r(o0r$2 zr_QRja{Q9^x1@K2wIyee)L5YAMv3h@$2MsNdJK?3kN+rPh>He6 zQDCg{+MhlDiz`qN5}W*&5_thW7#Ib?uKrO1mR;Br3_@Xr*Z(L{7vc#6LBQC^KT0I@ zLOh{p7!o`2mlAyeTPO;S#PV(cGg0sh20@W95Y}shP>2?E0U#&>0mgQ20PI7d5f_(* z1KecT1eS$BFGK|jL&C8(e<{%yq5_2iAXNWRqQDoG5EL3qvGqrZxacPY464-IB4i>U z03*RzuU$f}yP%6ef}!9G76C-)+b6t%{k}_hUmOX$urC;OkvU*cl?ZxaVHgCBLpWd@6t^@0?m;1c zjl=w3ssQk#z!xXOaTyGZ1FXlTJQxV$zhxjO2o1tr4+=sc|2Y@f9|Twfhx}9cM+N|o z##s;w1_t6p0rZETIB5;YAg~KEFz|I5TmxZHF*pu2K~TUuT==0V5E_SvK-^#mG}iD0 zh#>?Q4iFp(`Ino2c?1+coazLEq9GTpMla@}PjCi3>R#4u$-a_+)>)gMpDKBsTICfCI1rXEzuajeB$f6A-X}uJRW| zXxIfAFcfkQ{%^`rp$8 z1Vh2mxa2^9(a3*7^cSU(NB|Gq$w1oR#sP-9P(%Pv;4oaa!B9x}MHw3R6a@@KV=g2t z3=O)F4}k1~ConYl0z(1f;g(-UY3@$+cV&b+!AUIHFa6AD4hx&ghBw*Zw2XGDvy)YLB|C{{)|Nqlh0%?Lk;GzKZ zb6{K&0x}d13V%}UVx#q^RmCAS3JwDQyQ%^B1Hy|Q1NGqoTqw{5!og7&3pLOXUO)>D zz;$swBn*zhu@DIyfVd$6$}?`$3qnFrIG%t2*@bQ$gao=C+|))w5#ayqf%%`d9vFxO zUl<4jU5GJoN5JiTKu93rFN6++ys$9}4DkQo4r4%o5b^Jr03iXxFLV$nFzRB#hXa>} ze=i7x6MCUyqTr|tT?!zB!dU;RfB#nlPzcz+^Aku297J-RtS(?=*xia=&v*9FZ};g1OBv3xTO^+z(^#{#y~d! z!i5#+Z!e@BAVWfM909{2sDD&K%wLyU@SjqM11lH_Md4o1fa(p!%@W{h1;wp;fDEX6 z|3(W8L=N=7asNB)U^E5+!Sw_JgcFIwYCwhn;mROLAWi0!IMmT#)@wlm3tO5J)s`KL$Ynjo?2*@>l#JNZ>rc zfd-HvaoGx71TNgR0oldc4|EW?XhA_hX8&Vjz;B>!M_`b+76L+m|Nr(5xE!GV%S9FF zoN>khS3um;4!9IvI3oZVS*sfV+<1vzj^LJ%-=E-qO$ HIr9Gpm(2eo diff --git a/doc/doc.tex b/doc/doc.tex index 5f20f08..3d663ab 100644 --- a/doc/doc.tex +++ b/doc/doc.tex @@ -249,13 +249,12 @@ wobei $\Pi_h$ die $L_2$ Projektion auf $P^0(\T_h)$ ist. \mu_h &= \norm{\varrho^{1/2}(\phi_{h/2} - \phi_h)}_{L^2(\Gamma)}\\ \tilde\mu_h &= \norm{\varrho^{1/2}(\phi_{h/2} - \Pi_h\phi_{h/2})}_{L^2(\Gamma)} \end{align} -Dann gilt: +Dann gilt auf isotropen Netzen: \begin{itemize} - \item Schätzer sind equivalent\\ $\eta_h \sim \tilde \eta_h \sim \mu_h \sim \tilde \mu_h$ + \item Schätzer sind equivalent\\ $\tilde \mu_h \leq \mu_h \leq \sqrt{2}C_3 \eta_h$ \item sie sind effizient \item sie sind unter Saturationsannahme auch zuverlässig \end{itemize} - \end{sat} \begin{bew} Siehe S.F. Paper $\mapsto$ THM 3.2 \& 3.4 diff --git a/src/A_plot.m b/src/A_plot.m deleted file mode 100644 index dffc3a9..0000000 --- a/src/A_plot.m +++ /dev/null @@ -1,35 +0,0 @@ -function A_plot() - -type2str = ['Analytisch ' ; 'Quad Element' ; 'Quad Achse '; 'Quad Seite ']; - -%%G_D -global G_D; - -[m n] = size(G_D); - -step = round(n/3); - - -if step<1 - disp ('Error: No Data to show.') -else - -figure(4) -loglog(G_D(:,1),[G_D(:,[3+(0:step-1)*3]),20*G_D(:,1).^(-1/2),6*G_D(:,1).^(-1/4),5*G_D(:,1).^(-3/4)],... - G_D(1:end-1,1), sqrt(G_D(2:end,[4+(0:step-1)*3])-G_D(1:end-1,[4+(0:step-1)*3])),'--') - -title('Fehler') -xlabel('Elemente'); -ylabel('Schaetzer'); -legend([['m';'m';'m'] type2str(G_D(1,[2+(0:step-1)*3])',:); 'N^{-1/2} '; 'N^{-1/4} '; 'N^{-3/4} ';['g';'g';'g'] type2str(G_D(1,[2+(0:step-1)*3])',:)] ,'location','best'); - - -figure(5) -loglog(G_D(:,1),G_D(:,[4+(0:step-1)*3])) -title('Energie Norm') -xlabel('Elemente'); -ylabel('eNorm'); -legend(type2str(G_D(1,[2+(0:step-1)*3])',:)); -end - -end \ No newline at end of file diff --git a/src/A_plots.m b/src/A_plots.m index 5c694fe..0b603b7 100644 --- a/src/A_plots.m +++ b/src/A_plots.m @@ -53,7 +53,7 @@ for i = 1:length(files) leg2 = {leg2{:}... ['min hmin/max hmax ' l0 l1{i}]... ['min hmax/max hmax ' l0 l1{i}]... - ['hmin/hmax ' l0 l1{i}]... + ['min hmin/hmax ' l0 l1{i}]... }'; sym = {sym{:} type2sym{data(1,[2+(i-1)*rows])}}' end @@ -71,7 +71,7 @@ else % sol = interp1(1./X((round(1)):(end),4)',G_D((round(1)):(end),4)',0,'spline') % sol = 8.28466; % LShape - sol = 4.60957; % Quad + sol = 4.609193; % Quad % G_D diff --git a/src/A_run.m b/src/A_run.m index be5943d..1476eef 100644 --- a/src/A_run.m +++ b/src/A_run.m @@ -1,8 +1,8 @@ -function A_run(file,times,mu,type,eta,eps,out) +function A_run(file,times,zeta,type,eta,eps,out) % A_AnIso(file,times,mu,type,eta,eps,out) % file - starting mesh % times - how often -% mu - Zulaessigkeitsbed +% zeta - Zulaessigkeitsbed % type - art des Tests % eta - alle verfeinern oder nur wichtige? % eps - wie verfeinern iso or aniso @@ -11,7 +11,7 @@ function A_run(file,times,mu,type,eta,eps,out) % P. Schaefer %Parameter testen -assert(length(mu)==max(type)-1||length(mu)==1,... +assert(length(zeta)==max(type)-1||length(zeta)==1,... 'Pleas set right type and mu parameters'); %Netz laden @@ -36,7 +36,7 @@ for i = 1:times end % Ein kompletter Verfeinerungschritt - A_step(mu,type,eta,eps,out); + A_step(zeta,type,eta,eps,out); %Zeit Speichern usedTime = G_T(end,:); diff --git a/src/A_step.m b/src/A_step.m index 27ef6d5..71657c1 100644 --- a/src/A_step.m +++ b/src/A_step.m @@ -38,7 +38,7 @@ time = zeros(1,3); %Flaecheninhalte Berechnen (rhs) % b = sqrt(sum(quadNorm(coo_fine,ele_fine,'w').^2,2)); - b_fine = area(coo_fine,ele_fine); + b_fine = areaQuad(coo_fine,ele_fine); tic %data -> ErgebnisMatrix @@ -50,7 +50,7 @@ time = zeros(1,3); %Alle MatrixBrechenungsArten mit dem selben Netz berechnen for i = 1:length(type) %Matrix aufbauen -> MEX - A_fine = mex_build_AU(coo_fine,ele_fine,zeta,type(i)); + A_fine = mex_build_V(coo_fine,ele_fine,zeta,type(i)); %Testet auf Fehlerhafte Einträge (NaN +/-Inf) [r c] = find(isnan(A_fine)~=isinf(A_fine)); @@ -64,9 +64,7 @@ time = zeros(1,3); %Lösung Berechnen x_fine = A_fine\b_fine; - - - [b,tnormal,s] = area(G_C,G_E); + [b,tnormal,s] = areaQuad(G_C,G_E); hmin = min(s,[],2); hmax = max(s,[],2); @@ -75,17 +73,17 @@ time = zeros(1,3); %Fehlerschätzer 2 aufbauen - A = mex_build_AU(G_C,G_E,zeta,type(i)); + A = mex_build_V(G_C,G_E,zeta,type(i)); x = A\b; xo_fine(f2s) = repmat(x,1,4); xd_fine = xo_fine'-x_fine; % |||h/2 -h||| - tmu2 = xd_fine'*A_fine*xd_fine; + eta = xd_fine'*A_fine*xd_fine; % \tilde \mu ( h/2 -h + L_2 ) - tmu3 = hmin.*b.*sum((x_fine(f2s)'-repmat(x',4,1)).^2)'/4; + mu = hmin.*b.*sum((x_fine(f2s)'-repmat(x',4,1)).^2)'/4; %Energienorm^2 Berechnen |||h||| & |||h/2||| @@ -97,7 +95,7 @@ time = zeros(1,3); - data = [data type(i) sqrt(sum(tmu)) sqrt(tmu2) xe_fine sqrt(sum(tmu3))... + data = [data type(i) sqrt(sum(tmu)) sqrt(eta) xe_fine sqrt(sum(mu))... min(hmin)/max(hmax) min(hmax)/max(hmax) sqrt(sum((hmin./hmax).^2))]; end time(2) = toc; diff --git a/src/Aelements.m b/src/Aelements.m deleted file mode 100644 index aad35f5..0000000 --- a/src/Aelements.m +++ /dev/null @@ -1,19 +0,0 @@ -function A = Aelements(coordinates, elements) -% This function calculates the Sum of all Surfaces -% -% A=Aelements(coordinates,elements) -% -% P.Schaefer - -elles = size(elements,1); -A=0; -for i = 1:elles - tri = elements(i,:); - a = (coordinates(tri(1),:)-coordinates(tri(2),:)); - b = (coordinates(tri(1),:)-coordinates(tri(4),:)); - A1 = norm(cross(a',b')); - A = A+ A1; -end - -end - diff --git a/src/area.m b/src/areaQuad.m similarity index 92% rename from src/area.m rename to src/areaQuad.m index f4c657a..8473339 100644 --- a/src/area.m +++ b/src/areaQuad.m @@ -1,4 +1,4 @@ -function [area, vector, sites] = area(coordinates, elements,varargin) +function [area, vector, sites] = areaQuad(coordinates, elements,varargin) % % [area, vector, sites] = area(coordinates, elements) % diff --git a/src/computeEstSlpMuTilde.m b/src/computeEstSlpMuTilde.m deleted file mode 100644 index dce7564..0000000 --- a/src/computeEstSlpMuTilde.m +++ /dev/null @@ -1,34 +0,0 @@ -function ind=computeEstSlpMuTilde(x_fine,coo,ele,f2s) -% ind=computeEstSlpMuTilde(x_fine,coo,ele,f2s) -% Berechnet den Fehlerschätzer über die Elemente -% -% x_fine - Phi -% coo - coordinaten -% ele - elemente -% f2s - Vater Sohn beziehungen -% -% P. Schaefer - -xF2S = x_fine(f2s)'; - -if(size(xF2S,1)==1) - xF2S = xF2S'; -end - -tmp = abs([sum(coo(ele(:,2),:)-coo(ele(:,1),:),2)... - sum(coo(ele(:,4),:)-coo(ele(:,1),:),2)]); -hmin = min(tmp,[],2); -% hmax = max(tmp,[],2); - -% len = min(abs([sum(coo(ele(:,2),:)-coo(ele(:,1),:),2)... -% sum(coo(ele(:,4),:)-coo(ele(:,1),:),2)]),[],2); - -sizes = sqrt(sum(quadNorm(coo,ele,'w').^2,2)); - -% max(sizes ~= hmin.*hmax) - -x_var = sum((xF2S-repmat(sum(xF2S,1)/4,4,1)).^2)'; - -ind = hmin.* sizes .* x_var/4; - -end \ No newline at end of file diff --git a/src/expmlAA_2DQuad_error.eps b/src/expmlAA_2DQuad_error.eps deleted file mode 100644 index 7a86d86..0000000 --- a/src/expmlAA_2DQuad_error.eps +++ /dev/null @@ -1,824 +0,0 @@ -%!PS-Adobe-2.0 EPSF-1.2 -%%Creator: MATLAB, The MathWorks, Inc. Version 7.12.0.635 (R2011a). Operating System: Linux 3.0.0-16-server #28-Ubuntu SMP Fri Jan 27 18:03:45 UTC 2012 x86_64. -%%Title: ./expmlAA_2DQuad_error.eps -%%CreationDate: 03/26/2012 23:43:43 -%%DocumentNeededFonts: Helvetica -%%DocumentProcessColors: Cyan Magenta Yellow Black -%%Extensions: CMYK -%%Pages: 1 -%%BoundingBox: 53 195 550 602 -%%EndComments - -%%BeginProlog -% MathWorks dictionary -/MathWorks 160 dict begin -% definition operators -/bdef {bind def} bind def -/ldef {load def} bind def -/xdef {exch def} bdef -/xstore {exch store} bdef -% operator abbreviations -/c /clip ldef -/cc /concat ldef -/cp /closepath ldef -/gr /grestore ldef -/gs /gsave ldef -/mt /moveto ldef -/np /newpath ldef -/cm /currentmatrix ldef -/sm /setmatrix ldef -/rm /rmoveto ldef -/rl /rlineto ldef -/s {show newpath} bdef -/sc {setcmykcolor} bdef -/sr /setrgbcolor ldef -/sg /setgray ldef -/w /setlinewidth ldef -/j /setlinejoin ldef -/cap /setlinecap ldef -/rc {rectclip} bdef -/rf {rectfill} bdef -% page state control -/pgsv () def -/bpage {/pgsv save def} bdef -/epage {pgsv restore} bdef -/bplot /gsave ldef -/eplot {stroke grestore} bdef -% orientation switch -/portraitMode 0 def /landscapeMode 1 def /rotateMode 2 def -% coordinate system mappings -/dpi2point 0 def -% font control -/FontSize 0 def -/FMS {/FontSize xstore findfont [FontSize 0 0 FontSize neg 0 0] - makefont setfont} bdef -/reencode {exch dup where {pop load} {pop StandardEncoding} ifelse - exch dup 3 1 roll findfont dup length dict begin - { 1 index /FID ne {def}{pop pop} ifelse } forall - /Encoding exch def currentdict end definefont pop} bdef -/isroman {findfont /CharStrings get /Agrave known} bdef -/FMSR {3 1 roll 1 index dup isroman {reencode} {pop pop} ifelse - exch FMS} bdef -/csm {1 dpi2point div -1 dpi2point div scale neg translate - dup landscapeMode eq {pop -90 rotate} - {rotateMode eq {90 rotate} if} ifelse} bdef -% line types: solid, dotted, dashed, dotdash -/SO { [] 0 setdash } bdef -/DO { [.5 dpi2point mul 4 dpi2point mul] 0 setdash } bdef -/DA { [6 dpi2point mul] 0 setdash } bdef -/DD { [.5 dpi2point mul 4 dpi2point mul 6 dpi2point mul 4 - dpi2point mul] 0 setdash } bdef -% macros for lines and objects -/L {lineto stroke} bdef -/MP {3 1 roll moveto 1 sub {rlineto} repeat} bdef -/AP {{rlineto} repeat} bdef -/PDlw -1 def -/W {/PDlw currentlinewidth def setlinewidth} def -/PP {closepath eofill} bdef -/DP {closepath stroke} bdef -/MR {4 -2 roll moveto dup 0 exch rlineto exch 0 rlineto - neg 0 exch rlineto closepath} bdef -/FR {MR stroke} bdef -/PR {MR fill} bdef -/L1i {{currentfile picstr readhexstring pop} image} bdef -/tMatrix matrix def -/MakeOval {newpath tMatrix currentmatrix pop translate scale -0 0 1 0 360 arc tMatrix setmatrix} bdef -/FO {MakeOval stroke} bdef -/PO {MakeOval fill} bdef -/PD {currentlinewidth 2 div 0 360 arc fill - PDlw -1 eq not {PDlw w /PDlw -1 def} if} def -/FA {newpath tMatrix currentmatrix pop translate scale - 0 0 1 5 -2 roll arc tMatrix setmatrix stroke} bdef -/PA {newpath tMatrix currentmatrix pop translate 0 0 moveto scale - 0 0 1 5 -2 roll arc closepath tMatrix setmatrix fill} bdef -/FAn {newpath tMatrix currentmatrix pop translate scale - 0 0 1 5 -2 roll arcn tMatrix setmatrix stroke} bdef -/PAn {newpath tMatrix currentmatrix pop translate 0 0 moveto scale - 0 0 1 5 -2 roll arcn closepath tMatrix setmatrix fill} bdef -/vradius 0 def /hradius 0 def /lry 0 def -/lrx 0 def /uly 0 def /ulx 0 def /rad 0 def -/MRR {/vradius xdef /hradius xdef /lry xdef /lrx xdef /uly xdef - /ulx xdef newpath tMatrix currentmatrix pop ulx hradius add uly - vradius add translate hradius vradius scale 0 0 1 180 270 arc - tMatrix setmatrix lrx hradius sub uly vradius add translate - hradius vradius scale 0 0 1 270 360 arc tMatrix setmatrix - lrx hradius sub lry vradius sub translate hradius vradius scale - 0 0 1 0 90 arc tMatrix setmatrix ulx hradius add lry vradius sub - translate hradius vradius scale 0 0 1 90 180 arc tMatrix setmatrix - closepath} bdef -/FRR {MRR stroke } bdef -/PRR {MRR fill } bdef -/MlrRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lry uly sub 2 div def - newpath tMatrix currentmatrix pop ulx rad add uly rad add translate - rad rad scale 0 0 1 90 270 arc tMatrix setmatrix lrx rad sub lry rad - sub translate rad rad scale 0 0 1 270 90 arc tMatrix setmatrix - closepath} bdef -/FlrRR {MlrRR stroke } bdef -/PlrRR {MlrRR fill } bdef -/MtbRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lrx ulx sub 2 div def - newpath tMatrix currentmatrix pop ulx rad add uly rad add translate - rad rad scale 0 0 1 180 360 arc tMatrix setmatrix lrx rad sub lry rad - sub translate rad rad scale 0 0 1 0 180 arc tMatrix setmatrix - closepath} bdef -/FtbRR {MtbRR stroke } bdef -/PtbRR {MtbRR fill } bdef -/stri 6 array def /dtri 6 array def -/smat 6 array def /dmat 6 array def -/tmat1 6 array def /tmat2 6 array def /dif 3 array def -/asub {/ind2 exch def /ind1 exch def dup dup - ind1 get exch ind2 get sub exch } bdef -/tri_to_matrix { - 2 0 asub 3 1 asub 4 0 asub 5 1 asub - dup 0 get exch 1 get 7 -1 roll astore } bdef -/compute_transform { - dmat dtri tri_to_matrix tmat1 invertmatrix - smat stri tri_to_matrix tmat2 concatmatrix } bdef -/ds {stri astore pop} bdef -/dt {dtri astore pop} bdef -/db {2 copy /cols xdef /rows xdef mul dup 3 mul string - currentfile exch readhexstring pop - dup 0 3 index getinterval /rbmap xdef - dup 2 index dup getinterval /gbmap xdef - 1 index dup 2 mul exch getinterval /bbmap xdef pop pop}bdef -/it {gs np dtri aload pop moveto lineto lineto cp c - cols rows 8 compute_transform - rbmap gbmap bbmap true 3 colorimage gr}bdef -/il {newpath moveto lineto stroke}bdef -currentdict end def -%%EndProlog - -%%BeginSetup -MathWorks begin - -0 cap - -end -%%EndSetup - -%%Page: 1 1 -%%BeginPageSetup -%%PageBoundingBox: 53 195 550 602 -MathWorks begin -bpage -%%EndPageSetup - -%%BeginObject: obj1 -bplot - -/dpi2point 8.33333 def -portraitMode 0150 5100 csm - - 294 79 4147 3390 MR c np -76 dict begin %Colortable dictionary -/c0 { 0.000000 0.000000 0.000000 sr} bdef -/c1 { 1.000000 1.000000 1.000000 sr} bdef -/c2 { 0.900000 0.000000 0.000000 sr} bdef -/c3 { 0.000000 0.820000 0.000000 sr} bdef -/c4 { 0.000000 0.000000 0.800000 sr} bdef -/c5 { 0.910000 0.820000 0.320000 sr} bdef -/c6 { 1.000000 0.260000 0.820000 sr} bdef -/c7 { 0.000000 0.820000 0.820000 sr} bdef -c0 -1 j -1 sg - 0 0 4803 3605 PR -4.16667 w -0 2938 3722 0 0 -2938 624 3208 4 MP -PP --3722 0 0 2938 3722 0 0 -2938 624 3208 5 MP stroke -2.77778 w -DO -SO -4.16667 w -0 sg - 624 3208 mt 4346 3208 L - 624 270 mt 4346 270 L - 624 3208 mt 624 270 L -4346 3208 mt 4346 270 L - 624 3208 mt 4346 3208 L - 624 3208 mt 624 270 L - 624 3208 mt 624 3189 L - 624 270 mt 624 288 L - 624 3208 mt 624 3170 L - 624 270 mt 624 307 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 563 3334 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - - 655 3283 mt -(0) s - 904 3208 mt 904 3189 L - 904 270 mt 904 288 L -1067 3208 mt 1067 3189 L -1067 270 mt 1067 288 L -1184 3208 mt 1184 3189 L -1184 270 mt 1184 288 L -1274 3208 mt 1274 3189 L -1274 270 mt 1274 288 L -1348 3208 mt 1348 3189 L -1348 270 mt 1348 288 L -1410 3208 mt 1410 3189 L -1410 270 mt 1410 288 L -1464 3208 mt 1464 3189 L -1464 270 mt 1464 288 L -1511 3208 mt 1511 3189 L -1511 270 mt 1511 288 L -1554 3208 mt 1554 3189 L -1554 270 mt 1554 288 L -1554 3208 mt 1554 3170 L -1554 270 mt 1554 307 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - -1493 3334 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - -1585 3283 mt -(1) s -1834 3208 mt 1834 3189 L -1834 270 mt 1834 288 L -1998 3208 mt 1998 3189 L -1998 270 mt 1998 288 L -2114 3208 mt 2114 3189 L -2114 270 mt 2114 288 L -2204 3208 mt 2204 3189 L -2204 270 mt 2204 288 L -2278 3208 mt 2278 3189 L -2278 270 mt 2278 288 L -2340 3208 mt 2340 3189 L -2340 270 mt 2340 288 L -2394 3208 mt 2394 3189 L -2394 270 mt 2394 288 L -2442 3208 mt 2442 3189 L -2442 270 mt 2442 288 L -2485 3208 mt 2485 3189 L -2485 270 mt 2485 288 L -2485 3208 mt 2485 3170 L -2485 270 mt 2485 307 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - -2424 3334 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - -2516 3283 mt -(2) s -2765 3208 mt 2765 3189 L -2765 270 mt 2765 288 L -2928 3208 mt 2928 3189 L -2928 270 mt 2928 288 L -3045 3208 mt 3045 3189 L -3045 270 mt 3045 288 L -3135 3208 mt 3135 3189 L -3135 270 mt 3135 288 L -3209 3208 mt 3209 3189 L -3209 270 mt 3209 288 L -3271 3208 mt 3271 3189 L -3271 270 mt 3271 288 L -3325 3208 mt 3325 3189 L -3325 270 mt 3325 288 L -3372 3208 mt 3372 3189 L -3372 270 mt 3372 288 L -3415 3208 mt 3415 3189 L -3415 270 mt 3415 288 L -3415 3208 mt 3415 3170 L -3415 270 mt 3415 307 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - -3354 3334 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - -3446 3283 mt -(3) s -3695 3208 mt 3695 3189 L -3695 270 mt 3695 288 L -3859 3208 mt 3859 3189 L -3859 270 mt 3859 288 L -3975 3208 mt 3975 3189 L -3975 270 mt 3975 288 L -4065 3208 mt 4065 3189 L -4065 270 mt 4065 288 L -4139 3208 mt 4139 3189 L -4139 270 mt 4139 288 L -4201 3208 mt 4201 3189 L -4201 270 mt 4201 288 L -4255 3208 mt 4255 3189 L -4255 270 mt 4255 288 L -4303 3208 mt 4303 3189 L -4303 270 mt 4303 288 L -4346 3208 mt 4346 3189 L -4346 270 mt 4346 288 L -4346 3208 mt 4346 3170 L -4346 270 mt 4346 307 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - -4285 3334 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - -4377 3283 mt -(4) s - 624 3208 mt 642 3208 L -4346 3208 mt 4327 3208 L - 624 3208 mt 661 3208 L -4346 3208 mt 4308 3208 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 445 3238 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - - 537 3187 mt -(-3) s - 624 2986 mt 642 2986 L -4346 2986 mt 4327 2986 L - 624 2857 mt 642 2857 L -4346 2857 mt 4327 2857 L - 624 2765 mt 642 2765 L -4346 2765 mt 4327 2765 L - 624 2694 mt 642 2694 L -4346 2694 mt 4327 2694 L - 624 2636 mt 642 2636 L -4346 2636 mt 4327 2636 L - 624 2587 mt 642 2587 L -4346 2587 mt 4327 2587 L - 624 2544 mt 642 2544 L -4346 2544 mt 4327 2544 L - 624 2507 mt 642 2507 L -4346 2507 mt 4327 2507 L - 624 2473 mt 642 2473 L -4346 2473 mt 4327 2473 L - 624 2473 mt 661 2473 L -4346 2473 mt 4308 2473 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 445 2503 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - - 537 2452 mt -(-2) s - 624 2252 mt 642 2252 L -4346 2252 mt 4327 2252 L - 624 2123 mt 642 2123 L -4346 2123 mt 4327 2123 L - 624 2031 mt 642 2031 L -4346 2031 mt 4327 2031 L - 624 1960 mt 642 1960 L -4346 1960 mt 4327 1960 L - 624 1901 mt 642 1901 L -4346 1901 mt 4327 1901 L - 624 1852 mt 642 1852 L -4346 1852 mt 4327 1852 L - 624 1810 mt 642 1810 L -4346 1810 mt 4327 1810 L - 624 1772 mt 642 1772 L -4346 1772 mt 4327 1772 L - 624 1739 mt 642 1739 L -4346 1739 mt 4327 1739 L - 624 1739 mt 661 1739 L -4346 1739 mt 4308 1739 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 445 1769 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - - 537 1718 mt -(-1) s - 624 1517 mt 642 1517 L -4346 1517 mt 4327 1517 L - 624 1388 mt 642 1388 L -4346 1388 mt 4327 1388 L - 624 1296 mt 642 1296 L -4346 1296 mt 4327 1296 L - 624 1225 mt 642 1225 L -4346 1225 mt 4327 1225 L - 624 1167 mt 642 1167 L -4346 1167 mt 4327 1167 L - 624 1118 mt 642 1118 L -4346 1118 mt 4327 1118 L - 624 1075 mt 642 1075 L -4346 1075 mt 4327 1075 L - 624 1038 mt 642 1038 L -4346 1038 mt 4327 1038 L - 624 1004 mt 642 1004 L -4346 1004 mt 4327 1004 L - 624 1004 mt 661 1004 L -4346 1004 mt 4308 1004 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 445 1034 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - - 537 983 mt -(0) s - 624 783 mt 642 783 L -4346 783 mt 4327 783 L - 624 654 mt 642 654 L -4346 654 mt 4327 654 L - 624 562 mt 642 562 L -4346 562 mt 4327 562 L - 624 491 mt 642 491 L -4346 491 mt 4327 491 L - 624 432 mt 642 432 L -4346 432 mt 4327 432 L - 624 383 mt 642 383 L -4346 383 mt 4327 383 L - 624 341 mt 642 341 L -4346 341 mt 4327 341 L - 624 303 mt 642 303 L -4346 303 mt 4327 303 L - 624 270 mt 642 270 L -4346 270 mt 4327 270 L - 624 270 mt 661 270 L -4346 270 mt 4308 270 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 445 300 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - - 537 249 mt -(1) s - 624 3208 mt 4346 3208 L - 624 270 mt 4346 270 L - 624 3208 mt 624 270 L -4346 3208 mt 4346 270 L -gs 624 270 3723 2939 MR c np -/c8 { 0.000000 0.000000 1.000000 sr} bdef -c8 -68 -537 100 -204 65 7 95 -61 107 14 100 72 97 70 93 67 -87 68 72 63 55 60 84 58 55 56 64 55 60 55 60 57 -62 48 45 49 62 50 48 50 51 53 71 51 51 52 90 53 -71 56 104 77 122 51 190 65 254 54 560 71 1184 933 31 MP stroke -gs 1133 882 3146 1578 MR c np - 25 25 1184 933 FO - 25 25 1744 1004 FO - 25 25 1998 1058 FO - 25 25 2188 1123 FO - 25 25 2310 1174 FO - 25 25 2414 1251 FO - 25 25 2485 1307 FO - 25 25 2575 1360 FO - 25 25 2626 1412 FO - 25 25 2697 1463 FO - 25 25 2748 1516 FO - 25 25 2796 1566 FO - 25 25 2858 1616 FO - 25 25 2903 1665 FO - 25 25 2965 1713 FO - 25 25 3025 1770 FO - 25 25 3085 1825 FO - 25 25 3149 1880 FO - 25 25 3204 1936 FO - 25 25 3288 1994 FO - 25 25 3343 2054 FO - 25 25 3415 2117 FO - 25 25 3502 2185 FO - 25 25 3595 2252 FO - 25 25 3692 2322 FO - 25 25 3792 2394 FO - 25 25 3899 2408 FO - 25 25 3994 2347 FO - 25 25 4059 2354 FO - 25 25 4159 2150 FO - 25 25 4227 1613 FO -gr - -/c9 { 0.000000 0.500000 0.000000 sr} bdef -c9 -68 5 100 11 65 11 95 56 107 64 100 52 97 36 93 52 -87 44 72 29 55 33 84 52 55 40 64 43 60 46 60 50 -62 48 45 50 62 46 48 52 51 47 71 54 51 54 90 46 -71 61 104 69 122 59 190 59 254 59 560 92 1184 1285 31 MP stroke -gs 1133 1234 3146 1523 MR c np - 25 25 1184 1285 FO - 25 25 1744 1377 FO - 25 25 1998 1436 FO - 25 25 2188 1495 FO - 25 25 2310 1554 FO - 25 25 2414 1623 FO - 25 25 2485 1684 FO - 25 25 2575 1730 FO - 25 25 2626 1784 FO - 25 25 2697 1838 FO - 25 25 2748 1885 FO - 25 25 2796 1937 FO - 25 25 2858 1983 FO - 25 25 2903 2033 FO - 25 25 2965 2081 FO - 25 25 3025 2131 FO - 25 25 3085 2177 FO - 25 25 3149 2220 FO - 25 25 3204 2260 FO - 25 25 3288 2312 FO - 25 25 3343 2345 FO - 25 25 3415 2374 FO - 25 25 3502 2418 FO - 25 25 3595 2470 FO - 25 25 3692 2506 FO - 25 25 3792 2558 FO - 25 25 3899 2622 FO - 25 25 3994 2678 FO - 25 25 4059 2689 FO - 25 25 4159 2700 FO - 25 25 4227 2705 FO -gr - -/c10 { 1.000000 0.000000 0.000000 sr} bdef -c10 -68 1 100 0 65 0 95 2 107 2 100 2 97 3 93 5 -87 6 72 7 55 12 84 19 55 23 64 28 60 33 60 38 -62 49 45 49 62 51 48 51 51 51 71 53 51 52 90 53 -71 57 104 76 122 56 190 68 254 64 560 97 1184 1252 31 MP stroke -gs 1133 1201 3146 1111 MR c np - 25 25 1184 1252 FO - 25 25 1744 1349 FO - 25 25 1998 1413 FO - 25 25 2188 1481 FO - 25 25 2310 1537 FO - 25 25 2414 1613 FO - 25 25 2485 1670 FO - 25 25 2575 1723 FO - 25 25 2626 1775 FO - 25 25 2697 1828 FO - 25 25 2748 1879 FO - 25 25 2796 1930 FO - 25 25 2858 1981 FO - 25 25 2903 2030 FO - 25 25 2965 2079 FO - 25 25 3025 2117 FO - 25 25 3085 2150 FO - 25 25 3149 2178 FO - 25 25 3204 2201 FO - 25 25 3288 2220 FO - 25 25 3343 2232 FO - 25 25 3415 2239 FO - 25 25 3502 2245 FO - 25 25 3595 2250 FO - 25 25 3692 2253 FO - 25 25 3792 2255 FO - 25 25 3899 2257 FO - 25 25 3994 2259 FO - 25 25 4059 2259 FO - 25 25 4159 2259 FO - 25 25 4227 2260 FO -gr - -/c11 { 0.000000 0.750000 0.750000 sr} bdef -c11 -68 -525 100 118 65 -379 95 -65 107 14 100 70 97 70 93 63 -87 68 72 62 55 59 84 59 55 56 64 55 60 55 60 57 -62 49 45 48 62 51 48 50 51 54 71 51 51 51 90 54 -71 55 104 80 122 48 190 66 254 51 560 70 1184 931 31 MP stroke -gs 1133 880 3146 1569 MR c np - 25 25 1184 931 FO - 25 25 1744 1001 FO - 25 25 1998 1052 FO - 25 25 2188 1118 FO - 25 25 2310 1166 FO - 25 25 2414 1246 FO - 25 25 2485 1301 FO - 25 25 2575 1355 FO - 25 25 2626 1406 FO - 25 25 2697 1457 FO - 25 25 2748 1511 FO - 25 25 2796 1561 FO - 25 25 2858 1612 FO - 25 25 2903 1660 FO - 25 25 2965 1709 FO - 25 25 3025 1766 FO - 25 25 3085 1821 FO - 25 25 3149 1876 FO - 25 25 3204 1932 FO - 25 25 3288 1991 FO - 25 25 3343 2050 FO - 25 25 3415 2112 FO - 25 25 3502 2180 FO - 25 25 3595 2243 FO - 25 25 3692 2313 FO - 25 25 3792 2383 FO - 25 25 3899 2397 FO - 25 25 3994 2332 FO - 25 25 4059 1953 FO - 25 25 4159 2071 FO - 25 25 4227 1546 FO -gr - -DD -c8 -68 27 100 40 65 26 95 37 107 42 100 39 97 39 93 37 -87 34 72 29 55 21 84 33 55 22 64 25 60 24 60 24 -62 24 45 18 62 24 48 19 51 21 71 27 51 21 90 35 -71 28 104 41 122 48 190 75 254 101 560 221 1184 604 31 MP stroke -c9 -68 14 100 20 65 12 95 19 107 21 100 20 97 19 93 18 -87 18 72 14 55 11 84 16 55 11 64 13 60 12 60 11 -62 13 45 8 62 13 48 9 51 10 71 14 51 10 90 18 -71 14 104 21 122 24 190 37 254 50 560 111 1184 764 31 MP stroke -c10 -68 40 100 59 65 39 95 56 107 64 100 58 97 58 93 55 -87 52 72 42 55 33 84 50 55 32 64 38 60 35 60 36 -62 36 45 27 62 37 48 28 51 31 71 41 51 31 90 53 -71 42 104 61 122 73 190 112 254 151 560 331 1184 1115 31 MP stroke -gr - -c10 -DD -0 sg -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - -2367 206 mt -(Fehler) s -2308 3434 mt -(Elemente) s - 388 1927 mt -90 rotate -(Schaetzer) s -90 rotate - 612 3239 mt -( ) s -4335 300 mt -( ) s -SO -1 sg -0 854 1238 0 0 -854 665 3168 4 MP -PP --1238 0 0 854 1238 0 0 -854 665 3168 5 MP stroke -2.77778 w -DO -SO -4.16667 w -0 sg - 665 3168 mt 1903 3168 L - 665 3168 mt 665 2314 L - 998 2414 mt -(tilde ) s -%%IncludeResource: font Symbol -/Symbol /ISOLatin1Encoding 83.3333 FMSR - -1174 2414 mt -(m) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - -1222 2414 mt -( testAA Analytisch) s -gs 665 2314 1239 855 MR c np -c8 -256 0 716 2383 2 MP stroke -gs 793 2332 103 103 MR c np - 25 25 844 2383 FO -gr - -gr - -c8 -0 sg -%%IncludeResource: font Symbol -/Symbol /ISOLatin1Encoding 83.3333 FMSR - - 998 2527 mt -(h) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - -1048 2527 mt -( testAA Analytisch) s -gs 665 2314 1239 855 MR c np -c9 -256 0 716 2496 2 MP stroke -gs 793 2445 103 103 MR c np - 25 25 844 2496 FO -gr - -gr - -c9 -0 sg - 998 2635 mt -(error testAA Analytisch) s -gs 665 2314 1239 855 MR c np -c10 -256 0 716 2605 2 MP stroke -gs 793 2554 103 103 MR c np - 25 25 844 2605 FO -gr - -gr - -c10 -0 sg -%%IncludeResource: font Symbol -/Symbol /ISOLatin1Encoding 83.3333 FMSR - - 998 2744 mt -(m) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - -1046 2744 mt -( testAA Analytisch) s -gs 665 2314 1239 855 MR c np -c11 -256 0 716 2713 2 MP stroke -gs 793 2662 103 103 MR c np - 25 25 844 2713 FO -gr - -gr - -c11 -0 sg - 998 2875 mt -(N) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 66.6667 FMSR - -1058 2834 mt -(-1/2) s -gs 665 2314 1239 855 MR c np -DD -c8 -256 0 716 2834 2 MP stroke -SO -gr - -c8 -0 sg -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 998 3003 mt -(N) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 66.6667 FMSR - -1058 2962 mt -(-1/4) s -gs 665 2314 1239 855 MR c np -DD -c9 -256 0 716 2962 2 MP stroke -SO -gr - -c9 -0 sg -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 998 3131 mt -(N) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 66.6667 FMSR - -1058 3090 mt -(-3/4) s -gs 665 2314 1239 855 MR c np -DD -c10 -256 0 716 3090 2 MP stroke -SO -gr - -c10 - -end %%Color Dict - -eplot -%%EndObject - -epage -end - -showpage - -%%Trailer -%%EOF diff --git a/src/expmlAA_2DQuad_hminmax.eps b/src/expmlAA_2DQuad_hminmax.eps deleted file mode 100644 index 892747d..0000000 --- a/src/expmlAA_2DQuad_hminmax.eps +++ /dev/null @@ -1,865 +0,0 @@ -%!PS-Adobe-2.0 EPSF-1.2 -%%Creator: MATLAB, The MathWorks, Inc. Version 7.12.0.635 (R2011a). Operating System: Linux 3.0.0-16-server #28-Ubuntu SMP Fri Jan 27 18:03:45 UTC 2012 x86_64. -%%Title: ./expmlAA_2DQuad_hminmax.eps -%%CreationDate: 03/26/2012 23:43:44 -%%DocumentNeededFonts: Helvetica -%%DocumentProcessColors: Cyan Magenta Yellow Black -%%Extensions: CMYK -%%Pages: 1 -%%BoundingBox: 53 195 550 602 -%%EndComments - -%%BeginProlog -% MathWorks dictionary -/MathWorks 160 dict begin -% definition operators -/bdef {bind def} bind def -/ldef {load def} bind def -/xdef {exch def} bdef -/xstore {exch store} bdef -% operator abbreviations -/c /clip ldef -/cc /concat ldef -/cp /closepath ldef -/gr /grestore ldef -/gs /gsave ldef -/mt /moveto ldef -/np /newpath ldef -/cm /currentmatrix ldef -/sm /setmatrix ldef -/rm /rmoveto ldef -/rl /rlineto ldef -/s {show newpath} bdef -/sc {setcmykcolor} bdef -/sr /setrgbcolor ldef -/sg /setgray ldef -/w /setlinewidth ldef -/j /setlinejoin ldef -/cap /setlinecap ldef -/rc {rectclip} bdef -/rf {rectfill} bdef -% page state control -/pgsv () def -/bpage {/pgsv save def} bdef -/epage {pgsv restore} bdef -/bplot /gsave ldef -/eplot {stroke grestore} bdef -% orientation switch -/portraitMode 0 def /landscapeMode 1 def /rotateMode 2 def -% coordinate system mappings -/dpi2point 0 def -% font control -/FontSize 0 def -/FMS {/FontSize xstore findfont [FontSize 0 0 FontSize neg 0 0] - makefont setfont} bdef -/reencode {exch dup where {pop load} {pop StandardEncoding} ifelse - exch dup 3 1 roll findfont dup length dict begin - { 1 index /FID ne {def}{pop pop} ifelse } forall - /Encoding exch def currentdict end definefont pop} bdef -/isroman {findfont /CharStrings get /Agrave known} bdef -/FMSR {3 1 roll 1 index dup isroman {reencode} {pop pop} ifelse - exch FMS} bdef -/csm {1 dpi2point div -1 dpi2point div scale neg translate - dup landscapeMode eq {pop -90 rotate} - {rotateMode eq {90 rotate} if} ifelse} bdef -% line types: solid, dotted, dashed, dotdash -/SO { [] 0 setdash } bdef -/DO { [.5 dpi2point mul 4 dpi2point mul] 0 setdash } bdef -/DA { [6 dpi2point mul] 0 setdash } bdef -/DD { [.5 dpi2point mul 4 dpi2point mul 6 dpi2point mul 4 - dpi2point mul] 0 setdash } bdef -% macros for lines and objects -/L {lineto stroke} bdef -/MP {3 1 roll moveto 1 sub {rlineto} repeat} bdef -/AP {{rlineto} repeat} bdef -/PDlw -1 def -/W {/PDlw currentlinewidth def setlinewidth} def -/PP {closepath eofill} bdef -/DP {closepath stroke} bdef -/MR {4 -2 roll moveto dup 0 exch rlineto exch 0 rlineto - neg 0 exch rlineto closepath} bdef -/FR {MR stroke} bdef -/PR {MR fill} bdef -/L1i {{currentfile picstr readhexstring pop} image} bdef -/tMatrix matrix def -/MakeOval {newpath tMatrix currentmatrix pop translate scale -0 0 1 0 360 arc tMatrix setmatrix} bdef -/FO {MakeOval stroke} bdef -/PO {MakeOval fill} bdef -/PD {currentlinewidth 2 div 0 360 arc fill - PDlw -1 eq not {PDlw w /PDlw -1 def} if} def -/FA {newpath tMatrix currentmatrix pop translate scale - 0 0 1 5 -2 roll arc tMatrix setmatrix stroke} bdef -/PA {newpath tMatrix currentmatrix pop translate 0 0 moveto scale - 0 0 1 5 -2 roll arc closepath tMatrix setmatrix fill} bdef -/FAn {newpath tMatrix currentmatrix pop translate scale - 0 0 1 5 -2 roll arcn tMatrix setmatrix stroke} bdef -/PAn {newpath tMatrix currentmatrix pop translate 0 0 moveto scale - 0 0 1 5 -2 roll arcn closepath tMatrix setmatrix fill} bdef -/vradius 0 def /hradius 0 def /lry 0 def -/lrx 0 def /uly 0 def /ulx 0 def /rad 0 def -/MRR {/vradius xdef /hradius xdef /lry xdef /lrx xdef /uly xdef - /ulx xdef newpath tMatrix currentmatrix pop ulx hradius add uly - vradius add translate hradius vradius scale 0 0 1 180 270 arc - tMatrix setmatrix lrx hradius sub uly vradius add translate - hradius vradius scale 0 0 1 270 360 arc tMatrix setmatrix - lrx hradius sub lry vradius sub translate hradius vradius scale - 0 0 1 0 90 arc tMatrix setmatrix ulx hradius add lry vradius sub - translate hradius vradius scale 0 0 1 90 180 arc tMatrix setmatrix - closepath} bdef -/FRR {MRR stroke } bdef -/PRR {MRR fill } bdef -/MlrRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lry uly sub 2 div def - newpath tMatrix currentmatrix pop ulx rad add uly rad add translate - rad rad scale 0 0 1 90 270 arc tMatrix setmatrix lrx rad sub lry rad - sub translate rad rad scale 0 0 1 270 90 arc tMatrix setmatrix - closepath} bdef -/FlrRR {MlrRR stroke } bdef -/PlrRR {MlrRR fill } bdef -/MtbRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lrx ulx sub 2 div def - newpath tMatrix currentmatrix pop ulx rad add uly rad add translate - rad rad scale 0 0 1 180 360 arc tMatrix setmatrix lrx rad sub lry rad - sub translate rad rad scale 0 0 1 0 180 arc tMatrix setmatrix - closepath} bdef -/FtbRR {MtbRR stroke } bdef -/PtbRR {MtbRR fill } bdef -/stri 6 array def /dtri 6 array def -/smat 6 array def /dmat 6 array def -/tmat1 6 array def /tmat2 6 array def /dif 3 array def -/asub {/ind2 exch def /ind1 exch def dup dup - ind1 get exch ind2 get sub exch } bdef -/tri_to_matrix { - 2 0 asub 3 1 asub 4 0 asub 5 1 asub - dup 0 get exch 1 get 7 -1 roll astore } bdef -/compute_transform { - dmat dtri tri_to_matrix tmat1 invertmatrix - smat stri tri_to_matrix tmat2 concatmatrix } bdef -/ds {stri astore pop} bdef -/dt {dtri astore pop} bdef -/db {2 copy /cols xdef /rows xdef mul dup 3 mul string - currentfile exch readhexstring pop - dup 0 3 index getinterval /rbmap xdef - dup 2 index dup getinterval /gbmap xdef - 1 index dup 2 mul exch getinterval /bbmap xdef pop pop}bdef -/it {gs np dtri aload pop moveto lineto lineto cp c - cols rows 8 compute_transform - rbmap gbmap bbmap true 3 colorimage gr}bdef -/il {newpath moveto lineto stroke}bdef -currentdict end def -%%EndProlog - -%%BeginSetup -MathWorks begin - -0 cap - -end -%%EndSetup - -%%Page: 1 1 -%%BeginPageSetup -%%PageBoundingBox: 53 195 550 602 -MathWorks begin -bpage -%%EndPageSetup - -%%BeginObject: obj1 -bplot - -/dpi2point 8.33333 def -portraitMode 0150 5100 csm - - 294 78 4147 3391 MR c np -76 dict begin %Colortable dictionary -/c0 { 0.000000 0.000000 0.000000 sr} bdef -/c1 { 1.000000 1.000000 1.000000 sr} bdef -/c2 { 0.900000 0.000000 0.000000 sr} bdef -/c3 { 0.000000 0.820000 0.000000 sr} bdef -/c4 { 0.000000 0.000000 0.800000 sr} bdef -/c5 { 0.910000 0.820000 0.320000 sr} bdef -/c6 { 1.000000 0.260000 0.820000 sr} bdef -/c7 { 0.000000 0.820000 0.820000 sr} bdef -c0 -1 j -1 sg - 0 0 4803 3605 PR -4.16667 w -0 2938 3722 0 0 -2938 624 3208 4 MP -PP --3722 0 0 2938 3722 0 0 -2938 624 3208 5 MP stroke -2.77778 w -DO -SO -4.16667 w -0 sg - 624 3208 mt 4346 3208 L - 624 270 mt 4346 270 L - 624 3208 mt 624 270 L -4346 3208 mt 4346 270 L - 624 3208 mt 4346 3208 L - 624 3208 mt 624 270 L - 624 3208 mt 624 3189 L - 624 270 mt 624 288 L - 624 3208 mt 624 3170 L - 624 270 mt 624 307 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 563 3334 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - - 655 3283 mt -(0) s - 904 3208 mt 904 3189 L - 904 270 mt 904 288 L -1067 3208 mt 1067 3189 L -1067 270 mt 1067 288 L -1184 3208 mt 1184 3189 L -1184 270 mt 1184 288 L -1274 3208 mt 1274 3189 L -1274 270 mt 1274 288 L -1348 3208 mt 1348 3189 L -1348 270 mt 1348 288 L -1410 3208 mt 1410 3189 L -1410 270 mt 1410 288 L -1464 3208 mt 1464 3189 L -1464 270 mt 1464 288 L -1511 3208 mt 1511 3189 L -1511 270 mt 1511 288 L -1554 3208 mt 1554 3189 L -1554 270 mt 1554 288 L -1554 3208 mt 1554 3170 L -1554 270 mt 1554 307 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - -1493 3334 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - -1585 3283 mt -(1) s -1834 3208 mt 1834 3189 L -1834 270 mt 1834 288 L -1998 3208 mt 1998 3189 L -1998 270 mt 1998 288 L -2114 3208 mt 2114 3189 L -2114 270 mt 2114 288 L -2204 3208 mt 2204 3189 L -2204 270 mt 2204 288 L -2278 3208 mt 2278 3189 L -2278 270 mt 2278 288 L -2340 3208 mt 2340 3189 L -2340 270 mt 2340 288 L -2394 3208 mt 2394 3189 L -2394 270 mt 2394 288 L -2442 3208 mt 2442 3189 L -2442 270 mt 2442 288 L -2485 3208 mt 2485 3189 L -2485 270 mt 2485 288 L -2485 3208 mt 2485 3170 L -2485 270 mt 2485 307 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - -2424 3334 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - -2516 3283 mt -(2) s -2765 3208 mt 2765 3189 L -2765 270 mt 2765 288 L -2928 3208 mt 2928 3189 L -2928 270 mt 2928 288 L -3045 3208 mt 3045 3189 L -3045 270 mt 3045 288 L -3135 3208 mt 3135 3189 L -3135 270 mt 3135 288 L -3209 3208 mt 3209 3189 L -3209 270 mt 3209 288 L -3271 3208 mt 3271 3189 L -3271 270 mt 3271 288 L -3325 3208 mt 3325 3189 L -3325 270 mt 3325 288 L -3372 3208 mt 3372 3189 L -3372 270 mt 3372 288 L -3415 3208 mt 3415 3189 L -3415 270 mt 3415 288 L -3415 3208 mt 3415 3170 L -3415 270 mt 3415 307 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - -3354 3334 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - -3446 3283 mt -(3) s -3695 3208 mt 3695 3189 L -3695 270 mt 3695 288 L -3859 3208 mt 3859 3189 L -3859 270 mt 3859 288 L -3975 3208 mt 3975 3189 L -3975 270 mt 3975 288 L -4065 3208 mt 4065 3189 L -4065 270 mt 4065 288 L -4139 3208 mt 4139 3189 L -4139 270 mt 4139 288 L -4201 3208 mt 4201 3189 L -4201 270 mt 4201 288 L -4255 3208 mt 4255 3189 L -4255 270 mt 4255 288 L -4303 3208 mt 4303 3189 L -4303 270 mt 4303 288 L -4346 3208 mt 4346 3189 L -4346 270 mt 4346 288 L -4346 3208 mt 4346 3170 L -4346 270 mt 4346 307 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - -4285 3334 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - -4377 3283 mt -(4) s - 624 3208 mt 642 3208 L -4346 3208 mt 4327 3208 L - 624 3208 mt 661 3208 L -4346 3208 mt 4308 3208 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 445 3238 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - - 537 3187 mt -(-7) s - 624 3097 mt 642 3097 L -4346 3097 mt 4327 3097 L - 624 3032 mt 642 3032 L -4346 3032 mt 4327 3032 L - 624 2986 mt 642 2986 L -4346 2986 mt 4327 2986 L - 624 2951 mt 642 2951 L -4346 2951 mt 4327 2951 L - 624 2922 mt 642 2922 L -4346 2922 mt 4327 2922 L - 624 2897 mt 642 2897 L -4346 2897 mt 4327 2897 L - 624 2876 mt 642 2876 L -4346 2876 mt 4327 2876 L - 624 2857 mt 642 2857 L -4346 2857 mt 4327 2857 L - 624 2840 mt 642 2840 L -4346 2840 mt 4327 2840 L - 624 2840 mt 661 2840 L -4346 2840 mt 4308 2840 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 445 2870 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - - 537 2819 mt -(-6) s - 624 2730 mt 642 2730 L -4346 2730 mt 4327 2730 L - 624 2665 mt 642 2665 L -4346 2665 mt 4327 2665 L - 624 2619 mt 642 2619 L -4346 2619 mt 4327 2619 L - 624 2584 mt 642 2584 L -4346 2584 mt 4327 2584 L - 624 2554 mt 642 2554 L -4346 2554 mt 4327 2554 L - 624 2530 mt 642 2530 L -4346 2530 mt 4327 2530 L - 624 2509 mt 642 2509 L -4346 2509 mt 4327 2509 L - 624 2490 mt 642 2490 L -4346 2490 mt 4327 2490 L - 624 2473 mt 642 2473 L -4346 2473 mt 4327 2473 L - 624 2473 mt 661 2473 L -4346 2473 mt 4308 2473 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 445 2503 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - - 537 2452 mt -(-5) s - 624 2362 mt 642 2362 L -4346 2362 mt 4327 2362 L - 624 2298 mt 642 2298 L -4346 2298 mt 4327 2298 L - 624 2252 mt 642 2252 L -4346 2252 mt 4327 2252 L - 624 2216 mt 642 2216 L -4346 2216 mt 4327 2216 L - 624 2187 mt 642 2187 L -4346 2187 mt 4327 2187 L - 624 2163 mt 642 2163 L -4346 2163 mt 4327 2163 L - 624 2141 mt 642 2141 L -4346 2141 mt 4327 2141 L - 624 2123 mt 642 2123 L -4346 2123 mt 4327 2123 L - 624 2106 mt 642 2106 L -4346 2106 mt 4327 2106 L - 624 2106 mt 661 2106 L -4346 2106 mt 4308 2106 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 445 2136 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - - 537 2085 mt -(-4) s - 624 1995 mt 642 1995 L -4346 1995 mt 4327 1995 L - 624 1931 mt 642 1931 L -4346 1931 mt 4327 1931 L - 624 1885 mt 642 1885 L -4346 1885 mt 4327 1885 L - 624 1849 mt 642 1849 L -4346 1849 mt 4327 1849 L - 624 1820 mt 642 1820 L -4346 1820 mt 4327 1820 L - 624 1795 mt 642 1795 L -4346 1795 mt 4327 1795 L - 624 1774 mt 642 1774 L -4346 1774 mt 4327 1774 L - 624 1755 mt 642 1755 L -4346 1755 mt 4327 1755 L - 624 1739 mt 642 1739 L -4346 1739 mt 4327 1739 L - 624 1739 mt 661 1739 L -4346 1739 mt 4308 1739 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 445 1769 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - - 537 1718 mt -(-3) s - 624 1628 mt 642 1628 L -4346 1628 mt 4327 1628 L - 624 1563 mt 642 1563 L -4346 1563 mt 4327 1563 L - 624 1517 mt 642 1517 L -4346 1517 mt 4327 1517 L - 624 1482 mt 642 1482 L -4346 1482 mt 4327 1482 L - 624 1453 mt 642 1453 L -4346 1453 mt 4327 1453 L - 624 1428 mt 642 1428 L -4346 1428 mt 4327 1428 L - 624 1407 mt 642 1407 L -4346 1407 mt 4327 1407 L - 624 1388 mt 642 1388 L -4346 1388 mt 4327 1388 L - 624 1371 mt 642 1371 L -4346 1371 mt 4327 1371 L - 624 1371 mt 661 1371 L -4346 1371 mt 4308 1371 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 445 1401 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - - 537 1350 mt -(-2) s - 624 1261 mt 642 1261 L -4346 1261 mt 4327 1261 L - 624 1196 mt 642 1196 L -4346 1196 mt 4327 1196 L - 624 1150 mt 642 1150 L -4346 1150 mt 4327 1150 L - 624 1115 mt 642 1115 L -4346 1115 mt 4327 1115 L - 624 1085 mt 642 1085 L -4346 1085 mt 4327 1085 L - 624 1061 mt 642 1061 L -4346 1061 mt 4327 1061 L - 624 1040 mt 642 1040 L -4346 1040 mt 4327 1040 L - 624 1021 mt 642 1021 L -4346 1021 mt 4327 1021 L - 624 1004 mt 642 1004 L -4346 1004 mt 4327 1004 L - 624 1004 mt 661 1004 L -4346 1004 mt 4308 1004 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 445 1034 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - - 537 983 mt -(-1) s - 624 893 mt 642 893 L -4346 893 mt 4327 893 L - 624 829 mt 642 829 L -4346 829 mt 4327 829 L - 624 783 mt 642 783 L -4346 783 mt 4327 783 L - 624 747 mt 642 747 L -4346 747 mt 4327 747 L - 624 718 mt 642 718 L -4346 718 mt 4327 718 L - 624 694 mt 642 694 L -4346 694 mt 4327 694 L - 624 672 mt 642 672 L -4346 672 mt 4327 672 L - 624 654 mt 642 654 L -4346 654 mt 4327 654 L - 624 637 mt 642 637 L -4346 637 mt 4327 637 L - 624 637 mt 661 637 L -4346 637 mt 4308 637 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 445 667 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - - 537 616 mt -(0) s - 624 526 mt 642 526 L -4346 526 mt 4327 526 L - 624 462 mt 642 462 L -4346 462 mt 4327 462 L - 624 416 mt 642 416 L -4346 416 mt 4327 416 L - 624 380 mt 642 380 L -4346 380 mt 4327 380 L - 624 351 mt 642 351 L -4346 351 mt 4327 351 L - 624 326 mt 642 326 L -4346 326 mt 4327 326 L - 624 305 mt 642 305 L -4346 305 mt 4327 305 L - 624 286 mt 642 286 L -4346 286 mt 4327 286 L - 624 270 mt 642 270 L -4346 270 mt 4327 270 L - 624 270 mt 661 270 L -4346 270 mt 4308 270 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 445 300 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - - 537 249 mt -(1) s - 624 3208 mt 4346 3208 L - 624 270 mt 4346 270 L - 624 3208 mt 624 270 L -4346 3208 mt 4346 270 L -gs 624 269 3723 2939 MR c np -/c8 { 0.000000 0.000000 1.000000 sr} bdef -c8 -68 111 100 0 65 110 95 111 107 110 100 111 97 110 93 111 -87 111 72 0 55 110 84 111 55 110 64 111 60 0 60 110 -62 111 45 110 62 0 48 111 51 0 71 110 51 111 90 0 -71 111 104 0 122 110 190 111 254 110 560 0 1184 637 31 MP stroke -gs 1133 586 3146 2535 MR c np - 25 25 1184 637 FO - 25 25 1744 637 FO - 25 25 1998 747 FO - 25 25 2188 858 FO - 25 25 2310 968 FO - 25 25 2414 968 FO - 25 25 2485 1079 FO - 25 25 2575 1079 FO - 25 25 2626 1190 FO - 25 25 2697 1300 FO - 25 25 2748 1300 FO - 25 25 2796 1411 FO - 25 25 2858 1411 FO - 25 25 2903 1521 FO - 25 25 2965 1632 FO - 25 25 3025 1742 FO - 25 25 3085 1742 FO - 25 25 3149 1853 FO - 25 25 3204 1963 FO - 25 25 3288 2074 FO - 25 25 3343 2184 FO - 25 25 3415 2184 FO - 25 25 3502 2295 FO - 25 25 3595 2406 FO - 25 25 3692 2516 FO - 25 25 3792 2627 FO - 25 25 3899 2737 FO - 25 25 3994 2848 FO - 25 25 4059 2958 FO - 25 25 4159 2958 FO - 25 25 4227 3069 FO -gr - -/c9 { 0.000000 0.500000 0.000000 sr} bdef -c9 -68 0 100 0 65 111 95 110 107 111 100 0 97 110 93 0 -87 111 72 110 55 0 84 111 55 0 64 0 60 110 60 0 -62 111 45 0 62 110 48 0 51 0 71 111 51 0 90 111 -71 0 104 0 122 110 190 111 254 110 560 0 1184 637 31 MP stroke -gs 1133 586 3146 1761 MR c np - 25 25 1184 637 FO - 25 25 1744 637 FO - 25 25 1998 747 FO - 25 25 2188 858 FO - 25 25 2310 968 FO - 25 25 2414 968 FO - 25 25 2485 968 FO - 25 25 2575 1079 FO - 25 25 2626 1079 FO - 25 25 2697 1190 FO - 25 25 2748 1190 FO - 25 25 2796 1190 FO - 25 25 2858 1300 FO - 25 25 2903 1300 FO - 25 25 2965 1411 FO - 25 25 3025 1411 FO - 25 25 3085 1521 FO - 25 25 3149 1521 FO - 25 25 3204 1521 FO - 25 25 3288 1632 FO - 25 25 3343 1632 FO - 25 25 3415 1742 FO - 25 25 3502 1853 FO - 25 25 3595 1853 FO - 25 25 3692 1963 FO - 25 25 3792 1963 FO - 25 25 3899 2074 FO - 25 25 3994 2184 FO - 25 25 4059 2295 FO - 25 25 4159 2295 FO - 25 25 4227 2295 FO -gr - -/c10 { 1.000000 0.000000 0.000000 sr} bdef -c10 -68 110 100 0 65 111 95 0 107 110 100 111 97 110 93 111 -87 111 72 0 55 110 84 111 55 110 64 111 60 0 60 110 -62 111 45 110 62 0 48 111 51 0 71 110 51 111 90 0 -71 111 104 110 122 111 190 0 254 110 560 0 1184 637 31 MP stroke -gs 1133 586 3146 2424 MR c np - 25 25 1184 637 FO - 25 25 1744 637 FO - 25 25 1998 747 FO - 25 25 2188 747 FO - 25 25 2310 858 FO - 25 25 2414 968 FO - 25 25 2485 1079 FO - 25 25 2575 1079 FO - 25 25 2626 1190 FO - 25 25 2697 1300 FO - 25 25 2748 1300 FO - 25 25 2796 1411 FO - 25 25 2858 1411 FO - 25 25 2903 1521 FO - 25 25 2965 1632 FO - 25 25 3025 1742 FO - 25 25 3085 1742 FO - 25 25 3149 1853 FO - 25 25 3204 1963 FO - 25 25 3288 2074 FO - 25 25 3343 2184 FO - 25 25 3415 2184 FO - 25 25 3502 2295 FO - 25 25 3595 2406 FO - 25 25 3692 2516 FO - 25 25 3792 2627 FO - 25 25 3899 2737 FO - 25 25 3994 2737 FO - 25 25 4059 2848 FO - 25 25 4159 2848 FO - 25 25 4227 2958 FO -gr - -DD -c8 -68 14 100 20 65 13 95 18 107 21 100 20 97 19 93 19 -87 17 72 14 55 11 84 16 55 11 64 13 60 12 60 12 -62 12 45 9 62 12 48 9 51 11 71 13 51 11 90 17 -71 14 104 21 122 24 190 37 254 51 560 110 1184 437 31 MP stroke -c9 -68 7 100 10 65 6 95 9 107 11 100 10 97 9 93 9 -87 9 72 7 55 6 84 8 55 5 64 7 60 6 60 5 -62 7 45 4 62 6 48 5 51 5 71 7 51 5 90 9 -71 7 104 10 122 12 190 19 254 25 560 55 1184 517 31 MP stroke -c10 -68 20 100 30 65 19 95 28 107 32 100 29 97 29 93 28 -87 26 72 21 55 16 84 25 55 16 64 19 60 18 60 18 -62 18 45 13 62 19 48 14 51 15 71 21 51 15 90 27 -71 21 104 30 122 37 190 56 254 75 560 166 1184 692 31 MP stroke -gr - -c10 -DD -0 sg -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - -2280 205 mt -(hmin hmax) s -2308 3434 mt -(Elemente) s - 388 1927 mt -90 rotate -(Schaetzer) s -90 rotate - 612 3238 mt -( ) s -4335 299 mt -( ) s -SO -1 sg -0 714 1792 0 0 -714 665 3167 4 MP -PP --1792 0 0 714 1792 0 0 -714 665 3167 5 MP stroke -2.77778 w -DO -SO -4.16667 w -0 sg - 665 3167 mt 2457 3167 L - 665 3167 mt 665 2453 L - 989 2548 mt -(min hmin/max hmax testAA Analytisch) s -gs 665 2453 1793 715 MR c np -c8 -249 0 714 2517 2 MP stroke -gs 788 2466 103 103 MR c np - 25 25 839 2517 FO -gr - -gr - -c8 -0 sg - 989 2651 mt -(min hmax/max hmax testAA Analytisch) s -gs 665 2453 1793 715 MR c np -c9 -249 0 714 2621 2 MP stroke -gs 788 2570 103 103 MR c np - 25 25 839 2621 FO -gr - -gr - -c9 -0 sg - 989 2754 mt -(hmin/hmax testAA Analytisch) s -gs 665 2453 1793 715 MR c np -c10 -249 0 714 2724 2 MP stroke -gs 788 2673 103 103 MR c np - 25 25 839 2724 FO -gr - -gr - -c10 -0 sg - 989 2879 mt -(N) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 66.6667 FMSR - -1049 2838 mt -(-1/2) s -gs 665 2453 1793 715 MR c np -DD -c8 -249 0 714 2838 2 MP stroke -SO -gr - -c8 -0 sg -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 989 3005 mt -(N) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 66.6667 FMSR - -1049 2964 mt -(-1/4) s -gs 665 2453 1793 715 MR c np -DD -c9 -249 0 714 2964 2 MP stroke -SO -gr - -c9 -0 sg -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 989 3132 mt -(N) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 66.6667 FMSR - -1049 3091 mt -(-3/4) s -gs 665 2453 1793 715 MR c np -DD -c10 -249 0 714 3090 2 MP stroke -SO -gr - -c10 - -end %%Color Dict - -eplot -%%EndObject - -epage -end - -showpage - -%%Trailer -%%EOF diff --git a/src/expmlAA_2DQuad_norm.eps b/src/expmlAA_2DQuad_norm.eps deleted file mode 100644 index d50751f..0000000 --- a/src/expmlAA_2DQuad_norm.eps +++ /dev/null @@ -1,552 +0,0 @@ -%!PS-Adobe-2.0 EPSF-1.2 -%%Creator: MATLAB, The MathWorks, Inc. Version 7.12.0.635 (R2011a). Operating System: Linux 3.0.0-16-server #28-Ubuntu SMP Fri Jan 27 18:03:45 UTC 2012 x86_64. -%%Title: ./expmlAA_2DQuad_norm.eps -%%CreationDate: 03/26/2012 23:43:44 -%%DocumentNeededFonts: Helvetica -%%DocumentProcessColors: Cyan Magenta Yellow Black -%%Extensions: CMYK -%%Pages: 1 -%%BoundingBox: 44 195 550 602 -%%EndComments - -%%BeginProlog -% MathWorks dictionary -/MathWorks 160 dict begin -% definition operators -/bdef {bind def} bind def -/ldef {load def} bind def -/xdef {exch def} bdef -/xstore {exch store} bdef -% operator abbreviations -/c /clip ldef -/cc /concat ldef -/cp /closepath ldef -/gr /grestore ldef -/gs /gsave ldef -/mt /moveto ldef -/np /newpath ldef -/cm /currentmatrix ldef -/sm /setmatrix ldef -/rm /rmoveto ldef -/rl /rlineto ldef -/s {show newpath} bdef -/sc {setcmykcolor} bdef -/sr /setrgbcolor ldef -/sg /setgray ldef -/w /setlinewidth ldef -/j /setlinejoin ldef -/cap /setlinecap ldef -/rc {rectclip} bdef -/rf {rectfill} bdef -% page state control -/pgsv () def -/bpage {/pgsv save def} bdef -/epage {pgsv restore} bdef -/bplot /gsave ldef -/eplot {stroke grestore} bdef -% orientation switch -/portraitMode 0 def /landscapeMode 1 def /rotateMode 2 def -% coordinate system mappings -/dpi2point 0 def -% font control -/FontSize 0 def -/FMS {/FontSize xstore findfont [FontSize 0 0 FontSize neg 0 0] - makefont setfont} bdef -/reencode {exch dup where {pop load} {pop StandardEncoding} ifelse - exch dup 3 1 roll findfont dup length dict begin - { 1 index /FID ne {def}{pop pop} ifelse } forall - /Encoding exch def currentdict end definefont pop} bdef -/isroman {findfont /CharStrings get /Agrave known} bdef -/FMSR {3 1 roll 1 index dup isroman {reencode} {pop pop} ifelse - exch FMS} bdef -/csm {1 dpi2point div -1 dpi2point div scale neg translate - dup landscapeMode eq {pop -90 rotate} - {rotateMode eq {90 rotate} if} ifelse} bdef -% line types: solid, dotted, dashed, dotdash -/SO { [] 0 setdash } bdef -/DO { [.5 dpi2point mul 4 dpi2point mul] 0 setdash } bdef -/DA { [6 dpi2point mul] 0 setdash } bdef -/DD { [.5 dpi2point mul 4 dpi2point mul 6 dpi2point mul 4 - dpi2point mul] 0 setdash } bdef -% macros for lines and objects -/L {lineto stroke} bdef -/MP {3 1 roll moveto 1 sub {rlineto} repeat} bdef -/AP {{rlineto} repeat} bdef -/PDlw -1 def -/W {/PDlw currentlinewidth def setlinewidth} def -/PP {closepath eofill} bdef -/DP {closepath stroke} bdef -/MR {4 -2 roll moveto dup 0 exch rlineto exch 0 rlineto - neg 0 exch rlineto closepath} bdef -/FR {MR stroke} bdef -/PR {MR fill} bdef -/L1i {{currentfile picstr readhexstring pop} image} bdef -/tMatrix matrix def -/MakeOval {newpath tMatrix currentmatrix pop translate scale -0 0 1 0 360 arc tMatrix setmatrix} bdef -/FO {MakeOval stroke} bdef -/PO {MakeOval fill} bdef -/PD {currentlinewidth 2 div 0 360 arc fill - PDlw -1 eq not {PDlw w /PDlw -1 def} if} def -/FA {newpath tMatrix currentmatrix pop translate scale - 0 0 1 5 -2 roll arc tMatrix setmatrix stroke} bdef -/PA {newpath tMatrix currentmatrix pop translate 0 0 moveto scale - 0 0 1 5 -2 roll arc closepath tMatrix setmatrix fill} bdef -/FAn {newpath tMatrix currentmatrix pop translate scale - 0 0 1 5 -2 roll arcn tMatrix setmatrix stroke} bdef -/PAn {newpath tMatrix currentmatrix pop translate 0 0 moveto scale - 0 0 1 5 -2 roll arcn closepath tMatrix setmatrix fill} bdef -/vradius 0 def /hradius 0 def /lry 0 def -/lrx 0 def /uly 0 def /ulx 0 def /rad 0 def -/MRR {/vradius xdef /hradius xdef /lry xdef /lrx xdef /uly xdef - /ulx xdef newpath tMatrix currentmatrix pop ulx hradius add uly - vradius add translate hradius vradius scale 0 0 1 180 270 arc - tMatrix setmatrix lrx hradius sub uly vradius add translate - hradius vradius scale 0 0 1 270 360 arc tMatrix setmatrix - lrx hradius sub lry vradius sub translate hradius vradius scale - 0 0 1 0 90 arc tMatrix setmatrix ulx hradius add lry vradius sub - translate hradius vradius scale 0 0 1 90 180 arc tMatrix setmatrix - closepath} bdef -/FRR {MRR stroke } bdef -/PRR {MRR fill } bdef -/MlrRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lry uly sub 2 div def - newpath tMatrix currentmatrix pop ulx rad add uly rad add translate - rad rad scale 0 0 1 90 270 arc tMatrix setmatrix lrx rad sub lry rad - sub translate rad rad scale 0 0 1 270 90 arc tMatrix setmatrix - closepath} bdef -/FlrRR {MlrRR stroke } bdef -/PlrRR {MlrRR fill } bdef -/MtbRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lrx ulx sub 2 div def - newpath tMatrix currentmatrix pop ulx rad add uly rad add translate - rad rad scale 0 0 1 180 360 arc tMatrix setmatrix lrx rad sub lry rad - sub translate rad rad scale 0 0 1 0 180 arc tMatrix setmatrix - closepath} bdef -/FtbRR {MtbRR stroke } bdef -/PtbRR {MtbRR fill } bdef -/stri 6 array def /dtri 6 array def -/smat 6 array def /dmat 6 array def -/tmat1 6 array def /tmat2 6 array def /dif 3 array def -/asub {/ind2 exch def /ind1 exch def dup dup - ind1 get exch ind2 get sub exch } bdef -/tri_to_matrix { - 2 0 asub 3 1 asub 4 0 asub 5 1 asub - dup 0 get exch 1 get 7 -1 roll astore } bdef -/compute_transform { - dmat dtri tri_to_matrix tmat1 invertmatrix - smat stri tri_to_matrix tmat2 concatmatrix } bdef -/ds {stri astore pop} bdef -/dt {dtri astore pop} bdef -/db {2 copy /cols xdef /rows xdef mul dup 3 mul string - currentfile exch readhexstring pop - dup 0 3 index getinterval /rbmap xdef - dup 2 index dup getinterval /gbmap xdef - 1 index dup 2 mul exch getinterval /bbmap xdef pop pop}bdef -/it {gs np dtri aload pop moveto lineto lineto cp c - cols rows 8 compute_transform - rbmap gbmap bbmap true 3 colorimage gr}bdef -/il {newpath moveto lineto stroke}bdef -currentdict end def -%%EndProlog - -%%BeginSetup -MathWorks begin - -0 cap - -end -%%EndSetup - -%%Page: 1 1 -%%BeginPageSetup -%%PageBoundingBox: 44 195 550 602 -MathWorks begin -bpage -%%EndPageSetup - -%%BeginObject: obj1 -bplot - -/dpi2point 8.33333 def -portraitMode 0150 5100 csm - - 223 79 4218 3390 MR c np -76 dict begin %Colortable dictionary -/c0 { 0.000000 0.000000 0.000000 sr} bdef -/c1 { 1.000000 1.000000 1.000000 sr} bdef -/c2 { 0.900000 0.000000 0.000000 sr} bdef -/c3 { 0.000000 0.820000 0.000000 sr} bdef -/c4 { 0.000000 0.000000 0.800000 sr} bdef -/c5 { 0.910000 0.820000 0.320000 sr} bdef -/c6 { 1.000000 0.260000 0.820000 sr} bdef -/c7 { 0.000000 0.820000 0.820000 sr} bdef -c0 -1 j -1 sg - 0 0 4803 3605 PR -4.16667 w -0 2938 3722 0 0 -2938 624 3208 4 MP -PP --3722 0 0 2938 3722 0 0 -2938 624 3208 5 MP stroke -2.77778 w -DO -SO -4.16667 w -0 sg - 624 3208 mt 4346 3208 L - 624 270 mt 4346 270 L - 624 3208 mt 624 270 L -4346 3208 mt 4346 270 L - 624 3208 mt 4346 3208 L - 624 3208 mt 624 270 L - 624 3208 mt 624 3189 L - 624 270 mt 624 288 L - 624 3208 mt 624 3170 L - 624 270 mt 624 307 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 563 3334 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - - 655 3283 mt -(0) s - 904 3208 mt 904 3189 L - 904 270 mt 904 288 L -1067 3208 mt 1067 3189 L -1067 270 mt 1067 288 L -1184 3208 mt 1184 3189 L -1184 270 mt 1184 288 L -1274 3208 mt 1274 3189 L -1274 270 mt 1274 288 L -1348 3208 mt 1348 3189 L -1348 270 mt 1348 288 L -1410 3208 mt 1410 3189 L -1410 270 mt 1410 288 L -1464 3208 mt 1464 3189 L -1464 270 mt 1464 288 L -1511 3208 mt 1511 3189 L -1511 270 mt 1511 288 L -1554 3208 mt 1554 3189 L -1554 270 mt 1554 288 L -1554 3208 mt 1554 3170 L -1554 270 mt 1554 307 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - -1493 3334 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - -1585 3283 mt -(1) s -1834 3208 mt 1834 3189 L -1834 270 mt 1834 288 L -1998 3208 mt 1998 3189 L -1998 270 mt 1998 288 L -2114 3208 mt 2114 3189 L -2114 270 mt 2114 288 L -2204 3208 mt 2204 3189 L -2204 270 mt 2204 288 L -2278 3208 mt 2278 3189 L -2278 270 mt 2278 288 L -2340 3208 mt 2340 3189 L -2340 270 mt 2340 288 L -2394 3208 mt 2394 3189 L -2394 270 mt 2394 288 L -2442 3208 mt 2442 3189 L -2442 270 mt 2442 288 L -2485 3208 mt 2485 3189 L -2485 270 mt 2485 288 L -2485 3208 mt 2485 3170 L -2485 270 mt 2485 307 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - -2424 3334 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - -2516 3283 mt -(2) s -2765 3208 mt 2765 3189 L -2765 270 mt 2765 288 L -2928 3208 mt 2928 3189 L -2928 270 mt 2928 288 L -3045 3208 mt 3045 3189 L -3045 270 mt 3045 288 L -3135 3208 mt 3135 3189 L -3135 270 mt 3135 288 L -3209 3208 mt 3209 3189 L -3209 270 mt 3209 288 L -3271 3208 mt 3271 3189 L -3271 270 mt 3271 288 L -3325 3208 mt 3325 3189 L -3325 270 mt 3325 288 L -3372 3208 mt 3372 3189 L -3372 270 mt 3372 288 L -3415 3208 mt 3415 3189 L -3415 270 mt 3415 288 L -3415 3208 mt 3415 3170 L -3415 270 mt 3415 307 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - -3354 3334 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - -3446 3283 mt -(3) s -3695 3208 mt 3695 3189 L -3695 270 mt 3695 288 L -3859 3208 mt 3859 3189 L -3859 270 mt 3859 288 L -3975 3208 mt 3975 3189 L -3975 270 mt 3975 288 L -4065 3208 mt 4065 3189 L -4065 270 mt 4065 288 L -4139 3208 mt 4139 3189 L -4139 270 mt 4139 288 L -4201 3208 mt 4201 3189 L -4201 270 mt 4201 288 L -4255 3208 mt 4255 3189 L -4255 270 mt 4255 288 L -4303 3208 mt 4303 3189 L -4303 270 mt 4303 288 L -4346 3208 mt 4346 3189 L -4346 270 mt 4346 288 L -4346 3208 mt 4346 3170 L -4346 270 mt 4346 307 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - -4285 3334 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - -4377 3283 mt -(4) s - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 2328 mt 661 2328 L -4346 2328 mt 4308 2328 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 400 2358 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - - 492 2307 mt -(0.65) s - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 1024 mt 661 1024 L -4346 1024 mt 4308 1024 L -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 400 1054 mt -(10) s -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 55.5556 FMSR - - 492 1003 mt -(0.66) s - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 624 3208 L - 624 3208 mt 4346 3208 L - 624 270 mt 4346 270 L - 624 3208 mt 624 270 L -4346 3208 mt 4346 270 L -gs 624 270 3723 2939 MR c np -/c8 { 0.000000 0.000000 1.000000 sr} bdef -c8 -68 0 100 0 65 0 95 0 107 0 100 0 97 0 93 0 -87 0 72 -1 55 0 84 -1 55 -1 64 -1 60 -3 60 -3 -62 -5 45 -7 62 -10 48 -14 51 -19 71 -28 51 -38 90 -53 -71 -82 104 -164 122 -184 190 -333 254 -480 560 -1229 1184 3208 31 MP stroke -gr - -c8 -gs 1133 501 3146 2759 MR c np - 25 25 1184 3208 FO - 25 25 1744 1979 FO - 25 25 1998 1499 FO - 25 25 2188 1166 FO - 25 25 2310 982 FO - 25 25 2414 818 FO - 25 25 2485 736 FO - 25 25 2575 683 FO - 25 25 2626 645 FO - 25 25 2697 617 FO - 25 25 2748 598 FO - 25 25 2796 584 FO - 25 25 2858 574 FO - 25 25 2903 567 FO - 25 25 2965 562 FO - 25 25 3025 559 FO - 25 25 3085 556 FO - 25 25 3149 555 FO - 25 25 3204 554 FO - 25 25 3288 553 FO - 25 25 3343 553 FO - 25 25 3415 552 FO - 25 25 3502 552 FO - 25 25 3595 552 FO - 25 25 3692 552 FO - 25 25 3792 552 FO - 25 25 3899 552 FO - 25 25 3994 552 FO - 25 25 4059 552 FO - 25 25 4159 552 FO - 25 25 4227 552 FO -gr - -gs 624 270 3723 2939 MR c np -DD -/c9 { 1.000000 0.000000 0.000000 sr} bdef -c9 -68 0 100 0 65 0 95 0 107 0 100 0 97 0 93 0 -87 0 72 0 55 0 84 0 55 0 64 0 60 0 60 0 -62 0 45 0 62 0 48 0 51 0 71 0 51 0 90 0 -71 0 104 0 122 0 190 0 254 0 560 0 1184 547 31 MP stroke -gr - -c9 -DD -0 sg -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - -2226 206 mt -(Energie Norm) s -2308 3434 mt -(Elemente) s - 343 1883 mt -90 rotate -(eNorm) s -90 rotate -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 66.6667 FMSR - - 302 1631 mt -90 rotate -(2) s -90 rotate -%%IncludeResource: font Helvetica -/Helvetica /ISOLatin1Encoding 83.3333 FMSR - - 612 3239 mt -( ) s -4335 300 mt -( ) s -SO -1 sg -0 232 1000 0 0 -232 3305 3168 4 MP -PP --1000 0 0 232 1000 0 0 -232 3305 3168 5 MP stroke -2.77778 w -DO -SO -4.16667 w -0 sg -3305 3168 mt 4305 3168 L -3305 2936 mt 4305 2936 L -3305 3168 mt 3305 2936 L -4305 3168 mt 4305 2936 L -3305 3168 mt 4305 3168 L -3305 3168 mt 3305 2936 L -3305 3168 mt 4305 3168 L -3305 2936 mt 4305 2936 L -3305 3168 mt 3305 2936 L -4305 3168 mt 4305 2936 L -3628 3031 mt -(testAA Analytisch) s -gs 3305 2936 1001 233 MR c np -c8 -248 0 3354 3000 2 MP stroke -gs 3427 2949 103 103 MR c np - 25 25 3478 3000 FO -gr - -gr - -c8 -0 sg -3628 3134 mt -(extrapoliert) s -gs 3305 2936 1001 233 MR c np -DD -c9 -248 0 3354 3103 2 MP stroke -SO -gr - -c9 - -end %%Color Dict - -eplot -%%EndObject - -epage -end - -showpage - -%%Trailer -%%EOF diff --git a/src/mex_build_AU.cpp b/src/mex_build_V.cpp similarity index 88% rename from src/mex_build_AU.cpp rename to src/mex_build_V.cpp index 7bc237e..6f11e60 100644 --- a/src/mex_build_AU.cpp +++ b/src/mex_build_V.cpp @@ -12,7 +12,7 @@ //#include "gauss.hpp" */ -#define M_EPS 10e-8 +//#define M_EPS 10e-8 /* //#include "tbb/parallel_for.h" */ @@ -28,7 +28,9 @@ int dimOfVec(double* vec) { return 1; else if (vec[0] != 0) return 0; - else {/* + else { + mexErrMsgTxt("length of Site is zero"); + /* // cerr << "dimOfVec : (" << vec[0] << " " << vec[1] << " " << vec[2] // << ") alle Werte 0 \n";*/ return -1; @@ -36,6 +38,7 @@ int dimOfVec(double* vec) { } +// Gibt von [0 1 2] die Fehlende Zahl zurueck inline int dimOfThird(int a, int b) { return ((-(a + b) % 3) + 3) % 3; } @@ -52,7 +55,7 @@ void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) { //Sicherheitsabfragen zu Datengroessen if ((nrhs != 4)) mexErrMsgTxt( - "expected (coordinates(Nx3),elements(Mx4),mu(double),type(int))"); + "expected (coordinates(Nx3),elements(Mx4),zeta(double),type(int))"); if (nlhs > 1) mexErrMsgTxt("has only one output argument"); @@ -76,7 +79,7 @@ void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) { double * E = mxGetPr(prhs[1]); int type = (int) *(mxGetPr(prhs[3])); - double * mu = mxGetPr(prhs[2]); + double * zeta = mxGetPr(prhs[2]); //Initialisieren der Hilfsvektoren @@ -231,27 +234,27 @@ void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) { if (rx == ry) { if (rxa == rya) { tmp = ctypeP(fabs(xa[rxa]), fabs(xb[rxb]), fabs(ya[rxa]), - fabs(yb[rxb]), d[rxa], d[rxb], d[rx], mu); + fabs(yb[rxb]), d[rxa], d[rxb], d[rx], zeta); } else { tmp = ctypeP(fabs(xa[rxa]), fabs(xb[rxb]), fabs(yb[rxa]), - fabs(ya[rxb]), d[rxa], d[rxb], d[rx], mu); + fabs(ya[rxb]), d[rxa], d[rxb], d[rx], zeta); } } else { if (rxa == rya) { tmp = ctypeO(fabs(xb[rxb]), fabs(xa[rxa]), fabs(ya[rya]), - fabs(yb[ryb]), d[rxb], d[rxa], d[rx], mu); + fabs(yb[ryb]), d[rxb], d[rxa], d[rx], zeta); } else if (rxa == ryb) { tmp = ctypeO(fabs(xb[rxb]), fabs(xa[rxa]), fabs(yb[ryb]), - fabs(ya[rya]), d[rxb], d[rxa], d[rx], mu); + fabs(ya[rya]), d[rxb], d[rxa], d[rx], zeta); } else if (rxb == rya) { tmp = ctypeO(fabs(xa[rxa]), fabs(xb[rxb]), fabs(ya[rya]), - fabs(yb[ryb]), d[rxa], d[rxb], d[rx], mu); + fabs(yb[ryb]), d[rxa], d[rxb], d[rx], zeta); } else { tmp = ctypeO(fabs(xa[rxa]), fabs(xb[rxb]), fabs(yb[ryb]), - fabs(ya[rya]), d[rxa], d[rxb], d[rx], mu); + fabs(ya[rya]), d[rxa], d[rxb], d[rx], zeta); } } diff --git a/src/test_sol.m b/src/test_sol.m index 76bcb4c..23a1fa7 100644 --- a/src/test_sol.m +++ b/src/test_sol.m @@ -2,7 +2,7 @@ format longG % Matrix MEX Funktion neu Compilieren -mex mex_build_AU.cpp slpRectangle.cpp +mex mex_build_V.cpp slpRectangle.cpp % Test ausführen -- 2.47.3