<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Mar 7, 2015 at 6: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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">OK, thanks for letting me know.<div><br></div><div>I've tried GAMG and ML with MatNullSpaceCreateRigidBody and those both work well for me. (I tried ML after Jed pointed out that ML also uses the near nullspace.)</div><div><br></div><div>I have also tried hypre but the convergence hasn't been as good for the elasticity models I've been using, though I may not have been setting the hypre options in an optimal way for elasticity.</div></div></blockquote><div><br></div><div>As Jed said, SA should be better for elasticity but it is worth testing.<br><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc 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><div>Do you mean Poisson ratio?  Nearly incompressible is hard.  bigger smoothers, like ASM, can help.<br></div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>all iterative solvers I've tried fare poorly. I guess it's because the highly contrasting stiffnesses give a large condition number. Is there a standard way to get the preconditioner compensate for this?</div></div></blockquote><div><br> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><br></div><div>Thanks,</div><div>David</div><div><div class="h5"><div><br></div><div><br></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></div><br></div></div>