<p>Dear Barry:</p><p>I have a new question about the parameter "iy"(new index set of subvector). To parallel vector, how to provide "ix" and "iy"?<br></p><p>Just providing the local index subset for "ix" and "iy"? if it is, it is a little difficult to let the local vector know its global position in the new subvector? If the user needs to provide the global position of local vector for "iy". some MPI communication should be needed, it looks like not a good method. In exacting submatrix, the function hides this problem. could you give me any comments? thanks a lot.</p>
<p>Regards,</p><p>Yujie </p><br><div class="gmail_quote">On Fri, Jan 2, 2009 at 11:46 AM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<br>
VecScatter is for this purpose.<br>
<br>
Rational: extracting subparts of vectors for ghost points etc takes place many times in a simulation; maybe millions.<br>
Thus separating it into a set-up followed by many uses is a worthwhile optimization. Extracting submatrices occur<br>
much less often in a simulation, maybe tens, hundreds or thousands of times so it is not worth the extra complexity<br>
of having separate set-up followed by many uses. One could argue that uniformity of design means we should have<br>
handled matrices with a MatScatter concept to parallel the Vec approach, but it is too late now :-).<br><font color="#888888">
<br>
Barry</font><div><div class="Wj3C7c"><br>
<br>
On Jan 2, 2009, at 1:41 PM, Yujie wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Like MatGetSubMatrix(), whether is there a function to get a subvector in parallel mode? I have checked some scatter functios, they don't likely work. thanks a lot.<br>
<br>
Regards,<br>
<br>
Yujie<br>
<br>
</blockquote>
<br>
</div></div></blockquote></div><br>