Changeset 5570 for branches/eam_rel8_b1/psLib/src/math/psMinimize.c

Timestamp:

Nov 22, 2005, 9:54:55 AM (21 years ago)

Author:

magnier

Message:

fixed compilation bugs

File:

• branches/eam_rel8_b1/psLib/src/math/psMinimize.c (modified) (10 diffs)

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

@@ -10 +10 @@
  *  @author EAM, IfA
  *
- *  @version $Revision: 1.142.6.1 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-11-22 18:06:10 $
+ *  @version $Revision: 1.142.6.2 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-11-22 19:54:55 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -1441 +1441 @@
 }
 
-
-/******************************************************************************
-BuildSums4D(sums, x, y, z, t, nXterm, nYterm, nZterm, nTterm). equiv to 
-BuildSums2D(). The result is returned as a double ****
- *****************************************************************************/
-static psImage *BuildSums4D(
-    psF64 ****sums,
-    psF64 x,
-    psF64 y,
-    psF64 z,
-    psF64 t,
-    psS32 nXterm,
-    psS32 nYterm,
-    psS32 nZterm,
-    psS32 nTterm)
+// here is the definition for BuildSums4D.  ifdef'ed away until it is used
+// by psPolynomial4DFit..
+# if (0)
+    /******************************************************************************
+    BuildSums4D(sums, x, y, z, t, nXterm, nYterm, nZterm, nTterm). equiv to 
+    BuildSums2D(). The result is returned as a double ****
+     *****************************************************************************/
+    static double ****BuildSums4D(
+        psF64 ****sums,
+        psF64 x,
+        psF64 y,
+        psF64 z,
+        psF64 t,
+        psS32 nXterm,
+        psS32 nYterm,
+        psS32 nZterm,
+        psS32 nTterm)
 {
     psS32 nXsum = 0;
@@ -1471 +1473 @@
     nTsum = 2*nTterm;
     if (sums == NULL) {
-        sums = (psF64 ****) psAlloc (nXsum*PSELEMTYPE_SIZEOF(psF64));
+        sums = (psF64 ****) psAlloc (nXsum*sizeof(psF64));
         for (int i = 0; i < nXsum; i++) {
-            sums[i] = (psF64 ***) psAlloc (nYsum*PSELEMTYPE_SIZEOF(psF64));
+            sums[i] = (psF64 ***) psAlloc (nYsum*sizeof(psF64));
             for (int j = 0; j < nYsum; j++) {
-                sums[i][j] = (psF64 **) psAlloc (nZsum*PSELEMTYPE_SIZEOF(psF64));
+                sums[i][j] = (psF64 **) psAlloc (nZsum*sizeof(psF64));
                 for (int k = 0; k < nZsum; k++) {
-                    sums[i][j][k] = (psF64 *) psAlloc (nTsum*PSELEMTYPE_SIZEOF(psF64));
+                    sums[i][j][k] = (psF64 *) psAlloc (nTsum*sizeof(psF64));
                 }
             }
@@ -1503 +1505 @@
     return (sums);
 }
+# endif /* BuildSums4D */
 
 /******************************************************************************
@@ -1509 +1512 @@
 in the output vector.  This routine works on single-precision polynomials with
 double precision data.
+ 
+XXX EAM : this function is now deprecated: psPolynomial2DEvalVector handles F32 and F64
  *****************************************************************************/
 psVector *Polynomial2DEvalVectorD(
@@ -1899 +1904 @@
     // assign sequence vector if xIn is NULL
     if (xIn == NULL) {
-        x = psVectorCreate (NULL, 0, f->n, 1, f->type);
+        x = psVectorCreate (NULL, 0, f->n, 1, f->type.type);
     } else {
-        x = xIn;
+        x = (psVector *) xIn;
     }
 
@@ -1907 +1912 @@
     float minClipSigma;
     float maxClipSigma;
-    if (finite(stats->max)) {
+    if (isfinite(stats->max)) {
         maxClipSigma = fabs(stats->clipSigma);
     } else {
         maxClipSigma = fabs(stats->max);
     }
-    if (finite(stats->min)) {
+    if (isfinite(stats->min)) {
         minClipSigma = fabs(stats->clipSigma);
     } else {
@@ -2055 +2060 @@
         } else {
             // this filters fErr == 0 values
-            wt = (fErr->data.F64[k] == ) ? 0.0 : 1.0 / PS_SQR(fErr->data.F64[k]);
+            wt = (fErr->data.F64[k] == 0.0) ? 0.0 : 1.0 / PS_SQR(fErr->data.F64[k]);
         }
 
@@ -2097 +2102 @@
 
 
+// XXX EAM : I have implemented a single function to handle the mask/nomask cases
+//           this function can be deprecated
 psPolynomial2D* RobustFit2D_nomask(
     psPolynomial2D* poly,
@@ -2382 +2389 @@
     float minClipSigma;
     float maxClipSigma;
-    if (finite(stats->max)) {
+    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 (finite(stats->min)) {
         minClipSigma = fabs(stats->clipSigma);
-    } else {
-        minClipSigma = fabs(stats->min);
     }
     psVector *fit   = NULL;