Index: trunk/psLib/src/math/psStats.c
===================================================================
--- trunk/psLib/src/math/psStats.c	(revision 4962)
+++ trunk/psLib/src/math/psStats.c	(revision 4991)
@@ -17,6 +17,6 @@
  *
  *
- *  @version $Revision: 1.144 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-09-07 21:40:09 $
+ *  @version $Revision: 1.145 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-09-11 22:18:40 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -746,4 +746,7 @@
 XXX: Write a general routine which smoothes a psVector.  This routine should
 call that.  Is that possible?
+ 
+XXX: This is, or will be, prettier than the previous
+psVectorSmoothHistGaussian().  However, it is not being used, yet.
  *****************************************************************************/
 psVector *p_psVectorSmoothHistGaussianNEW(psHistogram *histogram,
@@ -1549,5 +1552,6 @@
 
         // Determine the coefficients of the polynomial.
-        myPoly = psVectorFitPolynomial1D(myPoly, x, y, yErr);
+        //        myPoly = psVectorFitPolynomial1D(myPoly, x, y, yErr);
+        myPoly = psVectorFitPolynomial1D(myPoly, NULL, 0, y, yErr, x);
         if (myPoly == NULL) {
             psError(PS_ERR_UNEXPECTED_NULL,
@@ -1825,5 +1829,5 @@
         psPolynomial1D *tmpPoly = psPolynomial1DAlloc(3, PS_POLYNOMIAL_ORD);
         // XXX: What about the NULL x argument?
-        tmpPoly = psVectorFitPolynomial1D(tmpPoly, NULL, y, NULL);
+        tmpPoly = psVectorFitPolynomial1D(tmpPoly, NULL, 0, y, NULL, NULL);
         if (tmpPoly == NULL) {
             psLogMsg(__func__, PS_LOG_WARN,
@@ -1936,5 +1940,5 @@
         }
     } else {
-        printf("XXX: Generate an error here.\n");
+        psError(PS_ERR_UNKNOWN, true, "Unallowable vector type.\n");
         return(NULL);
     }
@@ -1950,6 +1954,8 @@
 psVector *Fit1DGaussian(psVector *x, psVector*y)
 {
-    printf("XXX: Generate an error here.\n");
-    printf("XXX: Error: This function was previously part of psStats.c, was removed, was purged from CVS, and now needs to be retrieved from tape.\n");
+    psError(PS_ERR_UNKNOWN, true, "This code has not been implemented yet.\n");
+    // XXX: This function was previously part of psStats.c, was removed, was
+    // purged from CVS, and now needs to be retrieved from tape.
+
     return(NULL);
 }
@@ -2168,8 +2174,6 @@
             stats->robustLQ = binLo25F32;
             stats->robustUQ = binHi25F32;
-            // XXX: No idea how to calculate stats->stdev
 
             // PS_BIN_MIDPOINT(robustHistogram, modeBinNum);
-
             // XXX: I think sumNfit == sumN50 here.
             stats->robustNfit = -1;