[petsc-users] access to matnest block (0,1) ?

Jed Brown jedbrown at mcs.anl.gov
Fri Jan 27 06:14:51 CST 2012


On Fri, Jan 27, 2012 at 02:43, Klaij, Christiaan <C.Klaij at marin.nl> wrote:

> [0]PETSC ERROR: Unhandled case, must have at least two fields!


You can use PCFieldSplitSetIS(). The implementation could check whether a
MatNest is being used and automatically set the splits if it is, but you
would have to call PCFieldSplitSetIS() later when you wanted to assemble
into an AIJ format, so I'm hesitant to pick it up automatically. Just call
the function for now.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120127/e7ac67da/attachment.htm>


More information about the petsc-users mailing list