[petsc-dev] DMDA_*PERIODIC and DMDA_XYZGHOSTED

Barry Smith bsmith at mcs.anl.gov
Tue Dec 7 13:45:39 CST 2010


   DMDA_XYZGHOSTED does not exist for 2d and 3d it was added, I'm guessing, as an experiment and was never in the initial design of DMDA. To fully support it one needs to go back tot he design of DMDA and see how to have it properly done and not just bolt it on.  Some people like to use these types of ghost nodes so I agree it is a useful thing to have but who is going to properly add it?


On Dec 7, 2010, at 1:30 PM, Jed Brown wrote:

> On Tue, Dec 7, 2010 at 20:21, Ethan Coon <ecoon at lanl.gov> wrote:
> 'd like a DA where there are ghost cells on every boundary, and some of
> those ghost cells (but not all) are filled in with periodic values.
> 
> It would be useful to people doing explicit stuff if there was a way to get ghost nodes in the local vector without implying periodic communication (and weird coordinate management).
> 
> A related issue for purely explicit is to have a way to VecAXPY without needing to copy to and from a global vector.  (TSSSP has low-memory schemes, paying for an extra vector or two is actually significant in that context, and (less significant) I'm certain I can cook up a realistic benchmark where the memcpy costs more than 10%.)  I think I know how to implement this sharing transparently (more-or-less using VecNest) so we could make it non-default but be able to activate it at runtime.

   Why can you not use VecAXPY() on the local Vecs?

   Barry



> 
> Jed




More information about the petsc-dev mailing list