<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">Randy:</div><div class="gmail_quote">It could be a bug or a missing feature in our new MatCreateSubMatrix_MPIAIJ_<wbr>SameRowDist(). </div><div class="gmail_quote">It would be helpful if you can provide us a simple example that produces this example. </div><div class="gmail_quote">Hong</div><div class="gmail_quote"><br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">I'm running a Fortran code that was just changed over to using petsc 3.8 (previously petsc 3.7.6). An error was thrown during a KSPSetUp() call. The error is "unsorted iscol_local is not implemented yet" (see full error below). I tried to trace down the difference in the source files, but where the error occurs (MatCreateSubMatrix_MPIAIJ_<wbr>SameRowDist()) doesn't seem to have existed in v3.7.6, so I'm unsure how to compare. It seems the error is that the order of the columns locally are unsorted, though I don't think I specify a column order in the creation of the matrix:<div>     call MatCreate(this%comm,AA,ierr)</div><div>     call MatSetSizes(AA,npetscloc,<wbr>npetscloc,nreal,nreal,ierr)</div><div>     call MatSetType(AA,MATAIJ,ierr)</div><div>     call MatSetup(AA,ierr)</div><div>     call MatGetOwnershipRange(AA,low,<wbr>high,ierr)</div><div>     allocate(d_nnz(npetscloc),o_<wbr>nnz(npetscloc))</div><div>     call getNNZ(grid,npetscloc,low,<wbr>high,d_nnz,o_nnz,this%xgc_<wbr>petsc,nreal,ierr)</div><div>     call MatSeqAIJSetPreallocation(AA,<wbr>PETSC_NULL_INTEGER,d_nnz,ierr)</div><div>     call MatMPIAIJSetPreallocation(AA,<wbr>PETSC_NULL_INTEGER,d_nnz,<wbr>PETSC_NULL_INTEGER,o_nnz,ierr)</div><div>     deallocate(d_nnz,o_nnz)</div><div>     call MatSetOption(AA,MAT_IGNORE_<wbr>OFF_PROC_ENTRIES,PETSC_TRUE,<wbr>ierr)</div><div>     call MatSetOption(AA,MAT_KEEP_<wbr>NONZERO_PATTERN,PETSC_TRUE,<wbr>ierr)</div><div>     call MatSetup(AA,ierr) <div>

<div><br></div><div><br></div><div>[62]PETSC ERROR: --------------------- Error Message ------------------------------<wbr>------------------------------<wbr>--</div><div>[62]PETSC ERROR: No support for this operation for this object type</div><div>[62]PETSC ERROR: unsorted iscol_local is not implemented yet</div><div>[62]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" target="_blank">http://www.mcs.anl.gov/petsc/<wbr>documentation/faq.html</a> for trouble shooting.</div><div>[62]PETSC ERROR: Petsc Release Version 3.8.0, unknown[62]PETSC ERROR: #1 MatCreateSubMatrix_MPIAIJ_<wbr>SameRowDist() line 3418 in /global/u1/r/rchurchi/petsc/3.<wbr>8.0/src/mat/impls/aij/mpi/<wbr>mpiaij.c
</div><div>[62]PETSC ERROR: #2 MatCreateSubMatrix_MPIAIJ() line 3247 in /global/u1/r/rchurchi/petsc/3.<wbr>8.0/src/mat/impls/aij/mpi/<wbr>mpiaij.c
</div><div>[62]PETSC ERROR: #3 MatCreateSubMatrix() line 7872 in /global/u1/r/rchurchi/petsc/3.<wbr>8.0/src/mat/interface/matrix.c
</div><div>[62]PETSC ERROR: #4 PCGAMGCreateLevel_GAMG() line 383 in /global/u1/r/rchurchi/petsc/3.<wbr>8.0/src/ksp/pc/impls/gamg/<wbr>gamg.c
</div><div>[62]PETSC ERROR: #5 PCSetUp_GAMG() line 561 in /global/u1/r/rchurchi/petsc/3.<wbr>8.0/src/ksp/pc/impls/gamg/<wbr>gamg.c
</div><div>[62]PETSC ERROR: #6 PCSetUp() line 924 in /global/u1/r/rchurchi/petsc/3.<wbr>8.0/src/ksp/pc/interface/<wbr>precon.c
</div><div>[62]PETSC ERROR: #7 KSPSetUp() line 378 in /global/u1/r/rchurchi/petsc/3.<wbr>8.0/src/ksp/ksp/interface/<wbr>itfunc.c </div><span class="gmail-HOEnZb"><font color="#888888"><div><br></div>-- <br><div class="gmail-m_-4997150650000030250gmail_signature"><div dir="ltr">R. Michael Churchill</div></div>
</font></span></div></div></div>
</blockquote></div><br></div></div>