A conservative fully adaptive multiresolution algorithm for parabolic

and differences between adaptive mesh refinement and adaptive wavelet ap ...... by multiresolution ( circles ) at t = 0 5 for the viscous Burgers equation R = 1000 ...
878KB taille 1 téléchargements 665 vues
 

       !    

"# $ $ $%& ' (! ) * +,- .02 / 13

457698 67:=>?A@ BCBD:E5qY%=-"!$*,r%=(s-"t("*,'>-4!>*,u%QvUt*[,3RS"&(+*,'>[&5w!$r':%+J,x("U("Z-4(*/c-" #&*,P0-4%=-4!>*,2%'5*,%!$y]^&z("!$JP;(+-4%9-{2fJ,"!>5(&d | ("!7-T%U'>k-4**,'}-"*s5{2M%5%0-4!7[,I5!>("#"-"!$q&%=-4!>*,x("#T(k%="ob'7-4!$+&(+*,'$-"!$*, 4- #T!$;* ]w%sZ~#!$[,tP-"OR *O5(i2%Y[,kU&K!>-"(/L*,1("*,'7[!$Jy'$%"Jc(+#Y%'>b0"*UOR '$(dYƒ1!>(+-"*,"!>#Y%'>'>Q,V%52%=0-4!7[,aJ,+!$5FP-4*O5( '>!$W=D„@'7-4!7R_… [,' | 5%0-4!7[,a &#T!$;J%`("Z-D*/'$*#Y%'>'>Q "Z{2&5@J,"!>5()]^"f(+-4&0@J,4%5!$&9-"(p*!$J,n-""#&%=-4!>*,s&+"*,+(1%="F/L*,58d˜ƒ1*])R [&&VO-4F52%N-T%k#*,0"&(+("!>*,n4%=-"`!$(X!$J,]^&+`-4`(+*,'$O-4!$*y!$(X%'$P*,(+-X#&*,(+-T%=QV%52%=0-4!7[,1]X%Y[,'$-)Z-4*5(g-"*b("*'>[,iADCFED(X2%Y[,`U&n5[&'$*05™d š *,y%u*[&+[!>]›%=5r%œ#&'$%("(+!>{2#&%=-4!>*,*/-4K5!>6&+&9-sP-4*O5(]Xx+/L&MdžJd 2 -4*CF%P&Ÿˆ¡ ¢=‹BVh£X*,uˆ¡  ”N‹BV š2¤*,'$!$#Tt%=5œ¥#T&!>5&Mˆž”“=‹Bd™P*-"!>[=%=-"!$*,-4* ("x]X%Y[,&'>-P-4*O5(-4*.#&*,(+-4+#-5™d)š}!>"(3-YVg-" ("#&%'$mR_("0%#&g"&0"&(+&9-T%N-4!$*k*//L#-"!$*,(hml!$U!7-4!$J1=djJd '>*O#&%'$'7Qi(3-4&0bJ,"%5!$9-4( *, U*52%3QF':%YQ&"(V!†džd&/L#-"!$*,(™]^*,("g‡X&(+*["&J':%+!>-_Q1!$( '$%"J,8-4%!>-4( ¥*,U*'$[ "J,':%="!>-_QV2!$(1~P#&!>&9-`(+!$J*,'>!$&%p%00+*Yl!>M%N-4!$*™V2!†džd˜U9Qs-4+&("*,'$5!>J-" ])%Y[&'$Z-I#&*~P#!$&9-"(n?BC`Z¦D*,"zˆ¡ ‰=‹BV}CF%ˆ¡ ¢=‹:GZdhy+&(+'>-c!$(-42%N-k*,'>Qx/L] #&*~P#&!>&9-4(X%"1&#&(+(4%3Q-"*k"0"("*,/L*,g%cJ!>[,y%##&"%#Qd2¥O#&*,5™V %c'$%"J,#':%(+(g*=/ 5!>6&+&9-4!$%'˜%5!$9-4J,4%='2*,04%=-"*,"(H2%Y[,i%c(+02%"(+i+&0+&("9-T%=-"!$*, !$Ÿ%§]X%Y[,&'>-nU%("!>(K?‡XQ'$W ZŽk ˆj¨ ‹©Go%5#&%vU0"#&*,5!>-4!>*,5vU9Qª("!$P0'> "("#Y%='$!$Jn?†«,%N6™%+5xˆž”¬=‹BVCF%P&y%=5n­F*-"tˆ¡ ¬=‹:GZdš2*,gZ[,*,'>-4!>*,2%3Qo0"*U'$(V ])%Y[&'$Z-p("#T(1*6&)-"0*("("!>U!$'>!>-_Qo-4*P%5%0-)-"U2%("!>(^%-"*,M%N-4!$#&%'$'7Qo!$y-"!$P U9Qs("!$P0'7Qn(+]^!7-4#T!>J*]X%Y[,&'>-i#&*~P#&!>&9-4(1!$n-4f!$J,5&@%5wš˜%"J,ŸˆŠ “ N‹:GZd %U*[, ("#T(M%=!$'7Qo("i]X%Y[,'$-"(^%(X-""!:%='%= 5 =*,H-4(+-)/L#-4!>*,()!>yA -4+*[< R f%='$&+W(k("**] 5!>6&"9-I%5 #&%œUy("&.!$§-4M("0!>"!>-c*/ ƒi%3-4  j(^0!>*,&"!>JI]X*"Wtˆž”,¨OV•”¢=‹Bd˜¥-4%+-"!$JI0*,!$9-^!>(^%I{!>-"[*,'$`(+#T&P/L*, 9QO0&+U*'$!$#c#&*,("&3[=%=-4!>*,x':%Y]^(`*,t%y+&J,':%1J,"!>5™d™¥U(";"&(+*,'$O-4!$*x%2%'7Q("!$(i!$(`(+&5t-4*%Y[,*,!>5xZlO0&(+!>[,o\2Olx#&*,P0-4%=-4!>*,(i!>t("P*O*-" "J,!$*,(&V}{2+(+-I]^!7-4*,O-o"&5#&!>Jt*,3Q‚+&;"&P&9-"(&VadžJx/L*,y C 9Q0"U*,'$!># #&*,("&3[=%=-4!>*, '$%Y]^(P?Bƒi%=+-4ˆž”,¨N‹:GZVa YC #*,("+[=%=-"!$*,t':%Y]^(F]^!>-"x[O!>("#*,("!7-_QŸ?B‡X!>2%"! ˆž¢N‹©GmV8”=C 9Q0"U*,'$!>#f#*,(+&+[=%=-"!$*,(1'$%]^(k?B‡g!$2%+!h%5zƒi%3-4&œˆ  ‹©GmV8”=C #&*0+&(+(3R !$U'>fED'$&1&;*, ':%Y]^(k]^!7-4œ#+[!$'>!$&%f02%=-4#T&(?BCF%= ZŽ ˆ¡   N‹:Gc%5œ(3-4"#-4"&5ªP&("&( ? | UJ,4%='$'H%5 ƒ1%+-"&rˆ¡ &‹BV—£X*& ZŽ ˆ€ Y‘=‹©Gmd | /L'$'7Qt%52%0O-4!>[[&"(+!$*,8V}(+-"!$'>'a!$ -4z#&*,9-4ml-s*/9QO0&+U*'$!$#@#&*,("&3[=%=-4!>*,':%Y]^(VX2%(U&r5Z[,&'>*,0&5-4*‚"&5#& %'$(+*zP&P*,+Qt+&;"&P&9-"(s? F*-"-"("#T'>!$#TOR „Ÿ ¤ '$'>&%5.„Ÿ ¤ '$'>&ˆž”=‘‹BV—£X*,& ZŽ ˆ€ ¨N‹©Gmd™"]^!7-4x-4P(+*,'$O-4!$*K!$(`+&0+&(+&9-45 %5t#&*0O-4&5t*,t%n5Q2%=!>#Y%'>'>Q [*,'>[!$Jy%-"*,%=-4!>#Y%'>'>Q%52%0O-4&5@J,"!>5™d˜C`!>6&+&9-1(+-4"%=-4J,!$(1%[kU&K0"*0*,(+&5 -4*[=%'>2%=-")-4\2Olb]^!>-4*,-D";]w*,n%52%=0-4!7[,`b'>-4!>"("*,'>-4!>*,I-4&#T!$;*,u':%Y]^(&VH]Xx+/Ly-4*.„ª ¤ '>'$xˆŠ“9”N‹BdXš2*,yP*,+z5Z-T%!>'$(M*u("!$P!$'$%"!7-4!>&( %5 5!>6&+&#&(fUZ-_]X&.%52%=0-4!7[,P("‚+{2&P&9-b%=5œ%5%0-4!7[,I]X%[&'>-I%=0OR 0"*,%#T&(V2]Xf"Z/L&-4+Y%5-4*£X*,œˆ¡ “=‹Bd e_‚-4M#&""9-I02%0&f]X0"("&9-P%s/L'$'7Qt%52%0O-4!>[Mb'>-4!>"("*,'>-4!>*,x("#T/L*, ("*,'7[!$Jc02%4%=U*,'>!$#1ADCFED(X!$y*,VO-_]g*o%5y-""`("0%#&F5!>("!>*,(D]^!>-"y5!>6&+&9-_Q0(F*/DU*52%3Qz#*,5!7-4!>*,(& d r!7-4K+&("0&#Z-F-4*y0"[!>*,(i]g*,"W˜V]Xbml-4&5t-" %'$J*,"!7-4P(5[&'$*05œ/L*,k9Q0"U*,'$!>#;*,(sˆ¡  ¨N‹H/L*,f-4y#&%("M*/^02%"%U*'$!$# *,(&VX-"*,J,-"&y]^!7-4v-"x#&*""("0*,5!>J &+"*,y%2%'7QO(+!$(d | (y]gt%+K#&*,#&&+&5 ]^!>-"x0"*,U'$&P(`%="!$(+!$JP/L"*, 09Q("!>#Y%'€R_#T!>#Y%'#&*,9-4ml-YV™#*,("+[=%=-"!$*,x*/D0[,`\Ols[=%'>2%=-"!$*,™VO!$J*,!$Jb%5s*,-"J,*,!$Jb\2OlO&(VU*-"@%=0OR 0"*YlO!$%=-"&5.%=-I!>9-4&3/%#&(b/L"*, #&'$'€R %Y[,4%J,[%='$(b*/p5!>6&+&9-o'$Z[,'$(&V}%"n*&#&(+(4%+!$'>QU2%'$%#&5™dO&+/L*,+i]XF5[!$(+%k]r/L*,+b'$%=-4!>*,P/L*)%52%=0-4!7[,i\2l #&*,P0-4%=-4!>*,@%=-^-4c!$9-4+/%#&(iU-_]g&K-_]X*M5!>6&+&9-i'>[&'$(pU&!$J(3-4"!>#-"'>Qs#&*OR ("+[=%=-4!7[,=d p02%0&a!$(D*,+J9%!>q&5%=(D/L*'$'$*]^ ( !$o(+&#-"!$*,M”OV9%J,&"%'{2!7-4)[,*'$PZ-4*5 /L*,h#&*,("&3[=%=-4!>*,c'$%Y]^( !>( 0"("Qn%52%0-"!>[kb'7-4!>"&(+*,'$O-4!$*(+#T&Pb!>(15("#"!$U&58V8!$#&'$5!$J% ])%YQ-4*p&("+X#&*("+[=%=-4!7[!>-_Q*/O-4D\2Olc#&*0O-T%=-"!$*,(dYe_k(+&#Z-4!$*k‘VY]gX5("#"!$U -4X%'$J*,"!7-4w("5k/L*,h-4g"!$#&%'*,™dC`Z-T%!>'$(h#&*#&"!>J1-"X52%N-T% (+-""#Z-4+F%5P-41!$P0'$*,P%"1%'>("*kJ!>[,™d52%=-4!>*,œ*/1-4@%52%=0-4!7[,sZ-4*5™V—-"&(+-3R_0+*,U'>&P(I'$!$W= #&*,9[&#-"!$*,R_5!768(+!$*,z*,1[!$(+#&*,(‡X"J&"(`%"b(3-45!$&58d  b"0*,3-`J9%=!$(1!$t£)A  -4!>p%5M*,3Q,V'†dš!$2%'>'>QV]g#&*,#'$5c%=5@0"&(+&9-1("*,P0"("0&#Z-4!7[,&(^/L*,)/L-""F]g*,"W˜d



@ *,(&d e_y-4F0"&(+&9-^02%0&VO]gF"(+-4+!$#Z-*,+("&'7[,(^-"*k-4F#Y%(+i]^&+F-4`5!>6("!7[,i\2l !$(5Z{25zU9Q%IJ,4%5!$&9-^*,0&"%=-4*&V2%(+("P!$JP#*,(3-T%9-p5!>6(+!>[!>-_N Q MPO  V!†džd

EM? @ ,&A @ Q G GSR}? @  G TUMJA @WV 6&+&9-a-4(+-3R_#Y%=("&(H0+&("9-4&5!$I-4"!$#&%'+&("'>-4(a(+&#Z-4!$*™V]^!>#TPQ!$'$5I(+!$IR 0'$k*5&'>(i/L*,i[!>("#&*(F\2!$5x5OQO%!>#&(%5§#*,bU(+-4!>*,x0+*,U'$(d™š2*,i-"n C #&*,9[&#-"!$*,R_5!768(+!$*,s;%=-4!>*,™V]gf%[y?YXZO  G

R}? @ GJG[X @ I^? @ GJG  e_-"#Y%(+*/}-"o C

[O!>("#*,(1‡X"J&"(^;%=-4!>*,™V]gcJ-

‘

@

R}? @ GJG I^? @ GJG

” 

%5s/L*,)-4f+Y%#Z-4!>*,OR_5!>6("!>*,y&;*,™V}]g@5&*-"@U9Q@? @ , A @ G G TZA CJE? @ ,&A @ G  % ! Tc!  #  TFŽ  d2-+-4&!$y-"F/L*,"

6 I^? @

? @ G@? @ ,&A @ G V ? .  

Gc-"

?B”,G

  #4ZŽYN    ŽY"!$#

*o5!>("#"-"!$qM?B”,GZV]Xf(+c%I#':%(+("!$#&%'{!>-"`[,*,'>`/L*"b'$%=-4!>*,M!>n-4F(3-T%5%"5 #&*,("&3[=%=-4!7[,k/L*,+@d˜e_z-"kJ&4%'#Y%(+V8'>-1(i#*,(+!$5i-4c#&*,P0-4%=-4!>*,2%' 5*=R M%=!$ 2 !$§5!>("!>*,  @]^!>-"‚%‚%"U!>-4"%+Qz("2%=0=V}%5 '>-f(f02%+-"!>-"!$*,t!>-!>9-4* #&'$'$(i'? 2LG ™V UG"![&;%9-4!>-_Q - *, 2}%N-^!$(3-T%9- .mV

- + 3?.+G G

 

+? *-,/.+G3* . 2 . / 021

. 2 . G54 0 1 *K!>()-4F[,*'$PF*/ -4#&'$'Bde_9-4&J4%=-"!$J@?†”GX*, 2 Q!>&'$5(

]^&+

? @ ? *-,/.+G6*NG /  / ? .

021

021

?

@ +? *-, .+G, A @ +? *-, .+G"G7*

!†džd

? @+  +? .+G G 8  +  ?.+G ? .

?“9G

| 00'>Q!$JI-"5!>[&"J&#f-"&*,+&V]gcJ,Z-

9+ _?+.+G G T

 

. 2 . / : 021

Eˆ @ ? *-,/.+G, A @ +? *-,/.+G_‹JCQP!7/ V/L*,g-_]X*P% 5 "% R # !$5y-4*@?B“9Gmd



 M8!2Ž_ +Ž'"!

C`y-4*x-"s%52%0O-4!>[s("02%#s5!$(+#&"Z-4!>qY%=-"!$*,™Vh-4J"!$5 !>(I#T2%J!$Jx!>.-"!$PV}%5 -4"Z/L*,"i]gi{2+(+-5!$("#"Z-4!$qi!>-4!>1%5-4n!$M(+02%#d2ƒp&"1]X`("F%MZlO0'$!>#&!7("#&*,5R_*,+5&1%#&#4%N-4cŒpJ,ZR_­O-"-T%@?BŒp­f”,G)("#T=dC`&*-"!$JoU9 Q .X-4f-"!$P (+-"&0K%=5zU9Q @ +  G @ + _+? . GZV]^" . G .mV2-"kŒp­f”o(+#T&Pf("5K&+c%(p-" /L*,"

  

@+



@ +  

G@ +  6 .  +     @ +  " 6 @ +   G ”





6 .  +  

@ + -4[,#-"*,`?@ +  G C &*=-4!$JFU*,t*/X-"o*,0&4%N-4*,  + !>(  5&(+#&+!$U&5y!>-41/L*,'>'$*]^!$JF(+U(+&#-"!$*,8VO-4i(+-4%U!>'$!>-_Qk#&*,5!>-"!$*,(HU&!$Jc-4+Y%=-"&5!$ ("U("&#Z-4!>*,x”Od•¨Od

! "

M&"$# &%

# *,("!>5&o%{lO&5 -4!> . %58V}!$œ-"M/L*,'>'$*]^!$JV -4y(+0&"(+#&+!$0-'!>( *,P!>-"-"&5t[&3Q]^&"=d—š2*,`-"@ YC & # %("=V2g!$(f%s("J,9-yˆ * (  ,/*   ‹a]^!>-4§(+-"&0

‚M*]

("!>q&) *

 G *   T * * (  

9+  G T

 

* ,  +

E+  



d˜Ea;d?L‘[t%5u5!768("!>[@-"&"P(y%="x%00"*YlO!$%=-45v5!768"[,P02%3-YVhŒ^*O j(f%00"*YlO!$%=-4I(+*,'$O-4!$*K-4* -4Œ^!$&%x0"*,U'>& J!>[, -"o'$/©-?:R3GF%5§"!>J,9-?6Gi[=%'$(*/ @ d eS-4((+#Y%'$% [,"("!>*,!$(J,!7[,&U9Q

R  ?@ ( ,@

 

G

 G

 ”

(

R—? @

G

6 R—? @  G)T .  ? @ ( , @

 G .?@ 

T @

(

G 

?‰9G

]^&+

 ?@ ( ,@

 G

R}? @  G)TUR—? @ @  T @ (

G

RL? @ (

y'>/©-I%=5œ"!>J,9-k-"&"P(&V

@ + (  @ + 

  



G @ + H6  

G @ +   6

@ + ( 

 % 5  

!>Jx%

?3  

G

  

]^&+ !>(-" !   '$!>!7-4&V]^!>#Tœ#T**,("(k-4MP!$!>M%'("'>*,0PU-_]g&& -4'>/©-p%5s+!$J,9-^("!>5&(V2!†džd



 .  . .  .  ,ZG>G !7/ 

i !7/ .  . O .  . ?



x("*"#K-4&+ !>(y%00"*YlO!$%=-45U9Q I=+  I^? @ + LGmd)š2*,M%‚J,&&4%='`*,ORS'$!>Y% ("*,"#&-4"V-4!$(#T*,!>#&FQ!$'$5(^%I("#&*,5OR_*,+5&1%##&"%#Qd



E}lO-"&(+!$*,@-4*y!>J,15!>("!>*,!$x£)%=+-4("!$%KJ,*,Z-4+!$&(^!>(i0&+/L*"P&5@-""*J, %k-4("*,0"*5#-%00+*9%#T™d2š2*,g-4”=C #&%("=

V 2  !$()%b+&#Z-T%J'$i]^!7-4s%k[,*'$P */("!$q . 2  . G *    d˜Ea;d ?“9GX#Y%=@UF]^"!>-+-4@%=(

? @ +   +  ? .+G ?.+G G 9  ? .

?3   G

]^&+  

 +   G T 9

*

 +

E+  



T E +  ( 





 

T -





f(4%PF&'>&0!>05k]^!>-"M%`[,*'$P */("!$q . 2   . G *     =d˜ƒ1#&]XJ,Z-

? @ +    ?.+G G 8  +   =?+.+G ? .

?3 Y“9G

]^&+

 +   9

 

G T 

T

* ,  +  

  +

E+  



E +  



  

T E +  ( 



 

T E +   





T -

 



+ 

E +  

 



T E +   

(  

-

6 I=+    V

(  -

\2OlO&(1%"F!>n-4!>(#Y%(+

E+  E +  





 









G R   @ +  (  (

G R   @ +  





 

, @ +  





 

, @ +  



 

 

@ +  1T @ +  T M  *   



 

 T M ¬

@ +  





 

1T @ +   







 







?3 &‘^(+#T&P`*,M-"^{2&(3-J,+!$5™dš2*,a-41'$!>Y%D#*,9[,&#Z-4!>*,OR 5!>6(+!$*, ;%=-4!>*,™Va5*-4!>JxU9U Q X-4s[&'>*O#!>-_Q %=5ŸU9Q *Ÿ-4s("%'>'$&(3-I(+-"&0 ("!>q&V-"c£)š—…‚[,@U9Q

X .

*

; G

%5s-4P&(+zŒ^Q*,'>5(0'$!$(F-42%=- .G   ? * GZd `'>Q/L*, O8Ow ,V]X2%Y[ . G  ? *Gmd



c

 $  _ W ' U2 W    t W L LF & B  L Y" -H 

‚0"!$#&!$0'$x*/c-4§b'7-4!$+&(+*,'$-"!$*,%2%'7QO(+!$(s!$(n-"*v"0"("[&x*,t%P{2IJ,+!$5K%=(i[=%'$(`*,t%M#*9%+("&iJ,"!>5K0'$(F%M(""!>&(F*/D5!768"#&(%=-

‰

5 !>6&+&9-P'$Z[,&'>(*=/`(+-45v5Q,%5!>#@J,+!$5(dae_ª/%=#-YV—-"Qª#*,(+-"!>-"-4%ª&(+&bU'$  ]^&+IY%#TtJ,+!$5@!$(1-_]^!>#&k{&`-"2%z-"I0"Z[!$*,(i*,db5!768"#&(F#&*,9-T%=!$ -4P!$O/L*,"%=-4!>*,K*/g-4P("*'$-"!$*,t]^& J,*,!>Jn/L"*, %s#&*9%="("o-4*@%{&fJ,"!>5™d™e_ 02%3-4!$#':%=&V=-"&(+1#*OZ~P#&!$9-4(g%"^("%'>'!$o"J,!$*(—]^&+p-"p("*,'>-4!>*,I!$(a(+**-48d )52%N-T%1(3-4+#-"")&&5(—-4*iUX*"J9%!$q&5I%(}%15Q2%P!$# +    ŽY4T !>/2*g]X%9-4( -4*o#&*0+&(+(152%N-T%V]^!$'$F(3-4!>'$'8U&!>JP%U'>`-4*o2%Y[!$J9%=-"`-4+*,J,!>-&d





 " !  +    ŽYTT

e_c-"g]X%Y[,&'>--4&+!>*,'>*,JQ,VY%pJ,"%55c-4+&g(+-4+#-""g#&*,+"("0*,5( -4*1-"   ŽY%    % PŽ'"! dDeS-4(M5!>6&+&#z]^!>-"v-4z#&':%=("("!>#Y%' ! "! _!8"@   % PŽ'"! !>( -42%N-1-"c#*,#-"!>[!>-_Q!$s-4f-""&c(+-""#Z-4+k!$(1%'7])%YQ(1(""&58d8e_z*=-4p]g*,"5(&V *t*,'>n!$(I%5!7-"-45 !$(+!$5-4n-4+&=dDC`¦D*,+§ˆ€ ‰N‹)("* ]g&5Ÿ-42%=-k-4s5!768"&#& UZ-_]X&KU*-"@%=00"*YlO!$%=-"!$*,(X!$(J,'$!>J,!$U'>1!>n-4&+(^*/";JVH#Y%U #*,("!>5&+&5@%(p%I*O5 | T!2Ž *5y2%(I” =  ©  T! *O5(&V tU&!>JK-4y("02%=#&n5!>("!>*,‚*/-" 0+*,U'>& #T!>'$5+&s*5&(p*/ -4("%F0%"[&*5@2%(MY%+&(+-y&!>J,Y%=#T5!$+&#Z-4!$*™VD#Y%'>'$&5ª-4 !8"4 #ZŽ   #!$# dU+*-4"(^#Y%=z%'>("*oUf#*,(+!$5"&5@%(^&%"(+-p#&*,("!$(

F!>[&œ%#T!$'>5‚*5=V -"M&%"(+-k#&*,(+!$(c*/)-"M02%="&9-I*O5%+M#Y%='$'$5§-" !8"4 #ZŽ H! Z$ # | *5!$(#&%'$'>&5s% $" ]^&!>-^2%(^*o#T!>'$5+& e_o*,+5&a-4*#*,0-4-4p!>J,*,!>JF%5*,O-4J,*,!>J`\2OlO(g*/8%J,!>[&o'$&%=/ V]g1&5 !7-4(a&%"(+-g#&*,(+!$( d r&*,p*/-"& !>(D*-HZlO!$(3-4!>JV,!7-D!$(D#&+Y%=-"&5y%=( %mŽ  | d  [ $ ! 34 2   % 2 ' > ' Y  = % 8 / !$(D*-H#*,("!>5&+&5n%=(X%PZlO!>(+-4!>Jf*O5i%5!>(H*'>QI("&5/L*, $" \2ls#&*,P0-4%=-4!>*,(d | (p%I#&*,("&;(M%5F*,@!7-Yd

š}!$Jd K!$'>'$(3-44%N-4&(I-"zJ,"%55v-4+&K(3-4"#-4"K!>    CbdDK(3-T%52%"5v*O5&(n%+ "0"("&9-"&5ªUY%Y[,(IU9Qœ% N  >' !$=V -"M[!$3-42%='g'>Y%Y[,(bU9Q %  NŽBŽ_  '>!$=d   

 " ! MŽY4T

|

$! (b%z-""M]^!$#T #T%J,(I!$‚-"!$Pd r& &55™V—(+*,M*O5&( #Y%ŸU@%55&5*,o+&P* [&5™dH *‚"M%=!$ "    Va!7-Pb(+-"("0&#Z-M-4/L*,'>'$*]^!$J #&*,5!>-"!$*,(  rz%I#T!>'$5s!$(#"&%=-45™V%'$'!>-4(U"*-"&+(i%+%'$(+*o#&"&%=-45@!>n-4(4%=`-"!$P



 

|

J,!>[&Ÿ*O5@2%=(y%'7])%YQ(M!7-4( zY%="&(3-n#'$(!>vY%=#T5!$+&#Z-4!$*™Va5!:%=J,*,2%' !>#&'>5&58d r&*-pml!>(+-"!$JV#&"&%=-4F-"& %=(p*O5( | J,!7[,&*52%(p%'7])%YQ(^!>-4( Y%="&(3-i#*,(+!$(^!$s&%#Tz5!>"#-4!>*,™d r&z*ml!>(+-"!$JV#&"&%=-4F-"& %=(^[!>+-4%'™'$&%Y[,&(  

| i ( %M#&*,("&;( 505(1*,-"f%#&#4%=#Q*/ -4f("02%=#&f5!$(+#&"Z-4!>qY%=-"!$*,™dš2*, %K("#&*,5R_*,+5&o)¦FC %##&"%=-4n(+#T&PVD%(b-"n*,n5&("#"!>U5ª!>œ-4n0"[!>*,( ("#-4!>*,™Vg%§{[ZRS0*,!>9-y("0%#&@("#TK!>(M%=00'$!>&5ª/L*,&%#T5!>("!>*,™dD"/L*" ]XP2%Y[, G ”d e_‚%55!>-"!$*,™V-4P"&(+*,'$O-4!$*b"0"("&9-4%=-4!>*,t?Bƒ1%+-"&Kˆj”,¨ ‹©Gmd  P*O5&(f%"o#&'$'€R %Y[,4%Jb[%='$(%5t-_]X**,0&4%N-4*,+(c%="o5{&5§-4*s2%Y[!$J9%=-" -4+*,J,b-"X-""&=d | #*,0'$-"X5("#"!$0O-4!$*b*/2-4F Cb'7-4!$+&(+*,'$-"!$*,+&0+&("9-T%NR -4!>*,c!>(J!>[,k!>f-"!$( ("U("&#Z-4!>*,I%5k%1U+!$Z/5("#&+!$0-"!$*,cZlO0'$%!$( *]œ-"*iZl-"&5I!7-4*b!$J,)5!>("!>*,(X!>s£)%3-4&(+!:%MJ,&*,P-"+Q("!>JI%k-"&(+*,)0+*O5#-^%00"*9%=#T™d ‚5&*-4UY%Y[,(&Vh%5xU9Q  -"I"(+-""!$#Z-4!>*,x*/  -4*%y'>[,' _V



-4I+&(3-4"!>#-"!$*,x*/ d 

 



  





š2*,p-"y C #Y%(+V]gb5*-"bU9Q 2 G 2 K K -"b"**-1#&'$'†V 2 SV  V# 0  -4F5!>6&"9-*5`#&'$'>(&V - + ™-4F#&'>'7R %Y[&4%=J,1[=%'>`*/ -"F;[&' a/L"*, ("-4    O ?L*,"(+-4+!$#Z-4!$*˜Gg*,0&"%=-4*,



     +    + 

L6

-4I*,&(F*=/a-4I'>[,'

  —?Bš}!>Jd2”,Gmd

 ,V]g

?3  ,  G



!$(^*,0&"%=-4*p!>(  %[&s-42%N-p-"f02%"9-i#&'$'€R %Y[,4%JF!$(^*-"OR !$JU-1-4k]g&!>J,9-45x%Y[,&"%J,b*=/a-4c#T!$'>5"K#&'$'7R_%Y[,&"%J,(&dš˜*i%M+&J,':%pJ,"!>5

 , 

(+-""#Z-4+c!>§ YCbV!>-!>(("!>0'7QM5{25zU9Qn-4fPY%=s[=%'>

@+   G ?

     +   G G

  ”

u l+ 1, 2 i

?

@+

 

H6 @ +       G

?3 Y¬9G

u l + 1, 2 i + 1 l +1

l u l,

u l,

AÝ·æî(jî[&' _V8]g ?Bš}!>Jd“9Gmd

   

      +      

?3 Y‰9G



 



!$(X*,04%=-"*,XJ!>[,(%    % PŽY"! */ + 8%N-)-4i'$[&' 6Ÿ FU9QM!>#Y%U'$^!$P-415Q2%P!$#^J,4%=5&5 -4+&(+-""#Z-4+k%(^5Z{25K%U*[,V-"!$(*04%=-"*,b(3-1("%=-4!>(+/©Qy-_]g*P0"*,0&3-4!>&(   

eS-k2%(c-"*KU  " V !†dždh-4!$9-4"0*,':%=-"!$*,t/L*,c%z#T!>'$5‚!>(k%5P/L"*,›-4M#&'$'€R %Y[,4%J&(p*/}!7-4(02%="&9-1%5!>-"( f&%"(+-1#'$&(^!>&%#T5!$+&#-"!$*, eS-^2%=(^-"*oU  "!$# #ZŽ_!2Ž ( ©Ž sŽ 2  Y ZŽY"! V!Bdj=d  Ge_5

        

ul

+ 1, 2 i

l +1

l u l,

u l,

u l,

AÝ·æîF jîQ*,P!:%'}!$9-4"0*,':%=-"!$*,§*/5&J,+&t? T   Gmd ^5&J"&`!>(—-4"Z/L*,"^"':%=-"&5I-4*F-"#-45[=%'>de_s-"  C #Y%(+V2!>-^!>(5{&5z%(

+   G @ +   T @  

?†””,G



&(+#&*~P#&!>&9-4(—%"g"&552%9-&V=-4X(" */-4)5-T%=!$'$( /L*,}%'$'9-"XU+*-4&"(}*/˜% 02%+&9-H#&&'>'U&!$JF&;`*,y-"” = #T!$'>5"n!$(g&;[=%'>&9-)-"*b-4`W505&9-g5-4%!$'>(&d9!>(D#Y% Uml0"&(+("5K!$t C U9Q ?

@+

 

 , @+     G    



 +   ', @ +   +

š2*,gP*,+i5Z-T%!>'$(H*,-"!$(g&;[=%'>&#V]X`+/Lg-4*bƒi%3-4&xˆj”,¨ ‹†dOš2*,X%cJ,!>[&y'>[,' _V2!>-#&%@U("P%"!>q&&5sU9Q



 +   ?

 + + ,  +

( 

G



O 

Œp0&%=-4!>Jn-4!>(F*,0&4%N-4!$*x"&#"(+!>[&'>QK*, '$Z[,&'>(&V *,IJ,-"(-"o("*=RS#Y%'>'$5 Ž'4 #N Ž'"!xŽY+!$#  " *,n-"#&&'>'7R_%[&"%J,`[=%'$&(cˆj”¨N‹†d

 

 +  T 

?

 +  ,  + 

(



, V V V ,  +  , + K

G

?†”=“9G





š2*,I-4z”C #Y%(+VD]X5&*-4U9QB2 G 2 K K K -4@+*O*-o#&'>'†V 2   V  V ? # , G 0  ,-4g5!768"' 2   dhš2*,iY%(+*/H*=-T%=-"!$*,™V]Xo5*-4oU9Q + QG ?- +   G    -" &(+&bU'$*/ -4ZlO!>(+-4!>Jo#&'$'7R_%Y[,&"%J,`[=%'>&(p%=--"'$[&' _d







  



f0"*+#-"!$*,s*,0&"%=-4*,)/L*^%o+&J,':%p£)%3-4("!:%=J"!$5s!>()-4@5Z{25zU9Q

@ +    G ?

     +   G 

 

G ‘

 









K K  “

@+

 

 







% 5s-4f0+&5!>#-4!>*,s*,0&4%N-4*,U2%=("&5@*,@%P'>!$&%)0*,'>Q*,P!:%='8!$9-"&"0*,'$%=-4!>*,y!$(^5mR  {25@U9Q,V2/L*, , 0 !  ,Y  )505&9-15-4%!$'>(&d | 2  %='$*,J,*("'7Q§-"*K-4n”C #Y%(+V}/L*,f-4s“C #Y%(+y]Xs5*-4nU9Q 2 G 2 K K K K -" "**-c#&&'>'†VW2    V  ^V)?# , ,G 05 —-4y5!>6&"9-b*5M#&'>'$(&V - +   o-" #&'$'7R_%Y[,&"%J,[=%'$*/8-"1;Q*,P!:%='8!$9-"&"0*,'$%=-4!>*,y!$(^5mR  {25@U9Q,V2/L*, , , - 0 !  ,Y  )œ-4x C %=5Ÿ”C #Y%(+&(V—-4n5Z-T%!>'$(I%="  +   NG @ +    T @   %5œ-"  W505&9-15-4%!$'>(&d e_w#*,#'$(+!$*,™VH-4tW'7R %Y[&4%=J,K[=%'$(n*/c%'$'p-4t'>Y%Y[,(!>( &;*,U2%' &+"*,`U-_]g&§-"I#&'$'€R %Y[,4%J,f[=%'$(`*/D-4IZl%#-(+*,'$O-4!$*K%=-`-"  '$Z[,&' V25*-45zU9Q @ + ' V˜%5s-4*,(+*/}-"b'>-"!$+&("*'$-"!$*,y#*,P0-T%N-4!$*y]^!>-4@%





 Y‘



M% l!>b



@+ 

'$Z[,&' V25*-"&5zU9Q

.. @ +  ' T @ + 

 . .

V#Y%Uf5&#&*0*,("5@!>9-4*I-_]X*P&+"*,+(



.. @ +  ' T @ +   .. 6 . . @ +   T @ + 

 

..



?†”N‘5V#Y%'>'$5   #4ZŽYŽ'"!K& " !$(D-4p*,p*/™-"^{2!>-"p[*,'$^("#T&P`*, -4F{2&(+-1J,+!$5n*/'$Z[,&' d2eS-^#Y%=@UU*,5&5U9Q

.. @ +  ' T @ +   ..



( 

”

 V

 O

?†”¨,G 

]^&+ .!$(1-"k#*,9[,"J,#&I*,+5&1*/a-4c{2!7-4[,*,'>Pk("#T&Pd8e_@-4c0"("&9#Y%(+Vg%(]gz("z(+&#*,5ORS*,"5n%##&"%=-4@("#T(y!$.-4!>%5v("02%#VD]X@2%Y[, [G ”Od š2*,-4c("&#*,5K"+*,&V˜#Y%='$'$5 2mŽ H TŽY"!.& " V£X*,& ZŽ ˆ¡  ¨N‹ ("* ]g&5K-"2%=-&V !>/™-"i5-T%=!$'$(H*,n%c'$[&' %="i5&'$Z-45y]^&n(+M%'>'$D-"2%n%c0"("#&+!$U&5y-"*,'$4%# + !$(H#&*,9-"4%#Z-4!>[`!$o-4i#T*,(+&n*"V

LVO!7/8-"i5!>("#"-"i-4!>^Z[,*,'>-4!>*,*04%=-"*, w %5!>/h-4-"*,'$4%#& %N-p-"'$Z[,&' !$(("Z-^-4*





G

L

”



= (  



]^&+ F!$(}-")("0%#&5!>("!>*,™V-4&I-4X5!>6&+&#)UZ-_]X&o{2!7-4g[*,'$)(+*,'$R -4!>*,I*,o-4){2pJ"!$5P%5o-"^("*,'>-4!>*,I*,U-4%!$&5IU9Qob'7-4!$+&(+*,'$-"!$*,k%#&#b'$%=-4( !$n-4!>`%5s[,&+!>{2(

.. @ +  T @ + 

 ..

O

 V 

?†”=¢9G 

]^&+ 5&*-4(o-"y!$&%h#&*,9[&#-"!$*,R_5!768(+!$*,c&;#T]gc5&5#f-"2%=-

 ‘

*

M86 X 

*

G  ‘

”

M86 X

(



”

 (

  ¨

V 



 

V

eS/8]g1]X%9-g-"p0+-4"U2%=-"!$*,P&+"*,a-"*cU1*=/8-41("%^*,+5&g%(a-415!>("#"-"!$q&%=-4!>*, &+"*,V]X&5@-"2%=-



(  ”



.



!†džd

”

 +

‘

M 6 X ”

(



”

-

C {2!$J^-4gA#&'>-bU ` -4z%(

( 

”



( 



 G X [M

(

 V-4H0+[!$*,(#*,5!7-4!$*f#&%kUD"Z]^"!>-3R



 6v”  3

?†” ,  G

 

š *,H-"i!>9[O!>("#!$5M#&%("b?!Bdj=d$M G  *, 2 687uGZV?†” ,GD!>(g;[=%'>&9-g-"*k-4i"("'7(    de_M-"i'$&5 4 Y4! TŽ N$ +! T ]^!>'$'8U("Z-p-"*



 G 

”

( 





 6v”  3

?†”=¬9G

*o#T**,(+%%#&#&0-4%U'$1[%='$i/L*X-4`/%=#-4*  V%b("&+!$(*/}#*,P0-T%N-4!$*(g]^!7-4 5!>6&+&9-^-4*,'>&"%#&(]^!$'$'U#&("(4%=+Q,V%((+*]^z!$s(+U("#-"!$*,K¨Od> d

 "!$#&%NŽ'%

 !

# %

   OŽ†ŽY"!

*§!$'>'$(3-4"%=-4M-"s#&*,("&3[=%=-4!7[,s\2Olª#&*0O-T%=-"!$*,™V ]Xs{2+(+-#&*,(+!$5P%ª C '$&%=/ 2    ]^*,(+)#&*,("!$(}!$-4X"!>J,9-}5!$"#-"!$*, 2   3 %5 2    %="D[O!>+-"2%'†d &+/L*,+f-"&!>)/%=-4& 2  !>(^%I'$Y%N/ d













2 i +1

2i

2 i +2

i

2 i +3

i +1

l +1

l

AÝ·æî.î*,   E      3 !$(X*-U2%'$%#&5y]^!>-"y-4i*,-4J*,!$Jf\2OlM/L"*, 2   !$M-"F'$Z/©5!$+&#Z-4!$*3  E  Bd i/D#*,"(+V8]gI#&*'$5x5!>"#-4'7Q#&*,P0-"k-4I*-4J,*!$JP\2OlO&(  /L"*, 2  -"* 2  !>x CbVU-^("#Tz%I#&*0O-T%=-"!$*,y#&%*=-^U`ml-4&5&5-4*    !$J&5!>("!>*,(V2%(]X#&%@(+&f!$š}!>Jd2¨Od



      

 ¢



2 i +1, 2 j +1 2 i , 2 j+1 l +1 2i, 2 j

2 i+1, 2 j

i, j

i +1, j

l

AÝ·æî  -"*,-!>#&+Y%(+!$Js("!>J,!7{2#Y%9-"'>Qz-"o'$*]^!$JV -4s0+!$#!$0'>M*/^-4s%'$J,*"!>-" !$(b0"&(+&9-45™dHš}!$+(+-&V}505OR !$J*,z-"b!$!>-4!$%'™#*,5!7-4!$*™V% ! ©ŽY +    Ž'4T !$(p#&+Y%=-"&5™d8&™VJ,!>[&z-" J,4%=5&5 -4+&s(+-""#Z-4+VD% ŽY Mz%N Ž'"! !$(I%5M*, -4 $"% # d}&ª5-4%!$'>( %"#&*,P0-"&5œU9Q O¡ŽY4 #N OŽY"!ªŽ'+!$#  " Vh!$œ*,+5&f-4* T M # -"-""d} * Uy%U'$P-4*x2%Y[!$J,%=-4!$(+!$5P-4y-""&M(+-4+#-""=V}]gn0"*,0*,(+M-4*x(+s% 4  H #% %'$J*,"!7-4d™#T*,(".52%=-T%(+-""#Z-4+y#Y%‚2%5'>z CbV”C %5 “,C £)%+-"&("!$% J,&*Z-4"!>&(d



 

 







#& ŽY"!

Ÿ fD  © +  B  WF 9!>-I#&'$'^%=5Ÿ#*,P0-4n-"s#&&'>'7R_%[&"%J,y[=%'$&(o!$ -4#T!$'$5"&.#&'$'>( ]^!>-"s-"!$!7-4!$%'#*,5!7-4!>*, C1£X*,P0-45Z-T%!>'$()!>n-4#T!>'$5+&s#&'$'$(U9Q ¡Ž'4 #N Ž'"!zŽ'+!$#  " C^e š-4o5Z-T%!>'—!$§%n#T!$'>5t#&&'>'a!$(F!>J,`-42%x-"0"&(+#&+!$U&5t-4*,'>&4%=#&=V ^ƒiE  ("0'>!>-)-"!$(#T!>'$5

 

C^^/L*,"P&D#T!>'$5PU#&*,P&(g%f0%"#&'$%=-4H-4)[=%'>&(—!>b*5&(D%=5b[!$+-"2%'O'$Y%Y[&(aU9Qb0+*+&#Z-4!>*,b/L"*, -4'$&%Y[,&( nW

 hPF  W ˜" W vF  Y=  hF ?Bš!$Jd¢9G Cpš2*,™-4D]^*,'>a-4"H/L+*, -"D'$Y%Y[&(h-4*^-4H"**-g?  ,#>0  >G £X*,P0-"15Z-T%!>'!$o-4p*O5  ˜U9QP% O¡ŽY4 #N OŽY"!sŽ'+!$# 



"

 



e šª-"k5-T%=!$'$(^!>@-4!$(^*5b%5z!$@!>-"(pU"*-"&+(`%+k("%'>'$& -4%-"0"("#"!$U&5s-4*'$&"%# ^ƒiE  -4#&'$'™%5!>-4(U"*-"&+(1%="  Z$ZŽS&$ Cpš2*,™-4D]^*,'>a-4"H/L+*, -"D'$Y%Y[&(h-4*^-4H"**-g?  ,#>0  >G e š -4!>(*O5%5@%'$'™!7-4(#T!>'$5+&s*5&(p%+5&'>-T%=U'$ | 1C -4@#T!>'$5+&ª*5&(M%"("!>0'>n'$&%[&(K?!†džd—-"QŸ%+ '$&%Y[,&(^]^!>-"*,-[!>+-4%'™#T!>'$5+&˜G ^ƒiE  5&'>-"f#T!$'$5"& ?!$(^&%()-"2%=-^]XW=&&0z*,FP*,"F'>[,'©G e š -4!>(*O5F2%(^*P#T!$'>5" | 1C !7-^!$(*=-^5&'>-T%=U'$ | 1C !7-^!$(*=-p%=-^-4%NlO!$b '$Z[,' ^ƒiE  #&+Y%=-"-"#T!$'>5"n/L*,)-4!>(*5 ?!$(t&%(K-"2%=-t]XŸ%55 *,ªP*,+ '$[&'b/L*,x&%#T OR 5'$-4%U'>F'$Y%N/TG



 

  

    ‚+!>-"`-4P&("@%5s-4f#&'>'7R %Y[&4%=J,`[=%'$(!$9-4*P%b{2'>d

 ¬







E iC C    )   C`&%'$'>*O#&%=-4p-T%U'>&(





AÝ·æî  î^õ+݇à"ӜӂÞDåCâºßآݷåC×åDä|Ø¢ãÓ¥ØÃբӂÓ`ÖÃØÃՅß ?B¢9Gmd

 

Ž† Z# Ž' ZŽ 4

_e -"f!$P0'$*,y*/ -4k%'>J,*,+!>-4@V[,'>Q %sJ,!>[&‚#&'>' 5&0&5!$JP*,n-4f'>*O#&%'™+&J,'$%"!7-_QM*=/-"("*,'>-4!>*,™d e_K-4I0"*,J,"%V-"k-4+&I!>(`"0"("&9-"&5§U(+-4(f*/ %I("-^*=/   MZŽY&f!  #"  !2Ž'©ŽY† # V2%5 "!2Ž_ #FŽ PŽ 2 ©  4! V]^!$#T %"!>y-4mlO-p{2^("UJ"!$58d2  NŽ !>()-4F{2+(+-p#&'$'™*=/-"F-4"c(3-4+#-""=d



E %#T *O5I#&%§Uo%55"("("5‚U9Qx!7-4(`'>[,'aK-4!>( D '$Z[,&' #X?# , b!$M”CbV[,(g-"!$(DJ,*9%=' ?Bš}!>Jd ,G (+-4%+-"!$JI/L+*, -4c+*O*-^'$-{25(p-"02%=-"-"*P-"#&&'>'†V'>[&'hU9Q '$Z[,&'BVU9Qs*,0-"!$Jb/L*,)-"#T!$'>5n]^!$#T!$(^%@%#&(3-4*,^*/ -4F-4%"J,Z-^#&'$'†d







'$(&VN]X%[ G  ” = V!†džd G  '$*J V bU&!>Jc-"i5!>("!>*,™d&("(V/L*,f%s!$J'>Q #&*,P0+&("(+&5 (+*,'$-"!$*,8V8-4 (+-40(b%="P"&;"5œ*,'7Qt!$§("%'$'a+&J,!>*,(c%5 -" (("2%='$'>QMb#T("%'$'>&X-"2% 8$&%'Qp/L*, ':%="J,—(+#Y%'>D“,C.#&*,P0-4%=-4!>*,(&d 2%=-D!>(}]^9Qb]X5&#&!>5&5P")-4*F*0-—/L*,a%F(+*,'$-"!$*,c]^!$#To*,0-"!$P!$q&(-"P&P*,+Q ";!$"9-4(&V2]^!7-4%I("'$!>J,9-4'7QMP*,"F#*,P0'$mls%'>J,*,"!7-4d



@ ()("#-4!>*,™V]gF0+&("9-k CbV”CbVO“,C "!$#&%'8+&(+'>-"(("!>Jo%I("&#*,5ORS*,"5& %#&#4%N-4P("#TM]^!7-4 %@b'>-4!>"("*,'>-4!>*,x%##&"%#Q§*/ G “dh!$(fPY%(-4%=*,'7QM*,FY%+&(3-1#&'>F!$nY%=#T@5!>"#-4!>*,n!$(X"&;"5™d!>(#T*,!$#F&2%=U'$(p()-4* M%=!$9-T%!>o%f%"+* ]v-4"^(+-4+#-""i%"*5P"&J!$*,(a]^"1!$J,o'>[,'$(g%"^&55™d

   "! % ZŽY"!    #"!§ Ž'"!



+*-,/.+G10.ˆTf ,Y &‹/4iˆ  ,&687uGmV

‚a#*,(+!$5h%^'>!$&%3RS#&*,9[&#-"!$*,F5!>6("!>*,`&;#'>&J=-4.("#&%'$n-4(+!$q*=/i-45*,M%=!$ %5ª%( T# 2%"%#-"&"!>(+-"!$#i-"!$Pi("#&%'$ :G VO-4!$();%=-4!>*,n#Y%sUi]^"!7-"-"&n!$y-"i/L*,'>'$*]^!$J 5!$P&("!$*'$("()/L*,+



? @ ? @ 6 G ? . ? * ]^&+

 



? @ ? *

?B“O  G

F5*-"&()-4A&#'$Z-^bU ZG  

@LK +? *G G



 





!7/

* 

!7/

* O

 

” 

d 

`#T**,("%(!>!>-"!:%'2#*,5!7-4!>*,

%5@C`!$+!$#T'>-)#*,5!7-4!>*,(^%=--"'$/©-p%=5s"!$J