Iterative Methods Multi Grid.pdf






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

/

0214365879;:<>=)?

@ ACBEDGFIH

J$KML(NL(OQPSRTLVUWPXNSLVY[Z \]LOQ^_`ZbaScedSLf^gYh_jiklYh^gmnUoZprqVsutXvw R]x]y]OQz{_r\]LfUW|nY _rZ4Rn^gRTL(Y } z~Z\)R]Yh^S€!Zz{€(^SO YLf_rx]O~Z _V Jƒ‚[„…L(Yh_lZ\nL4R†PC__z~ynz~OQz‡Z`dePgˆ_PSOQNuz~mn‰Š^R]Y PSynO~LV‹ Z z~Pmn^gO,ZPZ\]L4mux]‹y†LVYjPgˆx]m]Œum]P } mn_V

} z‡Z \ } PYŒ8^SmnU_rZPY ^S‰SLR]Y PSRTPSYŽ

J‘’L(OQO6UoLVNSLVO~PR†LfU>ˆ“PSY”OQz~mnLV^gY”LVO~OQz~R]Zz{€ŠR]YPy]OQL(‹_[i^gR]R]OQz{€(^XZ z~Pm•Z PPSZ\]LVY”LV–Cxn^XŽ Z z~Pmn_jz{_j_`Z z~OQOG^Sm^€ZzQNSL4^SYLf^PSˆY LV_LV^SY €h\, —l˜X˜g™šœ›gŸž`˜g™S C¡ `˜(žŸ¢£¤ ¥Q¤ ž`¤h›¡h¤f¦§‘¨8©Sª «E¬‡­V®r¬“¯8°C©g«²±V®¬E³fª ´µ·¶ ¸,¶u¹Gž`º‡»V»f¼h´ ½ ¶ ¾;¶C¿¤h›C¼Ÿ˜V›†´CÀV›C™ŠÁ]¶ ¶uá!ĘVž`Å4ºÆ¡rÇo´oÁgšÉÈÊÃËØV›C˜V»VžrÀVÌ͆´oÎVÏfÏVÏS¶

Ð ÙÊÚrÛ

0214ÑÒ!3ÔÓÕ9Ò(=)3 ÒVÖ>×!Ø Ñ @

ÜlÝÞÕßáàrâ[ß,ã;ä

ACBEDGF’å

æ>ç>è;éEêEê ëuì(ç!íŠîXï…ðhçòñ(óCñXéEïníô“õ]çWîXï…çWögéE÷Šç(ï]ìé²îXï…ñXêø…ùîSú(ê~ç(÷ûôÉîüéEêEê‡ëuìhôýùñgôÉç•ôEõnçòúhñþìé²ð øTùéEï…ð!é ø…ê~ç!ìî`ÿ4÷ëoêÆôýé‡ógùé²öe÷ç!ô“õ]îföþì
…õ]çé²öç ñXìÊøTùç!ì(ç!ï†ôÉçhöðhñXï’úhçùçhñSögéEêçôÉç(ïTöç ö8ôÉî ÷ëoêÆôœé ø…ê~çòögéE÷Šç(ï]ìé²îgï]
ì 6éEùhìhôœí4õnîXè6
ç þç(ùhíè6ç•è;éEêEêjù
ç Vé²ç(è ì(îX÷ç)ìéE÷ ø…ê~çbÿ!ñð!ôœì•ñSúhîXëoô éEôŸç!ùñgôœé Xçe÷Šç(ôEõnîVöXìŠì!çhç(ï éEï ô“õ]çeê~ñXìhôê~çhð(ôœëoùçñXï…öòöS
ç Xç!ê~îhø$ì!îg÷Šçeé²öSçhñXì8ôÉî•ñðhð ç(ê~ç(ùñXôŸç ô“õ]çéEôÉç(ùñXôýé XçÊø…ùîVðhç!ìh ì

q •x]O‡Z z~‰Yz{U z{_Š^gm z~ZLVY ^gZzQNSL‹L(Z\]PoU ^ óîfîfö8éEï†éEôœé²ñgêóSënçìhì } zQOQOlY LVU]xn€!LZ\]L mux]‹y†LVY Pgˆz‡Z L(Yh^XZ z~Pmn_
 ZP_PSOQNSL   yud8z‡Z L(Yh^XZ z~Pm } Lb€(PSx]O{UeZ ^gŒL   !#" } \]L(Y L$%!& !'  (( (ç!ù ç  í ) ñgïTö  öç!ï…îXôŸçôEõnç)÷ñXôýùé*gí+ Xçhð!ôŸîXù>î`ÿëoï-,gïTîgè;ï]ìhíñXï…öŠÿ!îXùðhç 

Xçhð!ôŸîXùhíùçìýønçhð!ôýé þç(êþíô“õ]ñgôbùç!ìëoêÆô²ìÿhùîX÷ ögé{ì(ð!ùç(ôýé/.!éEïuó îgëoùe÷îVöç!êø…ùîSú(ê~ç(÷ îXï ñógùé²öWî`ÿìé/f. ç102 ) ÷Šçhñgï]ìŠôEõnñXô[ô“õ]çeéEïnéEôýé²ñXêógënç!ìhìÊÿ!îgùôEõnçéEôŸç!ùñgôœé Xç 3 4 ø…ùîfð ç!ìhìîgïòógùé²ö$ 0 é{ì6Ÿ5 ñ øøTùî7  éE÷ŠñgôÉçhö9 8 ÿhùîX÷ ôEõnçŠì(îXê‡ëoôœé²îgï î`ÿôEõnçø…ùîSú(ê~ç(÷ îXï ñóSùé²ö>î`ÿìé/f. ç a 02üæ>çëuì!çô“õ]çè6îXùö>ñhøSøTù î
  éE÷ŠñgôÉçúhç ðhñgëuì!ç   ñXï…ö   ñXùç X çhð!ôŸîXùhì8ô“õ]ñgô[õ]ñ:X ç>öXé M ; ç!ùç(ïnôê~ç(ïCóSôEõuì
æ>ç8ìõnñXêEêì!çhçeê~ñXôŸç(ùõ]îgè ô“õué{ìñ øøTùî7  é< ÷ñXôýé²îXï+é{ìŠð ñgùùé²çhö>îXëoô=$>ôjé{ìð(ê~çhñXùôEõnñXôÊô“õ]çMø…ùîVðhçhöXëoù çîXëoôýêÆéEïTçhöé{ìùçhð(ëoùhìhéX ç ? ô“õ]ñgôé{ìhí;ôEõnçøTùîú!êQç!÷ îXï>ô“õ]ç óSùé²ö a + 0 ðhñgï>ñgê ì(î8úhçì(îXêX

çhö8ú
  éEôŸç!ù ñXôýé²îXï]í;è;éEô“õ>ñXï éEï†éEôœé²ñgêCóSënçìhìøTù î:f é²öSçhöú
4  ì(îXê*V

éEïuó %@  )@  A  í ! ç Ÿ ô B7C78æ)ç;ønîXéEïnô;îXëoôôEõnñXô ð @  éEôŸç(ùñXôýé²îXïüîXïòðhîVñgùhì(ç!ù4÷çìõnçì”é{ìð õnç ñhø]ç(ùúhçhð ñgëuì(ç‘ D é{ìMì÷ñXêEê~ç(ùhílñXï…öŠô“õ]ç(ùç²ÿ!îgùç

q

EGFIH
JKH1LNM+O/H M!M+P Q:JKR$SIH
J+LE=H JT:EULVQ:WIX&T:WZY\[!H!]!T:^-M
H4_ H P H
J+LE=H JKT`EaLVQ`WIMbT`JKHJKH!c ^dLJKH!Y:e LVfHCf4ghLNMiMkjlT`OO/H J
f

m#ndo+p
p7p

oKqIr$s-nIt1m9r7uwv`x2y*or
u!z:oKy{v|s9}sIr
rC~dr7~\oKv }v€/‚riƒB„…4„1†‡ „ }ˆoy/€{€Š‰l‹ gŒŽ

 †

‘ ’`“ ‘

” LWZ]!HiM
Q:j•HM!j•QCQ:EGF–]!Q:j%SIQ`W—H
W#EaM$Q_$EGFIHbH JkJKQ`J$P˜LOOMkEULOOJKH
jlT`LW—f ™—šœ›Ÿž¡  ¢ p£
¤z:xor
u%zbxGr7¥y/oKr7uz`oKy{v|s#}7¦-oKq#r H JkJQ:J$LNMiMkjlQCQ:EGF ¦#¥)r+§
v|nI€N~3}v€/‚rxGvuoKqIrr
uuKvu vs¨z ]!QCT:J!M
H
J$jlH MKF ¦drp ©ªƒ „&« „ †¬ „ p Œ Œ Œ ­¯® QCQCY•LVY|H!T°[!H]!T`^-M HC± ²A³ tlv-v`oq\xGn#s#§koy/vs#}w§
z|s\m—r$ur
´Iur7}Kr
sœoKrC~ªvs°§
vz`u!}Kr
u©|uy{~#}
µ ²¶ vz|u}KrB©|uy{~°}v€/nIoKy{v|s#}z`ur$§!qIr7z|´9r7u7p

· F-LNMLVYH!T–LNMªLW3_ T|]
EiEF#H°]!H
W#EUJKT:O)LVY|H!T¡Qˆ_ªjb^dOEUL ® kJ LVY¸E=H]!F-W9LVc ^#H M
f¸¹ WºQ:JY|H
JE=Q¨EU^dJkW EGF-LNM$LVY|H!TLW#E=QTBSZJKT|]
EULV]T`O2T`O ® Q`JkLEFdje)M
H »:H JT:OLW ® KJ H!Y`LVH W9EVM$P˜LOO2[!HJKH!c ^dLJHYf

¼½a¼

¾¿ÁÀÀ˜ÂIÃwğÅ

™—šœ›Ÿž¸Æ

x oKq#r FdL ® F_!JKH!c ^#H
W—]
È § vtl´9vsIr
sœo}Šv`xdoqIr r7uKuv|uŠ~dr7§7zŽÉ%xVz}ˆor
uŠoKq#z|s$oKqIr O/Q:PÊ_!JKH!c ^#H
W—]
È Çd § vtl´9vsIr
sœo}7¦d¥4ri}zŽÉ•oq#z:ooKq#r$y*or
u!z:oKy{‚|ritlr oqIvd~°y{}wz Mkj•Q7Q`EF#H J p

ËÌNËÌVÍ

ÎÏIЁÑ#ÒÔÓ

™—šœ›Ÿž¡Õ LOW MODES

HIGH MODES

1 0.8 0.6

λ(R )

0.4

J

0.0 0.2 0.4 0.6 0.8 1.0

0.2 0 -0.2 -0.4 -0.6 -0.8 -1 0

2

4

v2 (mode k=2)

6

8

10

12

mode k n=19

Ç }&֜z|§ vmIyzl}tlv-v`oqIr
uk×

14

16

18

20

v15 (mode k=15)

¢

p
p7p!ØÚÙªÛ

ܨÝ+Þ
Ý!ÝbßGàIá:ßÔâZãVäå`æ4çéè ê9ëìãVäå:æ7èdçéè êZí|ãVäåŽæ7è˜á:î—ï:ð)ÞkñîZòÝóñNÞ+ßà#Ý+à-ñ*ôŽàIÝ Þ!ߊõ!öKÝ!÷ ø#Ý
î—ò
ù ï Ýkð ñß)ñNÞò ü/Ý!á`ö+ßà#á:ßýdá|ò û|þ ñ˜ñNÞ+î û ß4álÞ ú•û7û ßà#Ý ö
ÿ ú•û 





      !

"$# %'&)(+* 1

0.5

λ (Rω J)

ä-,-åBç/. äå10 ã32546.)æ7

ω =1/2

0

ω =1.1

-0.5

ω =2/3

(UNSTABLE)

ω =1 -1

ê$8dãVä9,-åŽæ ç:.)ê I 8 ãVäå:æ0Aã;2-46.)æç<2546.%ã3254

mode k

ê$8dãVäå`ææ>=

? @ ç 2=BACABAC=ó

Ü¨Ý ûþ Þ
Ý
öEDŽÝ°ßGàIá`ßÔõ û öF.HGI2|ðBýdá|ò û|þ ñBò!á`î6ñî°õ á|ò
ß þ Ý3á°ô ûCû ï¸Þ ú•û7û ßà#Ý ö
ÿKJVõML ÝÞ
Ý ß ßGàIÝ°ò û î—ï`ñßañ û îè êN{ëPOŠí3QSRTãaä ,-å æ
èŸç è ê9ëãVä ,-å æ7è ð>L Ý ûþ ß=á:ñîU. çWVYXPZ9ÿAÜ3Ý°á:ü Þ û î û ß Ý ßGàIá:ßZõ û ö[.:\]2ßGàIÝ ú Ý ßGà û ï þ Ý!ò û:ú Ý Þiødî#Þkß á þ ü/Ý^ˆï û ÝkÞ+î û ß&ò û îD:Ý
öaôÝ;_•Þkñî—ò!Ý)õ û ö+Þ û`ú Ý ? ðè ê 8#ãV9 ä ,-åŽæ7) è \W2œÿJ î Þ û:ú ÝlÞ
Ý
îIÞ
Ýkðßà#Ý>`—ökñVò!Ýß ûÊþ Ý>`#á:ñVï1õ û öªýdá|ò û|þ ñ&ß û–þ Ýáô ûCû ï Þ úlûCû ßGàIÝ
ö$ñNÞálÞkü û Lï û L˜î¸ñî¡ò û îD:Ý öaôœÝ
î—ò!Ý û õ$ßGàIÝü û L–õ!öKÝ!÷
ø#Ý îZò ù ú•û ï|Ý Þ ÿ

"$# %'&)(ba

cedfCgihjdlknmo pqglfrYsktgfBudlmvgfBuswCfFhPoxfCgglmYgzyvmufF{|h~}ShwdlmYg5mP}q2B€€ n=19

600

ω=1 ω=2/3

Number of iterations

500 400 300 200 100 0

0

2

4

6

8

10

mode k

12

14

16

18

20

ZàIÝ1ô`öKái`IàÊÞKà û LÔÞlßà#Ýlî#ø ú•þ Ý ö û õbñß Ý
öKá`ßañ û îIÞ1öÝ÷
ødñöKÝ!ï°ß û öKÝ!ï:ø#ò!ÝlßGàIÝlßGàIݪá ú `ZüñßUø#ï|Ý õ Ý!á|ò!à–Ý
ökö û ö ú•û ï|Ý þ ùªá&õ á|ò
ß û ö û õƒ‚Y„„|ÿ\ܨÝiÞ
Ý!Ý$ßGàIá`ßßà#ÝiÞkß=á:î—ïá:öKïªý-áò û|þ ñã'. ç@2Cæ û $ á:ü ô û ökñßà ú ðöÝ÷
ødñöKÝ Þ ú á:î9ùÊñß=Ý öá:ßUñ û îIÞ°ß û Ý
üñ ú ñî—á`ß Ý3ßGàIÝ\à-ñ*ôŽàIÝ Þ!ß4õ!öKÝ!÷ ø#Ý
î—ò
ù úlû ï|Ý Þ
ÿ … î¨ßà#Ý û ßGàIÝ
ö%àIá`î—ïŽðßGàIÝ+ødî—ïÝ öE†=öKÝ
ü/áB‡œÝ!ï•ýdáò ûþ ñ2Þ
ò!àIÝ ú Ý ð)Ý ü*ñ ú ñîZá:ß=ÝkÞ$ßGàIÝ%àdñ*ôCà úlû ï|Ý Þ D:Ý ökù3÷ ødñViò ˆ:ü*ùŽð þ ødß û î ßà#Ý û ßGàIÝ
öià#á:î—ï:ð&ßà#Ýü û L õ!öKÝ!÷ ø#Ý
î—ò
ù ú•û ïÝkÞ1ß áˆ|Ýbü û îœôœÝ ö!ðßGàIá:î L˜ñßGàºÞkß á`î—ïá:öKï ý-áò û|þ ñNðbß û ï:ñNÞ
iá ``#Ýá`ö
ÿ Ü3ݖÞKàIá:üüwÞ Ý!ݸßà#á:ß1ßàdñNÞ°Þkü û L ï û L˜î ñîAßGàIÝ ò û î DŽÝ
öaôœÝ îZòÝ û õ4ßGàIÝü û Lªõ!öKÝ!÷ ø#Ý
î—ò
ù úlû ï|Ý Þ4ñNÞî û ßöKÝ!á:üü*ùz á `—ö ûþ ü{Ý ú á:î—ï$ßà#á:ßað þ ùiø-Þkñîœô ò û á`ö!Þ
Ý ö ú ÝkÞKà#ÝkÞ!‰ ð L Š Ý L˜ñüü þ Ýbá þ ü/Ý1ß û ‹Þ `#Ý!Ýïª ø `¡ßGàIÝbò û  î D:Ý
öaôÝ
î—ò!Ý û õ$ßGàIÝ Þ
Ý ú•û ï|Ý Þ ÿ 

nŒ Ž  ‘ !S  “5fwChP”n”‹• Z

"$# %'&)(+’

1 0.8 0.6

λ(RGS)

0.4 0.2 0

-0.4 -0.6 -0.8 -1

0

2

4

6

8

10

mode k n=19

v2 (mode k=2)

12

14

16

18

20

v15 (mode k=15)

–;—[˜š™P› —l—3œžŸC ¡ŸC¢™¤£¥¥¡—3¦v¥¥P§¨ŸC©Eª

«C«C«

¬ '­ ®°¯i±ƒ²´³±v±C­¶µY±C®·j±i¯²ž¸j¹iº¤¸3»-¼>½¾U¿P®$À~ÁÂÀY¸Ã®$¸j²5¯i¸j­'®°¯­SÀY±EÄq²'³±¹±Š­ºƒÅ¶­'²²‹Åt±~Ɖ±M¯i¿j®KºC¿PÇ ¿Èi¸Pɲ1²´³±ŠºEÊM¸B¸P²'³­'® µŠË$¹¸Ë±¹E²‹­S±Eºv¸;»ÍÌ[¿PɺºEÎe¬±­SÀY±ÅÏÈÇÅt¸B¸Ðj­'®µK¿j²q²´³±v±C­¶µY±C®·j¿jÅÑɱº¤¸3» ²´³±Š­'²e±C¹l¿j²‹­S¸j®UÊv¿P²¹E­¶ÒÓ –e§ŸC©i™j§ t¥YÚ—1©lŸÛY› t©ŸB¡§l¥v©ŸB¡›ÜCŸƒ™PÚKÝ ±C¹E¹l¸j¹FÊM¸BÀY±5Þß ™~àS™Ü§l¥Y©5¥PàqáBââ

Ô$Õ Ö'×)Ø+Ù

n=19

200

Number of iterations

180 160 140 120 100 80 60 40 20 0

0

2

4

6

8

10

mode k

12

14

16

18

20

«B«C«i˜F —5™ µ ¸B¸ÀvºEÊM¸B¸P²'³±¹ «

Ç Åt¸¸BÐP­'® µK¿P²1²´³±~® ÉÊvÈi±C¹Š¸;»ƒ­'²ž±C¹l¿j²‹­S¸j®ºŠ¹l±iäÉ­'¹±À²ž¸Í¹l±iÀPɯ±~²´³±v¿jÊ-Ë$Ŷ­'²ÉÀY±M¸;»~±i¿Y¯³ ã Í ÊM¸BÀY±EÄM¸;»x²´³± Ý Êv¿P²¹E­¶ÒĤÈCÇå¿M»¿¯C²ž¸j¹6¸;»KæYççjĤƉ±U¯¿P®èÀ±C²ž±¹Eʤ­'®$±6²'³±xºEÊM¸B¸P²'³­'® µ Ë$¹¸Ë±¹E²‹­S±EºU¸3»é²´³±åÌ5¿PɺiºiÎe¬±­SÀY±Å-ºC¯³±ÊM±BÓëêE²Š²‹É¹E®ºK¸jɲ¤²´³¿j²MÈi¿ìºC±iÀå¸P®è²´³±+¿Èi¸j·j± µ¹l¿Ë³ Äí²´³±9³­¶µ³ƒ»i¹l±iäɱC®$¯CÇvÊv¸À±ºŠÀ¸v­'®~»¿¯C²qÀY±¯i¿PÇM¿P²1¿vʤɯi³š»¿ìºE²ž±¹š¹l¿P²e±ƒ²'³¿j®K²´³± Åt¸PÆK»i¹l±iäɱC®$¯CǸj®°±EºCÓ

îqïð

ñ:ò)óYôõö3÷$ôö;øzù

Ô$Õ Ö'×)ØûúPü

ý ±Kºl³¿jÅ'Å>¹l±iäÉ­'¹l±MË$¹¸B¯i±iÀjɹl±ºŠ»¸j¹U²¹¿j®º´»±C¹E¹E­'® µû­'®B»¸P¹EÊv¿P²­S¸P®þȱC²‹Æ‰±±C®µ¹E­SÀìºCÓ ÿ°³± Ë$¹¸B¯i±ºº-¸3»5²¹¿j®º´»±C¹E¹E­'® µ¤¿ƒ·j±i¯²ž¸j¹»i¹¸jÊ ¿
‰®$±5²ž¸Š¿¤¯i¸B¿P¹iºC±[ÊM±Eºl³Í­º-¯¿PÅ'Åt±iÀƒ¹±EºE²‹¹E­S¯²‹­S¸j®°Ó ˜> ŸCÚ 1Ÿš¥ Þ §i™P nÚ  Þ ß

 ! "     " $# ©ŸB—3§l© nܧl n¥Ú¥&%$ŸC©i™j§l¥Y©('S¦M™j§l© *)$+« ,

-/.0 1$2*3547681:9<;>=?35@>A$93B.C4EDGFHIJKLDNMF

O P"Q&R SUTVOQ&R PS

W ;X9ZY T\[X]^^^?]_/`bP a

c Iedgf:hiGihddkj>lmIknbo?MpedkDGl8qrisDNJDGKtun8Kvf$hKxwy [ DCdmhFzI{uIF|F:j>lm}"I?p~EBIkdkKLpkDNJ?KLDNMFz}t DGFHI"JKDNMF€p
—ƒ˜ 6
œ oˆKvf$Imd?Misj rel="nofollow">KLDNMFzDCdVd"lmM5MKvfUžkn8KGf:ImpFrJ?KLDNMFŸDCdmh …M5M¡hq&qrp<M¢XDGlmhKLDNMF£K‚M Kvf$I¤MpkD¥…DGFrhi¦…&pkDNmo"j>F‡JKLDNMF§MF§KGf:IE¨xFrIl I?dF‡JKDNMF$n¬h i*MK8M7oDGFo?MpklmhKLDNMF€DCdi*MdkKj>pkDGF…KGf:IepFrJ?KLDNMF€JhF¯}"Il hX&Il Mp
O Q T ¶_ ¸ ² ± ² ²k· a w¼ [ ¹ ˜ 2šm6<º:3E0ˆ;>@>3548» T\[X]^^^?] 9<3¿<™.C4<.*À$23ÁeÀ/š ´X9<.C@ X½ ÂÄà ½ ¾ ÅÆǎÈLǎÉ5ʎËg̬É͆ÎkÏLǎРPQ Õ ÊŽÇ Õ ÈË7Ð Ñ ÒkÓ w¼ [ Ô Ñ w y [ ÒkÓ Ô w¼ [X] w [X] ½ ]^^^?] ½ ]  ]  ^  ^  ^ ] ½

Ö

±U²>³‹3.´&3 ˜ ™X3=k6;&9"4; WŒµ

‡‘’G“U”–•&­

×ÄØÙUØNÚ

ÛÜNÝNÞ:ßÝàâá

ã‡äåGæUç–è&é 2nd Eigenvector (n=19)

1

18th Eigenvector (n=19)

1 0.5

0.5 0

0

-0.5

-0.5

-1

0

0.2

0.4

0.6

0.8

-1

1

2nd Eigenvector (n=9)

1

0

0.4

0.6

0.8

1

0.5

0.5 0

0

-0.5

-0.5

-1

0.2

2nd Eigenvector (n=9)

1

-1

ê‰ë>ì>íîïñðNòó£ô5õö÷eë&øù&úûCìüý‡íþë&ÿ í
EðNòózî$õ¦ÿ ë>ì>íEë&øù&úûCì÷Xü 0

0.2

0.4

0.6

0.8

1

0

0.2

0.4

0.6

0.8

1

    !#" %$'&()+*,- ./"0'1/ 23 "%/4-'056&/7 8 "% ")90:9%;<7=!"%90>  ?@ AB('C D>
EFG%¤üHDG");"AIFI"):"0$J%*KDGCE1:% ÷&üH"%E7 :L%(M3$J%*N )/9%!7 
O

×ÄØÙUØC×

PRQ,S.SzÞ(T0U

ã‡äåGæUç–è%V

ø YJZ[Y*ë%\!ÿ ë>ì>íCZû*ø¤üþ
]øýíEúí^$úí
<íø!_íì[\¦í
YJYbûø¤÷&üR WNX / ë \§ÿ ë>ì>í
ëXø¤üýí5þ?ë&ÿ íBb$ûù1b:íúZÿ ë>ì>í
ûø¤÷&üR Wa` % c d/"%9C!5"e"%fg A)/9%F7 g"%f
O

hjilk

ã‡äåGæUç–è1y

mon/p@qp@r sutwvxp@r

zI7C {4"%/l-["G9% f%gM*@ "z2@7g")7E|"[ 
")gD>}~)7%*€ "03I%$J)!!"%%O

 ƒû ‚&íø„<…†G\¬íEë&ý(_]ûø[„[†mý!Z[‡ T%ˆ‡Üˆ‡àâá‡Þ(‰5݈‡à …† „†gŠ+‹ † „[…† …A† ‹ †NŒ ^$ú<ëFY*ëXø$ùF])_<ûë&øë1^‡íú]%_<ëXú(ðNÿ>]%_<úûJ$õ Žô 

’‘!“/”ƒ•
–)—ƒ—J‘F˜(™’š4›/œšGž/Ÿf ¡¢/£%¤J¥)Ÿ£)ž¦ ²e³1´ ±µ ²e³1´ ±µƒ¶ ¯

¨

§w¨ª©F«¬
¬¬«!­-®°¯ ± ¨

² ± ³1´ µ ¯ ±g· ² ± ³1´ µ¹¸

² ± ³1´ µJ¶ ¯º

ÂÄðÅÆoÇ

»7¼!½¾R¿HÀ1Á

È>É(ÊÆ%ËCÌÍÎaπÃ(ËCÐÒѹÊÉ°ÅÍfÃ(Î

Ó=ÔwÕÖšB×Aš
œØÕ×A”ƒ•Õ”ƒÙ1Ú–)Ú/Û“×AÙ1—ƒÙ1ÚÜ!–%ÕA”JÙFÚ|Ý>–%Õ×A”}•š
œzœA–%ÕA”ƒœØÔl‘ Þ ³ ± ³ ¨Nß Þ ³ ³ · ± ºà ÔlÙ1×ᖠ•Ù1Ú/œØÕA–1Ú!Õwß%˜CÕÖ/š’ÝGšÕAÖÙ(Û{”}œáœ–1”ƒÛ4ÕAÙ Ö/–0âFš€ÕÖ/š’ãä(å%æä(ç
æè¹éwäê°ëwå%è/ë@ì!å0çíu¦€jÖ”}œR”}œ,– “×AÙ1“7š×Õ‘dÕÖ/–)ÕwÔ–F•”ƒ—J”JÕA–)ÕšCœ€œØÙFÝ:š3Ù)Ô¹ÕÖš3ÕAÖšÙFךՔ}•–)—(•Ù1Ú-â1š
×ÜFšÚ/•šu“×AÙ!Ù1Ôœ,ÔlÙ1×áÝg›—JÕ”ƒÜ1×A”}Û ÝGšÕAÖÙ(Ûœ¦>î Ù%™’š
â1š×C˜Ù)ÕAÖš×{œ”JÝG“—ƒš×{•Ö/Ù1”}•š
œ
˜°—ƒ”JïFšÕAÖš>Ù1Ú/š
œ4•ÙFÚ/œØ”}Û(š
ךCÛo–)ð7Ù%â1šF˜•–1Ú –)—}œÙ™3Ù1×Aï>™’š
—J—R”JÚ5“×–1•Õ”}•š1¦ ñ Ù1Ú¹œØ”}Û(š×C˜7ÔlÙ1×B”ƒ—J—ƒ›/œØÕ×–%ÕA”JÙFÚo“›×A“7ÙFœš
œ
˜¹ÕAÖ/–%Õ4ÕAÖš:ò/ÚšGÜ1×A”ƒÛÄó Ö/–1œ4ô.¨öõ<“¹ÙF”JÚ!ÕAœ
˜–)Ú¹Û ÕÖ/š{•ÙF–1×AœšdÜ1×A”ƒÛ5÷1ó5Ö¹–1œjôo¨+ø:“7Ù1”ƒÚ!ÕAœ
¦@jÖšB“×AÙ1—ƒÙ1ÚÜ!–%ÕA”JÙFÚ<Ù1“7š×–%ÕAÙ1×z•–1Ú5ð¹šB™z×A”'ÕÕš
Ú –1œ ûü © 
ûü ² 
ü
ü ¯
ü ÷
ü ²
±
ü
ü û ² ü ©þ©
ü 

ü ² 
¯ © ü ³ ± ü
ü ÷ ¨ ü  Þ ³úù ± ³ ¨ ² 

¨ ù ³ ý ² ± 
÷ ü ©ÿ© ²   ²   ± ³ ÷ ²  ý  ý  © ²

jÖš4•ÙF××Aš
œ“¹ÙFÚ/Û(”ƒÚÜg×Aš
œØÕ×A”}•Õ”ƒÙ1Ú[Ù1“7š×–%ÕAÙ1×@”}œj›/œ›/–)—ƒ—ƒ‘>ךÔlš×AךCÛGÕÙ>–1œj–gÔl›—ƒ—°™’š
”JÜFÖFÕA”JÚ/Ü Ù1“7š×–%ÕAÙ1×z–)Ú¹Û<Ö/–Fœ3ÕÖš4ÔlÙF×Ý

Þ ³e ¨ ± ³ ù ³ 

©

û ý

©

÷ ©

© ÷ ©

© ÷



©

ûü ² ü ¯ ü ² ± ü ü  ü ²  ü ü ²  ² ² ý  ²











¨

û ² ¯ ý ² ± ²

¨ 

ù ± ³

± ³



³

F›/œÕ@ÔlÙ1ג•Ù1ÝG“—ƒšÕAšÚšCœœ€™3šeÜF”JâFšeÖš×AšzÕAÖš4–)—ƒš×Aï-”JÚ>“/×ÙF—JÙFÚÜF–)Õ”ƒÙ1ÚG–)Ú/Û:Ôl›—ƒ—7™’š
”JÜFÖFÕA”JÚ/Ü ×Aš
œØÕ×A”ƒ•Õ”ƒÙ1ÚÄÙ1“7š×–%ÕÙF×Aœd”ƒÚÄÕ™’Ù Û(”ƒÝGšÚ/œ”ƒÙ1Ú/œ
¦|Ó=Ô ù ± ³ ¨öÞ ³ ³ ù ³ ˜ÕAÖš<•ÙFÝG“¹ÙFÚšÚ!ÕAœ4Ù1Ô ±

õ

 "!$#&%

'

(*),+ .-0/

D E

1)2+ '

E E E E E E E

'3(*) 4657+ .-

/

8 '()2+ - <='1)>415+ -@? 9 : ;

'3(*),+ A-*465

/

8 '()2+ - <='1)2+ -*415"? 9 : ;

EF

E

8IHKJLNMOHQP$R 9 E E E

8

E E E E E

'

()>4657+ A-*465B/

!TSVU2WX

'()2+ -

/



/ZY 

C



8



8 :

()2+ '

\[]

1)>415+ -

<='

_^`WXa\bcW!d"!

<='

WU1

[e

1)2+ -*415

<='

() 465+ -*415 ? G

*



'31*),ij57+ A-`ij5k<='3(*),ij57+ A-*465k<='36) 465+ A-`ij5l<m'31)>4657+ A-*465

8@fhg 9 < :

'

1)2+ A-`ij5

()2+ A-*465

<m'

<='

(*),ij57+ .-

<='

6) 465+ A- ?n<

C '

1*),+ A-"o

8pHKJLNMOHQPqR 9

stvu

8

L

r

wXxly{zn|~}Il€A1y~‚Akx„ƒh|~|~…|

†c‡‰ˆ‹ŠŒŽ‰

‘V’~“"”–•—‹˜*“™`š;—‹”œ›"“™`šN—,ž”–— Ÿ¢¡&Ÿ"“•¤£(š2’~“˜*“`Ÿ7šv˜7—,™"šv—,”¥A¦p§©¨>žªK«O@•“—‹”Tšv˜œ•¡d™“`Ÿ¢”¬\“"˜7˜*˜"­ ®¬’~§šn— Ÿ£6§¢¨ žª¯«O@•“…°¬£n”²±˜7—,•³n£´“™«O“`Ÿl«œ•X“…°$ž”²±˜7—,•

9 ³1­µ ”Oš2’~“3š‹’d“`˜…’~§”c•ž£

š2’~“¶c˜*¨>ž”&±X§šv—,”·¶~“"˜*§žš.ž˜²—‹”dšN˜*œ•¡d™“7Ÿ_“`˜7˜ž˜ŸIª6’&—,™’¯™”dš.§ž—‹”–’&—±¤’_¦˜*“¸"¡d“`”¬™`—,“`Ÿ"­

(n=19)

1

grid h grid 2h

0.5 0 -0.5 -1

0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

0.8

0.9

1

Interpolation error

1 0.5 0 -0.5 -1

¹ !

0

0.1

0.2

0.3

0.4

0.5

[ "bcWº  [  W!$ ! [ *W{S{»T^` e]

0.6

> ]

‘–“I´`Ÿ`“"˜7¾ž“…¦˜*«¿š‹’d“(À6±¡{˜*“š2’~§žšv£(š‹’d“p«O@•“¢° ¨>”‰±‰§žšN—,ž”~£;š2’~“p«@•“3°

/ÃÁ

0.7

U2œ¼‰»~œ!d^`%q"*WX eœ½ /ZÁ

ž”±˜7—,•

9 ³n£d¶¬˜*œ•¡d™“7ŸI§A¦š.“`˜k¶c˜*Â

”±˜7—,•I³$¶c¨Ä¡&ŸI§©Ÿ«§ž¨‹¨Ä£…´"¡{š;—‹«¢¶~˜7šÅ§”dšv£’{—±@’I¦˜*“¸`¡d“"”c™"Æ

“"˜7˜*ž˜"­

Ç

È

ÉËÊÍÌÏÎÑЅÒ"Ó

Ô~ÕQÌ©Ð;Ð;Ö3× Ø(Ò"Ì©ÙÛÚÝÜ×kÞhÖ¢ßàÖ

ácâ‰ã‹äåæç

èé~êIë`ì{ë"í>ê

î ï*ð1ñ ò ê 

üZýþ`ÿ
 ñ 

ü)(



 *,+

ü>?

í
@ê



4 G

ü)F

ê.-

 < 

9`H39

<

7

G 7

:\

n

4 þ

9 5

a

H

}~%€

\

ò

é

 

é3204

489;:

þ6587

í ~ê! î

 

-3< 7

ò

ò

òXû

òdùú

ïð6ñ%$& ò ' î

ò#"

-

>= < ò ò

ï

ï

'

' é32hë

1 î ò



:W: þ

H

êœë 

ò



ò ú

<

ï*ð ò î

é

 

$&

 

ï*ð6ñ $& ò î



í ~ê! î

ATò ' M <

ïð ò

$ &

îlïð6ñ òN'

"

GYX

GTJ489 : GL484 : 9;GLH :\ \ 9 HLacb ` 4 H Z7[S U þ 7 5 þ6]þ!^_ þ 7 5T d þV : L G 4 \ G b aib I 9 \ H ai a 4 G 9 9;GLH þ 7 ^ _ þ 7 7 5T7 þ þSj   ] L7 ú ò Z HWmGYX \ 4 9;GIH 4 489;: 9;GLH H ac 4 G GIHKJ 9;GLH 7 þ þ" ÿ hI7 þ8587 7 U ÿ L7 e e

ÿ>þ

X6GL4:69`HhJgR

489

Hc:

žÿ>þ

9`HKJ

a

b þ

: þ

b 5

9`R 7

ÿ‹ÿ

W

þ"ÿ

P : G

4

7 þ 5"þ 9 Jy\ 6 X GI4 4 7 Sw U

W4

X þ

484 þ

H

þ67[Skþþ

þ

\

L7wS 7 GL4 4 U] þV587

4 þ85

bp9`R

þÛþ o

< 4

5 e

\

7v7 H

U

9;GLH

L7

GIH þ

ï

< A

7

žÿ‹ÿB7

4

ÿ>þ

b

ÿ

GIb 5

]

 ñ

: þ

þ65 5

G

þV5

bp48H 7

î

ò

é

 


éDCE



ò ú

î

Ug ò e f 9;GLH$î ï*ð1ñ 4 ^lk²þ þ e ò P :489 P : H b þ6]þ65VUXþ85 þU e þ 5 þ U^ " žÿ

\

7

5

JIb

ò

>= ò

ê 

  :6489 P

þTSkþVUþ65

9

-3< 

ò

ATò

L7.þ

üZýþ`ÿ
#

OQP GLR

<



ò

ïð6ñ%$& ò 1 î

ò

ú0/

î ò



<

GIHKJ ÿ

 ò

òBA

é

 


óõô÷öø

G

R!9`HhJqX6GI4 ÿ

: GL484

U X6GI4

G

þUg7

5

JI489

b ]

\

x5.7

7 9;GIH

þ

ÿI7

H3b þ

 9;GIHD9`H

:

òr

7Åþ

4

GL4 ò

ú

]5

GYX 7

\

A s ò  484 GI4

 þ²þ

ò

^ut87 G

ÿ>þ Ux5Q7

^

P ]

î

þ

4VGYXzai4

þ6{%

H U

afG

587d{|58]

G!G 7

\p9`HKJq9

7Åþ

4 {

^

ácâ‰ã‹äåæ

‚lƒ„†…T‡Yˆ

‰ êV

Š

í@Xê



‹ /

‹cŽyŽ

– é

 

? 

í

™ 

? 

í’ ~ê!—(îz˜







é—‹

2që 

í@XêT



 



ê



 é

/

C 

 ø



é ø

L’‘ û

M”“

 

é

ø%

L’‘

ødŒ û



ø% û

 Œ

û ûT•

ø

Œ I’‘

ê —wE

*



‹



û

M  š C



ù



é

ø%

CE

éf,~é32{ê! %œ| ê"íLž{ê!2mŸ Xë

I’‘ û

 

!›  3



< 

 ¡

&

ácâ‰ã‹äå

¢p£ – é

 

ínë 

é32

 


é

“

Solution

Error

1.5

1.5

u(x)-u h (x)

1

0.5

u(x)

1

0.5

0

0.5

0

0.5

1

1

1.5

1.5 0

0.2

0.4

0.6

0.8

1

0

0.2

0.4

0.6

0.8

1

¤i¥f¦,¦§8§ ¨L§.©ª.«¬­¦®y¯°¨%±.²[³w¨q´ ¨L«µ¯f¨I¶c¦¶3²dª·z¸ ¹
º¼»d½I¾’¿cÀ.³Á¥p©;´¥©ª#¬Â
¨L³j«¨!­ ¦Z¨I¶²Ã¥f¦ Ä

¶i¦«,¦6ª ¥Kŵ¬L¶i­¬m¨I³ÇÆYÈ®p²¶c¨L²ª¨m¨I³iÉm«¨!­ ¦¨L¶Ê²Ã¥f¦´ ¨!¬I§ ª¦Z«,¦6ª ¥fËÍ̶β`¥3¦¨L²Ã¥f¦§ ¥f¬I¶c­LÅÁ¸%¹
º¼»%Ïy½L¾’¿iÀQ©ª#¬Â¨I³j«,¨y­¦Z¨L¶0²`¥3¦ «¦8ª ¥fË,¤i¥f¦T¶i¦ Ò ² Ä

Ä

¶c¦.«¦6ª%¥¬L¶i­q¬¥p© Ðy¥Ñ«¨!­ ¦¨I¶²Ã¥f¦#´ ¨y¬L§ ª¦

ÐI®p§¦8ª żª ¥3¨L³”²Ã¥f¦Q¦Óx¨IÂÔ®p²[©;¨L¶Ñ¨Y±l²Ã¥f¦Q¦§8§ ¨L§Q¬L¶i­V²Ã¥f¦µª¨LÂÔ®p²[©;¨L¶Í²Ã¥h§ ¨I®Ðy¥

²Ã¥f¦#«®pÂÔ²[© Ч8©;­T¯c§¨!´ ¦6ªªË ÕcÖK×`ؼÙÎÚ’Û ÜlÝ`Þßàá âlã)ä

¹ Þ ßà åLÞ ¹æ º3¸æº

Þ ç3ßQè º ß.é,ß ¸ ç Solution

Error

1.5

1.5

u(x)-u h (x)

1

0.5

u(x)

1

0.5

0

0.5

1

1

1.5

1.5 0

ê

0

0.5

0.2

0.4

0.6

0.8

1

0

0.2

0.4

0.6

0.8

1

± ²ë¦§#²Ã¥h§ ¦ ¦ª8«,¨y¨L²Ã¥h©`¶hÐÍ©`²ë¦6§ ¬I²d©;¨I¶fª#²Ã¥f¦¥h© Ðx¥°¦6§8§ ¨I§«¨y­¦#¥3¬xª,È ¦ ¦¶Ê¦ÂÔ©`«Z©`¶c¬I²|¦ ­ËΤi¥f¦ ¨I³)¦6§8§ ¨I§.«,¨y­¦.¥3¬xªQª8®6쵦6§¦­gÂÔ©`²[²d¦Z´ ¥f¬L¶hÐ ¦yË ÕcÖK×`ؼÙÎÚfÚ

ÜlÝ`Þßàµí æ åà ¸ ßî à ¹ï

í æ àà ßyí8Þ ¹æ º0»Ãð¹ Þ ß!àåIÞ ¹
æº3¸ À

Solution

Error

1

1

0.5

0.5

u(x)-u h (x)

1.5

u(x)

1.5

0

0.5

1

1

1.5

1.5 0

ê

0

0.5

0.2

0.4

0.6

0.8

1

0

0.2

0.4

0.6

0.8

1

± ²ë¦§g²`¥3¦Í´¨y¬L§ ª¦Ð§8©;­Ê´ ¨L§8§ ¦ ´²d©;¨I¶fÅ.²`¥3¦m¨L³Ç«,¨y­¦g¥3¬xªÈ ¦ ¦6¶ñª8®3È8ª8²ë¬L¶²d©;¬IÂ`ÂÔò0§¦­I®3´ ¦­LÅ È®p²B¶c¨I³ÅBª¨L«¦g¶i¦6³)¥h© Ðx¥ó«,¨y­¦6ªV¥3¬LÓL¦ÍȦ ¦¶ó©`¶3²[§ ¨y­L®3´ ¦ ­Í²Ã¥h§ ¨I®Ðy¥Î²Ã¥f¦q©`¶3²ë¦§[¯f¨I¬L²[©;¨L¶

¯c§¨!´ ¦6ªªË

Ï!ô

õcöK÷`ø¼ùÎúfû ülý`þÿ
 

%þ  Iþþ ÿIþ

!#";ÿ
$% Iý'&)(*&
+ÿ,.-

Solution

Error

1.5

1.5

u(x)-u h (x)

1

0.5

u(x)

1

0.5

0

0.5

0

0.5

1

1

1.5

1.5 0

0.2

0.4

0.6

0.8

1

0

0.2

0.4

0.6

0.8

1

/102#035460879.:!2);0<>=>2? 9@ A2B2
?DC2B:>EFB2B:!G õcöK÷`ø¼ùÎúH ülý`þÿ
IJ þ ÿ
BLþ  Solution

Error

1.5

1.5

u(x)-u h (x)

1

0.5

u(x)

1

0.5

0

0.5

0

0.5

1

1

1.5

1.5 0

0.2

0.4

0.6

0.8

1

0

0.2

0.4

0.6

0.8

1

õcöK÷`ø¼ùÎúK ülý`þÿ
&
! ÿL$M&
! ÿ.&8þ  Solution

Error

1

1

0.5

0.5

u(x)-u h (x)

1.5

u(x)

1.5

0

0.5

0

0.5

1

1

1.5

1.5 0

0.2

0.4

0.6

0.8

1

0

0.2

0.4

0.6

0.8

1

õcöK÷`ø¼ùÎúN ülý`þÿ
 þ ÿ
BLþ O"Ãÿ$% IýP&)(*&)+
ÿQ!Solution

Error

1.5

1.5

u(x)-u h (x)

1

0.5

u(x)

1

0.5

0

0.5

0

0.5

1

1

1.5

1.5 0

0.2

0.4

0.6

0.8

1

,!,

0

0.2

0.4

0.6

0.8

1

RTS#U
SBS#VXWY>VZVXWS#U
SB[\] ^_[
`![
aSQbLUQY>]1Y>a\cd\>eUOVf\VgWSihkjBUlVZ[
`![)a SlmZY] ^nVgWYVomiYpBVfS
jQ\>]a5` Vrqi\M[
`![)aS)UBmiVgWSQU)\ase*Vrbo\>]8bLUtYasu\>UlVv[B\>]xw6S
jycdSB^!z { |d}X~€‚xƒ „T…*†‡ˆ ‰!ŠˆL‹Œ
Ž!!‘Š‰‘
Œ
‘*’
“k’ˆ ‰!‡ ‘#‰!Šˆ ‹ 4 2 0

Weighted Jacobi, ω=2/3 2level multigrid: ν1=2, ν2=2

4

10

log ||r||

L

2

2

6 8 10 12 14 0

5

10

15

20

25

# Iterations

”1WbLUih•ce*jSQbXaXaseUlVyjY>V–SlUQVgWStueWOp)Y>UlVfS
j[B\>]xw6S
jycdS)]1[S\!—
V–Y>bX] SB^MeUlbX]dcMVXWSQue*asVrb5cjlbo^ ˜ j \[BSB^>e*jS.z™”1WSš[B\>]xw6S
jycdS)]1[SMaS
w6S
aibLUnUW\>q']›w6S
jBUleUVgWSM]e*un—BS)j_\pnbXVfS
jYVybo\]U\>] VgWSJhi] SunSlUWz‚œ]1S_[BY>]žVgWbX]*ŸD\pVgWSn]e*un—S
jn\ pbXVfS
jYVybo\]UbX]žVXWSJhk]1SunSlUWdm#Y6U —BS)bX]cY ˜!˜ j\
¡bXunY>V–S)a5` ˜ j\ ˜ \>jlVrbo\>]1Y>a1Vf\tVgWS#Y>un\>e*]xV'\p[B\>u ˜ e*V–Y>Vrbo\>]1Y>a qi\>jŸbX]w>\>a5w6SB^!z

¢

£

°²±³

¤T¥B¦'§
¨T¥)©«ª­¬k§)®T¯ ´Dµ!¶8·²¸º¹»

{ |d}X~€‚¼

”1WSMVyqi\Mcjlbo^TU)[BWS
unS ˜ jSlU
S
]V–S^8Y!—B\>w>S_a SYw6S)Une*]jS)U
\>asw>SB^šVgWSM½
eS)UBVrbo\>]›\ ptW\q¾Vf\ U
\>asw>Sp)\>j_¿ÀBÁÂbX]™VgWSÃ[B\YjBU)S8unSlUWzĔ1WS‚Y>]UlqiS
j8[B\>unSlUOpBj\>uÅjSBY>a5bÆ)bX]cžVgWYVtVgWS ˜ j\!—
aS)uUOÇÁÈZÁÉËÊ Á Y] ^nÇÀÁÌ¿ÀBÁÉÎÍxÀBÁÏa\.\ŸDS¡Y![
Vras`šVXWSU
Y>unS_Y>]1^MWS
] [BSVgWS U
Y>unS ˜ j\.[SB^e*jS[BY]‚—BSOeU
SB^nV–\U
\>a5w6StVXWS)unz Ð ‘OŒ)Ñ*Œ
‡‘

ÈZÒÁ ÓÕÔ ÖØ×QÙ Á Ú È Ò Á Û Ê ÁÜ

ÝÂÞ ‘
‡Lß>à%á Ô †ˆâ‘.’JŽ!šÇÁ²ÈÁtÉãÊ Áä ˆ†åTˆ ˆ5†ˆ ß!‡º‰!…‘.’’vÈ ÒÁtæ ÝêéfëkìMí ôi‡ ’‘

Œ)Žß!Š’‘’†v‰!Šˆ ‹ºî*‰Ž†ŽDïðñ }Xò ó

Ò Ó Ý ø ŽŠŠ‘Œ)†JÈ Á Ý

È ÒÁ ÓÕÔçè “

À

Í À Á ¿ ÀÁ Ö

Ö

Á õ À Á Ú Ê Á#ö Ç Á È ÒӕÔçè Ü ×QÙ BÀ ÁxÚr÷ Û Í BÀ Á Ü “

çè É™È ÒÁ ӕÔçè ùúõ Á ÀBÁ ¿ÀBÁ“

ïð
ñ }Xò ó Þ ‘‡ ßà%á À †ˆ â‘’vŽšÇ Á È Á ÉãÊ Á ä ˆ5†å%ˆ ˆ†ˆL߇û‰…‘’’È ÒÁ Ó ü.ý

À

çè æ

È ÒÁ ÓÕÔ “

þ1ÿ
 
    !"$#%g ÿ
"&(')*++,.-/*0213.465702#1$8:9gÿ
;<*8=,>$-/* ?A@B-/C#<8$'D8;.8&4613FEG +H
JIK< +L'M02gÿN"$# 02#020O"$4QP&-BM
RGISUTVW02Xÿ/02#?X + ') 'M02424Y02#X138=9Z[
\  ]8?.8&4613GE \  /\  IKX\  IS^ \ +  O [_ E    `B,=T2 T þ1ÿ/0ab<*8=, 'M02424`,8$#Xc02#X-!-/#024M')*++", ÿDX ÿL,8="$*.G!BÿH["& Õ ' ÿQ0OB, ÿW8$02#X(')N'M02424F8$#X465d>Z8 eZfFghe 8:8&X ÿ/02#HPi02].*"&c0O8&#
N"$#<>!<*8&4j8$#QP"&]W>Z8&'M#kg ÿ
*+-/4602#HPlB8&46-/0O8$# @+.4j8&'ipB8$*qpB8&-/*(P&*0O>$BT

rCsOt<sOt

uv3wCxQykz/{=|Ovz
}O}O~

<€H2‚ƒ „
…

rCsOt<sa†

†‡‰ˆ(ŠYzyD‹U}cx

<€H2‚ƒ Œ/

ŽQj‘’

šl›œŸž
 
¡j¢(£+¤žX¥&¦ ’

_ ”“–•=• g “–—B—$)I™˜

   ¢¨ ’©< ž
 Xª
¥&¦+«

“ IS§

<€H2‚ƒ Œ­¬

®M¥&¦¥&¯ ’ ¢ ’ ¦bª ’:°X’  ª ’  X± ’ 6

ν1=1 ν2=1 ν =2 ν =1 1 2 ν =1 ν =2 1 2 ν1=2 ν2=2

7

log10||r||L

2

8

9

10

11

12 0.5

0.55

0.6

0.65

0.7

0.75

ω

0.8

0.85

0.9

0.95

1

²#78$*>B*G]8?>Z.B*?02#<[gÿ
;8+–c02?-/³,8$@B02#–"$0O8$#d8´pM
"&*"$.].*µ.8$¶BB*02B#Xn "$0O8$#W0aG8´pB#7#<>Z> T þ1ÿ
A·UP&-/*(ÿ
8&'A(g ÿ
G*B0O>&-"$4<*>$-,.c0O8"op].*G¸¹?,8& –4jB mno,.5Z,.4jb-Q02#HPlº?P&*0O>$BD T þ1ÿb·q# ÿ
"$G»&¼L½7»Z¼¾02#X].*1$"$4¿.T¾V7%.%Xÿ"$µg ÿ
8+–c02?-/À13"&46-l8´p;ÁÂ>].#<>378$#Ãg ÿ
#X-/!@.*L8´p[–*Bn?"$#<>83cn]8:8&X ÿ/02#HP¾02B*"&n 0O8$#
.TÃV7"&4¿.8d8Z@B.B*1$gÿ
"&pB8$* eZf I e IĘ&µg ÿ
N"&4 PH8$*02g ÿQÀ>8:B!#<8$(,8&#1$B*cPH=T Å "$#¾5Z8$-i+¶–4j"$02# Õ ' ÿQ5–Æ ˜:Ç

Ï`ÐZÑQÒÓ.Ô+ÕZÓ.ÑXÖBÓG×ص×%ٔÚ
ÑÖÛ+ÜjÐÑ7Ð&ÙAÕZÔ+ÜÞÝ ßjÓ:ÒZÓ:ßÞØàcØ×Zá?Ó[âÑÓá?Ó=Øoã<ä

È<ÉHÊ2ËÌ Í
Î

0 1 2

log ||r||

L2

3 4

10

5 6

2 levels 3 levels 4 levels 5 levels 6 levels

7 8 9 10

1

2

3

4

5

6

7

8

9

10

# Iterations

åæ7ç”èQéaê)ë<ìjí&ç)î)ï[êè
í&îðç”è
ïJñ+í&æò$ïBóôï.æ<ñïMõBí&óö÷ùøúoñ.ûZñ.ìjïêí´õGç2èï;ü?ý/ì6çcéþôZóéOÿ
&ì ôí&óé2ç2è/ü  ï+ïë<é2æQô !ñ+í&æ
êç $æç qæ–ïGü!ïBêèLío õ   lï ê.ïïGç2è $çAç”è
ï[üí$óïbôZóéOÿ3ê(î)ïñí$æêéOÿZï.ó ç”è
ïõ $êç]ï.ó7ç2èïiñí$æXò3ï.ócôHïBæ–ñ+ï $êlï Bë
ïñ.ç]ïÿ 7 ï $óïdýQêé2æH ô !#" $%$æ–' ÿ &(!)*

&,+-/.0

Ï`ÐZÑQÒÓ.Ô+ÕZÓ.ÑXÖBÓG×ص×%ٔÚ
ÑÖÛ+ÜjÐÑ7Ð&ÙAÕZÔ+ÜÞÝ ßjÓ:ÒZÓ:ßÞØàcØ×Zá?ÓÖ.Ð×&ÔØÓGá!Ó:ØãXä

È<ÉHÊ2ËÌ Í
Í

0 0.5

4 levels 5 levels 6 levels 7 levels

1

log10||r||L

2

1.5 2 2.5 3 3.5 4

0

2

4

6

8

10

12

14

16

18

20

1–èQéa2ê CôZý/óïWé2ì2ì6ýQêç3ó $çïêWç2èïõ Zñ.çJç”#è &çç”è
ïdñí$æXò3ï.ócôHïBæ–ñï73ó $çï7éaêdÿZï.ç]ï.óü?é2æ<ï5 ÿ 4Bûùç”è
ï ñ í6&óê.ïêçüïêè#87`êê.ïBæXçcé9&ì2ì6ûùç2è/éaê7í4ê.ï.óò3ïÿ:4ïè#$ò=éOí$ý/ó%ë–óí$ò$ïêWç2è$çç”è
ïWæXý/ü4ï.óWíoõ é2çïBó3&çcéOí&æ
ê(éaêGé2æ<ÿZïëïBæ–ÿZï.æç`íoõ ç”è
ï[æý/ü4ïBó[í´õbô&óéOÿ(ëí$é2æXçOê;Xë<óí&ò:éOÿï+ÿ!êý ê $óï[ñí$æ
êéOÿïBóï,ÿ ?$æ–ÿ(è
ï.æ<ñï(ç2èï(í&ò3ï.3ó $ì2@ ì $ì ôHí$óé2ç”èQü  è 3ê-ñí3êçîCè/éOñ+èNéaB ê A Dà Cä  ECí$óç”èQéaêJ3ï &ü ë–ìjFï `î)Gï &óïýQêé2æQ2ô :%" $H&æLNM

# Iterations

O!PQHR>SUTWVYX

È<ÉHÊ2ËÌ ÍZ

.[

\^]`_ab2]dceaFfGgihjcj[klafm,bam knf[goegpf;fFq9r[kla6st@_#agr9u[a;j]dq`v,agWhxwzyij[cj[klaFf{q|f{]g}j;g ~j;a~y ]db3m,]a2€j_HgWho]D_#a‚ g,b3ƒ„g,~]`_aj;g=m,b;f[aG€af3_#afs … †ˆ‡ˆ‰[‰‹Š ‡ˆ‰FŒ?[ŽB’‘}“•”–-—˜– ^™ ši› œ^ž>ŸY z¡£¢#¤i@¡¦¥i¥W§ ¨©0ª`«Y¬®­#¯ °B±#²i²3³µ´#¶³l²-· ¸µ¸U²3¹¶º²»#º[¼o½6½[½ ¾@g}hm,b'‚ aˆ_m,v,a5g,~k¿c:f[aÀÁg,~@agih']D_#aH€m q`~Âq9Àam,f'q`~v,g,kÃv,a;ÀÁq`~€*kÃ]dq¿Ä bFq9ÀsÅt@_#a q9Àa;m„gWhofF]m,bF]dq`~0Ä}‚¦q`]`_mÄ0g=g6À}q`~q`]dq9m,k@Ä af;f-_mpf{~@g,]Æra;a[~Hq`€-ekla[€a~]Ça;À}ca]st_*q|fq|f mj;j;g,€-ek¿q|fi_a;À}‚¦q`]D_]D_#a•h;*k`kY€*kÃ]dq¿Ä bFq9Àf[j_a€a6s

È

Š ɏ¦–ËÊŐÌ2Í{ΓƖ>ÐϋÉш“ÆÒ Ô Õ ´³µÖ3¹ ²¼3¹×#³µØ¶¼3³|Ù#ÖFÚpº6Û0±· ²3³µ¹´Ö6½ Ô Õ ¸µ¶º6Ü#¼;· ³|ØxÝG±#¸l²3³µ¶¼³µÙU½ Ô Þ ¹´ßº[¼¶º[´Øºx²3»#º6¹¼à½ Ô á ¹,âI²3¹
Ù*º=· ¸UâƳl²3»¹ ²»#º[¼Æ¹×º6¼· ²3¹¼Ö6½

t@_#amr;g vpam,ba
f[g €a}]Çg;eq9jf
‚J_q9j;_ãm,b3aä]D_#a
f;r9u[a;j]2gWhämj[]Nq`v,aäb3aFf[a;m,bj_Ám,~@Àˆ‚J_q9j;_ ,m bara[cg,~@À]D_#aGf[j;gea
gWhG]`_*q|fGkla;j[]N*ba6s{å~äe#m bF]Nq9j[*klm,b;æ-m k Äa=r[b3m,q9j
€G*kÃ]dq¿ÄbFq9À„€a[]`_g6À,f m,ba-Ä0m,q`~q`~Ä
q`~@jbam,f[aÀoe#g;e*klm bFq`]Nc„ra;j;m,f[aGgWho]D_#a[q`b eg,]a[~]dq9m,kJm,f{mÄa[~a[b3m,k£aç[m,]dq9g,~ f[g,kÃv,a[b[s{èa[b3aæé‚ a_m,v,a
mpf;fF*€aÀä]D_#m ]-g,*b€m,]dbFq¿ê‚ m,fG]D_#ab3afF*kÃ]>gWhÀ,q|f[jba]dqlë[q`~0ÄHm ì¦íoî g,~ïm}Äq`v,a~5Ä bFq9À®gWh
fFqlë6að£sˆåF~‹]D_q|ff[a[]N]dq`~0Ä æ-q`]oq|ffF]Nb3m q¿Ä=_*] hg bF‚ m,b3À®]gHÀa9ñ ~a ]D_#a5j;g6m b;f[aˆÄ bFq9ÀHe@b3gr[kla[€f[sòå9hg,*bóa;çm ]Nq9g ~ÂfFcpfF]a€Gæô>õ÷öùøéæ„Àg=aFf~g,]}j;g €a h;b3g €ú]D_#aÀ q|fj[b3a[]Nqlë=m ]Nq9g ~5gihGm ì¦íoî æ^g b{]`_aÀ,q|f[jba]dqlë=m,]dq9g,~äeb3g=j;aFf;f{q|f{~@g,]>m,v,m,q`klmrkµaFæ ]D_#a-ebg6j;a;À,*b3afÄq`vpa[~5q`~5]D_q|fkla;j[]N*bajm ~ó~g,]-r;af[a;ÀÀ,q`b3a;j[]NkÃcs í a;m kÃq`~0Ä‚¦q`]D_HfFj_ Ä0a~@ab3m kUfFcpfF]Ça[€f;æ ‚¦q`]D_#g *]•baç[*q`bFq`~0Ä}q`~6hg bF€m ]Nq9g ~:mr;g,*]•]`_aÀ q|f[jb3a[]dqlë6m ]Nq9g ~
e@b3g=jaf;f;æ q|f2]`_a2fFr9u[a;j]^gWhm k Äa=r[b3m,q9j€*kÃ]dq¿Ä bFq9À
€a[]`_g6À,f[s û=ü

¨©0ª`«Y¬®­#Ó

ý@þ;þxÿ

#þ þ


    
   

þ þ 3þ ;þ

  "!  

ÿ

þ


#$%'&(*)+,- /.0. 1324*5,6

Bý

@L 1

ý

789

MNE$O<"N.P N Q7R16SRTUTPTP

V0W

:*;#-?@! AB1DCEF( *G/. @IH7 JK