<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Fri, Nov 6, 2015 at 9:29 AM, Denis Davydov <span dir="ltr"><<a href="mailto:davydden@gmail.com" target="_blank">davydden@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
> On 6 Nov 2015, at 16:22, Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> wrote:<br>
><br>
> Is it possible that the matrix is rank deficient? Jacobi will just chug along and sometimes work, but<br>
> AMG will fail spectacularly in that case.<br>
<br>
It should not. It is just a mass (overlap) matrix coming from linear FEs with zero Dirichlet BC assembled in deal.II.<br>
Due to elimination of some algebraic constraints on DoFs there are lines with only diagonal element,<br>
but it should still be SPD.<br>
<br>
More interestingly is that it does not fail immediately (i.e. the first time it’s used in SLEPc solvers),<br>
but only on the 4th step. So 3 times SLEPc worked just fine to solve GHEP with Gamg and zero shift.<br></blockquote><div><br></div><div>Then I think it is not doing what you suppose. I am not inclined to believe that it behaves differently on the</div><div>same matrix.</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">
Regards,<br>
Denis.<br>
<br>
<br>
</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>