[petsc-users] MatSetValues() for MATMPICUSP

[recrusader]

Dear Jed,

MATMPIAIJ works ;(.

Best,
Yujie

On Fri, Feb 10, 2012 at 2:26 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote:

> On Fri, Feb 10, 2012 at 14:24, recrusader <recrusader at 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.*
>
> http://www.mcs.anl.gov/petsc/documentation/changes/dev.html
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120210/8e6759d6/attachment.htm>