<p dir="ltr">What kind of VBR matrix? What are you partitioning using parmetis? A mesh? The blocks of the matrix? How do you create the entries in the matrix?</p>
<div class="gmail_quote">On May 14, 2013 4:36 PM, "Longxiang Chen" <<a href="mailto:suifengls@gmail.com">suifengls@gmail.com</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr"><div><div><div><div><div><div>To whom it may concern,<div><br></div>I use <span>parmetis</span> to partition a mesh for a sparse matrix.<br></div>Then I distribute the data to the  appropriate processors according to the result of partition.<br>


</div><br>The sparse matrix is stored in Variable Block Row(<span>VBR</span>) format.<br></div>After the distribution, I want to call <span>PETSc</span> <span>KSP</span> solver to solve Ax = b.<br>

</div>I tried to convert <span>VBR</span> to <span>AIJ</span> or <span>CSR</span> format, but the data would be re-distributed.<br><br>The ideal method is to keep the distribution result from <span>parmetis</span>.<br>

</div><div>For example, after <span>parmetis</span>, processor 0 has 0, 1, 4, and processor 1 has 2, 3, 5. I wish the <span>PETSc</span> would not change this distribution and solve Ax = b.<br>

</div><div><br>Are there any approaches to call <span>KSP</span> solver in <span>VBR</span> format from <span>PETSc</span>?<br></div>Or any suggestions for solving Ax = b?<br>
<br>
</div>Thanks in advance.<br><div><div><div><div><div><div><div><div><div><div><div><div><br></div><div>Regards,<div><span>Longxiang</span> Chen</div><div><br></div></div>
</div></div></div></div></div></div></div></div></div></div></div></div>
</blockquote></div>