[petsc-users] Global DA index using DAGetCorners

[khalid ashraf]

Sorry for the spam. Just to illustrate the problem. The following code produces 
the 
following output when run on 4 processors.

Code:
ierr=DAGetCorners(appctx->da,&xs,&ys,&zs,&xm,&ym,&zm);CHKERRQ(ierr);
 for (k=zs; k<zs+zm; k++) {
  for (j=ys; j<ys+ym; j++) {
  for (i=xs; i<xs+xm; i++) {
    {
    Vec_A[k][j][i] = (PetscReal)(k*my*mx+j*mx+i);
    }}}

Output: 
Vec_A (mx=16 X my=16 X mz=4):
0 1 2 34567
1617181920212223
3233343536373839
4849and so on

I don't understand why does the array value jump to 16 after 7. Meaning
why is the 8th element of the vector is 16. Thanks.


      
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20101125/3485910c/attachment.htm>