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

bug fixes in profiles

Browse source
This commit is contained in:
P.M. Sutter 2014-05-14 18:35:50 -05:00
parent f905ae92af
commit 9791ddcfa3
1 changed files with 2 additions and 1 deletions
Download patch file Download diff file Expand all files Collapse all files

3
python_tools/void_python_tools/voidUtil/profileUtil.py
View file

@ -67,7 +67,7 @@ def buildProfile(catalog, rMin, rMax):
shiftedPart = shiftPart(localPart, center, periodicLine, catalog.ranges) shiftedPart = shiftPart(localPart, center, periodicLine, catalog.ranges)
dist = np.sqrt(np.sum(shiftedPart[:,:]**2, axis=1)) dist = np.sqrt(np.sum(shiftedPart[:,:]**2, axis=1))
thisProfile, radii = np.histogram(dist, bins=nBins, range=(0,rMaxProfile)) thisProfile, radii = np.histogram(dist, bins=10, range=(0,rMaxProfile))
deltaV = 4*np.pi/3*(radii[1:]**3-radii[0:(radii.size-1)]**3) deltaV = 4*np.pi/3*(radii[1:]**3-radii[0:(radii.size-1)]**3)
thisProfile = np.float32(thisProfile) thisProfile = np.float32(thisProfile)
thisProfile /= deltaV thisProfile /= deltaV
@ -75,6 +75,7 @@ def buildProfile(catalog, rMin, rMax):
allProfiles.append(thisProfile) allProfiles.append(thisProfile)
binCenters = 0.5*(radii[1:] + radii[:-1]) binCenters = 0.5*(radii[1:] + radii[:-1])
nVoid = len(voidsToStack)
stackedProfile = np.std(allProfiles, axis=0) / np.sqrt(nVoids) stackedProfile = np.std(allProfiles, axis=0) / np.sqrt(nVoids)
sigmas = np.std(allProfiles, axis=0) / np.sqrt(nVoids) sigmas = np.std(allProfiles, axis=0) / np.sqrt(nVoids)