From 997ffa50b4cc27a363a062c05c4421226a3ded1e Mon Sep 17 00:00:00 2001 From: wolf Date: Thu, 11 Dec 2003 15:36:51 +0000 Subject: [PATCH] Fix one bug in the computation of the value of the solution in the neighbor cell. Move one variable out of the loop. git-svn-id: https://svn.dealii.org/trunk@8254 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/doc/news/c-4-0.html | 10 ++++++++ .../step-9.data/final-grid.gif | Bin 10221 -> 25571 bytes deal.II/examples/step-9/step-9.cc | 22 +++++++++++++----- 3 files changed, 26 insertions(+), 6 deletions(-) diff --git a/deal.II/doc/news/c-4-0.html b/deal.II/doc/news/c-4-0.html index 012c4cb7cc..e0d5dde0f1 100644 --- a/deal.II/doc/news/c-4-0.html +++ b/deal.II/doc/news/c-4-0.html @@ -79,6 +79,16 @@ contributor's names are abbreviated by WB (Wolfgang Bangerth), GK

General

    +
  1. + Fixed: step-9 had the computation of the value of the + solution in the mid-point of the neighbor cell wrong. This is now + fixed. Consequently, the resulting mesh also looks much nicer now (much + smoother). +
    + (Werner Scholz + 2003/12/11) +

    +
  2. New: The config.h file now declares a variable deal_II_numbers::invalid_unsigned_int. diff --git a/deal.II/doc/tutorial/chapter-2.step-by-step/step-9.data/final-grid.gif b/deal.II/doc/tutorial/chapter-2.step-by-step/step-9.data/final-grid.gif index 91c55c068e3e6a689fd053a53e1fd951120b808f..d46812483688bd633b244aecc5f89fb8f124ec51 100644 GIT binary patch literal 25571 zcmdRVWm6o$llKOKOA<7AfZ#5RYZBbu-QC?1+zIaP?(Xgmi!AQ2xa;!df6rajeTC~+ zQ$1Zh{p+rosp;+!lMv0f^mz-v|I<0018VKm-7!0e}wxK)p3Ud>nxI z5r8NJ(Dw$+{{XP!1BCufZvg-c0Khu{&6;(0KolE3gG@v9S8t``=^0%nH#Vo8&(3tvz){F z4lC4c1OUMN0|fy5)29GHa|PJL12BaEh(`a|`tJaecXbp1=ot_OK*a(q!UNu!O#w!0 z0GN6GoOe6eZXxpz?4bXG-VNwZ?jU5CUqpI;N9^ua0sx8t0Q7$_|7i#40|0ja{891& z00KOA10L~00RZI3983gQp??6&0D$tp%sCM5{&#$qZbYFLbp+5Npt~DD5ZVEln$P+8 zF2wx7gd2d34?si&7?%J5g8$=3K0CmV|M72uHK2DMK>5Fse*n0{|HnzS06^S-hXipM z0GI|yumXgH09Nm?EupX=WWYNnpj!w4AC1EOFC1`a0D!#&0MqApAdK+5pELmG1Hcbp z4FeeePcVAK0gFEX@A!aj#LuBQOUMiFnD5;}-B=|bmlgqEr{9*A{^9?}c_-Y0T@v8_ zFG0ZH-_DD@` zf5>NY$qez@yupaCY?f;ywfVy_-+qP@O4JpM##1YmX^++wjwdr4jb%vG7fq&fIi0MH z))!A_3k4t(Nj8)y4JC-TcT(#cresi)u-dw%e{TmjASh}TVyFV0@PIscEc6T^t1;-hswQher zh0SVXqP6~TI_GCVk~A~m_k4*`x$b0JBV?(@Xgo`%z3B|x;sn{4Y;V5U>L~pFd^2KHNB0mdCUUJUcz?_78{1 z=|R1o?3aqorUeDxpKn%h^?rZc!2R|co>G15_ao4q$n%pYRc^pDTFEzm6hrsi-%uBE z13z^CuiL(uVviSMN7Rhb{eo12CeCRsHd8>Nw` z9%l!`%z5Z%*x;8I2YG|c4uk#fX-YDT`<{x@` zB+KhY?`bXF(xZ`?Qtsxwm}-Z+=olL$J-xCUF26i~Z-P~_I?suA1T(k9LFrE0XIZUl z{llBu&)fEhAytkx*pQ0i8?te`(eKHpcEhbzktdBmH^ypvg~kZA+c7Srin97NT5E^6 zx=$TCNb@rsQ=bH^ZmY!(D{c}QA)fX<$oY^{rHdIXRvEC1;$- z!)bN|hQ@iavGm&-Cog=++>o!?4Rpi%;p}1aHv%J<{Wm#E);t+-1p8Vz=ec{^&j*61 zKtvQ;=0AT?&J6Y!^L!kZRg4H(R^+YFJPxy)J6n%WoA{opsBm6_GCj)FFu6Lc?x;)yBOjXAZTUG*I=QZtb zO=(^8WAzTROXl9O0rO4mk6^WKsO=ee+Vh?&dAjp-!;5e|?x20(nC|!TF41$t@b@W( zNc|^p>d9XOo*Ee@d82%i_gl9n1l6d3N?3!pv8JAPJaoG~g48p`k{>vE7eG|{?SWgI z?&X=9;EM3=3VXyB?6}XKP)hpAZFv)+_)=6iMrWljKOg(~j%$RUuQFDDxOZ-9#CQ-l zq#ZJrbf9K^3MEU}K}c24rJX{JlzCPEoXYl#f~+zmRhZAHbJj*Tj$DDol@Dxi_@Ls`;Z2xhrm5#YVbV z{C7}DNLy+)736V5>6+4C&3O+mbGiGA;&q(V#6AYfqv*iX9=i@Am+WC?L(L#dTjnSzjpznI1) zUMv*d1CFKj=UC&4p-QGWG+dKnrHW~{*55y}s7?NkMo)jBH&57s?~!!GNEWcBBHq=I z3g5XAxwS6h-XS4R&SDxek-1uYnowX}U5T7ciE^=|QqdTDy7FMT%4lMLXP0yde>t}P?!zy) zyE_K?!GfPIsY@njbv3|6h-~K2v=f6%3K7;ZnF7A5MV(7IW#<|&jB`AR^19xEfBSK7 zc_`y~t~zcIdLnnYLU7}d;G2gGk{NHTe0m+R`U%kTt5)dW&WZzIRtWLQfV|0?HM zS#n*jyvi(AoY(Aa+yR+MtTTx(((2y5xecHX60q()AX4mP?yz{xgem``l#NkL)k+PJqfk_9FPm`v= z=GGCrshg`k`eR`uoYX{tn%-#0F86clSQa|FcgRGAWPdQwQI~fLD79B)dy@u!J>YszA82QXn7ckV`J+GAb3;i*d7r1KY{AYk zQ+g*p+DC^n$|jKZ!P)b)g<(Ff3GOHtMZqQrSoYZk=(-$)$z8I!mR;ESVz0J(vS;8{ zxtw2hw{Y$;epc1I7UX-!jj7Wno!vG*Ty=ZO*m=A_VL6^nBvMwqcC4Yt=R1yeW$VO$ zE;M{sse$gw`p)02B*<6Ys^k8KFK`*~?6QfsnLU%8l@ojMx+kIg)HLk9nmo}y=-1(D zrZ)^MOR~Ro>cU@14}D~27YLOhvdYSrepsI3I%Kh%na*CAC?Vk+gok#S#>Ui!cTu&~ z6t_6?Rr}tH`*GK$XP2%mebpmK?&?SSWI}TV_UU`DyO0P|KI`{-?{ITHMzCpe^SEl3 z)D!a=<8%jcQk777J|HmfCflzMGmv5VCNvwumbjQxXw{0?pWr&2a$1=r6OFj}Q>XaT zw)mISI-M=)&xAO9C~+W8k!LyZfpfF%>DMtfaSjdn10!xo-u#CG!R6`Rk&xPd&ebgf zJ0Ok9a^CHC%tF8@yf%+}fZ?LKAmZ;^sz4YqFVdFZUShwQ5!6tFyqF(-#l8n0ni+G( z2T&%HJPo^@J1IB91skV$C@lZ}rS3J#A{Cgr)>>nHxGZ z`>26JzomqTE&qx$cX>xp`c3UIo%maEIoNT;+kenMSUsFR#apo?q~yS#zu2w`KB)L> zcw~wWwAk;HnlyGv1ni(S+p<+PLYT8!AiKNWL35a2XfO(8=%rdf5MqP_nR`V@MB-2c z)=_}>UW5^1WD$Z_0&Xz!vb}FgU^*zw`q3#lMI|fLHQzY`7tyG4#6I0DqdRwv>!Z8ok*8{Y)HJnxue$a_OL($*uy2d6)QD(lN-TbB zEa6J5NwRt#_1~?QC>KOC{1Mx_WQF#YxQ`OPPHJImCD6ZVh_>~>znkV^4qWjXLxDo# zac`KhN6q$B8agx@ia$r=K5;n@t9vjYxye09P#{@jYxu?v#bAcS)IG+Jd`;kY36XzN zQD8M79||#DGSo$iwG5F@MYLB+jW=(l;Khj>JW9%Hbj?2ur4CD^P5q0QnqVT4q+BW} zqmgLwl$bJpL{C+mR@|e^{{~A$Tp}fqtj1r1viMj2^1%9j zd(bm$O`@Wek^<#6OXS*y2$Ip}d^yRz)JQTr`1uX*?_yZK&exO()7%2H0#di^=;9#N z$9%RlX;Mv@w3U3$vJ@AytRMKvzuG(*En}GA|Ikq6+ej4oj>;*Vec5T;)TtP z7O85+XRV~f`axqIVtM?ir`IWz8Fe4|11X>10xiT)l@+Tpo#{kat@ zAy#RHY-!mi5+R4`d0O0sB58TuNE%QB(g{=lp)Iy2P)(mJDxyPQ^Vb;qvLZYsuG$!V;R&pTN4}UMTS&io7MBOq}+D{*+o_uF_47KM!y1qrIC!8Hz#H6q0Y8;Er*<5}*Wwbis)V#w*Qv^6qb^~#tv7G*W$6dv#>e@?h- zsHF067c0Mc`s5_nFFZCN4arnq6$8V=ozzp4G|B>2OEK$XbCJBsmupQ>a!kh>m#ON0 ziMFVL8(k!7A!4N%W$~nt@|d!;=sS@8J3N8sc zFp=G5!~1Ygd*?K3$ujZ-7g>Karn=D;|ES3MjLhV%HBg1$$@H4^9hHqzTD%+5(CFST zDcK=R8;#B#2;}X;%V?9*Y8m2fJZZ0h?wdD)b`v=7Pu{w%N1A8tZ zLve)7(i7QGlnQz2LdA&@J(OWp>p^wY(WjO{y_CwNmq9bDM&F6nljczA(wffej*hkQ zuGeuf5dZSipTDZ@lBt=q%bnr8*1Y3gGdQ(=(xXH#ohVwl??@dz9)D-NYurvK`lSOn zO(*+49|Ba~w4usv!X6DZgCiahq*k+$4E zg~^-ux?< zp>s}2XCaTSBV}UdC2UrhZ=tTe_oQfaIW6&H=S*K~q1@S0^L0*9`~1cHJZ|M8opgJn z^^ygB#Q{CYOJ)+;2E>zIU9mRsjlVmvW4z~eC5LF0jbYY{p`W&*{oS+3fp7STU<6WD z8G%+8abjH^6*PG}AZ)er0Iq+^DEF{V`kuZ@YM$|#udnhs7%{WhND&;(SI$$BJbBhy z0A&DoqeYhR&kQzAuCDi45|&?gRBBYLk@v0HPt0~@f;l=?gx0}i43kyXdVNvx;LNc@ zuhP9raK~h&7~e*S_J%8wo zhThY%HdCib5K~_MpaN9 zl_M%JQ`5P)M!%Usw9xx=FmsrL5PiRJVpL5s0|S3D%4&A_wC`^x2(rF&qV4AeYj^88 z_SnJ_QE~9=d{h+wc-=yfuKf>G+Xpjaw+*x9ly^cwoHQ`O36-6y<(D~_mFM5|w z1`^Jhwg*s2|4C1uZvk>C;k5-JU4EG=?mZu%w#h`OX8e?4MhTfHfk2P5iLU_P8=n}5 zU%DpVrgUCTAu9Bj!n!Nh(HkLA31m_|;yTypc4s4Ll^A|4M^&Z4RkP!f!?@k<=-(Hj zWiM-q&jYJ2rF6ptviH+Amf+6MwKwl{{hUOrcIFvdHg#^@yH>`b8TW*@;^;X7f_ro5 zH#QgBj_23-C%XnQOQh3R1-1qXpE~^b4-jAPiDY)ixgV^eSIP*!yS|pc9k(&^H<`c=)J?^VlzSYL7@U7 zq1V9Ju>ekgLhXLig#0e8xXW+tZ@Q@dz1}4sdG9jyQvP=N7yABi6t*lEQ{ANpVPtw` zqtBeI(y>NASTgKw@H_FvU%ryCu7*CPeR}I=f)?7|D0IaDs~$(}B*)}l!BeZ?>~11G z=truz^^F$*i_q4;Tss;vy}P~sk04Y6T5a>0?O*a=wL2nVMC86}=3k7(ZK`%)tg9We)#zM5g*VapC=MUERv>ttXL$IDG-FOa-v+O zQmH$brUFr|(&_YtBTzkEu7>x4d82VIkdC&t*jezH$<+?!(`aktR;rfn@OkL_IH5|h z-spxxGG+02w%Nf$M7|*rciG?OghuX*nwqEAoyH^vo#K0RIv6bc(9Cax@3On_mr5zv zTw`CO!4k9n95s30YJ*zheG;t$c)0yXF5*7&iOqa4Qf0n_h~If{W2BM|fm^5D2?WeH zLLh+x`!~HKs(l$dUUn`PtH&uNc&R_mDTr|Lf>g|pH(75c?st~)fjxS-{te~IzQ_%oL6D3%QkM3us^1|Tfee=es{$7gWSd6bm zRy9hFn;0+2X2s{sMrkxLJ3`k5EIM0>AjVlU>oCq;B}c70>ZJY zNaIzq*Z$Hf1S)cgsQ|@w$UPU^HPIN@tE92-8Wv@BKDsRG8nK~Ne#Wh zM9HnvExtV-A?idw?qvQ9XJ7j`PO)l{KP;Liid21&Y)VT6`%>~_W$5{5<>Hu!ep)9Z z{Yg~%sgpr1S}?P^Wbsd4+tw>ed9}jS1C6sr*d>h<3fQ4}%P}xWOS8-6zTIlXX42c< z8yom~JWqkp#%o7lw|b3YL7}xa`l0&GVO{Bv@No)nnCRAAgsHr)4WjGhfkx!FtbI5t z4R*6EMv(FV9DK2;q*0)a-vU%JG+3X%~ z!|O#i#t}`roGh%g?|v?wLAUO?KMfi8l-V@w)_~g+G{`|qI zp=&vp)*IzD8z%Jth+%7Q9_>Yf40oIVv@hlVB~n*RvxSHFI!61_=qtcV>8AlUeQ3?N z-sYIwFZ8?^H+pIHFBkWHh(njb?{uFojorHm$>_uO$O#vh5e6`k9R%=x8R0Lsh;lT| zMp5u?d^SlI=slK;Tua@0bHC*^WnppPaLzgV7}kr#U=vJ>oPtPd-YpnB8y!86xW>bk z$D?RzEV2EK;D$pSEAfm5fj%ked02w(8db-cF1#9D8C$+j-ZNL5oK!?SiEGL%T7FFj zsh6umRdf2MH~TitA82W^LS9Tl5&ic8^AvGG)vLs)(`~5hBfs&~oJt#NCLwM}v*h_) zdSSXT&33?me%YtLO6;aQnJLoZ3`(ZqdlF%an_m zn&ZdCt}%jW2h{nyQ)ILx8-!{H%n?=E#^hX`|NN1*z-3ayP<>ALCc2XPWsqjBx|bz; z3-j|GAGKKe6jr)|avrOTeIzShV3B;R=>q0R*MYM#Y3 z-P%0cMp_9cWawfeSzU&l4k^#ziM03FbVUgLoWw($@=-0YJch!Qjo%H_s+mwlYG`f> zY}c4_1=Sn*AAttcR{E2g?6uJ>w3(dL?C-CGC9%tX#Cr~dU_k!Hgccg@E2-J0SZUv^ ziAviLs{H(&E6VoeV0B)8HFDgynC)wFeGvz37@B)6?w7k%vubVoURpX*lu9=rYxgmN zzE_CYsyIQq%#v_h{{l~J7RP0`-;mm2rPKI_({|dR z!SOORxQvlWkZ-;4uD!hUCG)D=yB^ewu|7;4 zea7Ojo#SAjyByJFYH#9ZPjqdi+uKk77kn`J9#@Cpxu72`-||JYHQuJQ?&-kSEPhZ^ z8*t-1YCK~tuZXI5p+&nAyW$ljXWaXH*l&&5jm zQ&I(npDJmPYSKwMbVGvO_b3EC#?&@wVxUgjCPdv_F&QfYg!#J`(EaH{=v`r0^KIHX zzlFR@*T$ibXMz#fF-_L0&*7cy5O)1n&!!TjUhn?=ed=mHLSTbb?RlI$3b^#*bM?GH zb19Nyj{w=2)~0&dht}t*b>ol*u6S-A7K}Y6b{QKTT<^zBKZ~QcUB5&<^Ts^3 zi-CL}{V2Cv$1`7-^f-?oNw+^c)KI=n$8n7XutnAJaXPfK88mz}&O$hi{n$=F+tw-Y zTnF=!t$LS@44)5Xq3N8U>lI13REf2JyoYBQJ^ji`n6Tg*(JkSZz{oo=_w!RvrD4^Vh=Qm1-Yr;YPa*~v?SkM$Z1;v$ zH=Y(f?rBN;Y6ERX*S4lGAAa2s9Xm}yZM9P0t`<*%Qvp$d5RqgT(`la!(ogbRMApb& z24%q!tDYb113!NaO!RfIS1`m8(DOK#0vn2joEvS@g$V|Jm|`?YV*Om3|G8J*kUN&g zKhd6@Ud_hWa~NCEv@4L+_g&lRGfKgLV||;-ZGSCnpN4Y_x;JY^-1oC*jNY9d>I^h? z@H*XSDSU40#9E-0EF<}_lnjJ}* zJm@w~Bm-FM)Q(tB>JiZ0<|WEIJ}h1!{<%EOt>Kq<`j98qsC$KI-?->!=U;L3wT9Az zHh)H9P6x-mBs%NGtp-MlERl157OGWnJVSUh0&xw;2};NQos-oL14GFbB4MZE&bLC9 zy958SkGQi96szJ+ugns?qR0e%9`hZXSVO)U{j2>wA*Z8%VB1|FDdmLYH67y5<0C0f z0?`ZNwR|HL1O1hR67BR7121E_e@Q1ejJzSy@}V^Er%PmF{czQm zs5c&GYZr}KkQ7BtYc(F^Lt&kEl2%`9^SkBHENFrn@Ci<0@tJi<4ff|*k<0jwOBZ+0 zu*6E+e;iu$?wxHIWFsG+HJA{49rgB~sJH&v5GlF!y>jHVpKhU<@^oV3qr_*F$tBDm zAaBVeLbCeO?vQ#u{{@a#>R-Xq6S>UOiC=0GG^B$ZWKwR~K|f{vw8eJgev;>MPX7aG z&Ew_WTE{l{_u~YpV!EtiV(X(wLu4g--1PWI`RsxTvT?bS4kjqO#QU#NGRXL(aVG4f zEciPQ?5Qlwyv(CF|H|tm0NQzle~s8fiZdzx-=j(%~j3R6KVt*>8V79)M; zn~31%fUHxSEGnIb9n&M;?^PJ~;{G z7h}>^|m1442{Iw3NG~cuGr^#0zJ{CY1o~lMEtq z$jYkUE4V}mM&ONVLba!%p&-V_&V&XP7T}x)T(>3zNLpI{(_p3|xd?|>L6`ADm|+L9 z*SOxq;?e-94pmwvamIFCwFO3v6}mrD(yAZ?Y;sy0Hr*X|%bcg6RZnkNuz*V(3ew=5_^DA=f`)2N;}$z=1o zw5e68yC^Dui{%2R%IHny_O|SC3mVN>V6z0}0ax=E6({;xqoRKJ70ToWEu}_j@E2-$ zyy8rCtn_PD<1hd%QOV2+r>Q!o=YGj6nrKQpX?DVbY2r1-1Lw+!<^=df7hmY``N8To zORSn&Wb5ij6H4_AD^xp^hlHxj9bl43lkU6xI z@4z07iUlg$*rwZK^=nQUim!_+w?x9VQQecN+YKMNu`{#;Uel{6!~o~QOK^gT-fMM@ zUlow&$*wP*Wo{(LGL}$@=wH1900wbaU3`Fv)s5K=de3v`4;bC5vJog>unlal=Yh#|vG7 zn}039MTt33V~4`L^Tb@X(Ii%co2TpOd~k3eYI{dyF&Vhi8+G_Hc#wiB%6i6!YopOb zWI?rvuHmfzy^4Fj^3a-CV}J5EoI?_^%lv?0%dgfV!B8VYbt6;Na^_z@hlqvW{c4D= zrRVwD--E?O;@0nt<`r-Q8K(NK@D>eF%QD^7xAv3k%mb`oRS%Qul|-HVb>n+yin~+e zWRGpZh85~01EYN~NoG+5+;)+Saq8Jg*+z3Dj)i!zweS6!j;uwBD+I`>@}zSxpI{6n zCAR{jjinw~)2Ld6>MqU*C|;4zG7VY#B?;D|uMXbL#$gzMvKD83yW|M1TPZEjUFAma zCnu9EUajVgRF3vd&kwqGD&Q4sIBfDMPcyD{bcBTHf0&s zZ%=)1FZ5nw;oQWA-y}&?2L8RP$=Dnsb8*k80UL5*#|Ow91H{`Gb7gXnGZKZ~tY9_IJ_cZ&fK6Tj(h+)$u2^ zBw*$tl;Rq>cvukDKEbqiHEWiX%Dx1&M!soS#12u5h7b?VQpK1End+=3nMcVw+@Raf zMPGhnvY*a!)Wfsq5pzVmSZ0WJlo(r1ShNe?w_CF1vKVqCirBh4xURc*bOzlp+nGOu z#t&?+5fI3!)ts&?&+~i_IfpNGT+bD#?i3%aCeR(KbR|(h*ZjfG-w(By<2P7}d9x^; z-Pbg$56-|i?Q`o)t)gS}pRQGy3m zvAd6&T1D%0V_;puK@f7w=3(6w)!DcYsk6EsfKbUkv`s-27EKZPZt<#(m6!fyG~KY6 ztQE5Kf382ES2-8JKeD{q4BIU0e7bPIwQ&U9>rcDRpFh69KLog_#HzW*y-SsMD4S>B zu5<~54{j$Los!L5#e$q9?=6C=J?<&b>60YTqVI;#lyas|M;gxr?H=+k#2C69J>Kt1 zia~~Ib_YjpWyzA7*%y`2Cp}%oL_2pv_Q#y_2MQFg0FJxP!fQ8!ITv#Ehr@OhS6_9}t77IVV8Hm?&poo!@Y1az?4rU47tih9Ctt}2@~*c- z@>(k5er{{Jaxdymc(uFa@dy52n?Tbxm1QUbA^;X1gYaKfUN`^_6J7tvp7clP7b@Cc zh+AYsk$4pU(wsHl#DhtN6<^h_b48+P)#|K|uL~ukDfz?sT0c`xWiXoL-HzN8%Y^fr zVG!|PQOgvF1W=l^DVI#de)kb&z!L(1Dn4M~t3WmDja6z4N;w!_8m!fdVLnxx1*}@nr61WF|YHc&UJ!p@n9r^P&tkbOa3L*%U;5*yutXBSDvJH)=Uu>?l+GX3Gg1IjDM&Onx z1v<`emYpIvL*6voT)<_`&5YEZx2l;t!|_aZP`>BelcjpQ3n>5V87Vu?C^qVUT8TITR1WDD*pmhiL^$;brdo`Av((P5+9p-$$#H)hFe*A!f+LmKCgM@LsB+JN7i~%i;B&{(u1NHdrQv! zXbDFPA|vCSu)GjeoTnl&4P<%zSRG%EogA}2sie8TY0XX3m{rvc!)(IwO5-yOsEgA~ zR*wq{{gNz`^0@D5Xu=XtJj`;U4O=&JxrUy<=aa&qj`7>8P?$^W*M2=sil8%eta!+BrS5MEU-M_J_)7y@$j!2R?2%; zdaHi&uY+xMz1-xQALdR-R`#XeWA0@Q)z>gT9g(wTsBv!zB1)M zvl^RsU&i4xcM<#XR0UV2Y1WvSA*0f?j1Eo^mZO0a&dV<9q;?0^H?&w=`N39ZRZBQy z)qN6lo$T{6yl?Ca3Q})$4O@#=SG{yfNeqqxR{Ym4DV78cfkrnenziY<^o+?ikR2`iu@6(s>lC^$=px4)So4+`OkGyt`NqfHto@a zBk$J&D`n~fJewmPHulkHIQRUPJX#Jo`5{1v&?Xm9N)R^V>GY^IiU}0_iJ{~~w9~fc zh%f$JUP`4qrES0Fy$kxRqwW{6DXaQn$*j^<;wFZ>QQIy19ZAaXc675gp}6oocG~SSf*5=aQ+r5? zP7RsY-NnZhn9n?OJ#NRl-Ek0cRk{gw~IJxo{tfpWFH6&Qri5=lVNeA5y2r zhMP5Fg1bOE(azsMYvjBK#hq1JuTJE@n^pZVF3x_{!dObZLs0%iz3&mk<+Y6ow{x1W9Fme6pHaAm`%mX!X`ZOW)1>F?NQToOVCS4KB<b{t^GNN0^hjXE3Agh5SR3yvY&@jle z?VsB7fXu9|1K&7oMuqkJz6A2zQUE`wcdnzw$A@BVvd7?Rmr9%9= zCx~cWxP>PxVlzasjK5SaO>&PR*ikRwrCy8w)F{Ht!EP;#?(>k-^kq!{-9ho?X<1~@ z1T*7fgE2jormU_cW+thN;a7vuUvieUzEAYjbPjil>y4rb_tyEL{?7BM=g5ELj3{LV zxfBJFYoS3bLB zfyK=4?&3D4pJQchC_d)5)mlojW|0kAIy8+bi)FD#>jp{jB9N)tnocG@)>|RPrcGC& z-x;#|eAfYbf=U%97>91vMq)`rH{R5S*;U$nNDy*gz;YpQTGIU8*eLj*aa! zS2`@6O1+~ht*-Dhd#22jTK+{4=vKuv_;dJGSe!{ zY;&JP^T(X>|CK_;2sCseOlK`e&hkqqQoskEUlP;zsA^0jbS4F60<_Q0y%v!+#>LNC zppQf)ok;}ps;1Yi))Y*c13z^G!Td$)Psin&57;t=_ojs+S@QAPY zv2REF@?IU4wjO}r>c4RiKO5dJN~KvE!7-8W2~|zJTq80b{>o$U^)IVgnsx>hS`h10 z3DLP}TPjw*wnZt{pI}`Y^DU$P%vou+SeG=>x8X#=U#WLp3d~Ai zMd!AU>`NkH7uP^8qMa#%WB}_h$K=oR3(>B}iLGEJ5Bl{b8=TPj(L$el9qdF7kIn%v8;)b*%;ROeAB81c^*}k>?PWY!j>tI0=4@Y#u6kL&>3Zn0n;xpB z*D2mjo&Frl6CKO2wPgFiUJy-x(5$8B2M>lMxizjzF_^><-oF9Laf2QDpjmIe?p@un#WZ2Msb^l3{z5-%r ziEj@XYk~|qpz*G;ri5((YJfw^q*UD$_g_Y(e#_H z^wGw?eGuIo>m5gU*8AkDH5nW2inRUiq4vdd=F8AJ_E0w071rLU=A*(E3@8PrjpXqI z3I32l8XU?OMZF3_gRcj;?^luQ)EaX8oBZ-1~|leB)&MGcW}Sj=6*fd zemh&cFv+8|C8K)Qye!PPiQEp%jJ)sM3U|oe&dkN)P$!$(0yR!-YR}T3ubv3&ZpUJk zR&JrKW+V9KAuDE3-*1k*8P-YYyfxzMpCuBU8}6lMK+_GhS?anuVtl1!LeR5FNH(70lDI#7T$431Cei(&!0U9qv5i{zB=N&nW5tDOJ{biPUIX;LaJEG_-HgY|34@VBQrJ@x*P{Ao7 zeKP+GV=iM^n8ZO85nS>|ge=amB#s05#J1M|+~pXvW%C>Wpor z-Adg&w|6FNYM6Ii05xb3fXJzr)>~@n4FzQr(@MdWam$+L+WrcZERnuTR)$ZVHQCW& zRE&!*gM*&rMVtOuJ?K05=443Aq4~AGPQN0!{LdU!<|=Lb&~!~wd3+pQCq8|5Ce3ry0(3hBpR=0%WjI(xMiw8P;kqEX`@A2 zzpu4l8lB=0mGs;xk1f6P3eEy8S#MnF2#BuFmd+E&h7rT2ZE9~eNq4TO^l#tka$fo5 zkDWHk?HSop?Q_d%Q)tOFNV02#VdJ23W6XNah9SDCv_-&}V`KlQ)J$FzvMOM`K%BK* zN8f@|uw<)iOGde4yTxgn&tGTogl0QQ7W$Xyd|#4bcl~=pB*mDbsg`HKRJdhFU7@gP z`KAC(|GCAm>;7Kw{5jc04`|naini-oQxaWg zKPkuOS9Op;wLOEK3XZ`w$Yw-5u6@qdCam(pUmiBEYE4&HZ?5RW)G2ImR95hj0d4hK zooxw>?y0Wkqvzp5oZX@0>C>vobJ`)qQTYYls!r1Ms%Z69(-y+I`vzkPb&bS7v3=S( zvC}-c4`)1r8WxXbAtZO%xrkqQuG)P_zl0|_oo>>9JurLEB0Nr#1@Njo&>lWFS(UDo z;`Z73Fjk~0ltGFvQ#30--VARLvtUrMrbbp5Y@;!a3swLWJ$Fpo+_25#scGownporh{Gr6*bUwfUMv%_a(&{_f;x>RB#GphOft`5xT?&8V6;jb)TolP0c& zCFv4c^PwtT(UGuu=lH~dwfX9C-vLGMIA6w*DWZ(Rtm^X14f{_AQnN&=cX?Pn>Tj@D z)P5|dlPpFOPLy=Iz;;#&CEF=hhl#F*ru!8pHYc)(qZ7lcm!O;%fjYJ2YhIKqV7DbT zf+KVO>XPFPPBdFTDI4jDQzk9zF>BE$DRP*6UD|x7x>-Cvl7>3X+8-ey7}@NS#m*PA zH!z&eOe-t3j>Rmzj-7AjLqk_+u=c}E?A#Iy9Z(C@pY_W+jx^n8Qf`h!f+sGk#y`rB ze@*5Rg`1q&CA&DxBtJ5{%1v|(p1JC?Dp{eeQk%oyAc<9RJ-+X%RDwM5ax#^pQJAvee-Tg7><_0ce>O=3H$xB1fJA;Liwj-}oB~eR&Q;L5Gmru1)OUYxvCZnR+ zgFXV&L(^lowtaAfc_i**WJC=!+jdm8%^)lLFGbD<(H2&jmhey9s%rNhtXz(E?lv#1 zH2n8-h(v$kPA|?_BkWr9KCk??D{*|c|BYWR?d2itRRQ0e!As&1)6Bqbc{rTP#)SFc z)9&u8#A#F9VjVvlLf&Mn)sW-I9CqW5D7c$8YmnpTh>hbJJM3y(Zxa{bI$Fx)U+w5L zjvLK(hET|vDAt^~yXR=eL#Acr25wgBZt=xzDJdr`6)S@xxkLy*L8X$~cIBU=HQUO$ zYAZ;X4BU9q-MDM&t@B3gGLz?ZvAi}Ij{d}M>X#xJ|P$A|T1_RkOeM0(csiurG;dCnL7z8^hTFOQ?a zFv>%0I)gN)&iTz(wBW|PYiizVtG4k{^v5^5K_k4@k0{6YJlsT*fVQ+~yl_-D~I=2vIo zW6r0|usqK^%-Z?1X8qiAJ1(mG@NLo9>pj2sIphL4)qcL|!W@bpa>vJh!edSJ=kve! zE?gQZcv>^;JACTzlSh{P8$-W@>%LV#dhpZS=?}B=E3V+7H+Fab>m$F&w>nhQKFwsl z^ou068$3V&5IB%vL4yYoCM-CBVZ(+197a?)kzz%Q7cpj382<60MvossZsf=i;>ZCd zN1~Ky@+3=@D^t1*Ig@71nGOR8(77`JM4u#Ewp2OP;me{Ak0KQc^ygBiKT9@k3REf8 zm{+j|wW^gXN}^S}x)i(hE7_w_gFujI;8A#*h>+OT8DuH&wz z3tM#QwWmk_P3^UIZ^nr~pPg)acyZjnW3yIXn0fQ9vY-179{Trg=~{6U?`*p=!_1m{ zTD%V0s$H8E$0LS5ojq;ms%NKvd!C|v>(!I9?@Aw}{`CCO-y@Ga_ySZ+J@?A9@4m?R zy9g)EHY4pa{RUf%yW*D0PQmpCMDR5QJ4DPwqBPXdGPLkZQAGe349i3bU);{WDXFZ|$}4Ze%t+2q>v1k2 z!5s27wjOJ3MjsJ$&CC_OY!AdYYs4`t-o(rhLm2ODkwf)Z+Oba>ZM3n?1qZ!LyD$+o zj?4kR)v=kydV zS$jJ)w^I4)j?pa_0u?_j1-g#YS$);COj;rS4Ha2qM}@UKVC&nIMp0XJRmWy;T=v6a zcZ0UqA|d4R$2W1cwZujL8#mcv)m*mFGuMrmT2kudO!FdE-@A;Bn1W z^xg?~4cA+5Wm-6+h8>PL*nU-&xTaVC^>@8K0lv}Vk2QKwVQ%paa!eTm9=6_qffUwV zZHZ;szK%0S3ST?LeAn5H2jp1eo%bXe-hoLpIMEoX+p5tBNw%0ygoPQJk=H%Jh%P{7Cdgi@uHjE#S@e^Z_cp;eQ%x*JXr2q2e))u zt;7DA=2h)BJN3XHN4ayQc?DVb-Fffb_unxRyGDqCFWx4}9lf^bgLSo<>X&1OJyx41 zV?6qm_0*R~S_`g|+2<3N9@{?=-&bDIuZKI|&%t7wv+22}k8|+Po_%%YixTN>`y=-} zPW+oDzrGPparZMH`GhB`mR+xZkfO=_MB_g3Ew5zjs~`Oar?>Fcs(Sm=pyvJuKdE(3 zQs2|w_{6ud)@e+EQ-h$;)aM!grAC4P6x^v&#l8?;5OgP`VCFpdLJFM_eI!KP096Mx z1CDQn{%WE1>L>uWco0aa-`T?^W_oG z=f7?W%8>6}s8&_cQD z;2fv5;AYWtbUN1Y*$x>FH z`F!6I^Qg_fV$Y<0wJ1}+nby5|^p+gNA?g@e#q4R;sE_?8W-HlMj5bs#3yC9LotIdU z#xF~vO=^YykZM|j=1-)+gs+P36w5W^L8{gV-)2+7^CSR$$y*lo8NE`Vrd7xfFDo%zJ3!W$!}c8~e&4!4{x(p&IyNyuP#hc;KNqOEldOZYyWiB!czzA)?s)?k zF&|ym$3M2MPKVqriH#>A?p-iQT?|pdURc5}q_9&(i{(XXc*kPh*V=*^6t1IID$WVsI>NM_9d}V_OQJE1$dUH)mziB|T{;a-We)JA zTMTNM!g$kD_Vt{BZRBv{_}9fI@pzHV>O5O^&&?L`v;TdfUg?$EjpR1CCRb^RPn_Zv zw>a={O>5|y_u)zvvNt~+w~NkphX`+~ys0Qta~IQ?kTx@aOSABkpIWvr?X|aSsQz(U z8hpiF)-A~M&2TWkSWYHKt9%7bxRsZ4wBGLdsPQ!GnIrw++RCp4-K}fq>hRn{A5o)6 ze&?g2O?N8yFw<8}@29tQ%8O552O4%zZ@vEqqnh@M=-%54&H&~w5{5Hud6}D_?in_ z>^h3Aym7{6pC?dwT_^R=HK>edc0YO@AZx+_2Y$VuFaWWbSpvs?1Z2GgQcvu{`H4>>pc!W z_`QGpwS)NXsZYMnagXoD?#)n+{W5R!f(-r+ua}5R@l36Qe2?RRPkbQH^DOWJ$<1|Y z4+A;S%HU0(v@QVoZSFXV?(mJb-p;xr1@F(duF`J!X*|({S_- z@~vdB>+*4g&?6Gt@#WreM=Ed>HF6_4vLg=+*XXalC~>eh>l$RT3EC&^JiKzf2PS&{6F=?iA^77@ckjEiWgf5wQNAil?G+0uxaKmr^D@sTp|@ zDT^y2car*Qk{<_B97|C?cnkHe5*&}R)CiC!jqxH`Xe$%&!6YS3xhHk zvGNK<(%jmT$KbLe6|nyD3Y@0VDbKPgfl?jGQYURPL*~#PvxzDVaxTBJ2eA!G2y>=v8TDv|Fm+@mY|>&V7%CKVDq5tI|J(jR_h#quAV~!4w?F=ilC|_`^EK)-UCBMjzw2IR!S94e}PAN}x z4HvX8KPl^uRLi1Ee||In1g`FIv^SU259v-kL6qx6)Ihtl$V}8mK~WbmbVVDJD_zt* zbAv+f(Lw`LKARLP!EPUG>f(qeGtsU6z|@TN6P_T=vFb7iITReDDGa+0OA|9twUncf zFhPklMY;2B3@b0?W-;ATQhiZEPe&yM{_#3nl1+O+G`G%Wsxiqj)}5QPZwS4NuIrNLd}G zHnI&@YZOk|vrqdIQ~}FSQ!^ZYR8|cU@bofKHPay(l|fg~b|6(rlB_j_jyhLWS~K*} z@-$xMbzbQ;N%(UIdn^u7bw)|mNze5!Q|>+eGZbUc<&0EYe<((81<^#VP17?i`BE#X z6<W)?ul3 z2&DuW7d02N&?vnQI+;~EWp`R1Rv!-$9oIHHI~Outw{%xeb8Xh|*3oqPGy1l+c)PX> zTh~y0@di28a(PpAxiJ0}g_l0__I<0+i8c^^={HDm_p^4C2h}TJg?Dwg5-f9eSBo)!pHwQ3X9PoaW|Y@Q&lYtVcW*7&0o#K~zpe_O z*CIPOXhrupv9oC@_`^Jue=*p2KLvmrIE9DzdzYG79Y3jobYaik^ABz zs9so6_cz#75rS*D00nh@?+}Tv zjWq9gNQdGh?EbRghL#Bo~P;IC`F6*?KkB zT4j);{tY^fTA@35tGb!x zmK5z+#I+fWhgyac{1s){h&TDf8@zh?X*vK$QgwM{e_LVKd%aIv)mU7GXWX`ZSBYV~ zl-Fsv6`Yqh+?gL7W@FpNxwS9sb&AUx!&wr^x%;t~+-m1D#K)X^d;6gU+rP^>hu^!R z@jF^Rqi^T@W1{uR$-Dt6nT4G^%~h0z6IjZ1`l!*{%Kv-E2Ysk~9FYYmwBY`$Td_M= zuUyU-{Z^YCUQ>M2Ila^O^Tp58Zex51XphYv}jeX>iU6m*kIhw?)fEA@V-08V$ z+#OTnv)$%r-s6jW!L1(9bNlSOn;=ad(uoDu7xA%6D_ZYkUTUbF%r8wb7TG(Y)~w zpNJt>r%RvY-9FzX*_iDeV=G7VRd$^XKFebo^pQ&R(S6EWe~AAa^%;KKH(2rk9L_%& zn%Nv(v#i}aTEDkit3tn$>s%WFf52aI`NR6~YCQK{AGcw@{+KJ@3%u{g;xHEZ>(*0o#Fq74gI zV^@X_8$tvM7O&8pRHGtv={IoTt2Ezs-OH12Uc`EtHb%<0ui$}%(~7O^bn@TJY71k2 zsaf!6pQJgPj$7IAXwPRu8%(|Wbk(_a?f%4j8+YFRrAAXuExIyz)|Rq67hmi+sPg1` za|6wNo44@b#*?b9UR0`7-^Wc~U+WwA@YsnHe@AWn{CBhBr=x!lfA&tgb?-iA-h4OF z@YLCYZ_ibC){zGta`-hBS!2$XrdWAx@rPZ668`5Gcmm2dA%72!mtlgeT}WVt>Pgt% zgXx7hA#Jrm+=w6Bc zH3^r8OCDIDl2byq9F|&cdF6mJ0XPzfVm|p%h1|slTW3yc_~e*qat5HAa_zSvlpKQT z-k9e}d1XgN+Gb^!&*|6NkXiQ0V0nBd%4hzSiXQ4HXrZC_rbsnX8rE1lGRbLXX3AOU zkdt=0C!~4(DS&-SBAH}}bq2_0tXxLNsHs*u`Y3+5Hmatbl#U41WENJMDPe2kN!*Aj z(mAZMULor%tk9mh9(lYbiK?piv6?HIul6==m|fCp;H-h}YHPRSGAP}vdb+9=NR$=S zET(hzdS|rvf_mbYc@n#?zOf!Fu(PF7d1`H`a(n2YIf)u$t>%_1D8=aND)FokZ_IJW z9)ApS$dp-}@QMr{=&*yM)wz|MQeHgqsu{1_YP#l;Sud%^4lFFNI7?fwz#n@1Gq}<6 z8#K@?7Y%S(*4jLu!rDG|GR7jpi2muzGsB#+x-fGcscZIPrz3$Pzc#dCKaOqb+1nv~ zaLkNlTXn_d`pjm_#A;14#asVfGs0kf+jQE@mMI#nAQB$#bc7ooA>&EI&2+v?4_vrq z@LgLT-RB01D3h0q8+Yh|%S|=lr-RLI)S1gI^EB*80{@p)6&U!l=rSORlTT-iA-Ap&Y?+K7|3~Zo+%9kZvsVypI zgC5zcmnPYWs)EdWT>w!B{=TMYk9^!qpvoZEy{|zKWqnIw3yJo>mBo)t&Vvx_;I}^c zVXj}t>0jPlC%XL!&v!XP8Uh_?wSvn-W z+Is*PH`Ylpdtp?b7}a>l5L!`^T&toA2{l7%VP%e+vZ0r1!pBSsPiOWUA`tb5zx#Cu zkxx6LE5T^8M?Nr$=W`k)FG)*vSu$f1^>LK%D&_D_+09XsE`yod+%2uiNRx4Kp7gAz zJ@1LneCl&#N`xlpqB+KFIunxD6zCsi=*M;XYdm$lqa7V;MzcQyB}7ntvOxhWq7Jlp9T@A9zEzqzsk}r!LmLFT%IgndDf(+wVC^D ztYaSw+1mcCF0mr?=UCSWP_qWsZX(^R`U)9WJ?d3F!s@AA<=EDndNn4YP2WU&NKxCJ zHejbUUuOv@*2`kHtbkoC6+H`Fn5I&YyZfY_T+>U5K-U1e80I<0cT7O%_YVRa`P z+|F56w=4DS5mlPm-kMNT$K)N5{ul_^Q904{kYn@CQYFVOt4&|JJ;2KwxbTFtp{!U)RGZ4zq_TH^_~~m5Q})k zBpz{*@XKKq`ggw#W3k-ch*NIfG&`c9pIffl+8APTw&D!xubA~^ zr~bIRw{3ENooq80o06|Dr7?1m+n-3H7rKV*Wrcg(SEzy#%#Hi9R^2Gx4C7YGOMbGN zq1@sX=TcH#tm>uKi^w~_c*~xG;&tghsR*raUD}6z z>d=oq@v2#Ez$qGfQH=KTgzNWZDAeyc z<+b#xv^g4NiF=abq|SCL1xxOQpm@o+?zzu@4)maR%-0(S_F+3rPSrROwRu~oeyr7C zz9)aeSAX>PUlqt<7l>9xm3Gkff5g*+0Y+Nw1%QypQ8JfT=y!s{$AW^?J|ehK3r2!V z=zwgfghp5rPUwDZ=vt4XIS`m|QYdUWM~2*GhDgVM%Lj<*H7dqLOpZ7`!W188GjDju zfbBChi^zE}IEckYboeKDhq#0EH;Rx3L+CVnnl~I`b8U)8S6l~wPzQ@Bs9@&DhkEE~ zDQH}p*m{Eab)HCk3uAo82XNa5Yxu^ClIMe|xQ9;CP*~QAtM`huICFo6hwTT2$EJ!B z^?gFNiyO#yR!E3DXMtL%h{bqhmneqo$c~h_h2<4;z(`U>S9VI3gBjR=4tH!2XF2XA zTJM94+?acthCeIsQz1`n1;lNjP?R|TLm(VXhGI> zckcLl3pkRE_>CoJaSteiHAs#(XpY9Hed#ERw#bj%7%?dbkOE1Kng@~`IAuSXhYP8Y zw|F?pb&?4Nd@31{MpZ!V=#)6KpzmSHKD|CB9b zX_jY+mT9S$Ysr>v>6UNV7}*j4b4izVX_t42mwBm|d&!r5>6d>An1LyngGrc$X_$wJ zn2D*Fi^-Uc>6niRnUN`(lS!EhW0{wUnVG4Xo5`7->6xDinxQG0qe+^jX_}{rnyIOp ztI3+J>6)(zo3SaIvq_t^X`8o+o4Ki*yUClqDe0TP37o+xoWn_+#c7zDzNS5Y_rs~SJ?hD8AOxN~}=lUKT{0|HYhs2`sh;#vZ29bu@fG(ud2Q}oNY7^UK z^Xo-OyFM5L7kbao#BTEXqEK`Rwl+`NOssH_XruOC>Vxl zxnxqwhY@&*RR|gt+E_{YcDiU-3YXXFh06M>hoC;w)pEYjI8Zwk;C*8Jh<6_+!hH*qj{$dg>)*Dw$SH(2EJeq{~u$8?& z1`{4c2s7oxnc-xbEQfMr9*ygG=)4KECoN_{XA$j~b0#z&1r<_S_!O!qhe>P7oF{{- zP=Q<{X|>u;-AYa^V@5qY)#kjf^d5Fxi%6nIWE+if%q8nCU6LxrmV8U@r#G65Knb~{HrIi0F- zi-#qM_+{hgvnjTx-j-<$XXh-ZzMu7dxcRI8ExcC0-gMbPreA#Axg{M)e;wwXVankZ zUS`29Sd?Tq!B^XaZtnV8^B;Mph6DJxc}qIWz+D3>ebIb z=&|J7e^KBEWROCZwxfsp9Eqe$(Fw?;gY{5ooQQEf385Mi<`|`mB3)U+I#t3rU6Wug zM`lJ$rI%!gJ}z0JY*?~b(T!7@(`9By%?RV1Cv^A1fk}y}U2&Q8SHYT#fk>j9in91a zo22!F+GK34K_~@=stKEN5ui6E+5j%#j**CLc{uG^+Km{l~QV57Lnv9nd3&Eom1 zI&>-==)#>XTHZ)6k~yfXX=rz5ld+5ISrp#@>*0G$Y%fsA-%b=YeQ|hW43M{p|Mr-PA&*DMN z+`W*hkOjzMk4<*9B_nz((HLfJEUTcoOvHIH=A3Y`R6o47z*--Ck-b3+sJ74s_WZO` z;_kL|Zwp$@V43Mssr52=3;uq`tDipD_a+J(7lG@SYwrgY_VQtlXJ3ubW5UgJC8#}pWBl_9oT_{B>=FC*+^5Ecm!ysW4OmeBac3bEnIF>324) zVGFA#z56wddmgN!{uQq{5E0JLJ;|Gm&*0asmJN$$2Be(}=fbD_0B?p@3th_E7PW85 za40&Y;1>B+LN7WgfniMA7@cRpD01x1U0&FdXccNTZVF@jkWnY-kiotcRtZgmXKH0^E25^aZ%gtYlXv`xHmY}ht?PUKb%NphNX;pnH z#f<(K+592Jte73^i2{X#bguQZ;XJGZc_`eOZWghtEMs<6OF>GSbs(8^?K=u&-eYky zw&^3Hm6$4*RPpwUj-A@MhH2_BrYnaN?1cSZ-;SWE(*6>+vY0ai77p&JU@%v zZ))|4H*|3BoGZJVrgy^}9uj)k>r@I06v7VOjelkI-*OQZnfd%Dw?6tn5}!D<03IN2 z&jjML>i3oCqh|}$sUzFHEW3Jp@q}ACTnl$m$kEj>CJEEw(Qy(h;QJa9KAhMPGuF5M zgla2sn+#&AI9A7$uOK0dI2yaLn@H9!sg~S0z^-_4IL@+;!)i6>!caLR11k%^3;rxM z$MVQBt#DfNxYG3n2c=PNw4)_~WhpJyLX#%fLAg{~z7B`G^ew9*V~S}3+Bk8@5Hm8H zY33UPv}tvg^xei((&?2GJlh>9O{H5)c6u77@9eIJH9O>}0o8pOe029%*D+%^xz+e2 z>tbyT?G?Y*p}K`!rJ4HaUQhaRV7_l-r>GZQPqb5@&W)ixc56YEQqSCuwP@q@*$1E0 zS!-QwInjI`Co6m3{CvOJ+QzFp6eKk@!ePrIF=CXI-4r^#}NjvG8+!y z>ndo}rDmR|0-Wh@Z+p4=yR%T=ghYtOnbZ9=Q-obkYVaGwKuOAm>hZ`HFVG-=aO zh3!iJ$Vw=}mE;Pbr&r)ap&K;~v55M~J9jt%b z_t~q4_LveUYo0e_0Of85R&Fikd|q~63y4rE15_YaRcS^UZx$+{=Y3u!TgCQ0HdjT~ z6o0gYFshVv+W~YGQ-CGMP}(AZx<+q-=|@#wP#f~WcFtt6UTZnB7_X*dfXRi zAIL$#_IkT_XMab7bc1SKWbJ zj7&FpNT-mLHj&}8Do)`VUG)zfxhS$e292)2!)Ub zRfy%1gZOgLl$8@0QPMM&HD;2LrHAlEi|*2tixvc!6?F6GLf5r(TDU8f7mAeFlWWO) zxaN>9CVsGFI7Ru3NeP&ArH$R#P;AJJsuhRwV}qF}iz8=z%=eOHcWIg@IYs1($!A+l zU3*E4oVad`6N7@*mj`HgYKd#lID!oGmeL4`mFO^}l}j~2k;cZDOBtBtg^&~Z zEPVKYJJ^HnC1%UF6pI;)rb&|VRg0uDMxxo3V0L+~ct{FmfSf6bop^?S)`dnvokM7H zMR=F@_f@0VmVarSzDIBC2a~v&o9(Gh{us%U?)F#hMmP%zc4v8$nYWx{c6Q$Aoc?)8 zAUTlCX_Iy)GNJL8+gV+(Hk!d#j$k->nOBo><#V7&ibU9;%$O@aW1c&;q524zAbNKQ zIY}&*Z-xXeg-MH-DWEwvWg*FURhg2GQ;Sq#cq@+q2=1EKke zN;LU2IiygNu{j?pY>&y+Ln0dxua1^IY!EL9l1}y*?-kGXQUaA zSxKVoM@3^eRxzcFnn;>o8j`m5f@Au0;|GvTI-rIokDAG*!-$BAc#j|yKK%GW8{-a= z^i1>EoBRomS{j`)%2jjPBw&jEZ;={m+&Q1v_jpk%O%lndBjTt=nWY$7jXrsreTpHN znp)`8nCPdSXnKwOwWcpxn^mf-9Tk(22&?I4gK&9?$l0X~%9bYBph&82lc;$8hhap4 zpQ8#9gK30fI&E}S5TRC-1 z_^p#VU*ZT@u!y5w2YsXpsJP0W?n<(r)T`(jut(UPf!dg_CYxJjr>kkNo5-n2TC3^! zBD*QAp|r4Q=XoGyoCy9ara)!El&8j#U)gwrq@P|#ak9Aa*34tg|W@hTJ4-C1mXv0FQWmgw&R;#{AOeE|Z zp)jn1lUug%tGs!^tn#;G|Ju25tFA?ie(98}Bq&!`E23DdSTF3sR|de!Nml(kx|8ay z<}06amP0z+!o>%z+Oy1F<)Jch4*#m0M+vUq85Q2d$}DvB2h$VAM+qh!KD zC%sBc{>k^4#?WYFdK+*{OQ|6xx(&N>Z5vRbD!$Dmxmv7C+jF|cVy6H2pka$~AUMVM z%A0n%f;P)&U4^4LTe@cKYlTcl5{JZf;>aK~#^&|BL%g+IaiuC8!BdA{pq#_r z#4uW>hFi6$oRB`e%H7e*Rw%@!simn3&gN{(&A?X5s>qm1%gM~mQJcKLY`4N3!BD!clD`tp)1En`?(4`Ys9*%?dYffWC{5LoQiWFBcWKJd#7Ce9R9qRE zb_QC_+(|f=Gh~Fh zB`nqUOmP>yh-aMA8>~TgMa*bTWgB?a6TK{2{l=UD)?j&KG>wuUo0&DdnOJS0a*2ki z`@;sE*^~X(A0@?UE3<_Czd(t}u?xs$Ew=pw*(qp$p?Qqy%!+L_vQ;hIefF@#y188G zt$y)PbFFpL3wB$^*3eqUdR)E=msMLQu54G_<-FB9dBQ>()O?1XbQ>B+T^%~2$K%$mCl+PuXIO|MSHr7KI; zb_V$DeU}Ec-Olz@)yEwOnUv$i((&x0JUE@I17oDBd{CB2*d`~4# z-^h~52wtM>dg3|m;TWBarj^p+?Xnt<xe^a=9T3lPzB1c*){t zNK)bb;-$>HS~-o$Ylq=zVPeFm%r=3_X>9^m! zM)%z6->wL>UZ*!|*WZnEIG(EPR^_Xm+x*?3gUyg;u+H*JxvN&H|BQ|ejNPhxc7)ll z0iBqXjoVoqeb?N127WbTO@}w?K=%&rnnU9;NAUOU4;+}^)#8_64x|!A(&1@^y6)Sh zUGEev)XO5#=d6BWPR)MCr=jTZTTYpLk=fQhZD=%V*4MFg{yye9rqjBO@AZ7N{xsuO z$>;m}@(+*nLQslH7Rfgr>`foxXUoeIpWg5W^$CsPC!fL>e&9>L=PeiDH~qTNj&Znn z(?!Ky>ptR9cJWHt^Hr~bREww;-}S2P@f5!C^9=IwUdk#@zG+`s zE#Kr>&tqGEkh$&mpZw6A`rBBzsec&lLmi{+9`OhsF;6bB?obSR+ywR)2QtW*0D1P$(crX3YYp+@;V?PS%onO+b@BUeF z$lzW06`4T_MlZvlhtS4yLxlR!7tMi5;?20L)8pp z6$hT_+Fo5e-(TJ8M;%g#=qdV+NFFX_Ln(nnY4f&-IcHGeHrNFPMKTkXTsEK4Y4tiL zb#udU4=ibA-^)q-o?MrIjo!pARIZRDpeAw@WfTz$wU&N>UK}G!mI4L;HUlv_Gb^ni z92qH<-g*ohM;|d&Fa7GiMC3lQLg{i@HlFOvViyTlDRoZ=qu^vUx01V#E7gKOn^#(> zfF*-L1!s$VgCDj{Bc7}0A4 zfCR}-Yf2H}RfS2xHq2O&Yuc;SQci+M!6lQQFx|SVHFWAul-P*DYgf}3tX+-Of?3!W z@L;%zq0+r(jO|^2Xw>4ZY&9#;tyEE`R=t{aYuAe3b?sHR=iaegGfhM*TkKz^e?%7v z&J6hMR%7q|h^+e&K)HJ<3x(S;u>NEtrxUk?{1~)G=7(0|8a>@}b)3Y74?CkiXW#PT z1FngVp1k?rgw?-Jqts-5m~nJ*N-CKE)UQDO?(@&8|M0VeH{2k+4Ik*-`%arhAY4#7 z3qcv~t%)3R>7Uw$OW`Y6p0n$`4}DXRz9F!SMn;|Xh{-+8v}>`W;_UM&FyNMYtUVIb z3rC>+ydzLU7g+?yz$c-MQpzb`izlDiu2fM%j3^}X#Iqjikvh8O8xucOB#f&b#e!oxZW`MQ*GpTk4`qc zSIKP=;DXLI6rDu=VxHtMo`$};x~c^(8EbibMyW$djx1H}f7enuu6}*arRkgwvu*L+ ztPTn*7uV)FYnm}ljbm~)@7mJHClCFvb;f>~ZO>r#bGf#Kw(jp);mK=M)muHgk#O_C zH)P4Xpj*1=o1>~)5a%S$#D*L<@CUVjz3+WPQD{%YVo@80|G!KXu1 zM{Q=5v6bvrUT|uWhWpTQ-L9(kQ!xXyY=WIW<(7r5@WUbygPdOccEROEVueOr z-~gRRzy($hLlj(5;qq5R3z;Zs%3&Z7PnbqEuF*g)Oq8?EcEm8EZgr7bTdUs4LB{Nf za}}(j0Tpz|{-LoCOQc8R6lbTP&9OLXyPzPGgGD`pCwjtD+CKa!Mk$&Jl1;1L886ow z6%sO2_+w)^M#)HxO=)C|!lUlsb-%$m(vF0SCEDNzyZ(I9R-?2}Cq=fz33d^d(gGxi z>R2scYNlo@gIXz~b;;qS(wC!IrX?TPyilH!o89zg2!$EM$hlFNWo;bC; z{quS#n+CL`)*K+((t%XuXfCIy%h&-dJaNS3Hsi)abCQOl)|w|AbN45Gc#xkp4QS55 zmQS1hFL+uk*hi1YIx_YRdN*Py5BKG`?H$u2Ej?f{0}`Ld)xL! zQ%MhrHoCo(srC|`P($O&`VCYk8@;DKNeNf;>=cD%&FNVy!$yJ9#&})=s$bC>&4XQ& zAxKT!Qa9ODycTVH4uslQp~%yb+HZB`R25ekFYb}(}*{!Qtc zTGiTguT?EyFRggP9knk+8qF75jd)a+U~`%{v}=M8dOHziF|tf7(MeqE#;YDwNwIw` zq*{igBjJ`go8{wJWvR+QopYl?N>*H>sKe8O7b<@wZ-6w3Qlqd?l)g z`p|_Hb*RnMXe`k;pL7*eg#HWCx-yJJ`~A?sK9-{N=y}}`lD2{3OK}_iy6Ztb1j()} zE^lJBJ0>%GIg%7U?vFfL<&6v(f<>+^ka;Uj<1sbOvs`GN1;JF&)vainYX8zE?aeI_v&%xQ~|A<3M0FLXqtg)U<-9 zQ`KDU{nYoVm1C}yZ6ab$@8pwbt^wzwRldIlVdMV?$tGM(r+ zD>~Ncz3<-2aObo%Ru|tYbiMDrouVdrdh1E${gs=t1Z?-YQ$F(oM>nNH3cNUj>%Uy@ z7EgZfI*6^v?_k9VO~yu%vSi!b59^P%yjHN>>mK6o;+CYv%{wmm%wD8pe)M3JeH-WN zEh`0;h(HqjMP2pEsMBoGqF4M78Q*#e1V7PTZoEM&Z~g+_OE^88hpZ+2k^NrfIrNmS zz3bG4<3t)7SPjnT+zIFUFY}%K?FXmyi?6fC7k=|l`Zo0|t!T(l<()bZtw~(^5BNkOs@eC$n}P2<-`YZ0!^uUtM)cd|KKPAd(Z-}Z|wlk zWz?_Y+^-3pP+s8g0ey)1P7Xs*$lPr3?c7eZ1de`4=ivP9(jaVxMz8}v10RsE-*$t} z^sEMRX^yUr3x}q-tdIzAjrBh6;0SFbjtur@{*D8W4ckWN2@#|WJ#c_u_BDb_@?S2yys}0a5M{wGb5N zP3H`50_*Ik5{&IAjj(V;?V^yvVi6Hx%B$4t?=~(DlTa4NFf*)cka$iA#SpA$aOmdB z%ofj|(oNdd>;v1#5N&Yq#P6}7%Ep-N`;5_W;?4lcPREiX4ed@kbny%y%-G~mqz3Kz zRBafi&lMqYVpOjSwGO*}FCXc{)RM6(UQZMeQWSUZ!#aws1W*aF?-9Xo1rrbTbSJ6U zF%O+B47X_fb}Jqa(Xs0BeZURsdhGlDD)RA4X2ax%)GpBxDRIQWYa#z|jZRVuy^jeK zvY!r(Ayv&G?{FQ32lOOTeJHXT!4OV}@}o9V*vcr#d=R4g;rK)n220}Uq%a5_Fe#Jh z$(Zsbn<)SpF;qHnC(V-mXj0C+aRKen?1~a8QBsRgar@YCoW^nMzN*!}%Sntf*+}bi z-V*uDu;fGv_^{G?S`jf%M!FzpDh=$WfN(J3K-+H5d(u+7%7P8kku4D^GWD_N7Sq;P zk;stI{IX8LnuqQ>=pKBM9@~oqRZuIh1K6I9{z}utG%_Y_(`*v%EPWFp(=wjS(JwU9 z4b$%51kvXF5j90K^-k|N6{!9P87Oe}We@L>QaUbYhSLg3<}|PE8tZM{>hCa>a5<|> z1mh>_Zd0<1u06{#N~|%~017yBPp^7%ANg^f8mGU!Gh&j{lH$J0TQBr%mjT^Fo8p(~{1cDl_52G(W47>S&Z4 z|Fa*b(DQf#hF&n8fc_LjNwhG#bba=cJzo(P{p+w&a0P9YP9-o<)y}rCY(@>x7detH zhm<;f6eqrvMK!evT~t0Xv>j=WYHqSH7n4anG%VHAE_Jd`l?YBtG4rHvQ1cU^Ks7c` z6#o$PQjrVqSTR6(^&2M<)y~My$j>=n6+vSaKUXQUVh?*pbyFdeK-kPNleE_uRZ`(g zRIM!ezA}_R4k#;eC1dIfanTt&6-qxfnQF*cx3%=Tb@|FKG4yjORh7YZkhglZR%NU` zsP$I$3a6m;EYG#(j3H4^P)L`3{T@6SjH56j` zR2awe!phStQFUUisbXdIVqp`_SSoTlmSp+m0ACJYjnz#u@LIbvHDU7dILYT^Hf$;L zKDAR>owP7N~Kn2&| za<%MuEmH5zSC18MZB}m)mT%>>%s90^zezIhu3sS+T~ksrH#SN|)G&>9t?G1bKh!nj z*4<18IuA8lA9Qo~mh}3TDz6C2f`>T>A2eihf_8eRHz-E&*VfRJbOEAQC zc?D4Zb2pdfI(K$KciY5GOnDb|b=LC|5nuSnbcHRXnw2n}7jmH&E^$&419v*@i#wyn zDOYjN#ty?;L2vdmWz7wIowj(#S8vN#a?iJOJ}pgEOn{LTSNRo5kGFqicXGAIdpYqh z!<9ofaX0rAbvraZEN-;$xAz(Zfnk?EWyw@Fr6^o8aUdsWrKSoP$H_*9AaWpj97%T|GLHgaP# zgfV6-aTrQ@7i)w+xP%Awgb(xf zVzNwj`cFV9H2uhp2{Cxe+&FrjI1=UfdFfbrT~?8;xRGNAXAYTs6M2e}cU|ca&+_43 z?`$Bdwg}4tM`bv8?JX)TvT+TEfnS%4yV!)OxQJ#b6ZI)>S#_N(7-+?GLZj@4brxn* znTw_Phgms+9C>T77`+fxe?1tHS=K_ERfRY6S2CG1br6O}d9cthc)9irixm}h`EHrH zmv{J=pOk}<_>|n1HD=ITz1H@;nNB)*nBne=C3>PMx}qz(0000xxs6P{ diff --git a/deal.II/examples/step-9/step-9.cc b/deal.II/examples/step-9/step-9.cc index bedbda9b03..32da47f75f 100644 --- a/deal.II/examples/step-9/step-9.cc +++ b/deal.II/examples/step-9/step-9.cc @@ -1924,9 +1924,19 @@ GradientEstimation::estimate_interval (const DoFHandler &dof_handler, fe_midpoint_value.get_function_values (solution, this_midpoint_value); - // Now loop over all active - // neighbors and collect the - // data we need. + // Now loop over all active neighbors + // and collect the data we + // need. Allocate a vector just like + // ``this_midpoint_value'' which we + // will use to store the value of the + // solution in the midpoint of the + // neighbor cell. We allocate it here + // already, since that way we don't + // have to allocate memory repeatedly + // in each iteration of this inner loop + // (memory allocation is a rather + // expensive operation): + std::vector neighbor_midpoint_value(1); typename std::vector::active_cell_iterator>::const_iterator neighbor_ptr = active_neighbors.begin(); for (; neighbor_ptr!=active_neighbors.end(); ++neighbor_ptr) @@ -1943,15 +1953,15 @@ GradientEstimation::estimate_interval (const DoFHandler &dof_handler, // the value of the finite // element function // thereon. Note that for - // these information we + // this information we // have to reinitialize the // ``FEValues'' object for // the neighbor cell. fe_midpoint_value.reinit (neighbor); const Point neighbor_center = fe_midpoint_value.quadrature_point(0); - std::vector neighbor_midpoint_value(1); - fe_midpoint_value.get_function_values (solution, this_midpoint_value); + fe_midpoint_value.get_function_values (solution, + neighbor_midpoint_value); // Compute the vector ``y'' // connecting the centers -- 2.39.5