<pre><tt><tt>Dear all,<br><br>Does anyone have examples of using MatGetRow/MatRestoreRow?<br>I failed in using them.<br><br>My code is:<br><br>PetscInt </tt></tt><tt><tt><tt><tt>ncols_A</tt></tt></tt></tt>;<br><tt><tt>const PetscInt** cols_A_point;<br>const PetscScalar **vals_A_point;<br>for (i=0;i<row_A;i++){<br> ierr = MatGetRow(A,i,&ncols_A,cols_A_point,vals_A_point);CHKERRQ(ierr);<br> //do something<br> ierr = MatRestoreRow(A,i,&ncols_A,cols_A_point,vals_A_point);CHKERRQ(ierr);<br>}<br><br>and it gets errors as below:<br><br>Petsc Release Version 2.3.1, Patch 10, Thu Mar 9 22:48:00 CST 2006<br>BK revision: balay@asterix.mcs.anrank 0 in job 61 lab_43825 caused collective abort of all ranks<br> exit status of rank 0: killed by signal 9<br>[cli_0]: aborting job:<br>application called MPI_Abort(MPI_COMM_WORLD, 64) - process 0<br>l.gov|ChangeSet|20060310044535|22333<br>See docs/changes/index.html for recent updates.<br>See docs/faq.html for hints about
trouble shooting.<br>See docs/index.html for manual pages.<br>------------------------------------------------------------------------<br>./mytest1 on a linux-gnu named lab by root Wed Nov 14 17:57:00 2007<br>Libraries linked from /home/software/petsc-2.3.1-p10/lib/linux-gnu-cxx-complex-debug<br>Configure run at Thu Jun 15 13:08:29 2006<br>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<br>------------------------------------------------------------------------<br>[0]PETSC ERROR: PetscObjectDestroy() line 88 in src/sys/objects/destroy.c<br>[0]PETSC ERROR: Corrupt argument: see http://www.mcs.anl.gov/petsc/petsc-as/documentation/troubleshooting.html#Corrupt!<br>[0]PETSC ERROR: Invalid type of object: Parameter # 1!<br>[0]PETSC ERROR: PetscObjectRegisterDestroyAll() line 228 in src/sys/objects/destroy.c<br>[0]PETSC ERROR: PetscFinalize() line 599 in
src/sys/objects/pinit.c<br>[0]PETSC ERROR: main() line 329 in /home/myprogram/mypro/mytest1.c<br>make: [runmytest1] Error 137 (ignored)<br><br>Can anyone tell me where is wrong? THX<br></tt></tt></pre><p> 
<hr size=1><a href="https://member.cn.yahoo.com/cnreg/reginfo.html?id=89034" target=blank>ÑÅ»¢ÓÊÏ䣬ÖÕÉú»ï°é£¡</a>