methods on adjacency matrix

Barry Smith bsmith at mcs.anl.gov
Wed Aug 12 11:34:22 CDT 2009


On Aug 12, 2009, at 8:47 AM, Matthew Knepley wrote:

> All the matrices have the same interface. Adj is just a different  
> storage format.

    Yes, but very few methods are implemented for that type. It is  
there to support calls to partitioning packages.

    The function table is in src/mat/impls/adj/mpi/mpiadj.c it has only
MatGetRow_MPIAdj,
        MatRestoreRow_MPIAdj,
MatEqual_MPIAdj,
        MatSetOption_MPIAdj,
        MatGetRowIJ_MPIAdj,
        MatRestoreRowIJ_MPIAdj,
        MatDestroy_MPIAdj,
        MatView_MPIAdj,
        MatConvertFrom_MPIAdj,

    You can use AIJ matrix type instead.

    Barry

>
>    Matt
>
> On Wed, Aug 12, 2009 at 8:23 AM, Fredrik Bengzon <fredrik.bengzon at math.umu.se 
> > wrote:
> Hi,
> I'm trying to minimize my code managing a finite element mesh and  
> thought of using the adjacency matrix type. Is there a list of  
> methods working on this matrix type anywhere in the documentation? I  
> want to store the node-to-element information in such a matrix and  
> then use the MatGetSubMatrices and MatGetRow methods. Is this  
> possible? Thanks.
> /Fredrik Bengzon
> -- 
> 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



More information about the petsc-users mailing list