From 8345056d2ffc66adc5a730a69b9653c80c53a6b2 Mon Sep 17 00:00:00 2001 From: treecity Date: Sun, 28 Aug 2011 15:28:57 +0000 Subject: [PATCH] =?utf8?q?[src]=20eigenartig=20->=20Funktionen=20anders=20?= =?utf8?q?f=C3=BCttern=20geht=20aber=20Variablen=20tauschen=20nicht.=20.?= =?utf8?q?=20.?= MIME-Version: 1.0 Content-Type: text/plain; charset=utf8 Content-Transfer-Encoding: 8bit git-svn-id: https://drops.fb12.tu-berlin.de/svn/bacc/trunk@38 26120e32-c555-405d-b3e1-1f783fb42516 --- src/mex_build_As1.cpp | 29 ++++++--- src/save1.mat | Bin 0 -> 317 bytes src/save10.mat | Bin 0 -> 2823 bytes src/save11.mat | Bin 0 -> 3591 bytes src/save12.mat | Bin 0 -> 4485 bytes src/save13.mat | Bin 0 -> 5555 bytes src/save2.mat | Bin 0 -> 365 bytes src/save3.mat | Bin 0 -> 472 bytes src/save4.mat | Bin 0 -> 583 bytes src/save5.mat | Bin 0 -> 732 bytes src/save6.mat | Bin 0 -> 941 bytes src/save7.mat | Bin 0 -> 1235 bytes src/save8.mat | Bin 0 -> 1799 bytes src/save9.mat | Bin 0 -> 2225 bytes src/test_solveError_Grid.m | 57 ++++++++++++++++ src/test_solveError_Grid2.m | 126 ++++++++++++++++++++++++++++++++++++ 16 files changed, 204 insertions(+), 8 deletions(-) create mode 100644 src/save1.mat create mode 100644 src/save10.mat create mode 100644 src/save11.mat create mode 100644 src/save12.mat create mode 100644 src/save13.mat create mode 100644 src/save2.mat create mode 100644 src/save3.mat create mode 100644 src/save4.mat create mode 100644 src/save5.mat create mode 100644 src/save6.mat create mode 100644 src/save7.mat create mode 100644 src/save8.mat create mode 100644 src/save9.mat create mode 100755 src/test_solveError_Grid.m create mode 100755 src/test_solveError_Grid2.m diff --git a/src/mex_build_As1.cpp b/src/mex_build_As1.cpp index d54e512..485b74f 100644 --- a/src/mex_build_As1.cpp +++ b/src/mex_build_As1.cpp @@ -69,7 +69,7 @@ void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) { double * d = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); double * dt = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); - double tmp; + double tmp,tmp2; int itmp; @@ -206,11 +206,11 @@ void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) { if(xa[rxa]*xa[rxa]+xb[rxb]*xb[rxb]>ya[rya]*ya[rya]+yb[ryb]*yb[ryb]){ - /* printf("(%d,%d@%d,%d)",rxa,rxb,rya,ryb); + printf("+(%d,%d@%d,%d)",rxa,rxb,rya,ryb); printf("[%.2f,%.2f@%.2f,%.2f]",xa[rxa],xb[rxb],ya[rya],yb[ryb]); printf("{%d%d} %.1f %.1f %.1f",j,k,d[0],d[1],d[2]); - printf(" | "); -*/ + printf("\n"); + vtmp = xa; xa = ya; ya = vtmp; vtmp = xb; xb = yb; yb = vtmp; @@ -222,12 +222,12 @@ void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) { for (i = 0; i<3;++i) d[i] = -d[i]; - /* - printf("(%d,%d@%d,%d)",rxa,rxb,rya,ryb); + + printf("-(%d,%d@%d,%d)",rxa,rxb,rya,ryb); printf("[%.2f,%.2f@%.2f,%.2f]",xa[rxa],xb[rxb],ya[rya],yb[ryb]); printf("{%d%d} %.1f %.1f %.1f",j,k,d[0],d[1],d[2]); - printf("\n");*/ + printf("\n"); } @@ -243,10 +243,23 @@ void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) { //printf("%d%d %.1f %.1f | %.1f %.1f || %.1f %.1f | %.1f %.1f || %.1f %.1f %.1f\n",j,k,x0[0], x1[0], x0[1], x3[1], y0[0], // y1[0], y0[1], y3[1], d[0], d[1], d[2]); //printf("%d%d %.1f %.1f %.1f\n",j,k,d[0],d[1],d[2]); + + if(xa[rxa]*xa[rxa]+xb[rxb]*xb[rxb]<=ya[rya]*ya[rya]+yb[ryb]*yb[ryb]){ + tmp = apply0Int(F_par, fabs(xa[rxa]), fabs(xb[rxb]), - fabs(ya[rxa]), fabs(yb[rxb]), d[rxa], + fabs(ya[rya]), fabs(yb[ryb]), d[rxa], d[rxb], d[rx]); + + }else{ + + tmp + = apply0Int(F_par, fabs(ya[rya]), fabs(yb[ryb]), + fabs(xa[rxa]), fabs(xb[rxb]), -d[rxa], + -d[rxb], -d[rx]); + } + // if(fabs(tmp-tmp2)>10e-16) + // printf("wtf"); // printf("%d%d|%.2f\n",j,k,tmp); } else { //printf("%d%d %.1f %.1f | %.1f %.1f || %.1f %.1f | %.1f %.1f || %.1f %.1f %.1f\n",j,k,x0[0], x1[0], x0[1], x3[1], y0[0], diff --git a/src/save1.mat b/src/save1.mat new file mode 100644 index 0000000000000000000000000000000000000000..89b1459cd1e9bc50099033b6894e7a3cc6462e70 GIT binary patch literal 317 zcmeZu4DoSvQZUssQ1EpO(M`+DN!3vZ$Vn_o%P-2cQgHY2i*PhE(NS#k5il@%`tma{FeCzT#hk~<2?;w6o;YxT zZHc9Xf(ge4r=Lv@XLyb{sPf!#m@{X7{FF&Er$dBj*k^1=c3Y7<8GM@Qsu zBtJjQ(DR?U_gJPrU!jHk-B`XGZ=K8eckQ+lD7b!rk)iOT+$E4jK5&aV7?R8so;6%j zIKcI^*??QX*y-j<%auDwYpehO literal 0 HcmV?d00001 diff --git a/src/save10.mat b/src/save10.mat new file mode 100644 index 0000000000000000000000000000000000000000..d24664140dd6985b8d8da5aca63da159a18bc389 GIT binary patch literal 2823 zcmb7^`#%#3AIEu=j$6tlAxy*>x#qISWy$@LYa=xGgxSGnv2Zw{6bXljA=lZY+^-?I z%fj4aM`AOV85zdxai0I+`TqR=zCNGt>-+kAY{AYp;M)LgRds+Z*jXjuQHcL_fPIKJ zA^;u^(*u~>*tvprw5|h8!u`Dw{(b;>s2;!x843U+9|AN%01YiYt($rpAb_U2h6dpO z91c!PTR{#E4&WIMj>q}!Gapm@Ohx+Wjt(N-irYZO*+}Vbb5{qjFltPIpPySpT|Mwi zSC`iXF0Q}K_;u7NU0yvr{5l(5UJ2c6MX$k0D((eBREX1&?Sv0_w<45;CSaYQ3sn&h zl*{n6tex>RwVb#E2%6r?X6}}?e^`!aBjMUe?)$X3Z3ddwdAMC5$o@q8dx7@Tt{Aoz zzC?+fo<_m9nx>dJ_H7k^2yabqiF)*ie&1(Rco?nRWWk^OrvX<`vCXgQEs?yOJr~C> zwV$%CyfOL%SHMeH;z4-yflvMum&DBcF~sV>i_SK;k4p_6&k3_TA!7VFn{IWkbj&{| zjheh$W?s}Ay|(Q6qkm_0UwRU=qE7+k%dn1lq(L&cxw4h17|24v5w8sYaxF}YvJ4>sLCCBOQxD^o*5Zn{Z#l`8J#Q1+cVVz?&u zH?7~Sx}eciYCYq6A$Sw3A9e8=%l~vCiF#uV(Y<$e+4cKIf6*`(vKebIZ14tA8Om;6 z6=hj74N>zE6$ESq;Vje`N!?}Afg#l_wmxDphK*@8!Wj&op5`Goj{LJjvEkuiJDJ{F zck>JF8#1iulUgughuel3X;A6Q)`bS11G)J$p}9W zJ!LO6Hs=1yG-P-f&Ck`f^EJ?zc`=0hP=hxl%e~@``dD#N)G3~nR#Sa6*b<1Ohn$es z_p2CJ=4{MOI+iOt;fs6t z;yWkK@a^{Rb1KVPM3eDYkh|hm-C_Pjtup$b33;g&RS&f(cmm5i&N4di>u%A$8!Eie3G!Uf>v9fxUKXk~T3xnF+V~C{ zR+K0(N8&?!9&>wU5c-chH=8zNs56JQht^flQl|NM2i3>osrMDc$qWa_#WUl*?8j=3S2^d@MvV%bFNgoO?yxn zclB$T8kr2_S>59gF3EQ<^(}3qEwv3kAcMD1wQzr^l^pvm;lq@5@R|}0QLARCb2Q^z z#pXH{XWd$#rm>^)nq7#*5`0uR>#Z z>Lg77dJl#=N&5Sg?{yJ1C#=ut%qMb*hF-vrLnZwU8E33Y6p0x_L+_A&5Iui;{S}Dz zq$8)Oqhcx~>~H5sE6##0-<+ZyTz+K-;X2D7oogc+A5Gqo!L?wVlp7orkl#XAsZq#j zIcnKE=AFWz=Z@+YH9=}7ze-7TWa=0QId*1m5>YU5Ee40uvR?UB5oYl-bsLJ0|vrSZs-Sb)cP6L%GW4nuStpMxxk&B|K&9RYm1U4EL@CPr}km zdGB=4<^IPC+n0tB_dbG$0LXM53mHFt&2xONcVN`ods`rCQ?rKYl4!HH$95m7cwO7J*K0o= z@3z$?8a65=Nkm)DA3R%7KpR_vHU>9K(^-QdaSV0B)(9jNeb2_6h$&0a9!%RP!6*H$ zQUN2nAsftDTBhWNLYj!CN=@Of7q!ZAxh?Ov#`w)Db8=H1q!j`)QuB(_7o{clFqhW2 zM4y2et8VD~MB%og49xu~eYg9f*c#6E+b)JL%A3s;#YfFj{n``}_U3!JEIq7KE1vW# zxlc$_-#@t@D6iX_tuv-Ix~*N7YYOC!b{dU{Z}CX^{kTk~d5Ca4M0>`zke5p%u3m(b z!!?|HH(jq`R>5Bs6JkpgcgiHA}&*Uat2S2gT)Fal#+sYkT zDRpeFmo>I;cvj<7b<}YJ8^K51rmF(JS&OSu0lS8~Ze6g2X)f)nFC@kjh(TKm7I*Y8 z9~W@-qeqMi1`~mO|1G0*r~>3JRMN}8dzfKASf9}M#3I-zN+s9jb{X-W~0~-ve z1-7t&<)J}%m{M@;t4BDP74GRzhybkOJvF$T0mL4TjH&YCsswP4>RPLcFU?dXP<{FI z4j=+VG3p&FfxW(0M zsw#5bL7M#M;-vSHn4{_!;r0PwgT1+wcdII;vOp5A`8rJYXdTlhc13sWNTORI)hF33 zVxYV~+yiSl+&eJ*c<~?G8&cJ}TCXidyvp+C+J9vII&_}PA9m94@tn*&$)?Gr7zx3J zY~zE!|D({et^y90>8X_{Pal zW#VpBB^p>8`;5;|NZ{|;T+qv#%RVN3{bJ>1B$OMA$7MQVZEkH?EA%SAo&AZB*_lZ^ zR&@Z2Mkc42F|xvBhjO9UTpG_!Zh&Fam1)%3#6FDAjonAssYZOkO=P%fYi^w2UgTS5 zt2NgXGeL${(ro&Y1a?d!n0HxwY$16`%^!UGLr7tX9Ux!Y=-NBN(+{4-5#BAIlRPE_7RIF(CJhWkFfO<6 zDq}6n)rXx+=|4Yz@~0Ccngv=0tpj~xZ?jXHMXrrBhY&yTI^3u_`0^fHUCAxb4zu-q zn(n83e>^D?Q6yYo`|&N1QQsXyx$D{yfAMSaOD~ZUap9ZdbD&vMK7EZFO1)Q_d}Lqn zu*L9RWvD2X^_hWy##$*|%)-VB$I^MJrAljw&@LW`Cf;Y}^NZKrkPVuV)G=WZz2!-H VlmH0vq>$N?Iob!3B*Mk3{{@eqyLbQq literal 0 HcmV?d00001 diff --git a/src/save11.mat b/src/save11.mat new file mode 100644 index 0000000000000000000000000000000000000000..093cc3d1081841a1a733e96a8586e80d74b7cb52 GIT binary patch literal 3591 zcmb7^XEYm(!^LT}T3ajj))uv6hZwaw(5kIQBZyKWYOfk?^(mrSv5F|IsE|@4R>i1M zl29X7tXdU&HMalre0o2>_uS9--gE9bKMO+}bHiKW>Wa$Z7KSzoUQeES%8NgE3JLOp z2l|1;&CD$w3^ml`#f<_zAwixV;&6YExOK3;{Q)E zFy67?WME*hyTHH@mfJL1>1$)i)wN}7#nbv+*6j4$z24W{JjNyp^78Qx4i3WZ5D4P+ z%aCf66P+nWXt-4~c_s8&eb6c_*ai=w}k2hd$h2 zCMxz7>r;I(n3G>W(=FGwf2>a*OcLbq{*#VPF+1Dev7yhWL;W|VqY6`&>iTu4d&Bi7 zN1suryhH0O!pwBvP$66!cW&@80QYq>_mnT@)%Y-(!#Tl-pVS_HZ%xT)3H3=aL5vxqeO-LS zyzgVfaQ=$o@JQ^+x>Cx%{Cv(bdsJ%PS?)Bp-m8IK3jBredL2z$yPo_P4^wusUn7-Z zQO>I1SL4>9fl}gs6YISq7>8k>2K7dch(!c~n23w0RpOkpszS@{MB3l3^od>kpCc-y zFHa9c2J`8@qzx_Y&2h_@AD@aetRwY{)<*(KL1MJja0$}+)o8N%;QO3=Ld;cH2}O$? z+G-M44zm3px^?>3zhF+Yzmpp^V&Lb@rDi9WDbG|#%#-@VGVMB3 zw)>(2&rZ&{Y(OdW)PLLJt#qk<$%zCZ&gYsjhZ~_abo_69k)=pnw!gs={i^8-eDcnr zh)WO^1g~<#`3n>0D1DZ=u_xv9rx;NxcEdNX3!J&)hSmGCLkW#m0>yA8k|QY%*cI2D zt&!Aj9K8`u1aa>9xYQuX2Rg31lT+LGC&S%*(ZZ(_x;lF#Z&w8{s(Fun3F=&Uh&y!> zqv-dobxQwcxOf^w3xgh7htM*ta#KDWMPo%F!jsH-aGU84ON`#CuI|Vm%vQ>1V%&2| zDroP!cdXb9Lf z+PIp7XB3Yc_%cW;#oX4nn(T~#&A>6Pv_s>)v6AMb6U%6=fwXW)v%@sUYUq?wn&{n? zzLRqd3}j1x+)mE2`I7gu5+F2{ue001u= z)%XU|UtNFe|33UfZ)UH+58a(_A11z&$wXJ9{RP|(ray6GRM+ zHFZj8|55my}ns_8zrT6Axx%h$B~tONgmt zB`n?8rJ32GE28;DM@*9U35(`B=OI5Y0rQ3eoQ6Uc(ZQ{iwEK%I!n|9g!eBdn`wH4d z?MX=!XVQL8+o6sV?7+FD^tmko=++l78D6~ApR@HMG;76C?D%vue7J8_bBdYzH8FZL zRr=qJXqPCoeo=4+;WBZWegPX}Lr-der*qD{*gr$9Zv{)i%R3DBYgNq))r9@P>enJ=)9!`OUl8Lljk#{%lx`47F|b?pwESdd8Wg8QE0B0oz#p%|?)iWtraL93 z|W`vz)ju2)3iNjJ50+!1+o%o6HkEy! zaSiYT8PQ4h=~(w^ANtcxUhP>oC*>1+nu#Pvd6Gm~r{HlWR10~Tr&Z?0!z&iEYUCMt z|5m?UJGZ?N2KE6i_CY5cjG^m^oK*GuMs35R@VANK>4QPmz{&RVP8_glcv{r+piF(N|{kzs_V zHHY84*WzFZuMGrXrGn~5WLRnYwFw~2Y;6ZwhaQRyTPs{Fsf@3FM<07wK5Un||KfQS z7Nl$pUifTtwcP1c?&2AbjRSMUFdJ(K(TPD+3KODs$2`MZS`K7!>(qnAz6FS}TA1zQng zI#?yL{EanK?D65CA-=R#r49Rn1neTOa6>Npybu6isHktMV@~i%#~(D zX`Qc1IW!4)k*7h+gn~5+UnkZuKJ`;Xxrs$?y_ftPR zDnZ{;e~_SV(@r2pHMQ}4ktuR71ddXTR7?idjZ%IkRHaYK3lOyh*D9rQe>rq_IClba zmCQP>oWPvR`fqk1B&tH+*Jj1(RPasg3L?|M*PnsWw#6OfPACM~o%3-w!VpI6L z#jCJ@-tc5G$V4T^vH6wn0moUPKvkja1dx4P8N9O1L>|q+F`ESfd1CQtxm?-+|40z` zyLK~@<86^N#^#>*Z!)Q5nbIiJy390ZXYIFb{7NZP-pHa7k7Pg_H^mrqtx$NYV7cHpY8fo6m@S3;;V$1A*A!R# z2if2h8(9K(O*Xuc3PAQ{hNZ|o{gfJFnZ{;QfE1 zvywhGAzRA~U_;)w9%|9&lr3{K;VK+IQ4gPKqmixGMOo~}Y|MPBJAGK&YZ89>bVU;r zehS8Z0kvj>NL9Jz!mC}4S-u&_g)mpm^YAT>JQ!!)gQUGIbR|(ZlC188VZYmQjrdq` zc**iw`i^&Eh^(G(<+uTUrM5CQ1i{I<$(34JO_rqmO`9-YLrm2YJ)+t-o*FqX7)}V2 zSKJO7t;P-WhDM;ttU{kF# zEjDb`?obXbyXX+EcB^1(oO)xjhu2-?_Dm(h-BVJ{k^L&}n#YYKj_@l673wf#x z(tziN`?oY8yK1@B;utayUc&|TP*fGO%Vxy!X+T(Op_rN$(bX=)4@XwQ__40vWvC8q z=jt2;DR7Lfo@L}BKA=l*p?S@2M3&5IFN@=IodHb!2J6jS)k`pMG1mCyVfyIqoPzq` zc#uo4hbPCJqUMf@ex5>e5Hh{`ri`z!-S+NX+xuce--N#E35}Dcvn)obGXAZUr!d%b zzCUP^pUCys9WGT)aMb9&vA_aLPpnfm1?F7U{yCST>V+Im@{Q-e4=9`W`4<}vN6|v~ z!L8~Mqaki6O$!?S1{Ri2er7i4={jrcS|u5HVZLEWPsBErSScOlML%W*mF7awWq!H$ zF3&e~YV(|}_u+B_*OrTpZ(ef)*dkCivMAf|gP*96gi`D)A1scBjT4l8U@;GpBd`3W zm5^!0>w%VC@XuPBU9!&HZzxO_6dQ{{a}vdw^j=lTw09mPL)IBFcm&3Ez|jsToXyh( G&Hn>DqaZH; literal 0 HcmV?d00001 diff --git a/src/save12.mat b/src/save12.mat new file mode 100644 index 0000000000000000000000000000000000000000..095aa945bf087b57d4a3cc286d525948cacab20c GIT binary patch literal 4485 zcmb7F_g4~**VZg8XJ+KU0<|<(?k$j+S(-{p&eYuHBKOv;%p4zSxNzat6nBmsI7-c3 z4oFQLI1*D-0-X5v{S)4M?z!il`@_A@dCqgsW1@$+uXhWqc0&~paN;MGe-%(Jyz?~42bi!(9oAp@FDfU!??S(|h|hYNBY2AjND9(H6X#%DMY*-<+x67HBbaYNQuK@4?EJzWJ7Ya&MWEb%Q?y53^M#FQ zpM2wqkO(7_M=*VoC_m*a>`s5P&UJYbRF(J85Ec_#_6OznR9HWAWsW(tClE*+I9LxW z!&j=(3l)u`6<~B9wRC^?)SVy zs>Cgl;IFBL=GxcUY|$*ih5SRo;(idUTkbN}=2*`3cJjb~c08u%aZT|N^zy7IpBOEl zl4rluKlmMf`jnhqoBp&cZ{1~JM@3eiWKkZ~Z-swCgj60R>=* zJ%Jewg@W9vWz?#j2qHr~F9bRj!$`t~>_($$8Dkp9CV!h8ZTk#s(?tlf))=lvkUh#7 z*)p&EA0@139!xz9{6|XFyt(a-R3h~0Vgxb2qnZKjWyl;54soR8XiR{k%i+B%m+1GR zHK?EMllBRy;c7zLl}GBl1Yoq_#vh9lAB70Y4ApEd%pRGmbAa@_5s4XGBCmhMN2Lg% z!Uk|->oxO?zhjIWt5PKG`$u7$I>VTMk9>0SCsPc0@xl3ihQtm}`EBv@fH$wf=g;1f zIS^t;J*63R-;u9N@JmwmDC%)yuKg~Umw~)F3~Tx}OT`n4QqYd%tjII1fZYJkcjk;_f)~8q3_qf?t@G2WCJ~sw?J;wk zXziM~CVls^MP~^Wzor6vgxX(f5|4X(V%-FO+vo_IVg zXTBOC{_V&@P_**E6{oJZIGHVa@zh9!(2KK z+833@M0DO9wPMIXYshF?Uvw8OxR!TgN9(ubR!20{(bjbkL|Okee+oe>KS&~LByQwq zx5R{<;F`V_Qv}}bw6LV*?b}!VHH1LU#B~Z@vcYy4c~JvGO<(_a5AaBDjv{xIfLeJA0CJ?-0GC3CJBBuiwQvrXmzde?Td|E~@(a_)7*e@459{j?rum~An#nR zN+6%~+?#Rm-Fre{qYBl+=dpnm*UO)$7IR)O;(D<=Hzs_tJ%S0zwb{U`P;rN?c^g+b z;83l%T{Y+ULwI02ClMV2+Hrr%x~E8^lTP#4KpOccMQ0_Svt^J;zPOM3u#2^R?Xe*Q}>%e7BjBhMaHZ3AZV zFKg!uDCYx%1Ej-?6#Si7j23{ANto1Bg+;d<^J@=pdYX1uW~St2#)*Zp!AKXPyJxLm z8MV9yTBod`W2dI#H)>f0zGJfM0GLIt$b`$G0vxp*^3TpE)m;GPpYQBPWxf&0x#h>Z z>Xzz$#&QYSk~rx7Vo`M9K)Z({^bP7+&dI?DKS_yIKv@X|)?6Fi z^nmxW3+0beF#$l2CsM!y4h9w6xukKFUbGL1S=Kb|`9}(6+kTYb#vF`0$lcl0GInha z%V6vaf_C()@D}DEm#QtXcki#KmDB~M6P0b`Ov|`%d|G_%C?#9CU$?`L$8)A8=#It z6V;^EHrv9h>BHtfC8G-E4-goy>?YZ^l~2E#vi5x-Ty_s4x{ZFyq}KUccFRl$F$c(& z6Xw!IWIKJY;Zk*Dc^_dry8uZ;!wqlyO{EBHh3q8i^`XdeE;0T)GrpS>tk>>1h)xoi z%^^4aA}d8qr`>~ZIh7b%EgN1I-D&&44N;u9=eamEmcSfcSrh&9wa*2@ZMcUD(NU91 zed~BAGgMyQfc`Q4>xwg+0NFnK8Oi@&RTQ$de&Du^7 z{(23t(w?CYQ13mP;>>oL7N<=NykB=$o~GRzM+Arz{)9-WzD z!)#yu%xdezYDO&_^y44$|77=DRG~J`R$6;ES!b`W<{$u9P<^1^9BAjXl7P!(+z!ug+^WKn%i=|f#Y9RwC3J)(_AI45-@96* zjWV=ZCe$n%CWRnchBi4{Vz}ZAaRh43wSA!Zb0)@bOFUH#`me>f`k?@&3}yyV9ECLWRXfAKPVQeil>% z#ai+|Fdqx{%49wkegEY7o$H#9v9nIq#2*p16`JEPyt~ho$H6T}&-h?+RP~Zb^$NT1 zvLu73UVX;^=0b1V^xeLMY#T9`{h(Y4%~s5MsFXT1{Q~_x#|O>nlbIfIjKC({d@R09 z{bClay!9g#CPU78rRs$csS3_?v4~fdQ`SlPM_JjC(w}XN!?b2Q6hUq8FlVB8 zC<|U5I_D9o3eV1)70S03r4YN$JbISbpC9>ZwQI`$S3jS3aX&Lssy;z7K&N!~t=&~& zlo()Q+xgZ0aO8g~xfvfrc^dY$xdI@Cxn%S9;-Twn+y-?x-QHeQ+zKuw#?7POX<3ENIq@UtfRGmxW|C7q=8wV@~y?Gn~A;+N$ieDo{< zN%pZ{v5=$l5S`79eWgeCL?=G8k^Ll=q_;cuK+=luZT*s6&3T91!sIeIt)YNA=ccsY zs^YXsD6+DaQg2C=vE?!O;&jJ`?J}v%w|$|=aSsTrds%}^`X&;LYLwG->xwZIE?~!h zv-HDe)u9_#5kby^qX@OX?y`3IX_2oCJY8ZFP4T_0HqUysX`hI6IZ%T07 zz+4FsB9Z(3UalqJp#?CyW`jpm>u%*)_) zwpBRPuKyODAGYdoG(IO#Lf7n~h({+(vEl2&(+zljNeN8#yjJS-Rcw)f%+#~F$OGt0 zHb`BL0Qv<@n}PP*7T}VJi&-jc$`X|nXuO8Yl{M^J!c+c>RCRfNu?Ah4P@tur8X~06 zQ^$p;0}d=m2*8gEijZcJs1(y@f*xS}P(7+)ijS?sDn7I_W!o1>G6quqnsi;g->d4l z$309>&!u?GoT?d&Go5r^3;p;3*y<0!KZdA>t3)Tu?@cmVd1j*;+G@wR0%5#iwAW3L zC0CICle=XzcJxAwTY(4LSb@8~6|oy~6lb^xG~2#-B_Nq%e1AFhYLKb9a$}9MDpJxt z$iZJxgQ^U7coVz1BsVZVYSw#~V3`u=YJ#T&G`4Sf54(G>vwQdAz1J8WkS`UjO6icF z*-<%fnZEJeqJ8k5B8lQ;(rwdK1WsfT@)k;v)%Zl}xeK|ms9e3warCzGbIOLW`u4fz z176(rHSe_>TdNLN9&b-LIez75cDd+ot%w?F{oOHO7pZOmR#QtEm&*-adfKK-40f6|5Z|~_FfLR1JtR$>tESqvf}~SZbH{VZiw|-| zU^+6Q1~*Dxogz$FfwEdpVoV^q^~V_<%s#Idl3t0W>uMw4y*sMIvT7bD$Jt`xaejzb zmJUF_uMJTx#t|2BhL^{!w6T;x#GCY^!MG;^jV4m1bZETdsy|_6BpHAHKSeqESO5S3 literal 0 HcmV?d00001 diff --git a/src/save13.mat b/src/save13.mat new file mode 100644 index 0000000000000000000000000000000000000000..e028544fea94a4eefafab8c42d2f65dceda50716 GIT binary patch literal 5555 zcmb7mXEYm*_kN37t!hiDQIwBW#Hc-sYKz)4_NW;|5i61uMeQo3s&?%?W5?cVRg54; z&Ct{qL4N(d`@j7^FYb$b&$-XJ&$-VzkB+MG8&x$fDS#N4j;gVUvzv#LFqgiEy`MAK z$5Wo`)f-)NRcT3KE_EL#dq1bQTwpJGE+dE+mny`COG1`QTvA?APF_-)OF~Rsoa_G+ zH*RX`(BHUm!;$*Njll11b2YKn8jJ%6JbZc7*K`C1;XY?%b#;oJeG<1j6XSn0;t?Mh znQN-6&3n7M9g-LYg*hCOy4M~u2{LOYarZqi6s{eD(j+?7wf;Ge1Ct#MEeF&p1~CI{ zH-BczJ#&F?Ev6GqU)(`++hiP`oF|fOid;4V^fcP!H-x5%V>?$n1=?nS5ct7*8OT-{ zijUuc1&K8E3n`kznzt;42GVdQBT)U;4Dau4RZHqkx1#?3xB?%BB>N|z4gtR^{S#KC z**iTrepFuTeka*P2@5;M^p5NcQ481vqr1hu#X$!EGG_jXW!Uf$g%0U zvEaPbh%{Ha!TtO4Z~bggs#|h)9pwGuvK`o&DG4vNtWbjFpoGwiz8)yIKzPk&{y^Y&in)W>D#T#Wj){5*ne*1#`iZ6-Ft zLTSd<~N{uiI31GW+_bY(k(mv@!3l<=sO zciWBK4)5j0f^FD)c|-^!F~!H;I%zkv=%*5!MjqrORGhWx57fnHYU)8xFlW}q>)xDi z+k=Zxm;$E|x_y;s3{`=XRm2JFS;{lmcQH8Mzcyt9IOJrb&cc*F3;*L{!iEVo3NEUT z!%Yxon(*_?bgzz@_s*4sMI0YGAAs6gZP3sxI2jCuhXVH%KV)Us1~%-$W$VBlvvBk^ zt#M)kKX9D?AIbZRR+|Q_@O2#j5uSWISQLmQ{v)G!hWI*u9CdQh6MVIGu4G1pXOmp6 z_EJ+1Zu!Bf#4Hgc1|H(F(f|+)_d^-YuG)`Z zjs{$3S0i?a*_{LWN+mK^_LskW#`rH>ko!eP&{7U_Y|%eLYMjYcan>SY{~1a+>so6Y z5y6Wi;s33S0$xjE1&Eb_LG_DMg&L>d#JI z37ocsr|w*G@U~G6;nyca0G6nN`nYFeBiK#omD3ev{c?y`hqMZ@gLJFn455W3lNqO7-e_O#AKI<)v&nBTYfqbhRY>ii_B(?acIGkV2grF~9FSwdqx))a4#(nnb>B zy?!Qs`~fjWRu^YRJ6=XY_=j?CjrXWQlOuvR@h&Sj-V$T4^i8N%cn4(za!DmoyJ`4Y zJ4=PTs^lf27|bn#jiyTfV=NERPFn%WB)I>yFIfd~(~bTrOfmhqs3R>&6Q##>mES~1 z3)?(+_W_&84T?>H7h;veE-Qb(+EiU=C8vd&iAc;JaQw(Ru3F<0I9`_*uOCFiMG5KG z$6Lr8L_io^_xW+1bCTQlL6He3O%+eeA6E3*M;~y)SSa~ZrDmHpau&`QCH-SOB*WJu zsbXETj~+AM1+99X3&wI!pNn}7(oGbsRE25`%dWMpgapyVTu)!*t#w75aABczKADbI zg31&T)>5bR46$$E{st$uC0=h~@Vh?jECX?t_IfcRTU;I3xt|Ja~NV9GD zD4RSs<((d)B0Nt;4s7$jJhr93Y>ic5k0pWIg(>}V2u8X7^}x~sd`Uis5B$`30u|@gD)eD68b|rPo*dB zZmp9~X#8sLP|mh^%c_Fj`30seB@`dV4$H(KRhbPCLD$O2K}M`1eW6N%>cjCxv$AtX z#U=8H_5mi~cLZ8RCypqMw2HkD$IHdZYN%y9n-0c>76uHU~!Nc`VoohniZs&dd zVw!)|pW;6<@vfWHWIEJjH-VbeKs$R1bWDf+le*u~(F~OR}=a6T}AF z9J2*qx2=m&B7OqT8qO3oN6QU;>YA=uupF?+Dp>>~JKi}7v4+Phv`>|QioDvP5r|kk z2``STLZ;05SSYgz18Je%7J0Q!)rm=chCo&jh*0rurk2cfP=3K&$eFi1#i_xpyv?NE ztfEw>egTasw7Bnuq?9&|_e4n{nb5-v=JFu0g4noCeVTA9q>kHvLKz7@e{?aQA#I~q z6Iamb&E}S0Q5TQQk{;_o)}1D|yFA+yY;{fZ3@uQVH+E(&uQrjMYjLQkK6xAKo6AhIf48zWyZ5QOwM-Cert@v>lJMqcAcz0n+WY z*9T|s-P(?=@}?jc|5dnZJ=*K>W>xMoZs+o$_Cvjq6ROH}KCg8Eldm8r+2vsBQ6%q* zJFl5mfvH7-*~D%Md@gA=8{dYUIw8*!2*yiSl#ad>jL>p#$#>D|{5P*crb# z`eJf@xR70oV`EZ_QU8fjF1LjPud$ZXaSGrG#29Gua9GrNd|+YVsm5erw`jJJ1*i{<@AYi#ywJo`H(oQPadCJOyQ1Vm0e+_ zcjYDyCc?=0il*K$`YVtVl<6awT?sT9)zE+yE@_QK$`_}5lvMrF5G`vBE>5$VrkBt6 z9qvCG?$d6`E4xsig%}|rW@l3Ot(D&MeH{rj2y?y7@Jchk;5Cc2#)1GWB`ZBY^0k9p zo|QqCfUV>Zy!!K?=lCO$&@a&!yE=6GY9mhMQZLkJhO_EdOV@PEHQ0nY>(M_C%^_?h zZ&K&v=Vg&AK(I^(1Ok-Ki^^cgjcj+$hkQ_IBfg7CBwse6Ve+w_Y6U2WW1WC@a#p!{ z{(}pv6b#U<$6~4*vYm!&mI_13r(r_VBYYUQjN)t&M%$l&`G+G8Zb=*_=Z58B)wsUEf$xu<0rRDgGj1o#$iZU>O|qjVX#X+EzjOSHV6_YUu~%F!wgL*3hrZ$fLt zjhwe_RPwEe00M$+*M!=7%;uBZgT_UbtdfE})tmISCBsVEo4lmIt;=RXDxUF4Bolv4 zTW&M^y4{PMR|ur96QeR>3I5)N@#n*=iYFR=kaQlpnIN>nh**@7**P=Vw zJM#WQG|4qBhT?Q*C6Q`ih`!!Qdvl?{{kl_#LJqUHzRrpApSkSQ2#am)hEu+UlKfJ~ zvWN$XyG<=SA3CN5Gp}1e+4|WqLo^BJ#YqmrW3|HLxZr*rP=AJd55|3@b+?zu>WUc3 zQ@#^J_cx4BJZ4kOOHu_r`O0eKo3+0hgs~s5UGPSK&@2H*-IS+3%c}vzDY= zGy_F<9=4Ze?N%rUTc@XgnNJUIPUk~U2r!VqKH_4>^80_M7E)kR!w$o*YID|J--$VG1LXfHPfjq zWZ39x-ZJ?jgc+)&C|$+Q-IsDhG&(Cs_FW=uQ7>!6Ht%-gbUpx}R^wld)3^ zv<7(C1I(&3tb8tK{896-WNE%2VokI>Dt6G7d7h%qk1kVAFTk5V`V(X3!41*XcYZzZ zaM4rLsVCt{>>nVFla0;h$bAv*!2wwoqx-wr_i>0U+jVll1eHaks#S!BduYR;bN<{? zzqb&dH^2X?r*prMtY=b3Z^ond>t5IZ*$A#hU|Y{w;7xkZRO}c&@Qr4@mYU)Z|2}uS z28}ym6>Cv_iqlz{bTgH@c;E>zgLv=<|oJHDvwmcyqwY$k5nn8L3Mu@Iespl(okCtGP3rN>Xp zmQDo#Q7Pam)T?Paz-m)wE|KqLGRv%|PVuccPTsid3RfoH=+7F!bW5R6ww!A`Az^yf zA=z6q5Fq8&r_R(Vd92npDBt&53&04y2PZh~3b zej02+^m`tE5)q4$cHUFz?#LAe%!m#T#SV{jhfHP5ewpj`JdYQZc&(*ju((YQ@DWA& zCV_wnZg~l95pO1)D5$0vc+{^!c8bZk`?H>}_n>?wyuS3^KAJ8pB6VpQGik4M)3{OK z?_s8M%(_WmMG7k7MB`{phn}`Rn6Zjp=`G8>Il3z#~=`3r%y!jxyS`-wvG6N2AYnKh3C*dV=UkF(gz@5QW42tp3K`gYP#TV z%fUG>k6G#SY!tCs^2cl_I3xFu_8FO`?2$Xg5xi>v0=N7+p27a%t)y-~THNn%OVt?v z*+28-P}leO1}wE!GF3JU_|<9gPlb!_4DR*0Xx#m+=KWM9Ojj>YAg1+(lOwJ;#3SBV zLkeG=ijQk?egE>B(LK)DHPqBKdO`#0_kzFeBVaku?5u(+m}jq~qXtdoKpkmM@v^^{ z)>P-Pf{|(Pc#bkXmT-O}w)&BF4VlpW%5%r^GwtQy$hKRU8uDy=`mpM#%R>5f5PP3E z*47>z@dZMeyLD8F;kD{AvFf2)RIA<(6)nm6`pDFHLP$9N_Z@8>k?Cg)b?0h>{ohCX z+D6C2YHka6uL%!TKgc?KDO6L74S}HTbkGMQ;LhsZZ@~>Bb5qEkBM&?#x6T%1?rcr6 zajxpyg-y%!a@{>uj`U?5@D}F(>#4iRNO@G~%jDt36qFt-|U9KM1&d93>Xd23Ei zY|1931f~r9`;Raq#eJ!)_#Jo_b5OKd;FsCJOc@x@@&ym+e65{>wUoZPXmj*8A{nSedf68W(f=gs5nD;J zWd;{#ciY~lbgPbZ`|#2&uHU_F8U0pgMu8R-nCcA5>YjE?Vgi+jPMBEi<0J+VpoBi1*us>yZP(^Qq&?TNB_A4`Q@!~}Foxll0)s}^ i$ZvyE0SKej=p(1-Zl>FJDalE#J;N_a#k5il@%`tma{Fw_8X#hk~<2?;+A9yxG; zZHc9Xf(ge4r=Ob~&hQ-RaX7n?u&23`k+)`pyKsSm-=9X? zGCQ6FKkP4LlvL|4xG2xD!hQR9eyP6y2qor8WcoW0=DXCo0`tI{dTeK&uWgW|0T@qe8$Rx zPUaa5+JD=p={xck{QP+$ec!$I1NRuW$GPZeJBOIBKh|pJ$i33jclpGNaXPmjIAxW( T&UX9|=zHyW>TiZcQtW#HExLpA literal 0 HcmV?d00001 diff --git a/src/save3.mat b/src/save3.mat new file mode 100644 index 0000000000000000000000000000000000000000..cade101292786f36f54524be401c59cba3aaa7a7 GIT binary patch literal 472 zcmeZu4DoSvQZUssQ1EpO(M`+DN!3vZ$Vn_o%P-2cQgHY2i*PhE(NS#k5il@%`tma{Fl+|miaC#y6B0NQQWFxG zFUjx-NU=mLtT9$dV@q;Vc*geV$)7)qrkvSxWYU$S4OV^%VhO1Z%<^w8vdDA@*C+JP zye-{_?&h2`?ocxQ2O?h?;3mo~`zZ5v}r-P(-T#{T@ z(SGs8%Z>^5ZU+Mx8FG(iJA<4)73_4~bFO@c9e7wTT-cJR(3)9RSH$0O=DbYA;OAB)_7Z@%R8Ta$bYPXtFCc^-E4Kv7hWg4YrrYbWzLrc!fRJVQd5 Z{UZCGw-@$xugd(LdGr4-X0#k5il@%`tma{FgyX`iaBp58}c18;BoC0 z2xvJn^_F!36AO<{2ov8iL8G!WCvHq^T$s$xvZO(BUe?kao~RanYiD}4t94QSS7+8x&5nAZvR!2b&k)zCWziV zzsGz7i|mirf8}1wJ&)M`iN{F)^PPj-6->#s8}cXpOKEj{iZKxcNXo7{r$kkGXL7?-z%Ko zZGN@+OXkG?n%C0v_x_!3%Xn$avezKLUIzR1adJXJ1;ddXhck>z9UDYXU1v%XQ%Gz1 z6kHrs9Q@gv$Rv!XDhr&$$6!trML3kGu^|lU&{-qEMKu^X|BqVmkXbrJAUL) z*1_`!Po=8AUBDiaob+tbiZ!QJztcROn;;^glb^bH;bWCI%bq@b`}(}zi8H5;ojrY6 uUx8;!cvM)Nqo=F0cd)zv{R>CV9J+Mt+`)@SuU@`=`nbM<1jAby#YzD3kOr3k literal 0 HcmV?d00001 diff --git a/src/save5.mat b/src/save5.mat new file mode 100644 index 0000000000000000000000000000000000000000..e4d3227d0c3891622d98bbf38ad55d3cbc560360 GIT binary patch literal 732 zcmeZu4DoSvQZUssQ1EpO(M`+DN!3vZ$Vn_o%P-2cQgHY2i*PhE(NS#k5il@%`tma{FsLyyFjUNWJH?vokb^** zyGtO;`YAE4R}LL(xOFq4il;vzbLJtXR4#6BZ|?3XVv}~*g`JqLGkMFT#;3{CPM%Z# z{XW2DQKQH84eIjy)_neYJ#FjRteL5YQe;@VPv+P~`A&boV_)R$#sfn9(sf(^c89-< z-uzEC<(`Dkl!|Ay4dMnYujk1p+aD`mbN>^Mk^bi;4@B=YN6EjLo$`+(<-X!K=O?9S zr&k(%V0-8A?WpAh^$#4IdQ1c-;QJ{`a*zzxdAG zY`S&ZjUk`uYTJ5?w!2il()9L%P&+QlRpL_LfbJc>l2>tVy-_2faH+^61lDw>M zkvq0W#jpFFQu}J&-vYf^d8T(a+wC*HW$l{vJ}$Y}^qc=Tmb8g~4ZazV>zB;nJn@h5h|RRbJzBGB&N%&I>I*U|QcB{C zS{HX~muBX%WzV=Qw|!$h=Xy`ICsmZuQSq8!pkbkqr(K<7q*rBI);pg$p^9M-V^#?U z1sbw7`5jwUxad*mt2|GmJu6>II!1QwjH+yPIx4sH-P5U|Yo3{H`?l-cq5x*A(4|eExV*HaqGnlD{c2N76PtFe rjIU2`-Mo8oc6E&h7bZ@;xbY+BNzKfeH!Xj%x~is%US(+QaBKqr`vpS% literal 0 HcmV?d00001 diff --git a/src/save6.mat b/src/save6.mat new file mode 100644 index 0000000000000000000000000000000000000000..1a71948e9a5c54078711eccd79c696ae168d8688 GIT binary patch literal 941 zcmeZu4DoSvQZUssQ1EpO(M`+DN!3vZ$Vn_o%P-2cQgHY2i*PhE(NSy5Q}Nt($dr+wI>o>rULAn)A{- z%Vu%p*dBA(b-Tj!-Mi|~fB7fWc^BlKS$^Mt0~7a#_y1Ind9O+SGjF}Z-%ID&YM6^J zoNBI|{N+E7(R-%Htbek;E}hR-!=A42IZ5_}_y^7uo3L2J|9Qee(-a=3#;*N+)GGR; ziQ-eaH5>OG+qYH7{QNHoPd5x*tl;v(1I`nQ>qf2A~nVq`KnMjj_8 zB+NK?;=lp6C6*EjCL9}_ehNCAVLH+iP?zPGm6sLgX(&8lE`!LE%$#ST7oQ}HT@o*t zGHIf)p!bEMua~a*zBDs44L`AQtA4?|?|L(X%ybnrQv6&(eCD`JI^Wczv()CW>gljG zt5#*M+RT%k``+iZ?kuZ`8?6W?%C>??vfZRSzD;%suqz(x+2f7ykL`t2^=HN!y#kd(?gB z#^oAb&s%77!mTVV%q%ZXu8uKL(2z4TvXWI*lvmeR*0$DJ*qGZp+nU{7JlsEDUw;3A z2NyP+cyZ%L%abWv&b*oPr)kloO_4ZaN@Sb#5R*b!Q+lI3O`JIy5;tF)<)AFflP8ARr(hARr(hARr(hARr(hARr(h zARr(hARr(hARr(hARr(hARr(hARr(hARr(B0ZB~{000060ssJbob8ssN&`U5}@yXu#;o0|lo7J<^ z5!kGYiC(k$NwKt6)^Amg#>iuVWekTM1tFY9WZFV@+L$Ay#&u*CNz z)-9E|xepeZ*Dc6AGvu2gj;wy&o(n10AoJk z##4;1K*nv=4St6%q+EhLxAez-CwabxnCJbX)MwX)lo2?FpFjS;CoW`U)am88S3mZq z={I$p@oD!uzo#x!AGyAz8`MMUCC|M#O}`t~KYE`!+BOgB9raN8UzY1xdZNXr>9x_?mgw=-|J}X6Yz(*;~UpMLKWk zq4u9>`t51;CybMN!9J>OAI3h9{*xt3mMmGaWdFfF0XWJ1Uk?BP06hW#0C=3^V_;w~ zU}9ik1JWD}5E=wnq5L`^W&vVGC_gnPH8(Y{q!_4$;RjHVXTv656-_>7A$>twFKykFZevq(g&bS^0_T43k||{e0*;29jEIVflSuMekddX7XL$PTit|-9SL-aV zx3+ht^%hON`{>DxS=O)KEPD6h)0Z{t8k@Il+p%lUzG=q}1s$KWWZ9~Vmz%eHR8SQr0#)b8h8@i4ZU_F zGWp)r2algtTfBY?aw9u4D;qa255Iu0h?s=5th|DfiiVb^j-I}uiJ66^wVi{bi>s%% zmycgyaA;U$bZmS=a!Oiec20glQAt@vb!|gaV>8f}?!NxXQ>RRuIeX6h1&fy~Td`{G x`VE`5Zr{0k@BRZvjvYO5^7PsBm#$pBe(TQNhmRjUdH(X%8&Ldw2LQ2h`ur2YUV8ul literal 0 HcmV?d00001 diff --git a/src/save8.mat b/src/save8.mat new file mode 100644 index 0000000000000000000000000000000000000000..297ebba0e40f34f3beff4633a61f43da9d089268 GIT binary patch literal 1799 zcmb8ti9gc|003~QJ?S7fF@ByaWyvkq+S}7+Wt(G2j%em+%(cg2ybz_Qa_27R$dOzn z$CLY*VZ0nMx0b}rBIoNgv3iI*t z#qsm;g=Vj%LoiM-@mLJ=niq0MaH-_YTWG~;^&?`R@h#r#!gN?4Jz;>5O* zpK76IZAa3Pi=^if>=pT@W3t$#wkC_+BSGzTNyhCR-q#z@V@;rLfa?zKah(-@&Y?BV zAIwovANXZ7SIc&vknK1zDf;{?J4bTd!JK#Rh2jTRMe(_3H+|0`n9y~gur}y9#h)s@ zqPNq<^Ma#fqsqQJK|9w02DlC#nazQ8c>)V|sVI9}Zf5Buo)u7xlI|yQD7#TdC1!=_g?&^jHr=X$v!Afc z?g6dlHDgTeGX=;}HRBt8^jA|Z}ai8H=s>7GN!pj z2Zg(FT&%E_=HpPE=1!i~E}9K5(!?nQxcawZ6;2;b9OUmRW|k(cG(Xc^j7ZlY~gB20mMetE%oDdj^+y*~Hw;V~WbDJ1@O<^-D)a`@t2o^~cD|%Kd zqN1nt)5asSgF`K7!Tu<2G$l#Hd2C4xZL2ce)?Tz(YJ*gE%Scg74sB!9^-87%Y>4P$;4U0 z2=5Go!>`HFt_qa$Ai7)4}t|2aTn1_t_V^ zDN!XmEdM@rlvWwda821lfvb+=oo?j_)Zy%PykFkbvlgpRxecjJH~*%F0V@~#LleHh zue!4f7WUgCJBEi?nak^7*rkh1(~m^)F*i}4B=Hv!!akW>eDb4;j3yjI)zxS`7Quu=)+f@;Y+D z=QLc!MndKh5%pu|%Sx3-t$XZHRJ|c5?YcjW2O?w?S0rm%pCQ6c>Efb{J1-LI?QVl# zHgAsh{s{b%4-3YH@N7xE;a_qWz0nzhF7+Jw^{EcE@GixOrI5o#S`Nd#-kHREXf|lX z>XT;P{{86>Cd`IvQ@xyl&$Zg`v}~F1;+nu9g^40jdfo_u%i1+(g3!8*%`$kCMT46O zeEj#^Wc@i*5Ju7Qmbf+SxO7M03}DwD`e$YB%M$0Uy6AxHcKC9h?{4xVgYD9AbVL&0qldI=Q?t;o67St2}D9I#(OOPSPI+l=r z9cdh+&&xaM|wQ)nUW91&EvPx?v=67sj+6a4> zxGlu(mKcc%Z_J02zSDQ%dxl?(Pxp&_i*0$1f>s_~rc#VRixwtPPIB_3xj nfwOT)*9Aqg`>=43=dDx4A0ujQDsRPt*yTD?R_4ZaN@Sb#5R*b!Q+lI3O`JIy5;tF*G1DFflP8ARr(hARr(hARr(hARr(hARr(h zARr(hARr(hARr(hARr(hARr(hARr(hARr(B0ZB~{0000i0{{Scoa~pqZWBQiMsI|O zKuEks9>abVBcMn_5TT%SK|w)bAqol#3#CXwu|gCnD7K)Wbde&Z%M=lIVId1e8J5I@ zG-vkmd9;XXN3xD)y}S2k?(E#}zL`z;%yt{zu^av#w>R^n9kW{v8`!<(_weB0*`xg@ z51v0hyuR;l^Be6M#&UoQ!|krm2krAX;(QI#?hIUt`)T(QoP#x3`u-WN&+r`mk30iu zmv*kebFczc9E_W!{cj-Su5rDFT`#}16Tk~_3NkL{fjkAjg7in_759&?QrZbff0;@T%SS4Ss~sEepK!69fp*43^ER`5B>ZCc3h7-q_i=p z^3L^B`~>p5qt@xWGrz$62&CUCZ?HQJDeWD|x}`m?JIVT*;yS-CPJTuqrF{fXVCNt2 z?|DdRPeD~rtUpD5KUr@dK=pll3#xUcy~4Vhg_QORq&?Ohxe$c?Dx|dcAnQ}r=>*p& zh-1?3M^={|>u?!T+6Yv2$h!Tx;U(s4%=z55eJ7t|owF}^1F}x}JpFhH#u5Mj0dn2> z{=&LH4=L>hNPoC5NbUpJI}0i8b;|ML@67VobzgL>?{!F-jpe|Z58LN)c=gJP5<|-(;TSQ2n3RH+uRilIzGiX5ET-vbuyGkoqB> zxAA?fW3I2#FG~L?{iN!i`BM6i>!tP$^(gzIuJx)&eb0FM(T{%gqaXik`~m;~|NmoP z7zLwX6pVsVFbYP&U1ONbdoHdmNlvHH^M(;=0F5TU*bayvPvvkJ- zODwf?vvhYisDKiZDj}i*VqhQ=q96vg0ul!Jy#e;joINvV&%Ay2-v66$;C2%5gJeAY z9sDGI_WF{aiTu>F*U^KcheXH5jSSw~?>|mOq?#*?0u&`m2}wsOx~LA>$U!bLkcpl| z(14~CqXd1KPdO@3i98IZ2DKPQJr)y594i>jDO&L^QM^Y-I@6WzWT7|B*unq?l8;!* z@ID{1lT;edM4~mBDx}j)W;2(5EMx>Lxz0Cy%No{Eji301ZDdjw-62M=u$O5J(IJlT z8gFo%ldR$lr)k4kHn4*iiO?r}#&p_qh0+?XTa;H7Rnzy}riNNDN@G+{<8+^g>}QG| z8%_cxQW7Oo3SA&vpYsJ7^c#OLTRD_V@#Iqh6;cr`RdJQna;?;ARn=OFv20?!HfpoB zs*#%LZ~mpF+Nrg|X|Jy8s5Ep_5A{+XrBQ!n(;(&1P?gn4)lwZb&;(7=0j4QlGj)%- zny-agtleB=DPL1jg=jaP^S*;xc3Y3E_Sl@QcDRlF_Fa|_=%rnRnP?yN zQ&lFDURP}P8vX5aB+K}U{`!u9ZFPR;S2vScS#_6zCO1gaiB&$$rhpc*#Aa6$P&{q3 z8+n13C>G*5R2|$)8I8~lDws=kw^50jYRPDgReg;QF?yA$nx-Sf>o_xYmXG+@o-c8k zTl~(S%u&E@1~Xk!KyaCBiDHFTX^kprohoXai`}42+M;dRu7~`?e>|p*ve8lH>b_irX=-kx`_?7W^k1*eCZAbnSP{g4bWnJmHf>}D};=**1v`Vj=p+2%}o)(zxEatIDOH@pIXy*dT zse;dICG)PXnjz-(O?`(R5tRRg?6%cTQf4~&Y$a4urBrH{nNb6kTd}&t2vt^4&3HA| zJ{NPC*Ll-zJYW)OZEKmH7O_sxOGrvWdrD#pVP3=3q%rwm_pK zn|SYWZTm#dmN)~1eqG_a3IuVgb?*hyR8O?z|dLO0v)XI`S zW6f?7Q*CM*@yue5Y0PJl_jI}GuJE$0HT?~CzQvYpb%4Vh^}QYAgzxW^ z>p0`(J?G^O_Vt1*xJa;LSMC2gH+@EeO%2=!iyCa(9q#dfM+A%Y1es=H^Gt4{smwHu z(wSyH?K^c8C54I<9zdOQ=xq literal 0 HcmV?d00001 diff --git a/src/test_solveError_Grid.m b/src/test_solveError_Grid.m new file mode 100755 index 0000000..4c745b6 --- /dev/null +++ b/src/test_solveError_Grid.m @@ -0,0 +1,57 @@ + +dataIso =[]; +dataAniso =[]; +maxSize = 500; +Netz = 'exmpl_2DLShape'; %Energienorm sollte gegen 8.28466 konvergieren +% Netz = 'exmpl_3DFichCube'; %Energienorm sollte gegen 162265 konvergieren +% Netz = 'exmpl_2DQuad'; + +%% Anisotrope Verfeinerung +disp 'Anisotrop' +load(Netz) + +ref = 0; +while size(elements,1)