Changeset 4991 for trunk/psLib/src/math/psStats.c

Timestamp:

Sep 11, 2005, 12:18:40 PM (21 years ago)

Author:

gusciora

Message:

....

File:

• trunk/psLib/src/math/psStats.c (modified) (7 diffs)

trunk/psLib/src/math/psStats.c

@@ -17 +17 @@
  *
  *
- *  @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 +746 @@
 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 +1552 @@
 
         // 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 +1829 @@
         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 +1940 @@
         }
     } else {
-        printf("XXX: Generate an error here.\n");
+        psError(PS_ERR_UNKNOWN, true, "Unallowable vector type.\n");
         return(NULL);
     }
@@ -1950 +1954 @@
 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 +2174 @@
             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;