On Fri, Feb 10, 2012 at 14:24, recrusader <[email protected]> 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.* http://www.mcs.anl.gov/petsc/documentation/changes/dev.html