Index: /branches/eam_rel8_b2/psLib/src/math/psMinimize.c
===================================================================
--- /branches/eam_rel8_b2/psLib/src/math/psMinimize.c	(revision 5595)
+++ /branches/eam_rel8_b2/psLib/src/math/psMinimize.c	(revision 5596)
@@ -10,6 +10,6 @@
  *  @author EAM, IfA
  *
- *  @version $Revision: 1.142.6.2 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-11-22 19:54:55 $
+ *  @version $Revision: 1.142.6.2.2.1 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-11-25 23:25:16 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -1872,4 +1872,6 @@
 
 // This function accepts F32 and F64 input vectors.
+// XXX EAM : note that the value used to mask new points must be
+//           part of the maskValue
 psPolynomial1D *psVectorClipFitPolynomial1D(
     psPolynomial1D *poly,
@@ -1913,12 +1915,12 @@
     float maxClipSigma;
     if (isfinite(stats->max)) {
+        maxClipSigma = fabs(stats->max);
+    } else {
         maxClipSigma = fabs(stats->clipSigma);
+    }
+    if (isfinite(stats->min)) {
+        minClipSigma = fabs(stats->min);
     } else {
-        maxClipSigma = fabs(stats->max);
-    }
-    if (isfinite(stats->min)) {
         minClipSigma = fabs(stats->clipSigma);
-    } else {
-        minClipSigma = fabs(stats->min);
     }
     psVector *fit   = NULL;
@@ -1941,4 +1943,7 @@
         float minClipValue = -minClipSigma*stats->sampleStdev;
         float maxClipValue = +maxClipSigma*stats->sampleStdev;
+
+        psTrace (".psphot.VectorClipFit", 4, "residual stats for robust fit:  %g +/- %g\n",
+                 stats->sampleMedian, stats->sampleStdev);
 
         // set mask if pts are not valid