<div dir="ltr">Penalty methods are hell on solvers.  There is a good change that GAMG will do the right thing with a finite threshold. and keep this stiff section separate.  You might want to start wit a small number of levels and a large coarse grid and get it working.  Then add more levels and see if it breaks.<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Mar 8, 2015 at 11:38 AM, David Knezevic <span dir="ltr"><<a href="mailto:david.knezevic@akselos.com" target="_blank">david.knezevic@akselos.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span class="">On Sun, Mar 8, 2015 at 10:31 AM, Mark Adams <span dir="ltr"><<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>></span> wrote:<br><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"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div><br><br></div><span><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"><div dir="ltr"><div><br></div><div>Another situation I've been meaning to ask about: If the elasticity model includes some "near rigid" regions (very high Young's modulus) </div></div></blockquote><div><br></div></span><div>Do you mean Poisson ratio?  </div></div></div></div></blockquote><div><br></div><div><br></div></span><div>No, just high Young's modulus in subregions. Engineers sometimes implement a rotation boundary condition about a node by including rigid (or, in practice, very high Young's modulus) parts in a model. For example, imagine including a rigid pyramid in a model, in which the tip of the pyramid is is clamped so that the base of the pyramid (which is connected to a surface in the model) can rotate about the tip.<br></div><div><br></div><div>This works fine with a direct solver, so I was curious about getting a good iterative solver for this case too.</div><span class=""><div><br></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"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div>Nearly incompressible is hard.  bigger smoothers, like ASM, can help.</div></div></div></div></blockquote><div><br></div></span><div>Agreed. But the Poisson ratio is still just 0.3 or so in this context, so locking and near incompressibility aren't an issue here.</div><span class=""><div><br></div><div>Thanks,<br>David</div><div><br></div><div><br><br><br><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span><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"><div dir="ltr"><div><div><div><br></div><div class="gmail_extra"><div class="gmail_quote">On Sat, Mar 7, 2015 at 5:57 PM, Mark Adams <span dir="ltr"><<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>></span> wrote:<br><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"><div dir="ltr"><div>FYI, stüben used classical AMG for elasticity but he has articulated his code for elasticity more than Hypre as I understand it.  Hypre can work OK for elasticity in my experience.  Its worth a try.<span><font color="#888888"><br><br></font></span></div><span><font color="#888888">Mark<br></font></span></div><div><div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Mar 5, 2015 at 5:27 PM, David Knezevic <span dir="ltr"><<a href="mailto:david.knezevic@akselos.com" target="_blank">david.knezevic@akselos.com</a>></span> wrote:<br><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"><div dir="ltr">OK, got it, thanks!<span><font color="#888888"><div><br></div><div>David</div></font></span><div><div><div><br></div><div><br></div><div class="gmail_extra"><div class="gmail_quote">On Thu, Mar 5, 2015 at 5:08 PM, 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:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><span>David Knezevic <<a href="mailto:david.knezevic@akselos.com" target="_blank">david.knezevic@akselos.com</a>> writes:<br>
> I was just wondering if its possible to achieve the same sort of thing with<br>
> other AMG solvers (e.g. BoomerAMG)? I assume that MatSetNearNullSpace does<br>
> nothing for external solvers like hypre, right?<br>
<br>
</span>It is used by ML (smoothed aggregation), but not BoomerAMG (classical<br>
AMG) which uses an algorithm that doesn't have a natural place for such<br>
information.  To my knowledge, classical AMG is not widely used for<br>
elasticity.  It is very robust for M-matrices.<br>
</blockquote></div><br></div></div></div></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br></div></div></div></div>
</blockquote></span></div><br></div></div>
</div></span></div><br></div></div>
</blockquote></div><br></div>