petsc-3.3-p7 2013-05-11

DMDACoor2d

Structure for holding 2d (x and y) coordinates.

DMDACoor2d **coors; Vec vcoors; DM cda;

DMDAGetCoordinates(da,&vcoors); DMDAGetCoordinateDA(da,&cda); DMDAVecGetArray(cda,vcoors,&coors); DMDAGetCorners(cda,&mstart,&nstart,0,&m,&n,0) for (i=mstart; i<mstart+m; i++) { for (j=nstart; j<nstart+n; j++) { x = coors[j][i].x; y = coors[j][i].y; ...... } } DMDAVecRestoreArray(dac,vcoors,&coors);