On Dec 14, 2010, at 4:15 AM, Mohamad M. Nasr-Azadani wrote:
Hi guys,
Is it possible to update the ghost values from a global to a local vector for distributed arrays when global and local vectors are not from the same DA, but the global vectors are the same?
Yes
This is the the code that I have, (the only difference between the two DA's is the width. So, I am assuming that any global vector created based on those are going to be the same)
Yes
G_data is created based on DA_3D, whereas L_data2 is created based on DA_3D2.
Vec G_data, L_data; Vec G_data2, L_data2;
ierr = DACreate3d(PCW, DA_NONPERIODIC, DA_STENCIL_STAR, NX, NY, NZ, PETSC_DECIDE, PETSC_DECIDE, PETSC_DECIDE, 1, width, PETSC_NULL, PETSC_NULL, PETSC_NULL, &DA_3D); ierr = DACreate3d(PCW, DA_NONPERIODIC, DA_STENCIL_STAR, NX, NY, NZ, PETSC_DECIDE, PETSC_DECIDE, PETSC_DECIDE, 1, width+2, PETSC_NULL, PETSC_NULL, PETSC_NULL, &DA_3D2);
ierr = DACreateGlobalVector(DA_3D, &G_data); CHKERRQ(ierr); ierr = DACreateLocalVector(DA_3D, &L_data); CHKERRQ(ierr);
ierr = DACreateGlobalVector(DA_3D2, &G_data2); CHKERRQ(ierr); ierr = DACreateLocalVector(DA_3D2, &L_data2); CHKERRQ(ierr);
/* =====> Is this possible? */ ierr = DAGlobalToLocalBegin(DA_3D2, G_data, INSERT_VALUES, L_data2);CHKERRQ(ierr); ierr = DAGlobalToLocalEnd(DA_3D2, G_data, INSERT_VALUES, L_data2);CHKERRQ(ierr);
Thanks, Mohamad