<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Fri, Aug 29, 2014 at 8:55 AM, 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 trying PCFieldSplitSetSchurPre with<br>
PC_FIELDSPLIT_SCHUR_PRE_SELFP in petsc-3.5.1 using fortran.<br>
<br>
The first problem is that PC_FIELDSPLIT_SCHUR_PRE_SELFP seems to<br>
be missing in fortran, I get the compile error:<br>
<br>
This name does not have a type, and must have an explicit type.   [PC_FIELDSPLIT_SCHUR_PRE_SELFP]<br>
<br>
while compilation works fine with _A11, _SELF and _USER.<br></blockquote><div><br></div><div>Mark Adams has just fixed this.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

The second problem is that the call doesn't seem to have any<br>
effect. For example, I have<br>
<br>
CALL PCFieldSplitSetSchurPre(pc,PC_FIELDSPLIT_SCHUR_PRE_USER,aa,ierr)<br>
CALL PCFieldSplitSetSchurFactType(pc,PC_FIELDSPLIT_SCHUR_FACT_LOWER,ierr)<br>
<br>
This compiles and runs, but ksp_view tells me<br>
<br>
PC Object:(sys_) 3 MPI processes<br>
  type: fieldsplit<br>
    FieldSplit with Schur preconditioner, factorization LOWER<br>
    Preconditioner for the Schur complement formed from A11<br>
<br>
So changing the factorization from the default FULL to LOWER did<br>
work, but changing the preconditioner from A11 to USER didn't.<br>
<br>
I've also tried to run directly from the command line using<br>
<br>
-sys_pc_fieldsplit_schur_precondition user -sys_ksp_view<br>
<br>
This works in the sense that I don't get the "WARNING! There are<br>
options you set that were not used!" message, but still ksp_view<br>
reports A11 instead of user provided matrix.<br></blockquote><div><br></div><div>Can you send a small test code, since I use this everyday here and it works.</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<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>
<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
</div></div>