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

higher range for interpolating distances in prunevoids

Browse source
This commit is contained in:
P.M. Sutter 2014-02-19 16:33:20 -06:00
parent 10b2fdf575
commit b8ad7f8752
1 changed files with 4 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

8
c_tools/stacking/pruneVoids.cpp
View file

@ -170,8 +170,8 @@ int main(int argc, char **argv) {
double result, error; double result, error;
size_t nEval; size_t nEval;
int iZ, numZ = 4000; int iZ, numZ = 10000;
double maxZ = 2.0, z, *dL, *redshifts; double maxZ = 5.0, z, *dL, *redshifts;
dL = (double *) malloc(numZ * sizeof(double)); dL = (double *) malloc(numZ * sizeof(double));
redshifts = (double *) malloc(numZ * sizeof(double)); redshifts = (double *) malloc(numZ * sizeof(double));
for (iZ = 0; iZ < numZ; iZ++) { for (iZ = 0; iZ < numZ; iZ++) {
@ -460,7 +460,7 @@ int main(int argc, char **argv) {
// load voids *again* using Guilhem's code so we can get tree // load voids *again* using Guilhem's code so we can get tree
clock3 = clock(); clock3 = clock();
if (!args.isObservation_flag) { //if (!args.isObservation_flag) {
printf(" Re-loading voids and building tree..\n"); printf(" Re-loading voids and building tree..\n");
ZobovRep zobovCat; ZobovRep zobovCat;
if (!loadZobov(args.voidDesc_arg, args.zone2Part_arg, if (!loadZobov(args.voidDesc_arg, args.zone2Part_arg,
@ -488,7 +488,7 @@ int main(int argc, char **argv) {
} }
voids[iVoid].level = level; voids[iVoid].level = level;
} }
} // end re-load //} // end re-load
clock4 = clock(); clock4 = clock();
interval = 1.*(clock4 - clock3)/CLOCKS_PER_SEC; interval = 1.*(clock4 - clock3)/CLOCKS_PER_SEC;
printf(" Re-read voids (%.2f sec)...\n", interval); printf(" Re-read voids (%.2f sec)...\n", interval);