<div dir="ltr">Hello guys,<div><br></div><div>I am getting this error while using fieldsplit:</div><div><br></div><div>[3]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------<br><b>[3]PETSC ERROR: Nonconforming object sizes<br>[3]PETSC ERROR: Local column sizes 6132 do not add up to total number of columns 9200</b><br>[3]PETSC ERROR: See <a href="https://petsc.org/release/faq/">https://petsc.org/release/faq/</a> for trouble shooting.<br>[3]PETSC ERROR: Petsc Development GIT revision: v3.18.1-191-g32ed6ae2ff2  GIT Date: 2022-11-08 12:22:17 -0500<br>[3]PETSC ERROR: flubio_coupled on a gnu named alienware by edo Wed Nov  9 08:16:29 2022<br>[3]PETSC ERROR: Configure options PETSC_ARCH=gnu FOPTFLAGS=-O3 COPTFLAGS=-O3 CXXOPTFLAGS=-O3 -with-debugging=no -download-fblaslapack=1 -download-superlu_dist -download-mumps -download-hypre -download-metis -download-parmetis -download-scalapack -download-ml -download-slepc -download-hpddm -download-cmake -with-mpi-dir=/home/edo/software/openmpi-4.1.1/build/<br>[3]PETSC ERROR: #1 MatCreateSubMatrix_MPIBAIJ_Private() at /home/edo/software/petsc/src/mat/impls/baij/mpi/mpibaij.c:1987<br>[3]PETSC ERROR: #2 MatCreateSubMatrix_MPIBAIJ() at /home/edo/software/petsc/src/mat/impls/baij/mpi/mpibaij.c:1911<br>[3]PETSC ERROR: #3 MatCreateSubMatrix() at /home/edo/software/petsc/src/mat/interface/matrix.c:8340<br>[3]PETSC ERROR: #4 PCSetUp_FieldSplit() at /home/edo/software/petsc/src/ksp/pc/impls/fieldsplit/fieldsplit.c:657<br>[3]PETSC ERROR: #5 PCSetUp() at /home/edo/software/petsc/src/ksp/pc/interface/precon.c:994<br>[3]PETSC ERROR: #6 KSPSetUp() at /home/edo/software/petsc/src/ksp/ksp/interface/itfunc.c:406<br>[3]PETSC ERROR: #7 KSPSolve_Private() at /home/edo/software/petsc/src/ksp/ksp/interface/itfunc.c:825<br>[3]PETSC ERROR: #8 KSPSolve() at /home/edo/software/petsc/src/ksp/ksp/interface/itfunc.c:1071<br><br></div><div>Do you have any ideas? Probably something missing in my brief implementation here:</div><div><br></div><div><i>            call PCSetType(mypc, PCFIELDSPLIT, ierr)  <br><br>            call PCFieldSplitSetBlockSize(mypc, 4-bdim, ierr)<br>           </i></div><div><i>            !2D, 3x3 block<br>            if(bdim==1) then <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><br>             ! 3D 4x4 block<br>            else <br><br>                ufields(1) = 0<br>                ufields(2) = 1<br>                ufields(3) = 2<br>                pfields(1) = 3<br>                call PCFieldSplitSetFields(mypc, "u", 3, ufields, ufields, ierr)<br>                call PCFieldSplitSetFields(mypc, "p", 1, pfields, pfields, ierr)<br><br>            endif <br>            <br>            ! Field split type ADDITIVE, MULTIPLICATIVE (default), SYMMETRIC_MULTIPLICATIVE, SPECIAL, SCHUR<br>            call PCFieldSplitSetType(mypc, PC_COMPOSITE_SCHUR, ierr)</i><br></div><div><i><br></i></div><div>Thanks for the help!</div></div>