<div dir="ltr">Hi Jed<div>That's a great idea. Do you have any ex* to start with? Or suggesting me a starting point.</div><div class="gmail_extra"><br></div><div class="gmail_extra">Thanks<br clear="all"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Giang</div></div></div>
<br><div class="gmail_quote">On Tue, Oct 18, 2016 at 3:49 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="">Hoang Giang Bui <<a href="mailto:hgbk2008@gmail.com">hgbk2008@gmail.com</a>> writes:<br>
<br>
> Dear PETSc folks,<br>
><br>
> While searching literature on the preconditioner for contact/mesh tying<br>
> problem, I saw the paper by Dr. Adams "Algebraic multigrid methods for<br>
> constrained linear systems with applications to contact problems in solid<br>
> mechanics, NLAA, 2004". Given the promising aspects the paper has shown for<br>
> constrained linear system, I wonder if some code's also available in PETSc<br>
> for testing/further extension?<br>
<br>
</span>This particular algorithm is not available within the PETSc library.<br>
Mark might have some code around, but I doubt it's maintained or easy to<br>
tinker with. You should be able to implement the algorithm fairly<br>
quickly using PETSc. It would make a great example if you're willing to<br>
contribute. And we can advise if you decide to write such an example.<br>
</blockquote></div><br></div></div>