Changeset 5596

Timestamp:

Nov 25, 2005, 1:25:16 PM (21 years ago)

Author:

magnier

Message:

fixed bug in psVectorClipFitPolynomial1D test for stats->min,max

File:

• branches/eam_rel8_b2/psLib/src/math/psMinimize.c (modified) (4 diffs)

branches/eam_rel8_b2/psLib/src/math/psMinimize.c

@@ -10 +10 @@
  *  @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 +1872 @@
 
 // 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 +1915 @@
     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 +1943 @@
         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