On Tue, 20 Apr 2021 at 01:06, Constantine Khrulev <[email protected]> wrote:
> Hi, > > I would like to transfer values from one DMDA-managed Vec (i.e. created > using DMCreateGlobalVector() or equivalent) to a Vec managed using a > different DMDA instance (same number of elements, same number of degrees > of freedom, *different* domain decomposition). > > What approach would you recommend? > > VecScatter. > Bonus question: what about DMDAs using different MPI communicators? > VecScatter. :D We do exactly this in PCTELESCOPE using VecScatter. It might be worth snooping through telescope_dmda.c which you can find here https://www.mcs.anl.gov/petsc/petsc-current/src/ksp/pc/impls/telescope/telescope_dmda.c.html Thanks, Dave > Thanks! > > -- > Constantine > >
