MPIAIJ and MatSetValuesBlocked

[Ryan Yan]

Hi, all,
I am using MPIAIJ for my matrix A, but when I call the function
MatSetValuesBlocked, I got the error:

PetscPrintf(PETSC_COMM_WORLD, "breakpoint 1\n");
               ierr =
MatSetValuesBlocked(A,1,&irow,1,(col_ind+icol),temp_vector,INSERT_VALUES);
 PetscPrintf(PETSC_COMM_WORLD, "breakpoint 2\n");


breakpoint 1
[0]PETSC ERROR: --------------------- Error Message
------------------------------------
[0]PETSC ERROR: No support for this operation for this object type!
[0]PETSC ERROR: Mat type mpiaij!
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 2.3.3, Patch 15, Tue Sep 23 10:02:49
CDT 2008 HG revision: 31306062cd1a6f6a2496fccb4878f485c9b91760
[0]PETSC ERROR: See docs/changes/index.html for recent updates.
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
[0]PETSC ERROR: See docs/index.html for manual pages.
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR:
/home/vyan2000/local/PPETSc/petsc-2.3.3-p15/src/ksp/ksp/examples/tutorials/ttt2/kspex1reader_binmpiaij
on a linux-gnu named vyan2000-linux by vyan2000 Fri May  1 17:58:48 2009
[0]PETSC ERROR: Libraries linked from
/home/vyan2000/local/PPETSc/petsc-2.3.3-p15/lib/linux-gnu-c-debug
[0]PETSC ERROR: Configure run at Thu Feb  5 21:10:10 2009
[0]PETSC ERROR: Configure options --with-mpi-dir=/usr/lib/
--with-debugger=gdb --with-shared=0 --download-hypre=1 --download-parmetis=1
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: MatSetValuesBlocked() line 1289 in
src/mat/interface/matrix.c
breakpoint 2

Is it the reason that  MatSetValuesBlocked only works for MPIBAIJ or
SeqBAIJ?

Thanks,

Yan
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20090501/f9735e83/attachment.htm>