<div dir="auto">Can you send a minimal working example so that I can fix the code?<div dir="auto"><br></div><div dir="auto">Thanks</div><div dir="auto">Stefano</div></div><div class="gmail_extra"><br><div class="gmail_quote">Il 19 Giu 2017 15:20, "Franck Houssen" <<a href="mailto:franck.houssen@inria.fr">franck.houssen@inria.fr</a>> ha scritto:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="font-family:times new roman,new york,times,serif;font-size:12pt;color:#000000"><div>Hi,<br></div><div><br></div><div>I try to call MatISGetMPIXAIJ on a MatIS (A) that has been feed locally by sequential (Aloc) dense matrix.</div><div>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 ?)<br></div><div><br></div><div>This (pseudo code) works fine:<br></div><div>MatCreateIS(..., A)</div><div>MatCreateSeqAIJ(..., Aloc)<br></div><div>MatISSetLocalMat(pcA, pcALoc) <br></div><div>MatISGetMPIXAIJ(A, ...) // OK !<br></div><div><br></div><div>When I try to replace MatCreateSeqAIJ(..., Aloc) with MatCreateSeqDense(..., Aloc), it does no more work. <br></div><div><br></div><div>Franck<br></div><div><br></div><div>PS: running debian/testing with gcc-6.3 + petsc-3.7.6<br></div></div></div></blockquote></div></div>