Index: trunk/psLib/src/math/psSpline.c
===================================================================
--- trunk/psLib/src/math/psSpline.c	(revision 2218)
+++ trunk/psLib/src/math/psSpline.c	(revision 2221)
@@ -7,6 +7,6 @@
  *  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,5 +32,4 @@
 #include "psTrace.h"
 #include "psError.h"
-#include "psAbort.h"
 #include "psLogMsg.h"
 #include "psFunctions.h"
@@ -145,5 +144,5 @@
  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,5 +166,5 @@
     }
 
-    // NOTE: Should I free r as well?
+    // XXX: Should I free r, T as well?  This is a memory leak.
     return(gauss);
 }
@@ -426,18 +425,4 @@
     }
 
-
-    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,5 +493,5 @@
         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,5 +589,5 @@
         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,5 +712,5 @@
         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,5 +861,5 @@
         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,9 +1170,4 @@
     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,5 +1216,5 @@
         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,5 +1310,5 @@
         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,5 +1432,5 @@
         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,5 +1583,5 @@
         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);