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

[Stefano Zampini]

Can you send a minimal working example so that I can fix the code?

Thanks
Stefano

Il 19 Giu 2017 15:20, "Franck Houssen" <franck.houssen at inria.fr> ha scritto:

> 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/9b44e113/attachment.html>