yeah, I meant that when "stencilwise" access is required, then <div>VecGetDM(v,&dm);</div><div>DMVecGetStencilIterator(dm,v,&itor);</div><div><something to dereference the itor></div><div>etc</div><div>

In a way, this is similar to DAVecGetArray(), which generalizes VecGetArray </div><div>to use the geometric information about the function space (DA) the Vec belongs to.<br><br><div class="gmail_quote">On Tue, Dec 7, 2010 at 3:42 PM, Jed Brown <span dir="ltr"><<a href="mailto:jed@59a2.org">jed@59a2.org</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div class="gmail_quote"><div class="im">On Tue, Dec 7, 2010 at 22:40, Dmitry Karpeev <span dir="ltr"><<a href="mailto:karpeev@mcs.anl.gov" target="_blank">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>And yet DMDA already defines this local fine-granularity "stencil" access.<div>Where would it go, if not in the DM?  In my opinion this is part of the function space definition:</div><div>it generalizes VecGetArray to something more geometric.</div>


</div></blockquote><div><br></div></div><div>I just think it should be a DM-level call.  Not VecGetUberCustomIterator.  Perhaps I'm wrong, but I think it belongs in the DM.</div><div class="im"><div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">


<div>

<div>Not all DMs have to implement this.</div></div></blockquote><div><br></div></div><div>Maybe that's what you were saying too.</div><div><br></div><font color="#888888"><div>Jed </div></font></div>
</blockquote></div><br></div>