Index: /trunk/Ohana/src/relphot/src/liststats.c
===================================================================
--- /trunk/Ohana/src/relphot/src/liststats.c	(revision 41669)
+++ /trunk/Ohana/src/relphot/src/liststats.c	(revision 41670)
@@ -278,5 +278,4 @@
 
   // 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,8 +287,4 @@
   // 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,10 +320,4 @@
   // 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,5 +332,5 @@
       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,11 +348,4 @@
   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++) {