<div dir="ltr">For this particular case (one subdoanin), it may be easy to fix in petsc. We could create a partitioning index filled with zeros. <div><br></div><div>Fande,</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Aug 17, 2020 at 5:04 PM Fande Kong <<a href="mailto:fdkong.jd@gmail.com">fdkong.jd@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">IIRC, Chaco does not produce an arbitrary number of subdomains. The number needs to be like 2^n.</div><div dir="ltr"><br></div><div dir="ltr">ParMETIS and PTScotch are much better, and they are production-level code. If there is no particular reason, I would like to suggest staying with ParMETIS and PTScotch.</div></div><div dir="ltr"><br></div><div dir="ltr">Thanks,</div><div dir="ltr"><br></div><div dir="ltr">Fande,<br><div><br></div><div><br></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Aug 14, 2020 at 10:07 AM Eda Oktay <<a href="mailto:eda.oktay@metu.edu.tr" target="_blank">eda.oktay@metu.edu.tr</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr">Dear Barry,<div><br></div><div>Thank you for answering. I am sending a sample code and a binary file. </div><div><br></div><div>Thanks!</div><div><br></div><div>Eda</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Barry Smith <<a href="mailto:bsmith@petsc.dev" target="_blank">bsmith@petsc.dev</a>>, 14 Ağu 2020 Cum, 18:49 tarihinde şunu yazdı:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><br>
   Could be a bug in Chaco or its call from PETSc for the special case of one process. Could you send a sample code that demonstrates the problem?<br>
<br>
  Barry<br>
<br>
<br>
> On Aug 14, 2020, at 8:53 AM, Eda Oktay <<a href="mailto:eda.oktay@metu.edu.tr" target="_blank">eda.oktay@metu.edu.tr</a>> wrote:<br>
> <br>
> Hi all,<br>
> <br>
> I am trying to try something. I am using the same MatPartitioning codes for both CHACO and ParMETIS:<br>
> <br>
> ierr = MatConvert(SymmA,MATMPIADJ,MAT_INITIAL_MATRIX,&AL);CHKERRQ(ierr);  <br>
>   ierr = MatPartitioningCreate(MPI_COMM_WORLD,&part);CHKERRQ(ierr);<br>
>   ierr = MatPartitioningSetAdjacency(part,AL);CHKERRQ(ierr); <br>
>   <br>
>   ierr = MatPartitioningSetFromOptions(part);CHKERRQ(ierr);<br>
>   ierr = MatPartitioningApply(part,&partitioning);CHKERRQ(ierr);<br>
> <br>
> After obtaining the IS, I apply this to my original nonsymmetric matrix and try to get an approximate edge cut. <br>
> <br>
> Except for 1 partitioning, my program completely works for 2,4 and 16 partitionings. However, for 1, ParMETIS gives results where CHACO I guess doesn't since I am getting errors about the index set.<br>
> <br>
> What is the difference between CHACO and ParMETIS that one works for 1 partitioning and one doesn't?<br>
> <br>
> Thanks!<br>
> <br>
> Eda<br>
<br>
</blockquote></div>
</blockquote></div>
</blockquote></div>