On Fri, Feb 10, 2012 at 14:24, recrusader
<recrusader@gmail.com> wrote:
I use Libmesh with PETSc for my FEM simulation.
The program works when using MATSEQCUSP. However, When I test
MATMPICUSP with 2 GPU cards, I met the following errors:
"
[1]PETSC ERROR: --------------------- Error Message
------------------------------------
[1]PETSC ERROR: Argument out of range!
[1]PETSC ERROR: New nonzero at (1,5) caused a malloc!
Can you try with plain MATMPIAIJ? I suspect you are not preallocating correctly.
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.