Since MATMPIAIJ works (I didn't change anything. Just set the vec and mat types to mpicusp and mpiaijcusp for GPU), I think the problem is likely from MatSetValues_MPIAIJ().
Which PETSc examples can test this function? I can try to use it to test. Thanks a lot, Jed. Best, Yujie On Fri, Feb 10, 2012 at 2:44 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote: > On Fri, Feb 10, 2012 at 14:35, recrusader <recrusader at gmail.com> wrote: > >> MATMPIAIJ works ;(. > > > Can you send a reproducible test case? (Instructions for how to run a > Libmesh example that would cause this problem is okay, producing it with a > PETSc example would be better.) > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120210/dfbc86e9/attachment.htm>
