<div dir="ltr">Thanks, <div><br></div><div>the stuff I am doing is within my code, so I am not sure you can reproduce it. </div><div><br></div><div>I am just doing this: </div><div><br></div><div>            call PCSetType(mypc, PCFIELDSPLIT, ierr)  <br><br>            call PCFieldSplitSetBlockSize(mypc, 4-bdim, ierr)<br>           <br>                ufields(1) = 0<br>                ufields(2) = 1<br>                pfields(1) = 2<br>                call PCFieldSplitSetFields(mypc, "u", 2, ufields, ufields, ierr)<br>                call PCFieldSplitSetFields(mypc, "p", 1, pfields, pfields, ierr)<br></div><div><br></div><div>On an MPIBAIJ matrix with bs = 3. </div><div><br></div><div><br></div><div><br></div><div><br></div></div>