Changeset 6185 for trunk/psLib/src/math/psMinimizePolyFit.c

Timestamp:

Jan 23, 2006, 10:44:29 AM (20 years ago)

Author:

gusciora

Message:

psMinimize is now split into 3 files.

File:

• trunk/psLib/src/math/psMinimizePolyFit.c (modified) (5 diffs)

trunk/psLib/src/math/psMinimizePolyFit.c

@@ -10 +10 @@
  *  @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 +862 @@
     psF32 minClipSigma;
     psF32 maxClipSigma;
-
     if (isnan(stats->max) || isfinite(stats->max)) {
         maxClipSigma = fabs(stats->clipSigma);
@@ -873 +872 @@
         minClipSigma = fabs(stats->min);
     }
+
     psVector *fit   = NULL;
     psVector *resid = psVectorAlloc(f->n, PS_TYPE_F64);
@@ -881 +881 @@
     // 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 +901 @@
         }
 
+        // XXX: Check error codes
         fit = psPolynomial1DEvalVector(poly, x);
         for (psS32 i = 0 ; i < f->n ; i++) {