impls/redundant/redundant.c

11   Vec                xdup, ydup;            /* parallel vector that congregates xsub or ysub facilitating vector scattering */
115       /* create working vectors xdup and ydup.
117        ydup concatenates all ysub and has empty local arrays because ysub's arrays will be place into it.
121       PetscCall(VecCreateMPIWithArray(PetscSubcommContiguousParent(red->psubcomm), 1, mloc_sub, PETSC_DECIDE, NULL, &red->ydup));
206   /* place ysub's local array into ydup */
208   PetscCall(VecPlaceArray(red->ydup, (const PetscScalar *)array));
210   /* scatter ydup to y */
211   PetscCall(VecScatterBegin(red->scatterout, red->ydup, y, INSERT_VALUES, SCATTER_FORWARD));
212   PetscCall(VecScatterEnd(red->scatterout, red->ydup, y, INSERT_VALUES, SCATTER_FORWARD));
213   PetscCall(VecResetArray(red->ydup));
244   /* place ysub's local array into ydup */
246   PetscCall(VecPlaceArray(red->ydup, (const PetscScalar *)array));
248   /* scatter ydup to y */
249   PetscCall(VecScatterBegin(red->scatterout, red->ydup, y, INSERT_VALUES, SCATTER_FORWARD));
250   PetscCall(VecScatterEnd(red->scatterout, red->ydup, y, INSERT_VALUES, SCATTER_FORWARD));
251   PetscCall(VecResetArray(red->ydup));
267     PetscCall(VecDestroy(&red->ydup));