[petsc-users] Building MatIS with dense local matrix ?

[Franck Houssen]

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 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20170619/ac5e05d7/attachment.html>