<div class="gmail_quote">On Thu, Jan 26, 2012 at 12:24, Dominik Szczerba <span dir="ltr"><<a href="mailto:dominik@itis.ethz.ch">dominik@itis.ethz.ch</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div id=":1kd">I am having problems with the following piece of code in petsc 3.2 but<br>
not 3.1 (with syntax changes where appropriate):<br>
<br>
Mat XBI;<br>
Mat adj = 0;<br>
MatPartitioning part = 0;<br>
IS is = 0;<br>
<br>
// Now declare XBI as MPIAIJ and fill it.<br>
<br>
MatConvert(XBI, MATMPIADJ, MAT_INITIAL_MATRIX, &adj);<br>
MatPartitioningCreate(PETSC_COMM_WORLD, &part);<br>
MatPartitioningSetAdjacency(part, adj);<br>
MatPartitioningSetNParts(part, npMax);<br>
MatPartitioningSetFromOptions(part);<br>
MatPartitioningApply(part, &is);<br>
<br>
It crashes in parmetis at the last line above and I have already<br>
reported that in their bug system.<br>
I want to eliminate the changes between Petsc 3.1/3.2 API as a<br>
possible cause - would there be any changes possibly affecting<br>
behavior this code?</div></blockquote></div><br><div>PETSc-3.1 with --download-parmetis uses ParMetis-3.1.1 with some custom patches.</div><div><br></div><div>PETSc-3.2 uses ParMetis-3.2 with other custom patches.</div><div>
<br></div><div>PETSc-dev uses ParMetis-4.0.2 with other custom patches.</div><div><br></div><div><br></div><div>I suggest that you try building your code with petsc-dev to see if ParMetis-4.0.2 fixed the bug you are hitting.</div>