<div dir="ltr"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Thanks for the reply.</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div>1) This is block-Jacobi, why not use PCBJACOBI? Is it because you want to select rows?</div></div></div></div></blockquote><div><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">I'm only using it to understand the performance behavior of PCFieldSplit since I'm also</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">having the same issue in a large and more complex problem.</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">​</div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div>2) We cannot tell anything without knowing how many iterates were used</div></div></div></div></blockquote><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div>  -ksp_monitor_true_residual -ksp_converged_reason -pc_fieldsplit_[0,1]_ksp_monit<wbr>or_true_residual</div><div><br></div><div>3) We cannot say anything about performance without seeing the log for both runs</div><div>  -log_view</div></div></div></div></blockquote><div><br></div><div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">I'm sending to you the log files with the recommended command line arguments for the three cases.</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">1-scalar case</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">2-PCFieldSplit (as we were initially running)</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">3-PCFieldSplit with Preonly/Jacobi in each block, as suggested by Patrick.</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default"><span style="font-size:12.8px">As Patrick pointed out, with Preonly/Jacobi the behavior is closer to what I expected.</span></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Please note that the log was taken for 100 calls to KSPSolve, I just simplified it.</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">What would be the proper way of creating this block preconditioner <br><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">As you can see, the timing with PCFieldSplit is bigger for case 3.<br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">For case 2 it is nearly 2x, as I expected (I don't know if this idea makes sense).</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">So for the case 2, the reason for the large timing is due to the inner/outer solver?</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div></div><div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">​Does the "machinery" behind the PCFieldSplit for a block preconditioner results</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">in some performance overhead? (neglecting the efficiency of the PC itself)</div><br></div><div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Best regards,</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Bernardo​</div><br></div></div></div></div>