guilhem_lavaux/vide_public
0
0
Fork 0
mirror of https://bitbucket.org/cosmicvoids/vide_public.git synced 2025-07-04 15:21:11 +00:00
Code Issues Projects Releases Packages Wiki Activity Actions

Merge branch 'master' of https://bitbucket.org/cosmicvoids/vide_public

Browse source
This commit is contained in:
P.M. Sutter 2014-05-15 18:34:25 -05:00
commit a2239750bb
2 changed files with 39 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

38
python_tools/void_python_tools/voidUtil/catalogUtil.py
View file

@ -492,3 +492,41 @@ def filterVoidsOnCentralDen(catalog, maxCentralDen):
return catalog
# -----------------------------------------------------------------------------
def stackVoids(catalog, stackMode = "ball"):
# builds a stack of voids from the given catalog
# catalog: void catalog
# stackMode:
# "ball": spherical cut
# "voronoi": only void member particles
#
# returns:
# stackedPart: array of relative particle positions in the stack
rMax = 100.
periodicLine = getPeriodic(catalog.sampleInfo)
if stackMode == "ball":
partTree = getPartTree(catalog)
stackedPart = []
for void in catalog.voids:
center = void.barycenter
if stackMode == "ball":
localPart = catalog.partPos[ getBall(partTree, center, rMax) ]
else:
voidPart = getVoidPart(catalog, void.voidID)
localPart = np.zeros((3,len(voidPart)))
localPart[0,:] = getArray(localPart, 'x')
localPart[1,:] = getArray(localPart, 'y')
localPart[2,:] = getArray(localPart, 'z')
shiftedPart = shiftPart(localPart, center, periodicLine, catalog.ranges)
stackedPart.extend(shiftedPart)
return stackedPart