<div class="gmail_quote">On Tue, Dec 7, 2010 at 20:53, 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;">
BTW: I consider the DMDA Vec just a very special case that gets way to much attention and can sometimes distort the overall design of PETSc in issues that are not important in that non-special case.</blockquote></div><br>
<div>Along these lines, I'm curious about using an API like VecGhost for the DA.  That is, what if there was a level where a Vec could know about a local form.  When you get the local form, one implementation would return the contiguous local form (no copy) and another implementation would return a separate vector.  Then you could VecUpdateVector(Xglobal,Xlocal,imode,scattermode) or so (need an unambiguous name).  Then DMComposite could operate without forcing separate local and global vectors.</div>
<div><br></div><div>Jed</div>