<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Mar 22, 2017 at 2:58 PM, Klaij, Christiaan <span dir="ltr"><<a href="mailto:C.Klaij@marin.nl" target="_blank">C.Klaij@marin.nl</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I'm solving the Navier-Stokes equations using PCFieldSplit type<br>
Schur and Selfp. This particular case has only Neumann conditions<br>
for the pressure field.<br>
<br>
With left preconditioning and no nullspace, I see that the KSP<br>
solver for S does not converge (attachment "left_nonullsp") in<br>
either norm.<br>
<br>
When I attach the constant null space to A11, it gets passed on<br>
to S and the KSP solver for S does converge in the preconditioned<br>
norm only (attachment "left").<br>
<br>
However, right preconditioning uses the unpreconditioned norm and<br>
therefore doesn't converge (attachment "right"), regardless of<br>
whether the nullspace is attached or not. Should I conclude that<br>
right preconditioning cannot be used in combination with a null<br>
space?<br></blockquote><div><br></div><div>No, neither of your solves is working. The left preconditioned version is just hiding that fact.</div><div><br></div><div>You should start by checking that the exact solve works. Namely full Schur factorization</div><div>with exact solves for A and S. Since you do not have a matrix for S (unless you tell it do</div><div>use "full"), just use a very low (1e-10) tolerance. My guess is that something is off with</div><div>your null space specification.</div><div><br></div><div>  Thanks,</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">
Chris<br>
<br>
<br>
dr. ir. Christiaan Klaij  | Senior Researcher | Research & Development<br>
MARIN | T <a href="tel:%2B31%20317%2049%2033%2044" value="+31317493344">+31 317 49 33 44</a> | mailto:<a href="mailto:C.Klaij@marin.nl">C.Klaij@marin.nl</a> | <a href="http://www.marin.nl" rel="noreferrer" target="_blank">http://www.marin.nl</a><br>
<br>
MARIN news: <a href="http://www.marin.nl/web/News/News-items/H2020-marinergi-project-launched.htm" rel="noreferrer" target="_blank">http://www.marin.nl/web/News/<wbr>News-items/H2020-marinergi-<wbr>project-launched.htm</a><br>
<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="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>