On Thu, Nov 25, 2010 at 11:24, khalid ashraf <khalid_eee at yahoo.com> wrote:
> Here is the call to DA creation: > > ierr = > DACreate3d(PETSC_COMM_WORLD,DA_NONPERIODIC,DA_STENCIL_BOX,appctx.l,appctx.m,appctx.n, > > PETSC_DECIDE,PETSC_DECIDE,PETSC_DECIDE,1,1,PETSC_NULL,PETSC_NULL,PETSC_NULL, > &appctx.da);CHKERRQ(ierr); > If you created the DA as above and then call DAGetCorners(appctx.da,&xs,&ys,&zs,&xm,&ym,&zm); then only on rank=0 will xs,ys,zs all be equal to zero. Your other message about PETSc versus natural ordering indicates that this is working correctly. Is it working now? If not, please send more context, including whatever code you are using to print output. Jed -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20101125/6da538a3/attachment.htm>
