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_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: call MatCreate(this%comm,AA,ierr) call MatSetSizes(AA,npetscloc,npetscloc,nreal,nreal,ierr) call MatSetType(AA,MATAIJ,ierr) call MatSetup(AA,ierr) call MatGetOwnershipRange(AA,low,high,ierr) allocate(d_nnz(npetscloc),o_nnz(npetscloc)) call getNNZ(grid,npetscloc,low,high,d_nnz,o_nnz,this%xgc_petsc,nreal,ierr) call MatSeqAIJSetPreallocation(AA,PETSC_NULL_INTEGER,d_nnz,ierr) call MatMPIAIJSetPreallocation(AA,PETSC_NULL_INTEGER,d_nnz,PETSC_NULL_INTEGER,o_nnz,ierr) deallocate(d_nnz,o_nnz) call MatSetOption(AA,MAT_IGNORE_OFF_PROC_ENTRIES,PETSC_TRUE,ierr) call MatSetOption(AA,MAT_KEEP_NONZERO_PATTERN,PETSC_TRUE,ierr) call MatSetup(AA,ierr)
[62]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- [62]PETSC ERROR: No support for this operation for this object type [62]PETSC ERROR: unsorted iscol_local is not implemented yet [62]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting. [62]PETSC ERROR: Petsc Release Version 3.8.0, unknown[62]PETSC ERROR: #1 MatCreateSubMatrix_MPIAIJ_SameRowDist() line 3418 in /global/u1/r/rchurchi/petsc/3.8.0/src/mat/impls/aij/mpi/mpiaij.c [62]PETSC ERROR: #2 MatCreateSubMatrix_MPIAIJ() line 3247 in /global/u1/r/rchurchi/petsc/3.8.0/src/mat/impls/aij/mpi/mpiaij.c [62]PETSC ERROR: #3 MatCreateSubMatrix() line 7872 in /global/u1/r/rchurchi/petsc/3.8.0/src/mat/interface/matrix.c [62]PETSC ERROR: #4 PCGAMGCreateLevel_GAMG() line 383 in /global/u1/r/rchurchi/petsc/3.8.0/src/ksp/pc/impls/gamg/gamg.c [62]PETSC ERROR: #5 PCSetUp_GAMG() line 561 in /global/u1/r/rchurchi/petsc/3.8.0/src/ksp/pc/impls/gamg/gamg.c [62]PETSC ERROR: #6 PCSetUp() line 924 in /global/u1/r/rchurchi/petsc/3.8.0/src/ksp/pc/interface/precon.c [62]PETSC ERROR: #7 KSPSetUp() line 378 in /global/u1/r/rchurchi/petsc/3.8.0/src/ksp/ksp/interface/itfunc.c -- R. Michael Churchill