<div dir="ltr">I think there is something wrong with your matrix. Use any solver and verify that you like the solution first.</div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Oct 28, 2016 at 9:13 AM, Jeremy Theler <span dir="ltr"><<a href="mailto:jeremy@seamplex.com" target="_blank">jeremy@seamplex.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi! I want to use PCGAMG as a preconditioner for a 3D linear elasticity<br>
problem (displacement-based FEM formulation) over an unstructured grid.<br>
I am not using DMPlex, I just build the stiffness matrix myself and pass<br>
it to PETSc.<br>
<br>
I set MatSetBlockSize() to 3 and pass the node coordinates through<br>
PCSetCoordinates(). But using gamg and gmres I get:<br>
<br>
PETSc error 77-0 'Eigen estimator failed: DIVERGED_NANORINF at iteration<br>
0' in /home/gtheler/libs/petsc-3.7.<wbr>4/src/ksp/ksp/impls/cheby/<wbr>cheby.c<br>
KSPSolve_Chebyshev:440<br>
<br>
Any suggestion? Another PC/KSP combination to try?<br>
<br>
Thanks<br>
--<br>
jeremy<br>
<br>
<br>
</blockquote></div><br></div>