<div class="gmail_quote">On Wed, Aug 22, 2012 at 8:35 AM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">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=":1l4"> Should MatCreateSchurComplement() be done properly the  PETSc way:   that is normal usage is<br>
<br>
     MatSetType(mat,MATSCHURCOMPLEMENT);<br>
     /* now if you want you can muck with the inner KSP */<br>
           MatSchurComplementGetKSP(mat,&ksp);<br>
           …….<br>
     MatSchurComplementSetSubMats(mat,……)       note: could merge with MatSchurComplementUpdate()<br>
      …..<br>
<br>
       and MatCreateSchurComplement() is just a convenience routine that calls the pieces in order and is NOT used directly by field split</div></blockquote></div><br><div>Yes</div>