Dear Jed,<br><br>MATMPIAIJ works ;(.<br><br>Best,<br>Yujie<br><br><div class="gmail_quote">On Fri, Feb 10, 2012 at 2:26 PM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov">jedbrown@mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im"><div class="gmail_quote">On Fri, Feb 10, 2012 at 14:24, recrusader <span dir="ltr"><<a href="mailto:recrusader@gmail.com" target="_blank">recrusader@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>I use Libmesh with PETSc for my FEM simulation.<br>
The program works when using MATSEQCUSP. However, When I test<br>
MATMPICUSP with 2 GPU cards, I met the following errors:<br>
"<br>
[1]PETSC ERROR: --------------------- Error Message<br>
------------------------------------<br>
[1]PETSC ERROR: Argument out of range!<br>
[1]PETSC ERROR: New nonzero at (1,5) caused a malloc!</div></blockquote></div><br></div><div>Can you try with plain MATMPIAIJ? I suspect you are not preallocating correctly.</div><div><br></div><div><i>Preallocation routines now automatically set MAT_NEW_NONZERO_ALLOCATION_ERR, if you intentionally preallocate less than necessary then use MatSetOption(mat,MAT_NEW_NONZERO_ALLOCATION_ERR,PETSC_FALSE) to disable the error generation.</i></div>
<div><br></div><div><a href="http://www.mcs.anl.gov/petsc/documentation/changes/dev.html" target="_blank">http://www.mcs.anl.gov/petsc/documentation/changes/dev.html</a></div>
</blockquote></div><br>