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>
