Changeset 6028 for branches/eam_rel9_b1/psLib/src/math/psMinimize.c

Timestamp:

Jan 16, 2006, 11:26:50 PM (21 years ago)

Author:

magnier

Message:

fixed asserts in clip fit functions

File:

• branches/eam_rel9_b1/psLib/src/math/psMinimize.c (modified) (8 diffs)

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

@@ -10 +10 @@
  *  @author EAM, IfA
  *
- *  @version $Revision: 1.147.4.1.2.6 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2006-01-10 04:47:27 $
+ *  @version $Revision: 1.147.4.1.2.7 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-01-17 09:26:50 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -206 +206 @@
 }
 
+# include "psAbort.h"
+
 // XXX EAM: this needs to respect the mask on params
 // alpha, beta, params are already allocated
@@ -246 +248 @@
         delta = ymodel - y->data.F32[i];
         chisq += PS_SQR(delta) * dy->data.F32[i];
+        // fprintf (stderr, "chi: %f : %f %f %f\n", chisq, ymodel, y->data.F32[i], dy->data.F32[i]);
+        // if (!isfinite(chisq)) psAbort ("min", "chi");
 
         for (int j = 0; j < params->n; j++) {
@@ -2951 +2955 @@
     }
     PS_ASSERT_VECTOR_NON_NULL(x, NULL);
-    PS_ASSERT_VECTORS_SIZE_EQUAL(f, y, NULL);
+    PS_ASSERT_VECTORS_SIZE_EQUAL(f, x, NULL);
     PS_ASSERT_VECTOR_TYPE_F32_OR_F64(x, NULL);
     PS_ASSERT_VECTOR_NON_NULL(y, NULL);
     PS_ASSERT_VECTORS_SIZE_EQUAL(f, y, NULL);
     PS_ASSERT_VECTOR_TYPE_F32_OR_F64(y, NULL);
-    PS_ASSERT_VECTOR_NON_NULL(f, NULL);
-    // PS_ASSERT_VECTORS_SIZE_EQUAL(f, f, NULL);
-    // PS_ASSERT_VECTOR_TYPE_F32_OR_F64(f, NULL);
+    PS_ASSERT_VECTOR_NON_NULL(z, NULL);
+    PS_ASSERT_VECTORS_SIZE_EQUAL(f, z, NULL);
+    PS_ASSERT_VECTOR_TYPE_F32_OR_F64(z, NULL);
     if (fErr != NULL) {
-        PS_ASSERT_VECTORS_SIZE_EQUAL(fErr, mask, NULL);
+        PS_ASSERT_VECTORS_SIZE_EQUAL(f, fErr, NULL);
         PS_ASSERT_VECTOR_TYPE_F32_OR_F64(fErr, NULL);
     }
@@ -3211 +3215 @@
     psFree(B);
 
+    if (Sums == NULL) {
+        psLogMsg (".psLib.VectorFitPolynomial3DOrd", 4, "all data has been clipped away\n");
+        return (myPoly);
+    }
+
     for (int ix = 0; ix < nXterm; ix++) {
         for (int iy = 0; iy < nYterm; iy++) {
@@ -3438 +3447 @@
     PS_ASSERT_PTR_NON_NULL(stats, NULL);
     PS_ASSERT_VECTOR_NON_NULL(f, NULL);
-    PS_ASSERT_VECTOR_TYPE(f, PS_TYPE_F32, NULL);
+    PS_ASSERT_VECTOR_TYPE_F32_OR_F64(f, NULL);
     if (mask != NULL) {
         PS_ASSERT_VECTORS_SIZE_EQUAL(f, mask, NULL);
@@ -3444 +3453 @@
     }
     PS_ASSERT_VECTOR_NON_NULL(x, NULL);
-    PS_ASSERT_VECTORS_SIZE_EQUAL(f, y, NULL);
-    PS_ASSERT_VECTOR_TYPE(x, PS_TYPE_F32, NULL);
+    PS_ASSERT_VECTORS_SIZE_EQUAL(f, x, NULL);
+    PS_ASSERT_VECTOR_TYPE_F32_OR_F64(x, NULL);
     PS_ASSERT_VECTOR_NON_NULL(y, NULL);
     PS_ASSERT_VECTORS_SIZE_EQUAL(f, y, NULL);
-    PS_ASSERT_VECTOR_TYPE(y, PS_TYPE_F32, NULL);
-    PS_ASSERT_VECTOR_NON_NULL(f, NULL);
-    PS_ASSERT_VECTORS_SIZE_EQUAL(f, f, NULL);
-    PS_ASSERT_VECTOR_TYPE(f, PS_TYPE_F32, NULL);
+    PS_ASSERT_VECTOR_TYPE_F32_OR_F64(y, NULL);
+    PS_ASSERT_VECTOR_NON_NULL(z, NULL);
+    PS_ASSERT_VECTORS_SIZE_EQUAL(f, z, NULL);
+    PS_ASSERT_VECTOR_TYPE_F32_OR_F64(z, NULL);
     PS_ASSERT_VECTOR_NON_NULL(t, NULL);
     PS_ASSERT_VECTORS_SIZE_EQUAL(f, t, NULL);
-    PS_ASSERT_VECTOR_TYPE(t, PS_TYPE_F32, NULL);
+    PS_ASSERT_VECTOR_TYPE_F32_OR_F64(t, NULL);
     if (fErr != NULL) {
         PS_ASSERT_VECTORS_SIZE_EQUAL(fErr, mask, NULL);
-        PS_ASSERT_VECTOR_TYPE(fErr, PS_TYPE_F32, NULL);
+        PS_ASSERT_VECTOR_TYPE_F32_OR_F64(fErr, NULL);
     }
 
@@ -3520 +3529 @@
                  Nkeep, x->n);
 
+        stats->clippedNvalues = Nkeep;
         psFree (fit);
     }