<p>Dear Matthew:</p><p>Two processors are used. The matrix dimension is 105*108;<br></p><p>MatGetOwnershipRange() </p><p>proc1: 0->56</p><p>proc2: 56->105</p><p>Mat->camp.rstart; Mat->cmap.rend</p><p>proc1: 0->54</p>
<p>proc2: 54->108</p><p>MatGetLocalSize()</p><p>proc1: row 56 col 54</p><p>porc2: row 49 col 54</p><p>MatDenseGetLocalMatrix()</p><p>proc1: 56*108</p><p>proc2: 49*108</p><p>thanks.</p><p>Yujie</p><br><div class="gmail_quote">
On Tue, Jan 6, 2009 at 11:33 AM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div><div class="Wj3C7c">On Tue, Jan 6, 2009 at 11:52 AM, Yujie <span dir="ltr"><<a href="mailto:recrusader@gmail.com" target="_blank">recrusader@gmail.com</a>></span> wrote:<br></div></div><div class="gmail_quote">
<div><div class="Wj3C7c"><blockquote class="gmail_quote" style="border-left:1px solid rgb(204, 204, 204);margin:0pt 0pt 0pt 0.8ex;padding-left:1ex">
<p>Dear PETSc developers:</p><p>I am trying to use MatGetArray() and MatSetValues() to combine several MPIDense matrices into one matrix. At the beginning, I use </p><p>MatGetOwnershipRange() and Mat->camp.rstart; Mat->cmap.rend (2.3.3-p8 version) to get the start, end row and column. I can calculate the local rows and columns.</p>
<p>I also use MatGetLocalSize() to confirm the accuracy. However, I always find some data loses in the combined matrix.</p><p>And then, I try to use MatDenseGetLocalMatrix() to get the lcoal matrix and output it. I find column information by MatGetLocalSize() is not consistent with by MatDenseGetLocalMatrix(), is it bug? could you give me some advice? thanks a lot.</p>
</blockquote></div></div><div>What information?<br><br> Matt<br> </div><blockquote class="gmail_quote" style="border-left:1px solid rgb(204, 204, 204);margin:0pt 0pt 0pt 0.8ex;padding-left:1ex">
<p>Regards,</p><p>Yujie</p>
</blockquote></div>-- <br><font color="#888888">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<br>
</font></blockquote></div><br>