harmonized array ordering

This commit is contained in:
P.M. Sutter 2014-05-30 11:05:37 -05:00
parent fe885ea008
commit 42ff05eab1
2 changed files with 8 additions and 8 deletions

View file

@ -542,10 +542,10 @@ def stackVoids(catalog, stackMode = "ball"):
localPart = catalog.partPos[ getBall(partTree, center, rMax) ]
else:
voidPart = getVoidPart(catalog, void.voidID)
localPart = np.zeros((3,len(voidPart)))
localPart[0,:] = getArray(voidPart, 'x')
localPart[1,:] = getArray(voidPart, 'y')
localPart[2,:] = getArray(voidPart, 'z')
localPart = np.zeros((len(voidPart),3))
localPart[:,0] = getArray(voidPart, 'x')
localPart[:,1] = getArray(voidPart, 'y')
localPart[:,2] = getArray(voidPart, 'z')
shiftedPart = shiftPart(localPart, center, periodicLine, catalog.ranges)

View file

@ -253,13 +253,13 @@ def plotVoidCells(catalog,
voidPart = getVoidPart(catalog, voidID)
newpart = np.zeros((3,len(voidPart)))
newpart = np.zeros((len(voidPart),3))
volume=np.zeros(len(voidPart))
radius=np.zeros(len(voidPart))
newpart[0,:] = getArray(voidPart, 'x')
newpart[1,:] = getArray(voidPart, 'y')
newpart[2,:] = getArray(voidPart, 'z')
newpart[:,0] = getArray(voidPart, 'x')
newpart[:,1] = getArray(voidPart, 'y')
newpart[:,2] = getArray(voidPart, 'z')
volume = getArray(voidPart, 'volume')
radius = (3.*volume/(4.*np.pi))**(1./3.)