[petsc-users] partition an MPIAIJ matrix

[Jed Brown]

On Mon, Feb 20, 2012 at 20:22, recrusader <recrusader at gmail.com> wrote:

> What types of AIJ matrices can directly use partitioning implementation?
> The matrix I am using doesn't work?
> I have to use ParMETIS since it can provide fill-reducing matrix ordering.
> From my understanding, MatConvert can only reorder the matrix and doesn't
> provide the reordering information.
>

We've exchanged a lot of emails now, but I'm still having trouble
extracting context from your confusion. You can

1. Make a normal matrix (e.g. AIJ), create a MatPartitioning object and
call MatPartitioningApply().

2. Create an MPIAdj matrix (perhaps by creating an AIJ matrix and using
MatConvert()), then use either MatPartitioning or (if you insist) calling
ParMETIS directly (you'll write the same code as is in
src/mat/partition/impls/pmetis/pmetis.c).
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120220/cdcd2c8d/attachment.htm>