I am sorry for these emails. I will further try them. Thanks a lot.<br><br>Best,<br>Yujie <br><br><div class="gmail_quote">On Mon, Feb 20, 2012 at 8:26 PM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov">jedbrown@mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im"><div class="gmail_quote">On Mon, Feb 20, 2012 at 20:22, recrusader <span dir="ltr"><<a href="mailto:recrusader@gmail.com" target="_blank">recrusader@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>What types of AIJ matrices can directly use partitioning implementation? The matrix I am using doesn't work?<br>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.</div>
</blockquote></div><br></div><div>We've exchanged a lot of emails now, but I'm still having trouble extracting context from your confusion. You can</div><div><br></div><div>1. Make a normal matrix (e.g. AIJ), create a MatPartitioning object and call MatPartitioningApply().</div>
<div><br></div><div>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).</div>
</blockquote></div><br>