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

estimate halo density taken from highest redshiftfor universal minVoidRadius

Browse source
This commit is contained in:
P.M. Sutter 2014-02-10 03:23:41 -06:00
parent a9e040540b
commit 5ab0a34129
1 changed files with 4 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

7
python_tools/pipeline_source/prepareCatalogs.in.py
View file

@ -604,11 +604,12 @@ if (args.script or args.all) and haloFileBase != "":
sys.stdout.flush() sys.stdout.flush()
# estimate number of halos to get density # estimate number of halos to get density
# take from *highest* redshift to get lowest density
if haloFileDummy == '': if haloFileDummy == '':
dataFile = catalogDir+haloFileBase+fileNums[0] dataFile = catalogDir+haloFileBase+fileNums[-1]
else: else:
dataFile = catalogDir+haloFileBase.replace(haloFileDummy, dataFile = catalogDir+haloFileBase.replace(haloFileDummy,
fileNums[0]) fileNums[-1])
numPart = 0 numPart = 0
if dataFormat == "sdf": if dataFormat == "sdf":
SDFcvt_PATH = "@CMAKE_BINARY_DIR@/external/libsdf/apps/SDFcvt/SDFcvt.x86_64" SDFcvt_PATH = "@CMAKE_BINARY_DIR@/external/libsdf/apps/SDFcvt/SDFcvt.x86_64"
@ -865,7 +866,7 @@ if (args.halos or args.all) and haloFileBase != "":
exit(-1) exit(-1)
numKept = 0 numKept = 0
inFile = open(dataFile, 'r') inFile = open(tempFile, 'r')
outFile = open(outFileName, 'a') outFile = open(outFileName, 'a')
for (iLine,line) in enumerate(inFile): for (iLine,line) in enumerate(inFile):
if np.random.uniform() > keepFraction: continue if np.random.uniform() > keepFraction: continue