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

fixed handling of high-z boxes in shifting particles

Browse source
This commit is contained in:
P.M. Sutter 2014-02-07 10:25:03 -06:00
parent 03d28e8e03
commit edeb23cd82
1 changed files with 20 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

22
python_tools/void_python_tools/partUtil/partUtil.py
View file

@ -107,7 +107,7 @@ def loadPart(sampleDir):
isObservationData = isObservation == 1 isObservationData = isObservation == 1
return partData, boxLen, volNorm, isObservationData return partData, boxLen, volNorm, isObservationData, ranges
# ----------------------------------------------------------------------------- # -----------------------------------------------------------------------------
def loadPartVel(sampleDir): def loadPartVel(sampleDir):
@ -137,9 +137,22 @@ def loadPartVel(sampleDir):
# ----------------------------------------------------------------------------- # -----------------------------------------------------------------------------
def shiftPart(inPart, newCenter, periodicLine, boxLen): def shiftPart(inPart, center, periodicLine, ranges):
part = inPart.copy() part = inPart.copy()
newCenter = 1.*center;
boxLen = np.zeros((3))
boxLen[0] = ranges[0][1] - ranges[0][0]
boxLen[1] = ranges[1][1] - ranges[1][0]
boxLen[2] = ranges[2][1] - ranges[2][0]
# shift to box coordinates
part[:,0] -= ranges[0][0]
part[:,1] -= ranges[1][0]
part[:,2] -= ranges[2][0]
newCenter[:] -= ranges[:,0]
part[:,0] -= newCenter[0] part[:,0] -= newCenter[0]
part[:,1] -= newCenter[1] part[:,1] -= newCenter[1]
@ -155,6 +168,11 @@ def shiftPart(inPart, newCenter, periodicLine, boxLen):
if ("z" in periodicLine): part[shiftUs,2] -= \ if ("z" in periodicLine): part[shiftUs,2] -= \
np.copysign(boxLen[2], part[shiftUs,2]) np.copysign(boxLen[2], part[shiftUs,2])
#part[:,0] += ranges[0][0]
#part[:,1] += ranges[1][0]
#part[:,2] += ranges[2][0]
return part return part
# ----------------------------------------------------------------------------- # -----------------------------------------------------------------------------