Grid3D getFaceCenters incorrect when mesh is offset

[guyer]

(as reported in http://article.gmane.org/gmane.comp.python.fipy/1418)

The Grid3D getFaceCenters doesn't appear to give correct
printout result when the mesh is offset. Below is my script:

from fipy import *

region_3D =  Grid3D(dx=1,dy=1,dz=1,nx=2,ny=2,nz=2) + ((0,),(0,),(2,))
x,y,z = region_3D.getFaceCenters()

Then, if I print out z vector, I don't see offset being applied on z-axis:

z = [ 0. ,  0. ,  0. ,  0. ,  1. ,  1. ,  1. ,  1. ,  2. ,  2. ,  2. ,
2. ,  0.5,  0.5,  0.5,  0.5,  0.5,  0.5,  1.5,  1.5,  1.5,  1.5,
1.5,  1.5,  0.5,  0.5,  0.5,  0.5,  0.5,  0.5,  1.5,  1.5,  1.5,
1.5,  1.5,  1.5]

The getCellCenters works just fine. The z vector from getCellCenters gives:

x,y,z =  region_3D.getCellCenters()
array([ 2.5,  2.5,  2.5,  2.5,  3.5,  3.5,  3.5,  3.5])

Imported from trac ticket #190, created by guyer on 06-09-2009 at 13:45, last modified: 06-09-2009 at 14:04

[guyer]

Fixed in r3052 and merged to trunk in r3056

Trac comment by guyer on 06-09-2009 at 14:04