[petsc-users] MatSetValues() for MATMPICUSP

[recrusader]

Dear Matt,

I added the print codes in libmesh after creating the matrix as follows:
    "    ierr = MatCreateMPIAIJ (libMesh::COMM_WORLD,
                                m_local, n_local,
                                m_global, n_global,
                                PETSC_NULL, (int*) &n_nz[0],
                                PETSC_NULL, (int*) &n_oz[0], &_mat);
             CHKERRABORT(libMesh::COMM_WORLD,ierr);

      MatSetOption(_mat,MAT_NEW_NONZERO_ALLOCATION_ERR,PETSC_FALSE); //by
Yujie
      std::cout<<"MatSetOption"<<std::endl;"

I run the same codes in CPU and GPU modes (the same parameters except that
GPU uses '-vec_type mpicusp -mat_type mpiaijcusp'). I can find
"MatSetOption" output from both the modes. Does that mean that the codes
set the options for both the modes?
Thank you very much.

Best,
Yujie


On Fri, Feb 10, 2012 at 3:14 PM, Matthew Knepley <knepley at gmail.com> wrote:

> On Fri, Feb 10, 2012 at 3:12 PM, recrusader <recrusader at gmail.com> wrote:
>
>> Dear Jed,
>>
>> The first example works. However, the example uses MatSetValuesStencil()
>> not MatSetValues(). Are they same?
>>
>
> MatSetValuesStencil() calls MatSetValues(). I suspect that your MPIAIJ
> matrix does not have the option set to
> throw an error when inserting a new nonzero, and your MPICUSP matrix does.
>
>    Matt
>
>
>> Thanks a lot,
>> Yujie
>>
>> On Fri, Feb 10, 2012 at 2:51 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote:
>>
>>> On Fri, Feb 10, 2012 at 14:49, recrusader <recrusader at gmail.com> wrote:
>>>
>>>> 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?
>>>
>>>
>>> Try src/ksp/ksp/examples/tutorials/ex43.c and
>>> src/snes/examples/tutorials/ex48.c
>>>
>>
>>
>
>
> --
> What most experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to which their
> experiments lead.
> -- Norbert Wiener
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120210/9534ce61/attachment.htm>