Re: [petsc-dev] FieldSplit fixes in 3.3

22 Aug 2012


      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