Richardson Cascadic Multigrid Method for 2D Poisson Equation Based on a Fourth Order Compact Scheme

<table>Comparison of the maximum error <svg height="20.237499" id="M163" style="vertical-align:-3.39066pt" version="1.1" viewbox="0 0 70.099998 20.237499" width="70.099998" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,15.95)"><path d="M162 -163h-60v866h60v-866zM309 -163h-60v866h60v-866z" id="x2016"></path></g>
<g transform="matrix(1.25,0,0,-1.25,90,90)">
<path d="m -66.36,67.71 7.32,0" style="fill:none;stroke:#000000;stroke-width:0.68000001;stroke-linecap:butt;stroke-linejoin:miter;stroke-miterlimit:10;stroke-opacity:1;stroke-dasharray:none"></path>
</g>
<g transform="matrix(.017,-0,0,-.017,7.05,15.95)"><path d="M515 96q-37 -44 -84.5 -76t-71.5 -32q-12 0 -18 7t-6 34t10 75l19 89h-2q-87 -113 -181 -176q-43 -29 -83 -29q-46 0 -46 63q0 31 9 67l52 222q7 36 -1 36q-10 0 -76 -50l-13 24q47 45 92 71.5t67 26.5q18 0 21 -16.5t-8 -65.5l-56 -242q-16 -67 13 -67q41 0 114 74
t114 146l32 145l74 26h11q-52 -199 -80 -347q-8 -39 6 -39q7 0 32 17.5t47 39.5z" id="x1D462"></path></g>
<g transform="matrix(.012,-0,0,-.012,16.2,7.8)"><path d="M559 163q-23 -66 -68 -163h-474l6 26q62 4 79.5 19.5t28.5 75.5l78 409q7 35 8.5 49t-8 25t-24 13t-51.5 5l5 28h266l-6 -28q-65 -5 -79.5 -18t-25.5 -74l-76 -406q-10 -57 14 -75q12 -13 96 -13q93 0 126 29q41 40 76 109z" id="x1D43F"></path></g>
<g transform="matrix(.017,-0,0,-.017,27.387,15.95)"><path d="M535 230h-483v50h483v-50z" id="x2212"></path></g><g transform="matrix(.017,-0,0,-.017,41.139,15.95)"><path d="M515 96q-37 -44 -84.5 -76t-71.5 -32q-12 0 -18 7t-6 34t10 75l19 89h-2q-87 -113 -181 -176q-43 -29 -83 -29q-46 0 -46 63q0 31 9 67l52 222q7 36 -1 36q-10 0 -76 -50l-13 24q47 45 92 71.5t67 26.5q18 0 21 -16.5t-8 -65.5l-56 -242q-16 -67 13 -67q41 0 114 74
t114 146l32 145l74 26h11q-52 -199 -80 -347q-8 -39 6 -39q7 0 32 17.5t47 39.5z" id="x1D462"></path></g><g transform="matrix(.017,-0,0,-.017,50.285,15.95)"><path d="M162 -163h-60v866h60v-866zM309 -163h-60v866h60v-866z" id="x2016"></path></g>
<g transform="matrix(.012,-0,0,-.012,57.275,20.038)"><path d="M983 225q0 -112 -67 -174.5t-150 -62.5q-91 0 -154.5 43.5t-113.5 129.5q-49 -85 -104 -129t-138 -44q-98 0 -158.5 66t-60.5 154q0 59 21 106.5t54.5 75.5t70.5 43t73 15q90 0 152.5 -43.5t112.5 -128.5q48 84 104.5 128t140.5 44q93 0 155 -65t62 -158zM478 196
q-27 49 -47 80t-50 67t-64 54t-73 18q-48 0 -81.5 -47t-33.5 -128q0 -96 37.5 -157.5t99.5 -61.5q68 0 117.5 47t94.5 128zM889 204q0 91 -35.5 151t-99.5 60q-68 0 -119 -47t-95 -127q27 -49 47 -80.5t50 -67.5t65 -54t74 -18q113 0 113 183z" id="x221E"></path></g>
</svg> and cpu time for Example <a href="https://static.hindawi.com/articles/jam/volume-2014/490540/figures/#ex2">5</a> with <svg height="11.375" id="M164" style="vertical-align:-0.1638pt" version="1.1" viewbox="0 0 37.412498 11.375" width="37.412498" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,11.112)"><path d="M559 163q-23 -66 -68 -163h-474l6 26q62 4 79.5 19.5t28.5 75.5l78 409q7 35 8.5 49t-8 25t-24 13t-51.5 5l5 28h266l-6 -28q-65 -5 -79.5 -18t-25.5 -74l-76 -406q-10 -57 14 -75q12 -13 96 -13q93 0 126 29q41 40 76 109z" id="x1D43F"></path></g><g transform="matrix(.017,-0,0,-.017,14.477,11.112)"><path d="M535 323h-483v50h483v-50zM535 138h-483v50h483v-50z" id="x3D"></path></g><g transform="matrix(.017,-0,0,-.017,29.181,11.112)"><path d="M285 378v-2q65 -13 102 -54.5t37 -97.5q0 -57 -30.5 -104.5t-74 -75t-85.5 -42t-72 -14.5q-31 0 -59.5 11t-40.5 23q-19 18 -16 36q1 16 23 33q13 10 24 0q58 -51 124 -51q55 0 88 40t33 112q0 64 -39 96.5t-88 32.5q-29 0 -64 -11l-6 29q77 25 118 57.5t41 84.5
q0 45 -26.5 69.5t-68.5 24.5q-67 0 -120 -79l-20 20l43 63q51 56 127 56h1q66 0 107 -37t41 -95q0 -42 -31 -71q-22 -23 -68 -54z" id="x33"></path></g>
</svg>, taking step lengths <svg height="16.3125" id="M165" style="vertical-align:-3.276pt" version="1.1" viewbox="0 0 64.3125 16.3125" width="64.3125" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,12.162)"><path d="M499 88q-41 -46 -83.5 -73t-68.5 -27q-41 0 -19 98l52 225q10 43 8 59.5t-16 16.5q-37 0 -112 -65t-117 -135q-14 -42 -38 -179q-36 -7 -75 -20l-7 8l127 594q8 40 2.5 49t-35.5 9h-35l5 25q36 3 71.5 13t56.5 18t25 8q11 0 3 -38l-88 -417h2q77 102 195 167q44 24 79 24
q61 0 30 -130l-52 -220q-8 -33 4 -33q8 0 31 14.5t43 33.5z" id="x210E"></path></g>
<g transform="matrix(.012,-0,0,-.012,8.837,16.25)"><path d="M559 163q-23 -66 -68 -163h-474l6 26q62 4 79.5 19.5t28.5 75.5l78 409q7 35 8.5 49t-8 25t-24 13t-51.5 5l5 28h266l-6 -28q-65 -5 -79.5 -18t-25.5 -74l-76 -406q-10 -57 14 -75q12 -13 96 -13q93 0 126 29q41 40 76 109z" id="x1D43F"></path></g>
<g transform="matrix(.017,-0,0,-.017,20.975,12.162)"><path d="M535 323h-483v50h483v-50zM535 138h-483v50h483v-50z" id="x3D"></path></g><g transform="matrix(.017,-0,0,-.017,35.662,12.162)"><path d="M384 0h-275v27q67 5 81.5 18.5t14.5 68.5v385q0 38 -7.5 47.5t-40.5 10.5l-48 2v24q85 15 178 52v-521q0 -55 14.5 -68.5t82.5 -18.5v-27z" id="x31"></path></g><g transform="matrix(.017,-0,0,-.017,43.821,12.162)"><path d="M368 703l-264 -866h-60l265 866h59z" id="x2F"></path></g><g transform="matrix(.017,-0,0,-.017,50.825,12.162)"><path d="M766 88q-40 -45 -83.5 -72.5t-63.5 -27.5q-39 0 -16 103l49 224q15 68 -10 68q-36 0 -112.5 -78.5t-105.5 -133.5q-27 -101 -38 -165q-35 -4 -77 -18l-6 6q42 153 70 288q12 55 10 78t-19 23q-38 0 -114 -80.5t-101 -131.5q-24 -68 -41 -165q-36 -3 -76 -18l-8 6
q56 202 87 347q7 33 -3 33q-8 0 -31.5 -17.5t-41.5 -35.5l-12 28q42 45 84 72t63 27q44 0 10 -124l-20 -75h2q92 121 193 178q36 21 64 21q62 0 28 -159l-8 -37h2q50 67 103.5 112t94.5 65q39 19 62 19q60 0 23 -156l-45 -189q-9 -39 2 -39q17 0 71 49z" id="x1D45A"></path></g>
</svg>, <svg height="11.225" id="M166" style="vertical-align:-0.1638pt" version="1.1" viewbox="0 0 49.275002 11.225" width="49.275002" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,10.963)"><path d="M766 88q-40 -45 -83.5 -72.5t-63.5 -27.5q-39 0 -16 103l49 224q15 68 -10 68q-36 0 -112.5 -78.5t-105.5 -133.5q-27 -101 -38 -165q-35 -4 -77 -18l-6 6q42 153 70 288q12 55 10 78t-19 23q-38 0 -114 -80.5t-101 -131.5q-24 -68 -41 -165q-36 -3 -76 -18l-8 6
q56 202 87 347q7 33 -3 33q-8 0 -31.5 -17.5t-41.5 -35.5l-12 28q42 45 84 72t63 27q44 0 10 -124l-20 -75h2q92 121 193 178q36 21 64 21q62 0 28 -159l-8 -37h2q50 67 103.5 112t94.5 65q39 19 62 19q60 0 23 -156l-45 -189q-9 -39 2 -39q17 0 71 49z" id="x1D45A"></path></g><g transform="matrix(.017,-0,0,-.017,18.183,10.963)"><path d="M535 323h-483v50h483v-50zM535 138h-483v50h483v-50z" id="x3D"></path></g><g transform="matrix(.017,-0,0,-.017,32.887,10.963)"><path d="M137 343l67 33q37 17 63 17q79 0 129.5 -53t50.5 -131q0 -92 -58 -156.5t-147 -64.5t-147 68t-58 182q0 63 17 119t43 95.5t61.5 72t69 52t67.5 31.5q62 22 128 33l6 -32q-56 -11 -108 -35q-149 -71 -184 -231zM227 337q-47 0 -95 -27q-6 -23 -6 -70q0 -93 36 -155.5
t96 -62.5q53 0 78 45.5t25 105.5q0 68 -35 116t-99 48z" id="x36"></path></g><g transform="matrix(.017,-0,0,-.017,41.046,10.963)"><path d="M456 178h-96v-72q0 -51 12.5 -62.5t72.5 -16.5v-27h-256v27q65 5 78 17t13 62v72h-260v28q182 271 300 426h40v-407h96v-47zM280 225v295h-2q-107 -148 -196 -295h198z" id="x34"></path></g>
</svg>, <svg height="11.125" id="M167" style="vertical-align:-0.1638pt" version="1.1" viewbox="0 0 24.6 11.125" width="24.6" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,10.862)"><path d="M384 0h-275v27q67 5 81.5 18.5t14.5 68.5v385q0 38 -7.5 47.5t-40.5 10.5l-48 2v24q85 15 178 52v-521q0 -55 14.5 -68.5t82.5 -18.5v-27z" id="x31"></path></g><g transform="matrix(.017,-0,0,-.017,8.222,10.862)"><path d="M412 140l28 -9q0 -2 -35 -131h-373v23q112 112 161 170q59 70 92 127t33 115q0 63 -31 98t-86 35q-75 0 -137 -93l-22 20l57 81q55 59 135 59q69 0 118.5 -46.5t49.5 -122.5q0 -62 -29.5 -114t-102.5 -130l-141 -149h186q42 0 58.5 10.5t38.5 56.5z" id="x32"></path></g><g transform="matrix(.017,-0,0,-.017,16.381,10.862)"><path d="M249 635q70 0 116 -43t46 -105q0 -46 -28 -80q-22 -25 -80 -64q62 -35 97 -75t35 -99q0 -81 -63 -131t-135 -50q-83 0 -137.5 45.5t-54.5 123.5q0 52 45 95q29 28 89 64q-109 62 -109 155q0 66 50.5 115t128.5 49zM238 603q-42 0 -67.5 -31t-25.5 -72q0 -50 32.5 -79.5
t98.5 -62.5q61 48 61 124q0 59 -29.5 90t-69.5 31zM248 20q46 0 76.5 33.5t30.5 89.5q0 50 -39 85.5t-110 71.5q-81 -54 -81 -137q0 -67 35.5 -105t87.5 -38z" id="x38"></path></g>
</svg>, <svg height="11.225" id="M168" style="vertical-align:-0.1638pt" version="1.1" viewbox="0 0 24.6 11.225" width="24.6" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,10.963)"><path d="M412 140l28 -9q0 -2 -35 -131h-373v23q112 112 161 170q59 70 92 127t33 115q0 63 -31 98t-86 35q-75 0 -137 -93l-22 20l57 81q55 59 135 59q69 0 118.5 -46.5t49.5 -122.5q0 -62 -29.5 -114t-102.5 -130l-141 -149h186q42 0 58.5 10.5t38.5 56.5z" id="x32"></path></g><g transform="matrix(.017,-0,0,-.017,8.222,10.963)"><path d="M153 550l-26 -186q79 31 111 31q90 0 141.5 -51t51.5 -119q0 -93 -89 -166q-85 -69 -173 -71q-32 0 -61.5 11.5t-41.5 23.5q-18 17 -17 34q2 16 22 33q14 9 26 -1q61 -50 124 -50q60 0 93 43.5t33 104.5q0 69 -41.5 110t-121.5 41q-53 0 -102 -20l38 305h286l6 -8
l-26 -65h-233z" id="x35"></path></g><g transform="matrix(.017,-0,0,-.017,16.381,10.963)"><path d="M137 343l67 33q37 17 63 17q79 0 129.5 -53t50.5 -131q0 -92 -58 -156.5t-147 -64.5t-147 68t-58 182q0 63 17 119t43 95.5t61.5 72t69 52t67.5 31.5q62 22 128 33l6 -32q-56 -11 -108 -35q-149 -71 -184 -231zM227 337q-47 0 -95 -27q-6 -23 -6 -70q0 -93 36 -155.5
t96 -62.5q53 0 78 45.5t25 105.5q0 68 -35 116t-99 48z" id="x36"></path></g>
</svg>, and <svg height="11.125" id="M169" style="vertical-align:-0.1638pt" version="1.1" viewbox="0 0 24.6 11.125" width="24.6" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,10.862)"><path d="M153 550l-26 -186q79 31 111 31q90 0 141.5 -51t51.5 -119q0 -93 -89 -166q-85 -69 -173 -71q-32 0 -61.5 11.5t-41.5 23.5q-18 17 -17 34q2 16 22 33q14 9 26 -1q61 -50 124 -50q60 0 93 43.5t33 104.5q0 69 -41.5 110t-121.5 41q-53 0 -102 -20l38 305h286l6 -8
l-26 -65h-233z" id="x35"></path></g><g transform="matrix(.017,-0,0,-.017,8.222,10.862)"><path d="M384 0h-275v27q67 5 81.5 18.5t14.5 68.5v385q0 38 -7.5 47.5t-40.5 10.5l-48 2v24q85 15 178 52v-521q0 -55 14.5 -68.5t82.5 -18.5v-27z" id="x31"></path></g><g transform="matrix(.017,-0,0,-.017,16.381,10.862)"><path d="M412 140l28 -9q0 -2 -35 -131h-373v23q112 112 161 170q59 70 92 127t33 115q0 63 -31 98t-86 35q-75 0 -137 -93l-22 20l57 81q55 59 135 59q69 0 118.5 -46.5t49.5 -122.5q0 -62 -29.5 -114t-102.5 -130l-141 -149h186q42 0 58.5 10.5t38.5 56.5z" id="x32"></path></g>
</svg>, respectively.</table>

Journal of Applied Mathematics

fig4

Figure 4

Figure 4: Richardson Cascadic Multigrid Method for 2D Poisson Equation Based on a Fourth Order Compact Scheme 