<div dir="ltr">On Tue, Aug 27, 2013 at 2:45 PM, 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>
What is the proper way to use GAMG on a vanilla 3D linear elasticity problem. Should I use<br>
<br>
-pc_type gamg -pc_gamg_type agg -pc_gamg_agg_nsmooths 1<br></blockquote><div><br></div><div>This is fine.</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">

or<br>
<br>
-pc_type fieldsplit -pc_fieldsplit_block_size 3 -fieldsplit_pc_type gamg -fieldsplit_pc_gamg_type agg -fieldsplit_pc_gamg_agg_<u></u>nsmooths 1<br>
<br>
Do these options even make sense? With the second set of options the % increase in number of iterations with increasing problem size is lower than the first but not optimal.<br>
<br>
Also, ksp/ksp/examples/ex56 performs much better in that the number of iterations remain more or less constant unlike what I see with my own problem. What am I doing wrong?<br></blockquote><div><br></div><div>You need to give it a good near null space. Usually, we use the 3 translational and 3 rotations modes that are the</div>
<div>null space of the elastic operator. MatNullSpaceCreateRigidBody() is an example of making them. PyLith does this by default :)</div><div><br></div><div>  Thanks,</div><div><br></div><div>      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">

The output of -ksp_view for the two set of options used is attached.<br>
<br>
Thanks in advance.<span class=""><font color="#888888"><br>
<br>
Tabrez<br>
<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>