[petsc-users] MatSetValues() for MATMPICUSP

[recrusader]

Dear Jed,

I am sorry that I cannot generate the errors using PETSc by itself.
However, it easy to generate the errors with libmesh.
There is an example from libmesh, that is
http://libmesh.sourceforge.net/ex4.php
I can run it using the command "ibrun -n 2 -o 0 ./introduction_ex4-dbg -d 3
-pc_type none -ksp_type gmres -vec_type mpicusp -mat_type mpiaijcusp
-ksp_view -ksp_monitor -log_summary -malloc_debug -cuda_show_devices" (You
should replace ibrun -n 2 -o 0 using mpiexec -n 2)
I get the same errors. However, when using "ibrun -n 2 -o 0
./introduction_ex4-dbg -d 3 -pc_type none -ksp_type gmres -ksp_view
-ksp_monitor -log_summary". I can run it successfully.

Thank you very much.

Best,
Yujie

On Sat, Feb 11, 2012 at 4:43 AM, Jed Brown <jedbrown at mcs.anl.gov> wrote:

> On Sat, Feb 11, 2012 at 00:00, recrusader <recrusader at gmail.com> wrote:
>
>> in
>> http://www.mcs.anl.gov/petsc/petsc-dev/src/ksp/ksp/examples/tutorials/ex43.c.html
>>
>> Why are two matrices set to MATAIJ? if I set MATAIJCUSP, is them changed?
>>
>> 1468: DMCreateMatrix<http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/DM/DMCreateMatrix.html#DMCreateMatrix>
>> (da_Stokes,MATAIJ<http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Mat/MATAIJ.html#MATAIJ>
>> ,&A);
>>
>> 1469:   DMCreateMatrix <http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/DM/DMCreateMatrix.html#DMCreateMatrix>(da_Stokes,MATAIJ <http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Mat/MATAIJ.html#MATAIJ>,&B);
>>
>>
> You have to follow the code down to where those matrices are used.
>
> 1478:   AssembleA_Stokes(A,da_Stokes,da_prop,properties);
> 1479:   AssembleA_PCStokes(B,da_Stokes,da_prop,properties);
>
> This uses a "block preconditioner", the B matrix is intentionally
> different from the A matrix in order to handle the saddle point.
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120211/c316ab69/attachment.htm>