[petsc-dev] FieldSplit fixes in 3.3

Barry Smith bsmith at mcs.anl.gov
Wed Aug 22 20:38:55 CDT 2012


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





More information about the petsc-dev mailing list