<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">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:1px solid 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>