<div dir="ltr">On Wed, Jul 3, 2013 at 9:57 AM, Tabrez Ali <span dir="ltr"><<a href="mailto:stali@geology.wisc.edu" target="_blank">stali@geology.wisc.edu</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Hello<br>
<br>
Is there a PETSc example that really demonstrates the advantage of using a multigrid solver (in that that the number of iterations are relatively constant as the problem size increases)?<br>
<br>
The command line options used to run it would also be helpful.<br></blockquote><div><br></div><div style>You can run both GMG and GAMG on any structured example to compare. For example,</div><div style><br></div><div style>
./ex5 -da_grid_x 21 -da_grid_y 21 -da_refine 6 -ksp_rtol 1.0e-9 -pc_type mg -pc_mg_levels 4 -snes_monitor -snes_view</div><div style><br></div><div style>or</div><div style><br></div><div style> ./ex5 -da_grid_x 21 -da_grid_y 21 -da_refine 6 -ksp_rtol 1.0e-9 -pc_type gamg -snes_monitor -snes_view</div>
<div style><br></div><div style>and just adjust the da_refine to see the constant number of iterates. Same thing in 3D</div><div style><br></div><div style> ./ex48 -M 5 -N 5 -da_refine 5 -ksp_rtol 1.0e-9 -thi_mat_type baij -pc_type mg -pc_mg_levels 4 -snes_monitor -snes_view -log_summary</div>
<div style><br></div><div style>I use these examples in my latest tutorial, which I will put up today.</div><div style><br></div><div style> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
Thanks in advance.<span class=""><font color="#888888"><br>
<br>
Tabrez<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>