<div dir="ltr"><div dir="ltr"><br><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Feb 8, 2024 at 3:15 AM Maruthi NH <<a href="mailto:maruthinh@gmail.com">maruthinh@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi Barry, <br></div><div>Thanks. Yes, the global column index was wrong. I have one more question regarding 
MatCreateMPIAIJWithArrays. If I have 100 elements in rank 0 and 50 in rank 1, does PETSc redistribute equally among procs before solving?</div></div></blockquote><div><br></div><div>No, the m  you provided in MatCreateMPIAIJWithArrays(comm, m, n, ..., i, j, ..) determines the number of local rows and petsc won't change that. </div><div>If you want petsc to redistribute it, first create a balanced matrix and then call MatSetValues() to insert values.  Petsc will ship these values around to match the matrix layout you prescribed. </div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div> <br></div><div><br></div><div>Regards,</div><div>Maruthi <br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Feb 5, 2024 at 2:18 AM Barry Smith <<a href="mailto:bsmith@petsc.dev" target="_blank">bsmith@petsc.dev</a>> wrote:<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>
   Is each rank trying to create its own sequential matrix with MatCreateSeqAIJWithArrays() or did you mean MatCreateMPIAIJWithArrays()?<br>
<br>
   If the latter, then possibly one of your size arguments is wrong or the indices are incorrect for the given sizes.<br>
<br>
   Barry<br>
<br>
<br>
> On Feb 4, 2024, at 3:15 PM, Maruthi NH <<a href="mailto:maruthinh@gmail.com" target="_blank">maruthinh@gmail.com</a>> wrote:<br>
> <br>
> Hi all,<br>
> <br>
> I have a row, col, and A values in CSR format; let's say rank 0 has 200 unknowns and rank 1 has 100 unknowns. If I use MatCreateSeqAIJWithArrays to create a Matrix, it crashes. However, if each rank has an equal number of unknowns, it works fine. Please let me know how to proceed<br>
> <br>
> <br>
> Regards,<br>
> Maruthi<br>
<br>
</blockquote></div>
</blockquote></div></div>