[petsc-users] On PCFIELDSPLIT and its implementation

[Edoardo alinovi]

Thanks,

the stuff I am doing is within my code, so I am not sure you can reproduce
it.

I am just doing this:

            call PCSetType(mypc, PCFIELDSPLIT, ierr)

            call PCFieldSplitSetBlockSize(mypc, 4-bdim, ierr)

                ufields(1) = 0
                ufields(2) = 1
                pfields(1) = 2
                call PCFieldSplitSetFields(mypc, "u", 2, ufields, ufields,
ierr)
                call PCFieldSplitSetFields(mypc, "p", 1, pfields, pfields,
ierr)

On an MPIBAIJ matrix with bs = 3.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20221109/797287a1/attachment.html>