cleanup from last round

This commit is contained in:
Paul M. Sutter 2025-05-24 08:56:59 -04:00
parent 05030b24e8
commit be797c05dc
3 changed files with 5 additions and 5 deletions

View file

@ -49,7 +49,7 @@ logDir = os.path.join(workDir,"logs","example_observation")
numZobovThreads = 1 numZobovThreads = 1
# optimization: number of subdivisions of the volume # optimization: number of subdivisions of the volume
numZobovDivisions = 2 numZobovDivisions = 1
# buffer size (in fraction of the boz size) when sub-diving volumes # buffer size (in fraction of the boz size) when sub-diving volumes
zobovBuffer = 0.1 zobovBuffer = 0.1

View file

@ -131,8 +131,8 @@ def launchPrep(sample, binPath, workDir=None, inputDataDir=None,
log.close() log.close()
# add tracers along the bounding box to contain the tessellation # add tracers along the bounding box to contain the tessellation
log.write("\nAdding guards along box faces...\n") #log.write("\nAdding guards along box faces...\n")
nGuard = 42 #nGuard = 42
# load in file # load in file

View file

@ -263,9 +263,10 @@ def findEdgeGalaxies(galFile, maskFile, edgeGalFile, contourFile,
contourPixels = np.nonzero(contourMap)[0] contourPixels = np.nonzero(contourMap)[0]
log.write(" We have %d rays to work with\n" % (len(contourPixels))) log.write(" We have %d rays to work with\n" % (len(contourPixels)))
for pixel in contourPixels: for pixel in contourPixels:
#print("Working with pixel %d" % pixel)
vec = healpy.pix2vec(nside,pixel) vec = healpy.pix2vec(nside,pixel)
x = raySteps * vec[0] x = raySteps * vec[0]
y = raySteps * vec[1] y = raySteps * vec[1]
z = raySteps * vec[2] z = raySteps * vec[2]
@ -275,7 +276,6 @@ def findEdgeGalaxies(galFile, maskFile, edgeGalFile, contourFile,
flagList[nearest] = 1 flagList[nearest] = 1
# flag galaxies near redsfhit boundaries # flag galaxies near redsfhit boundaries
log.write(" Flagging galaxies near redshift bounds...\n") log.write(" Flagging galaxies near redshift bounds...\n")
sphereIndices = np.arange(len(contourMap)) sphereIndices = np.arange(len(contourMap))
vec = healpy.pix2vec(nside, sphereIndices) vec = healpy.pix2vec(nside, sphereIndices)