<div class="gmail_quote">On Fri, Feb 10, 2012 at 14:24, recrusader <span dir="ltr">&lt;<a href="mailto:recrusader@gmail.com">recrusader@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div id=":399">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>
&quot;<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>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">http://www.mcs.anl.gov/petsc/documentation/changes/dev.html</a></div>