[petsc-users] Fieldsplit with redistribute

[Carl-Johan Thore via petsc-users]

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