#include "petscdm.h" #include "petscdmlabel.h" #include "petscds.h" PetscErrorCode DMGetCoordinatesLocalTuple(DM dm, IS p, PetscSection *pCoordSection, Vec *pCoord)Not collective
dm | - the DM | |
p | - the IS of points whose coordinates will be returned |
pCoordSection | - the PetscSection describing the layout of pCoord, i.e. each point corresponds to one point in p, and DOFs correspond to coordinates | |
pCoord | - the Vec with coordinates of points in p |
This creates a new vector, so the user SHOULD destroy this vector
Each process has the local and ghost coordinates
For DMDA, in two and three dimensions coordinates are interlaced (x_0,y_0,x_1,y_1,...) and (x_0,y_0,z_0,x_1,y_1,z_1...)