Changeset 41670

Timestamp:

Jun 29, 2021, 9:59:28 AM (5 years ago)

Author:

eugene

Message:

remove some commented-out blocks

File:

• trunk/Ohana/src/relphot/src/liststats.c (modified) (5 diffs)

trunk/Ohana/src/relphot/src/liststats.c

@@ -278 +278 @@
 
   // make the initial guess based on the median (not weighted mean)
-  // ALLOCATE_PTR (values, double, Npoints);
   for (int i = 0; i < Npoints; i++) {
     dataset->values[i] = dataset->flxlist[i];
@@ -288 +287 @@
   // if (!fit_least_squares (&value, dataset->flxlist, dataset->errlist, dataset->wgtlist, NULL, Npoints)) return FALSE;
   
-  // XXX add to dataset elements?
-  // ALLOCATE_PTR (wtvals, double, Npoints);
-  // ALLOCATE_PTR (wtlist, double, Npoints);
-
   int converged = FALSE;
   for (int iterations = 0; !converged && (iterations < MAX_ITERATIONS); iterations++) {
@@ -325 +320 @@
   // double WtThreshold = WEIGHT_THRESHOLD * Sum_W / (1.0 * Npoints);
   double WtThreshold = WEIGHT_THRESHOLD * WtMedian;
-
-  // generate the unmasked subset
-  // XXX add these to the dataset elements?
-  // ALLOCATE_PTR ( ykeep, double, Npoints);
-  // ALLOCATE_PTR (dykeep, double, Npoints);
-  // ALLOCATE_PTR (wtkeep, double, Npoints);
 
   // save unmasked points
@@ -343 +332 @@
       continue;
     }
-     dataset->ykeep[Nkeep] = dataset->flxlist[i];
+    dataset-> ykeep[Nkeep] = dataset->flxlist[i];
     dataset->dykeep[Nkeep] = dataset->errlist[i];
     dataset->wtkeep[Nkeep] = dataset->wgtlist[i]; // externally-supplied weight
@@ -359 +348 @@
   stats->sigma = sqrt (dSig / (Nkeep - 1));
 
-  // bootstrap resampling to generate the errorbars
-  // XXX add these to the dataset elements?
-  // ALLOCATE_PTR (ysample,  double, Nkeep);
-  // ALLOCATE_PTR (dysample, double, Nkeep);
-  // ALLOCATE_PTR (wtsample, double, Nkeep);
-  // ALLOCATE_PTR (bvalue,   double, NBOOTSTRAP); // vector to save the bootstrap values
-
   int Nboot = 0;
   for (int iboot = 0; iboot < NBOOTSTRAP; iboot++) {