<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Jun 1, 2015 at 11:24 PM, Eduardo <span dir="ltr"><<a href="mailto:erocha.ssa@gmail.com" target="_blank">erocha.ssa@gmail.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">I am solving a FEM solid mechanics linear elasticity model, for now the only problem is the mesh that has needle-shaped and very flat elements. Have you any suggestion of preconditioner (and references).<div></div></div></blockquote><div><br></div><div>The problem here is your discretization. WIth quasi-regular elements, -pc_type gamg works fine. However, with flat elements, your FEM</div><div>basis becomes very ill-conditioned since the normal basis functions are almost linearly dependent. I think the best use of time here is</div><div>to investigate better discretization strategies for this problem, since no solver is really going to help you.</div><div><br></div><div>  Thanks,</div><div><br></div><div>    Matt</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>Thanks,</div><div>Eduardo</div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jun 2, 2015 at 4:11 AM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span><br>
> On Jun 1, 2015, at 4:06 PM, Eduardo <<a href="mailto:erocha.ssa@gmail.com" target="_blank">erocha.ssa@gmail.com</a>> wrote:<br>
><br>
> I am solving a linear system for which the preconditioned residual decreases, but the true residual increases (or have an erratic behavior). According to Petsc FAQ, this is due to a preconditioner that is singular or close to singular. What can I do in this case? I used GMRES with ILU preconditioner.<br>
</span>                                              ^^^^^^^^^^^^^^^<br>
<span>><br>
> Incidentally, I tried to solve a smaller system with a direct solver (superlu_dist) and it ran, so the system apparently is not singular.<br>
<br>
</span>    ILU can produce very badly conditioned (one could say singular PRECONDITIONER) from not singular sparse matrices. So it doesn't have anything to do with the system itself being singular.<br>
<br>
   What type of problem are you solving? Different problems need different preconditioners.<br>
<br>
  Barry<br>
<br>
<br>
<br>
><br>
> Thanks in advance,<br>
> Eduardo<br>
<br>
</blockquote></div><br></div>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">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></div>