Changeset 2221 for trunk/psLib/src/math/psPolynomial.c

Timestamp:

Oct 27, 2004, 1:31:44 PM (22 years ago)

Author:

gusciora

Message:

Parameter checking for the psCOord.c file.

File:

• trunk/psLib/src/math/psPolynomial.c (modified) (14 diffs)

trunk/psLib/src/math/psPolynomial.c

@@ -7 +7 @@
  *  polynomials.  It also contains a Gaussian functions.
  *
- *  @version $Revision: 1.56 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2004-10-27 21:25:52 $
+ *  @version $Revision: 1.57 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2004-10-27 23:31:43 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -32 +32 @@
 #include "psTrace.h"
 #include "psError.h"
-#include "psAbort.h"
 #include "psLogMsg.h"
 #include "psFunctions.h"
@@ -145 +144 @@
  Gaussian distribution as well.
  
- NOTE: XXX: There is no way to seed the random generator.
+XXX: There is no way to seed the random generator.
  *****************************************************************************/
 psVector* psGaussianDev(float mean, float sigma, psS32 Npts)
@@ -167 +166 @@
     }
 
-    // NOTE: Should I free r as well?
+    // XXX: Should I free r, T as well?  This is a memory leak.
     return(gauss);
 }
@@ -426 +425 @@
     }
 
-
-    if (NULL == myPoly) {
-        psAbort(__func__, "psPolynomial1DEval(): myPoly is NULL\n");
-    }
-
-    // NOTE: Do we want to flag this case?
-    if (myPoly->n == 0) {
-        return(1.0);
-    }
-
-    if (NULL == myPoly->coeff) {
-        psAbort(__func__, "psPolynomial1DEval(): myPoly->coeff is NULL\n");
-    }
-
     for (loop_x = 0; loop_x < myPoly->n; loop_x++) {
         if (myPoly->mask[loop_x] == 0) {
@@ -508 +493 @@
         return(p_psChebPolynomial1DEval(x, myPoly));
     } else {
-        psAbort(__func__, "Unknown polynomial type 0x%x\n", myPoly->type);
+        psError(__func__, "Unknown polynomial type 0x%x\n", myPoly->type);
     }
     return(0.0);
@@ -604 +589 @@
         return(p_psChebPolynomial2DEval(x, y, myPoly));
     } else {
-        psAbort(__func__, "Unknown polynomial type 0x%x\n", myPoly->type);
+        psError(__func__, "Unknown polynomial type 0x%x\n", myPoly->type);
     }
     return(0.0);
@@ -727 +712 @@
         return(p_psChebPolynomial3DEval(x, y, z, myPoly));
     } else {
-        psAbort(__func__, "Unknown polynomial type 0x%x\n", myPoly->type);
+        psError(__func__, "Unknown polynomial type 0x%x\n", myPoly->type);
     }
     return(0.0);
@@ -876 +861 @@
         return(p_psChebPolynomial4DEval(w,x,y,z, myPoly));
     } else {
-        psAbort(__func__, "Unknown polynomial type 0x%x\n", myPoly->type);
+        psError(__func__, "Unknown polynomial type 0x%x\n", myPoly->type);
     }
     return(0.0);
@@ -1185 +1170 @@
     double xSum = 1.0;
 
-    // NOTE: Do we want to flag this case?
-    if (myPoly->n == 0) {
-        return(1.0);
-    }
-
     for (loop_x = 0; loop_x < myPoly->n; loop_x++) {
         if (myPoly->mask[loop_x] == 0) {
@@ -1236 +1216 @@
         return(p_psDChebPolynomial1DEval(x, myPoly));
     } else {
-        psAbort(__func__, "Unknown polynomial type 0x%x\n", myPoly->type);
+        psError(__func__, "Unknown polynomial type 0x%x\n", myPoly->type);
     }
     return(0.0);
@@ -1330 +1310 @@
         return(p_psDChebPolynomial2DEval(x, y, myPoly));
     } else {
-        psAbort(__func__, "Unknown polynomial type 0x%x\n", myPoly->type);
+        psError(__func__, "Unknown polynomial type 0x%x\n", myPoly->type);
     }
     return(0.0);
@@ -1452 +1432 @@
         return(p_psDChebPolynomial3DEval(x, y, z, myPoly));
     } else {
-        psAbort(__func__, "Unknown polynomial type 0x%x\n", myPoly->type);
+        psError(__func__, "Unknown polynomial type 0x%x\n", myPoly->type);
     }
     return(0.0);
@@ -1603 +1583 @@
         return(p_psDChebPolynomial4DEval(w,x,y,z, myPoly));
     } else {
-        psAbort(__func__, "Unknown polynomial type 0x%x\n", myPoly->type);
+        psError(__func__, "Unknown polynomial type 0x%x\n", myPoly->type);
     }
     return(0.0);