[petsc-users] partition an MPIAIJ matrix

[recrusader]

I am sorry for these emails. I will further try them. Thanks a lot.

Best,
Yujie

On Mon, Feb 20, 2012 at 8:26 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote:

> 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/9ce38757/attachment.htm>