Hi, I try to call MatISGetMPIXAIJ on a MatIS (A) that has been feed locally by sequential (Aloc) dense matrix. Seems this ends up with this error: [0]PETSC ERROR: New nonzero at (0,1) caused a malloc. Is this a known error / limitation ? (not supposed to work with dense matrix ?)
This (pseudo code) works fine: MatCreateIS(..., A) MatCreateSeqAIJ(..., Aloc) MatISSetLocalMat(pcA, pcALoc) MatISGetMPIXAIJ(A, ...) // OK ! When I try to replace MatCreateSeqAIJ(..., Aloc) with MatCreateSeqDense(..., Aloc), it does no more work. Franck PS: running debian/testing with gcc-6.3 + petsc-3.7.6