about MatGetRow/MatRestoreRow

[jiaxun hou <jiaxun_hou@xxxxxxxxxxxx>]

Dear all,

Does anyone have examples of using MatGetRow/MatRestoreRow?
I failed in using them.

My code is:

PetscInt ncols_A;
const PetscInt** cols_A_point;
const PetscScalar **vals_A_point;
for (i=0;i<row_A;i++){
 ierr = MatGetRow(A,i,&ncols_A,cols_A_point,vals_A_point);CHKERRQ(ierr);
 //do something
 ierr = MatRestoreRow(A,i,&ncols_A,cols_A_point,vals_A_point);CHKERRQ(ierr);
}

and it gets errors as below:

Petsc Release Version 2.3.1, Patch 10, Thu Mar  9 22:48:00 CST 2006
BK revision: balay@xxxxxxxxxxxxxxxxxx 0 in job 61  lab_43825   caused collective abort of all ranks
  exit status of rank 0: killed by signal 9
[cli_0]: aborting job:
application called MPI_Abort(MPI_COMM_WORLD, 64) - process 0
l.gov|ChangeSet|20060310044535|22333
See docs/changes/index.html for recent updates.
See docs/faq.html for hints about
 trouble shooting.
See docs/index.html for manual pages.
------------------------------------------------------------------------
./mytest1 on a linux-gnu named lab by root Wed Nov 14 17:57:00 2007
Libraries linked from /home/software/petsc-2.3.1-p10/lib/linux-gnu-cxx-complex-debug
Configure run at Thu Jun 15 13:08:29 2006
Configure options --with-cc=gcc --with-fc=gfortran --download-f-blas-lapack=1 --with-mpi-dir=/home/software/mpich2 --with-scalar-type=complex --with-shared=0
------------------------------------------------------------------------
[0]PETSC ERROR: PetscObjectDestroy() line 88 in src/sys/objects/destroy.c
[0]PETSC ERROR: Corrupt argument: see http://www.mcs.anl.gov/petsc/petsc-as/documentation/troubleshooting.html#Corrupt!
[0]PETSC ERROR: Invalid type of object: Parameter # 1!
[0]PETSC ERROR: PetscObjectRegisterDestroyAll() line 228 in src/sys/objects/destroy.c
[0]PETSC ERROR: PetscFinalize() line 599 in
 src/sys/objects/pinit.c
[0]PETSC ERROR: main() line 329 in /home/myprogram/mypro/mytest1.c
make: [runmytest1] Error 137 (ignored)

Can anyone tell me where is wrong? THX

---

雅虎邮箱，终生伙伴！