[petsc-users] Fieldsplit with redistribute

[Carl-Johan Thore]

Hello,
I'm solving a blocksystem
[A C;
C' D],
where D is not zero, using the PCFIELDSPLIT preconditioner and set the split using PetscFieldSplitSetIS. This works very well until I try PCREDISTRIBUTE (which is attractive as I have many locked DOFs). I suspect something goes wrong when constructing the IS for the split (I've tried various things using the IS-routines). Can PETSc do this automatically? Or else, any hints?
Kind regards,
Carl-Johan