Index: trunk/psLib/src/math/psMinimizePolyFit.c
===================================================================
--- trunk/psLib/src/math/psMinimizePolyFit.c	(revision 6101)
+++ trunk/psLib/src/math/psMinimizePolyFit.c	(revision 6185)
@@ -10,6 +10,6 @@
  *  @author EAM, IfA
  *
- *  @version $Revision: 1.1 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2006-01-21 02:43:31 $
+ *  @version $Revision: 1.2 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-01-23 20:44:29 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -862,5 +862,4 @@
     psF32 minClipSigma;
     psF32 maxClipSigma;
-
     if (isnan(stats->max) || isfinite(stats->max)) {
         maxClipSigma = fabs(stats->clipSigma);
@@ -873,4 +872,5 @@
         minClipSigma = fabs(stats->min);
     }
+
     psVector *fit   = NULL;
     psVector *resid = psVectorAlloc(f->n, PS_TYPE_F64);
@@ -881,12 +881,11 @@
     // for now, for the SAMPLE_MEDIAN and SAMPLE_STDEV to be used
     stats->options |= (PS_STAT_SAMPLE_MEDIAN | PS_STAT_SAMPLE_STDEV);
-
     psTrace(__func__, 4, "stats->clipIter is %d\n", stats->clipIter);
     psTrace(__func__, 4, "(minClipSigma, maxClipSigma) is (%.2f, %.2f)\n", minClipSigma, maxClipSigma);
+
     //
     for (int N = 0; N < stats->clipIter; N++) {
         psTrace(__func__, 6, "Loop iteration %d.  Calling psVectorFitPolynomial1D()\n");
         int Nkeep = 0;
-        // XXX: Check error codes
         if (psTraceGetLevel(__func__) >= 6) {
             if (mask != NULL) {
@@ -902,4 +901,5 @@
         }
 
+        // XXX: Check error codes
         fit = psPolynomial1DEvalVector(poly, x);
         for (psS32 i = 0 ; i < f->n ; i++) {