<div class="gmail_quote">On Mon, Jan 9, 2012 at 22:47, Mohamad M. Nasr-Azadani <span dir="ltr">&lt;<a href="mailto:mmnasr@gmail.com">mmnasr@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>It might be a naive question, but I am wondering how KSPSetNullSpace() works when it is passed to a linear system ksp context?</div><div>Say, for instance, we have the simple case of Poisson equation solved in a square domain and Neumann boundary condition applied to all boundaries. </div>


<div>Does it take the integral of the solution and set it to zero as an extra constraint? </div><div></div></blockquote></div><br><div>It just projects out whatever you provide as a null space, so the Krylov method effectively runs in the remaining subspace.</div>