<div dir="ltr"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Mar 20, 2013 at 3:13 PM, Mehrdad H Farahani <span dir="ltr"><<a href="mailto:mh.farahani@gmail.com" target="_blank">mh.farahani@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 id=":5fb">I am solving a Poisson equation with Neumann boundary conditions applied along an arbitrary curve using a ghost fluid treatment. The matrix is non-symmetric and the null space contains more than the constant vector. I have implemented a solve for the null space itself and have previously supplied it using MatSetNullSpace. Unfortunately the null space solve is extremely expensive and does not play well with GAMG or ML. We've checked that the answers we get without supplying this null space are very close to the answers we get when we do include it. The matrix for the two solves in the previous attachment are identical. The problem also seems to be converging and then suddenly diverges. Is this common behavior for a singular preconditioner?<br>
</div></blockquote></div><br>If the RHS is not consistent or the preconditioner has a different null space than the operator, you will have problems.</div></div>