<div dir="ltr">Dear all, I have a four field problem with dof structure stored in a DM shell. The matrix is a MatNest. This works fine except for when I try to group fields. For example the following commands<div><br></div><div><div>-ksp_type gmres </div><div>   -pc_type fieldsplit </div><div>   -pc_fieldsplit_type multiplicative </div><div>      -pc_fieldsplit_0_fields 0,1 </div><div>          -fieldsplit_0_ksp_type preonly </div><div>          -fieldsplit_0_pc_type lu </div><div>       -pc_fieldsplit_1_fields 2,3 </div><div>          -fieldsplit_1_ksp_type preonly </div><div>          -fieldsplit_1_pc_type ilu </div></div><div><br></div><div>results in the following error</div><div><br></div><div><div>[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</div><div>[0]PETSC ERROR: Arguments are incompatible</div><div>[0]PETSC ERROR: Could not find index set</div><div>[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.</div><div>[0]PETSC ERROR: Petsc Release Version 3.5.2, Sep, 08, 2014 </div><div>[0]PETSC ERROR: /Users/Christina/Desktop/COLIN/Work/feapver83/ver83/parfeap/feap on a gnu-dbg named new-host.home by christina Mon Jan 12 00:36:52 2015</div><div>[0]PETSC ERROR: Configure options --with-debugging=1 --with-shared-libraries=0 --with-mpi-dir=/usr/local --with-parmetis-dir=/usr/local --with-metis-dir=/usr/local --with-hypre-dir=/usr/local --download-suitesparse --with-external-packages-dir=/usr/local</div><div>[0]PETSC ERROR: #1 MatNestFindIS() line 310 in /Users/Christina/Desktop/COLIN/Work/petsc-3.5.2/src/mat/impls/nest/matnest.c</div><div>[0]PETSC ERROR: #2 MatNestFindSubMat() line 371 in /Users/Christina/Desktop/COLIN/Work/petsc-3.5.2/src/mat/impls/nest/matnest.c</div><div>[0]PETSC ERROR: #3 MatGetSubMatrix_Nest() line 387 in /Users/Christina/Desktop/COLIN/Work/petsc-3.5.2/src/mat/impls/nest/matnest.c</div><div>[0]PETSC ERROR: #4 MatGetSubMatrix() line 7345 in /Users/Christina/Desktop/COLIN/Work/petsc-3.5.2/src/mat/interface/matrix.c</div><div>[0]PETSC ERROR: #5 PCSetUp_FieldSplit() line 559 in /Users/Christina/Desktop/COLIN/Work/petsc-3.5.2/src/ksp/pc/impls/fieldsplit/fieldsplit.c</div><div>[0]PETSC ERROR: #6 PCSetUp() line 902 in /Users/Christina/Desktop/COLIN/Work/petsc-3.5.2/src/ksp/pc/interface/precon.c</div><div>[0]PETSC ERROR: #7 KSPSetUp() line 305 in /Users/Christina/Desktop/COLIN/Work/petsc-3.5.2/src/ksp/ksp/interface/itfunc.c</div><div>[0]PETSC ERROR: #8 KSPSolve() line 417 in /Users/Christina/Desktop/COLIN/Work/petsc-3.5.2/src/ksp/ksp/interface/itfunc.c</div></div><div><br></div><div>My guess is that petsc is expecting a MatNest corresponding to the outer field groupings instead of individual submatrices. Is it possible to resolve this, preferably from the command line? </div><div><br></div><div>I should also mention that grouping and recursive field split works fine with the same DM shell and a MATAIJ, however a significant amount of time is spent copying the submatrices and I'd prefer to avoid that.<br></div><div><br></div><div>Thanks and all the best,</div><div>Colin</div></div>