<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Sep 4, 2014 at 7:06 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 playing with the selfp option in fieldsplit using<br>
snes/examples/tutorials/ex70.c. For example:<br>
<br>
mpiexec -n 2 ./ex70 -nx 4 -ny 6 \<br>
-ksp_type fgmres \<br>
-pc_type fieldsplit \<br>
-pc_fieldsplit_type schur \<br>
-pc_fieldsplit_schur_fact_type lower \<br>
-pc_fieldsplit_schur_precondition selfp \<br>
-fieldsplit_1_inner_ksp_type preonly \<br>
-fieldsplit_1_inner_pc_type jacobi \<br>
-fieldsplit_0_ksp_monitor -fieldsplit_0_ksp_max_it 1 \<br>
-fieldsplit_1_ksp_monitor -fieldsplit_1_ksp_max_it 1 \<br>
-ksp_monitor -ksp_max_it 1<br>
<br>
gives the following output<br>
<br>
  0 KSP Residual norm 1.229687498638e+00<br>
    Residual norms for fieldsplit_1_ solve.<br>
    0 KSP Residual norm 2.330138480101e+01<br>
    1 KSP Residual norm 1.609000846751e+01<br>
  1 KSP Residual norm 1.180287268335e+00<br>
<br>
To my suprise I don't see anything for the fieldsplit_0_ solve,<br>
why?<br></blockquote><div><br></div><div>Always run with -ksp_view for any solver question.</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">

Furthermore, if I understand correctly the above should be<br>
exactly equivalent with<br>
<br>
mpiexec -n 2 ./ex70 -nx 4 -ny 6 \<br>
-ksp_type fgmres \<br>
-pc_type fieldsplit \<br>
-pc_fieldsplit_type schur \<br>
-pc_fieldsplit_schur_fact_type lower \<br>
-user_ksp \<br>
-fieldsplit_0_ksp_monitor -fieldsplit_0_ksp_max_it 1 \<br>
-fieldsplit_1_ksp_monitor -fieldsplit_1_ksp_max_it 1 \<br>
-ksp_monitor -ksp_max_it 1<br>
<br>
  0 KSP Residual norm 1.229687498638e+00<br>
    Residual norms for fieldsplit_0_ solve.<br>
    0 KSP Residual norm 5.486639587672e-01<br>
    1 KSP Residual norm 6.348354253703e-02<br>
    Residual norms for fieldsplit_1_ solve.<br>
    0 KSP Residual norm 2.321938107977e+01<br>
    1 KSP Residual norm 1.605484031258e+01<br>
  1 KSP Residual norm 1.183225251166e+00<br>
<br>
because -user_ksp replaces the Schur complement by the simple<br>
approximation A11 - A10 inv(diag(A00)) A01. Beside the missing<br>
fielsplit_0_ part, the numbers are pretty close but not exactly<br>
the same. Any explanation?<br>
<br>
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>