How to get information about the off-diagonal submatrix

Satish Balay balay at mcs.anl.gov
Mon Aug 10 09:07:38 CDT 2009


alternatively - just access the private datastructe and print..

#include "../src/mat/impls/aij/mpi/mpiaij.h"
Mat_MPIAIJ     *mpiaij = (Mat_MPIAIJ*)mat->data;
MatView(mpiaij->B,...);


Satish


On Mon, 10 Aug 2009, Matthew Knepley wrote:

> If this is just for debugging, the easiest way I think is to call
> MatGetSubmatrix() with
> a column IS that omits the local columns.
> 
>   Matt
> 
> On Mon, Aug 10, 2009 at 8:28 AM, Thomas Witkowski <
> thomas.witkowski at tu-dresden.de> wrote:
> 
> > Hi,
> >
> > how is it possible to print just the off-diagonal submatrix of one process
> > of a parallel mpiaij matrix?
> >
> > Thomas
> >
> 
> 
> 
> 



More information about the petsc-users mailing list