<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">Florian:<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<br>
just to make sure I don't get something wrong:<br>
<br>
An SBAIJ (especially MPISBAIJ) matrix just stores the uper triangle. When I invoke MatSetSizes with the global dimension set (and local to PETSC_DECIDE) petsc tries to evenly distribute rows amoung ranks. Same when I set local dimensions to the same on every processor.<br></blockquote><div>Yes. </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<br>
This leads to unbalanced behavior regarding the number of stored elements on each rank. The first row (on the first rank) is completely stored, the last row (on the last rank) contains just one element.<br></blockquote><div> </div><div>Not necessarily, e.g., for a sparse banded matrix, this distribution is well-balanced. User can set their own distribution by input local rows and set global rows as PETSC_DECIDE.</div><div><br></div><div>Hong</div></div><br></div></div>