<div class="gmail_quote">On Tue, Aug 28, 2012 at 10:33 AM, shu guo <span dir="ltr"><<a href="mailto:greatgs2008@gmail.com" target="_blank">greatgs2008@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear All,<br>
<br>
I want to extract a sequential matrix from a MPIBAIJ matrix A with a<br>
dimension of 4n*4n. After assembly A, I use ISCreateStride to define<br>
the index by<br>
<br>
          call ISCreateStride(PETSC_COMM_WORLD,n,0,1,isrow,ierr)<br></blockquote><div><br></div><div>Fortran does not do type conversion without interfaces (which are about as much hassle as benefit). If you configured PETSc with non-default integer sizes, be sure to declare PetscInt zero,one instead of using numeric literals.</div>
<div><br></div><div>Otherwise, run in valgrind to check that there isn't earlier corruption and run in a debugger to get a stack trace.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<br>
And extract the matrix as<br>
          call MatGetSubMatrices(Kmatem,1,isrow,isrow,<br>
     &                           MAT_INITIAL_MATRIX,Kmatemseq)<br>
<br>
But it gives me a segmentation fault when I call ISCreateStride. Is<br>
there anything I did wrong? Thanks for any advice.<br>
<span class="HOEnZb"><font color="#888888"><br>
Shu<br>
</font></span></blockquote></div><br>