<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>Thanks,</div><div>Eduardo</div></div><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 class=""><br>
> On Jun 1, 2015, at 4:06 PM, Eduardo <<a href="mailto:erocha.ssa@gmail.com">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 class="">><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>