Hi, all,<br>I am using MPIAIJ for my matrix A, but when I call the function MatSetValuesBlocked, I got the error:<br><br>PetscPrintf(PETSC_COMM_WORLD, "breakpoint 1\n");<br> ierr = MatSetValuesBlocked(A,1,&irow,1,(col_ind+icol),temp_vector,INSERT_VALUES);<br>
PetscPrintf(PETSC_COMM_WORLD, "breakpoint 2\n"); <br><br><br>breakpoint 1<br>[0]PETSC ERROR: --------------------- Error Message ------------------------------------<br>[0]PETSC ERROR: No support for this operation for this object type!<br>
[0]PETSC ERROR: Mat type mpiaij!<br>[0]PETSC ERROR: ------------------------------------------------------------------------<br>[0]PETSC ERROR: Petsc Release Version 2.3.3, Patch 15, Tue Sep 23 10:02:49 CDT 2008 HG revision: 31306062cd1a6f6a2496fccb4878f485c9b91760<br>
[0]PETSC ERROR: See docs/changes/index.html for recent updates.<br>[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<br>[0]PETSC ERROR: See docs/index.html for manual pages.<br>[0]PETSC ERROR: ------------------------------------------------------------------------<br>
[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<br>[0]PETSC ERROR: Libraries linked from /home/vyan2000/local/PPETSc/petsc-2.3.3-p15/lib/linux-gnu-c-debug<br>
[0]PETSC ERROR: Configure run at Thu Feb 5 21:10:10 2009<br>[0]PETSC ERROR: Configure options --with-mpi-dir=/usr/lib/ --with-debugger=gdb --with-shared=0 --download-hypre=1 --download-parmetis=1<br>[0]PETSC ERROR: ------------------------------------------------------------------------<br>
[0]PETSC ERROR: MatSetValuesBlocked() line 1289 in src/mat/interface/matrix.c<br>breakpoint 2<br><br>Is it the reason that MatSetValuesBlocked only works for MPIBAIJ or SeqBAIJ? <br><br>Thanks,<br><br>Yan<br>