On Wed, Jun 20, 2012 at 9:00 AM, Klaij, Christiaan <span dir="ltr"><<a href="mailto:C.Klaij@marin.nl" target="_blank">C.Klaij@marin.nl</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
In order to implement SIMPLE-type preconditioners for the<br>
incompressible Navier-Stokes equations (Elman e.a. JCP 227, 2008)<br>
using the PCFieldSplit framework, it looks like I need to replace<br>
inv(A00) by some cheap approximation<br>
<br>
1) in the Schur complement<br></blockquote><div><br></div><div>When you have a Schur FS, the '0' solver is this approximation.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

2) in the L and/or U of the LDU factorization<br></blockquote><div><br></div><div>You can use whatever PC you want for the solver mentioned above.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

3) while keeping A00 in the D<br></blockquote><div><br></div><div>I think what you want here is -pc_fieldsplit_real_diagonal.</div><div><br></div><div>I would be really nice to make a PETSc example that did SIMPLE. I was</div>
<div>going to do this, but if you do it first, I will put it in.</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">
What would be the best way?<br>
<br>
<br>
dr. ir. Christiaan Klaij<br>
CFD Researcher<br>
Research & Development<br>
E mailto:<a href="mailto:C.Klaij@marin.nl">C.Klaij@marin.nl</a><br>
T <a href="tel:%2B31%20317%2049%2033%2044" value="+31317493344">+31 317 49 33 44</a><br>
<br>
MARIN<br>
2, Haagsteeg, P.O. Box 28, 6700 AA Wageningen, The Netherlands<br>
T <a href="tel:%2B31%20317%2049%2039%2011" value="+31317493911">+31 317 49 39 11</a>, F <a href="tel:%2B31%20317%2049%2032%2045" value="+31317493245">+31 317 49 32 45</a>, I <a href="http://www.marin.nl" target="_blank">www.marin.nl</a><br>

<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>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<br>