On Aug 22, 2012, at 8:46 AM, Jed Brown <jedbrown at mcs.anl.gov> wrote:
> On Wed, Aug 22, 2012 at 8:35 AM, Barry Smith <bsmith at mcs.anl.gov> 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