[petsc-users] ISCreateStride and MatGetSubmatrices

Matthew Knepley knepley at gmail.com
Tue Aug 28 10:36:26 CDT 2012

On Tue, Aug 28, 2012 at 10:33 AM, shu guo <greatgs2008 at gmail.com> wrote:

> Dear All,
> I want to extract a sequential matrix from a MPIBAIJ matrix A with a
> dimension of 4n*4n. After assembly A, I use ISCreateStride to define
> the index by
>           call ISCreateStride(PETSC_COMM_WORLD,n,0,1,isrow,ierr)
> And extract the matrix as
>           call MatGetSubMatrices(Kmatem,1,isrow,isrow,
>      &                           MAT_INITIAL_MATRIX,Kmatemseq)
> But it gives me a segmentation fault when I call ISCreateStride. Is
> there anything I did wrong? Thanks for any advice.

Run in the debugger, -start_in_debugger, and get a stack trace.


> Shu

What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120828/32798704/attachment.html>

More information about the petsc-users mailing list