<div class="gmail_quote">On Thu, Mar 8, 2012 at 15:40, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div id=":2uh">But the LSC_L reflects properties of the true operator while LSC_Lp reflects properties of the matrix from which you build the preconditioner.</div></blockquote><div><br></div><div>Only if you used -pc_fieldsplit_real_diagonal. My concern is that pc->mat is a MatSchurComplement that is not likely to be made by the user (the user never has access to it in PCFieldSplit so there is no way it can be anything else). If we actually used the MatGetSchurComplement() interface, the user would have an opportunity to see it, but it's not practical to do that currently because it would duplicate the submatrices.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div id=":2uh"> Hence LSC_L should be attached to pc->mat<br>
<br>
  So why not check in pc->mat first and if not there check in pc->pmat</div></blockquote></div><br>