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 

Reply via email to