<div dir="ltr"><div dir="ltr">Dear Matt,<div><br></div><div>I am trying to distribute the matrix after loading it. So I tried something like this:</div><div><br></div><div><div>  Mat B;</div><div>  PetscInt bm,bn;</div><div>  MatSetSizes(B,kk,kk,PETSC_DETERMINE,PETSC_DETERMINE);</div><div>  MatDuplicate(A,MAT_COPY_VALUES,&B);</div><div>  MatGetLocalSize(B,&bm,&bn);</div></div><div><br></div><div>where A is the original matrix (10*10) and kk is one the local sizes of A (kk=4 so I want to divide A into 4*4 and 6*6). However, I get error in MatSetSizes part and when I printed bm and bn, I get 5. In other words, B is divided equally even though I tried to divide it unequally. Am I using MatSetSizes wrong?</div><div><br></div><div>Thanks,</div><div><br></div><div>Eda</div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>>, 15 May 2019 Çar, 14:51 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"><div dir="ltr"><div dir="ltr">On Wed, May 15, 2019 at 7:35 AM Eda Oktay via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>> wrote:<br></div><div class="gmail_quote"><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">Hello,<div><br></div><div>I am trying to divide a matrix into unequal sized parts into different processors (for example I want to divide 10*10 matrix into 4*4 and 6*6 submatrix in two processors). When my program reads a matrix from file, it automatically divides it into equal parts and then I can't change local sizes.</div><div><br></div><div>How can I decompose a matrix that is read from a file?</div></div></blockquote><div><br></div><div>MatLoad() takes a matrix argument. I believe you can use MatSetSizes() before loading to get the distribution you want.</div><div><br></div><div>   Matt</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>Thanks,</div><div><br></div><div>Eda</div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail-m_-8784747921586512880gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>
</blockquote></div>