<div dir="ltr">Thanks Jed.<div><br></div><div>So the nullspace would be initialized with an array of 3 petsc vectors: (u,v)= (0,1), (1,0), and (-y,x), correct?</div><div><br></div><div style>And also to be sure -- this is usefuil only for multigrid preconditioners, yes? </div><div style><br></div><div style>Thanks</div><div style>Dan</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Feb 23, 2015 at 2:11 AM, Jed Brown <span dir="ltr"><<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">Barry Smith <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>> writes:<br>
>> I will try gamg as I know it can reduce the number of CG iterations required. (I'm guessing you mean algebraic, not geometric?)<br>
><br>
> By default GAMG is an algebraic multigrid preconditioner. Look at its documentation at <a href="http://www.mcs.anl.gov/petsc/petsc-master/docs/index.html" target="_blank">http://www.mcs.anl.gov/petsc/petsc-master/docs/index.html</a> it will be a bit better than the older documentation. The documentation for GAMG is still pretty thin so feel free to ask questions.<br>
<br>
</span>You might want to call MatSetBlockSize and MatSetNearNullSpace to define<br>
both translation and rotation modes.  This is most relevant for large<br>
ice shelves.<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><br>Daniel Goldberg, PhD<br>Lecturer in Glaciology<div>School of Geosciences, University of Edinburgh<br><span style="color:rgb(0,0,0);font-family:arial,helvetica,clean,sans-serif;line-height:15px">Geography Building, Drummond Street, Edinburgh EH8 9XP</span><br><br><br>em: <a href="mailto:dgoldber@mit.edu" target="_blank">D</a><a href="mailto:an.Goldberg@ed.ac.uk" target="_blank">an.Goldberg@ed.ac.uk</a><br>web: <a href="http://ocean.mit.edu/~dgoldberg" target="_blank">http://ocean.mit.edu/~dgoldberg</a></div></div></div>
</div>