On Aug 22, 2012, at 8:46 AM, Jed Brown <[email protected]> wrote:

On Wed, Aug 22, 2012 at 8:35 AM, Barry Smith <[email protected]> wrote: Should MatCreateSchurComplement() be done properly the PETSc way: that is normal usage is

MatSetType(mat,MATSCHURCOMPLEMENT); /* now if you want you can muck with the inner KSP */ MatSchurComplementGetKSP(mat,&ksp); ……. MatSchurComplementSetSubMats(mat,……) note: could merge with MatSchurComplementUpdate() …..

and MatCreateSchurComplement() is just a convenience routine that calls the pieces in order and is NOT used directly by field split

Yes

Working on this now. Will push after tested. Barry