<div class="gmail_quote">On Thu, Nov 25, 2010 at 11:24, khalid ashraf <span dir="ltr"><<a href="mailto:khalid_eee@yahoo.com">khalid_eee@yahoo.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div>Here is the call to DA creation:</div><div><br></div><div>ierr = DACreate3d(PETSC_COMM_WORLD,DA_NONPERIODIC,DA_STENCIL_BOX,appctx.l,appctx.m,appctx.n,</div><div> PETSC_DECIDE,PETSC_DECIDE,PETSC_DECIDE,1,1,PETSC_NULL,PETSC_NULL,PETSC_NULL,</div>
<div> &appctx.da);CHKERRQ(ierr);</div></blockquote></div><br><div>If you created the DA as above and then call</div><div><br></div><div> DAGetCorners(appctx.da,&xs,&ys,&zs,&xm,&ym,&zm);</div>
<div><br></div><div>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.</div>
<div><br></div><div>Jed</div>