Zhang <[email protected]> writes: > Dear All, > Now I am confused with the way to access a global vector defined from DMDA. > > Here is the codes. When I switch on _DEBUG1_ the results get wrong. While if > it's off, > > the results fine. I just wonder why I gave to use the local form of ul,vl,wl > to access the > values such as ul[k][j][i+1], and etc.
The global vector doesn't contain space for those ghost points, so you have to scatter to the local vector if you want to access ghost points.
pgpI8Tyuf4IzW.pgp
Description: PGP signature
