[petsc-users] Fieldsplit with redistribute

Carl-Johan Thore carl-johan.thore at liu.se
Sun Apr 16 13:50:58 CDT 2023

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,

More information about the petsc-users mailing list