<div class="gmail_quote">On Tue, Dec 7, 2010 at 22:52, Dmitry Karpeev <span dir="ltr"><<a href="mailto:karpeev@mcs.anl.gov">karpeev@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 id=":ev">Actually, I don't think this would require DM to know that much than it does now:</div></blockquote><div><br></div><div>I don't see how this is DM-generic, it seems to be specialized to certain DMs.  Yes, sometimes I would like to have "one DM to rule them all", but I think there will always be value in specialization, for performance or for special applications.  Maybe, if we had a proven system for traversal that we wanted to promote, it could be moved up from DMOne to DM, but I'm skeptical about this.</div>
<div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div id=":ev"><div>DMDA has GlobalToLocal as well as VecGetArray, which provides both "restriction to local" (a coarse-granularity, heavy-weight operation), and "restrict to (a bunch of) stencil(s)" (a fine granularity lightweight "scatter").  </div>


<div>This encapsulates your notions of a "local space" and an "extended space", which, clearly, require</div><div>different access patterns.</div></div></blockquote></div><br><div>A distinction is that we need iterators on the expanded space to be non-collective.</div>
<div><br></div><div>Jed</div>