<div dir="ltr">Hi,<br><br>I'm trying to use PCFIELDSPLIT with MATSBAIJ in PETSc 3.19.4. According to the manual<br>"[t]he fieldsplit preconditioner cannot currently be used with the MATBAIJ or MATSBAIJ data <br>formats if the blocksize is larger than 1". Since my blocksize is exactly 1 it would seem that I can<br>use PCFIELDSPLIT. But this fails with "PETSC ERROR: For symmetric format, iscol must equal isrow"<br>from MatCreateSubMatrix_MPISBAIJ. Tracing backwards one ends up in fieldsplit.c at<br><br>/* extract the A01 and A10 matrices */<br>ilink = jac->head;<br>PetscCall(ISComplement(ilink->is_col, rstart, rend, &ccis));<br>if (jac->offdiag_use_amat) {<br> PetscCall(MatCreateSubMatrix(pc->mat, ilink->is, ccis, MAT_INITIAL_MATRIX, &jac->B));<br>} else {<br>        PetscCall(MatCreateSubMatrix(pc->pmat, ilink->is, ccis, MAT_INITIAL_MATRIX, &jac->B));<br>}<br><br>This, since my A01 and A10 are not square, seems to explain why iscol is not equal to isrow.<br>From this I gather that it is in general NOT possible to use PCFIELDSPLIT with MATSBAIJ even<br>with block size 1?<br><br>Kind regards,<br>Carl-Johan</div>