[petsc-users] DAGetCorners; xs,ys,zs return 0 for all the processors

Jed Brown jed at 59A2.org
Thu Nov 25 04:35:43 CST 2010


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>


More information about the petsc-users mailing list