<br><br><div class="gmail_quote">On Tue, Jan 4, 2011 at 4:30 PM, 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;">

<div class="im"><br>
On Jan 4, 2011, at 3:45 PM, Dmitry Karpeev wrote:<br>
<br>
> VecScatter, at least the PtoP and PtoS versions, computes what amounts to a "localToGlobal" map, which is used<br>
> to set up the communication datastructures, but is later discarded.  This information appears, for example,<br>
> as the the inidx and inidy index array pair that VecScatterCreate_PtoP passes to VecScatterCreate_PtoS<br>
> and later frees.  It would be very useful to retain this information for my work on FETI-adapted matrices.<br>
> Would it be okay if I stashed this info in _p_VecScatter?<br>
><br>
</div>  Fine, whatever. What about saving it only if some option is set, so we don't waste the memory otherwise?<br></blockquote><div>I'll do that.  </div><div>Dmitry.</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">


<font color="#888888"><br>
   Barry<br>
</font><div><div></div><div class="h5"><br>
> I realize that this cuts a bit too close to PETSc's critical functionality, but it seems to me that this is only<br>
> minimally intrusive. The alternative would be to reconstruct this information from VecScatter's to/from->indices,<br>
> but these are (re)arranged by the receiving proc and searching through those isn't as convenient.  Of course,<br>
> for the special case scatters (SStoSG, MPItoAll, etc) the reconstruction has to take place in any event,<br>
> but I thought it would be useful to avoid it where the equivalent info is already computed as part of the<br>
> scatter's constructor.<br>
><br>
> Any thoughts/objections?<br>
><br>
> Dmitry.<br>
<br>
</div></div></blockquote></div><br>