Dear Matthew:<br><br>thank you very much. Regarding my problem about MatAXPY() for MPIDense, I have checked SEQDense has this function. To my knowledge, in each node, MPIDense's format is SEQDense, right? if it is, I should ba able to write a MPIDense-based function using SEQDense? thanks a lot.<br>
<br>Regards,<br>Yujie<br><br><div class="gmail_quote">On Wed, Dec 31, 2008 at 8:56 AM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Yes, you can have 0 length IS arguments.<br><br> Matt<div><div></div><div class="Wj3C7c"><br><br><div class="gmail_quote">On Wed, Dec 31, 2008 at 10:21 AM, Yujie <span dir="ltr"><<a href="mailto:recrusader@gmail.com" target="_blank">recrusader@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hi, PETSc Developers<br><br>In parallel mode, I have 3 MPIDense matrices, A1, A2, A3. Now, I need to exact submatrices B1, B2, B3 in parallel mode from A1, A2, A3 respectively. I know I should use MatGetSubMatrix(). My problem is in <br>
<pre>MatGetSubMatrix(Mat mat,IS isrow,IS iscol,PetscInt csize,MatReuse cll,Mat *newmat)<br></pre>, the parameter "isrow" means "rows this processor should obtain". if Bi in some nodes (cpus) of the cluster (Ai have rows in them) don't have rows, is it work?<br>
<br>After getting Bi, I will create MPIDense-based C and combine them into C like<br><br> B1<br>C= B2<br> B3<br>MatGetRow();<br>MatSetValues() should work, right? <br><br>Could you give me some comments about these operations? thanks a lot. Happy new year!<br>
<br>Regards,<br><font color="#888888">Yujie<br>
</font></blockquote></div><br><br clear="all"><br></div></div><font color="#888888">-- <br>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>