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

optional periodic in prunevoid barycenters, updated multidark script

Browse source
This commit is contained in:
P.M. Sutter 2012-12-18 12:14:25 -06:00
parent fddaf69984
commit c578ddf1f9
2 changed files with 20 additions and 8 deletions
Download patch file Download diff file Expand all files Collapse all files

12
c_tools/stacking/pruneVoids.cpp
View file

@ -331,6 +331,13 @@ int main(int argc, char **argv) {
voids[iVoid].barycenter[1] += voids[iVoid].center[1]; voids[iVoid].barycenter[1] += voids[iVoid].center[1];
voids[iVoid].barycenter[2] += voids[iVoid].center[2]; voids[iVoid].barycenter[2] += voids[iVoid].center[2];
if (periodicX)
voids[iVoid].barycenter[0] = fmod(voids[iVoid].barycenter[0], boxLen[0]);
if (periodicY)
voids[iVoid].barycenter[1] = fmod(voids[iVoid].barycenter[1], boxLen[1]);
if (periodicZ)
voids[iVoid].barycenter[2] = fmod(voids[iVoid].barycenter[2], boxLen[2]);
// compute central density // compute central density
centralRad = voids[iVoid].radius/args_info.centralRadFrac_arg; centralRad = voids[iVoid].radius/args_info.centralRadFrac_arg;
centralRad *= centralRad; centralRad *= centralRad;
@ -476,13 +483,14 @@ int main(int argc, char **argv) {
voids[iVoid].accepted = -2; voids[iVoid].accepted = -2;
} }
// *alwas* clean out near edges since there are no mocks there // *always* clean out near edges since there are no mocks there
if (tolerance*voids[iVoid].maxRadius > voids[iVoid].nearestEdge) { if (tolerance*voids[iVoid].maxRadius > voids[iVoid].nearestEdge) {
voids[iVoid].accepted = -3; voids[iVoid].accepted = -3;
} }
// assume the lower z-boundary is "soft" in observations // assume the lower z-boundary is "soft" in observations
if (voids[iVoid].redshift < args_info.zMin_arg) { if (args_info.isObservation_flag &&
voids[iVoid].redshift < args_info.zMin_arg) {
voids[iVoid].accepted = -3; voids[iVoid].accepted = -3;
} }
} }

16
pipeline/datasets/multidark.py
View file

@ -33,13 +33,16 @@ particleFileBase = "mdr1_particles_z"
# list of file numbers for the particle files # list of file numbers for the particle files
# to get particle file name, we take particleFileBase+fileNum # to get particle file name, we take particleFileBase+fileNum
fileNums = (("0.0", "0.53", "1.0")) #fileNums = ["0.53"]
fileNums = ["0.0", "0.53", "1.0"]
# redshift of each file in the above list # redshift of each file in the above list
redshifts = (("0.0", "0.53", "1.0")) #redshifts = ["0.53"]
redshifts = ["0.0", "0.53", "1.0"]
# how many independent slices along the z-axis? # how many independent slices along the z-axis?
numSlices = 4 numSlices = 1
#numSlices = 4
# how many subdivisions along the x- and y- axis? # how many subdivisions along the x- and y- axis?
# ( = 2 will make 4 subvolumes for each slice, = 3 will make 9, etc.) # ( = 2 will make 4 subvolumes for each slice, = 3 will make 9, etc.)
@ -57,7 +60,8 @@ haloFileBase = "mdr1_halos_z"
# minimum halo mass cuts to apply for the halo catalog # minimum halo mass cuts to apply for the halo catalog
# use "none" to get all halos # use "none" to get all halos
minHaloMasses = ["none", 1.2e13] minHaloMasses = [1.2e13]
#minHaloMasses = ["none", 1.2e13]
# locations of data in the halo catalog # locations of data in the halo catalog
haloFileMCol = 6 haloFileMCol = 6
@ -73,8 +77,8 @@ haloFileNumComLines = 0
# adjust these two parameters given the memory contraints on your system: # adjust these two parameters given the memory contraints on your system:
# numZobovDivisions: how many sub-volumes per dimension will zobov process # numZobovDivisions: how many sub-volumes per dimension will zobov process
# numZobovThreads: how many sub-volumes to process at once? # numZobovThreads: how many sub-volumes to process at once?
numZobovDivisions = 2 numZobovDivisions = 4
numZobovThreads = 2 numZobovThreads = 4
# simulation information # simulation information
numPart = 1024*1024*1024 numPart = 1024*1024*1024