<div dir="ltr">Thanks Satish. That all makes sense. <div><br></div><div>Truth is I am trying to move data from a parallel vector to a serial one on rank 0. I make an index set, in parallel, with stride 1 and use it for both source and destination and i believe that should do the job ... but it does not. Somehow it seems that the values from other processors do not get inserted in the serial vector on rank 0!</div>
<div><br></div><div>This seems it should be so simple and I'm not sure if there is a bug somewhere in my code or I'm still confused about IS and VecScatter!</div><div><br></div><div>Mohammad<br><div><br></div><div>
<br><div class="gmail_quote">On Thu, Feb 9, 2012 at 3:42 PM, Satish Balay <span dir="ltr"><<a href="mailto:balay@mcs.anl.gov">balay@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">
<div class="im">On Thu, 9 Feb 2012, Mohammad Mirzadeh wrote:<br>
<br>
</div>> So this actually means it somehow matches the local indices?<br>
<br>
Hm - no local indices here.<br>
<div class="im"><br>
> What i mean by that is,<br>
><br>
> y[ iy[i] ] = x[ ix[i] ].<br>
<br>
</div>Sure - with global indices.<br>
<div class="im"><br>
><br>
> Is that why ix and iy should have the same size?<br>
<br>
</div>When you specify data movement - you specify both source and<br>
destination for each element that is to be moved. If you are moving n<br>
elements - you have n sources, and n destination values - hence ix[n],<br>
iy[n]<br>
<span class="HOEnZb"><font color="#888888"><br>
Satish<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
><br>
> Mohammad<br>
><br>
> On Thu, Feb 9, 2012 at 3:17 PM, Satish Balay <<a href="mailto:balay@mcs.anl.gov">balay@mcs.anl.gov</a>> wrote:<br>
><br>
> > On Thu, 9 Feb 2012, Mohammad Mirzadeh wrote:<br>
> ><br>
> > > Hi guys,<br>
> > ><br>
> > > I'm just wondering if I understand how the VecScatter works. Considering<br>
> > > (petsc 3.2-p6 manual page 53):<br>
> > ><br>
> > > VecScatterCreate(Vec x,IS ix,Vec y,IS iy,VecScatter *ctx);<br>
> > > VecScatterBegin(VecScatter ctx,Vec x,Vec y,INSERT VALUES,SCATTER<br>
> > FORWARD);<br>
> > > VecScatterEnd(VecScatter ctx,Vec x,Vec y,INSERT VALUES,SCATTER FORWARD);<br>
> > > VecScatterDestroy(VecScatter *ctx);<br>
> > ><br>
> > > is the following statement correct?<br>
> > ><br>
> > > VecScatter looks into "ix" and "iy" index sets and `matches' the global<br>
> > > indecies between the two to copy data from vector "x" to vector "y". For<br>
> > > example, if "ix" maps local index "1" to global index "10", VecScatter<br>
> > > looks inside "iy" to find a local index that is mapped to global index<br>
> > "10"<br>
> > > and sends the data accordingly to the correct processor.<br>
> ><br>
> > nope - it means - if you have x[10],y[10]:<br>
> ><br>
> > ix = {1,5,9}<br>
> > iy = {0,2,1}<br>
> ><br>
> ><br>
> > Then you get:<br>
> > y[0] = x[1]<br>
> > y[2] = x[5]<br>
> > y[1] = x[9]<br>
> ><br>
> > [all numbers above are global indices]<br>
> ><br>
> > Satish<br>
> ><br>
> > ><br>
> > > Thanks,<br>
> > > Mohammad<br>
> > ><br>
> ><br>
> ><br>
><br>
<br>
</div></div></blockquote></div><br></div></div></div>