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

Timestamp:

Feb 12, 2007, 5:08:17 PM (19 years ago)

Author:

Paul Price

Message:

Fixing compilation problems following change to histogram API.

File:

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

trunk/psLib/src/math/psStats.c

@@ -13 +13 @@
  * use ->min and ->max (PS_STAT_USE_RANGE)
  *
- *  @version $Revision: 1.200 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2007-02-07 23:52:54 $
+ *  @version $Revision: 1.201 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2007-02-13 03:08:17 $
  *
  *  Copyright 2006 IfA, University of Hawaii
@@ -87 +87 @@
 /******************************************************************************
 MISC PRIVATE STATISTICAL FUNCTIONS
- 
+
 NOTE: it is assumed that any call to these statistical functions will have
 been preceded by a call to the psVectorStats() function.  Various sanity tests
 will only be performed in psVectorStats().
- 
+
 For many of these private stats routines, it is possible that there are no acceptable elements
 in the input vector (if no elements lie within range, or there are no unmasked elements, or the
@@ -106 +106 @@
 mean of the input vector.  If there was a problem with the mean calculation,
 this routine sets stats->sampleMean to NAN.
- 
+
 using the method below, with a single loop for various options costs only a small amount and is
 much easier to debug. running 10000 tests of 1000 point vectors for the two methods gives:
- 
+
                      separate    single          w/errors
 (mask: 0, range: 0): 0.067 sec  0.073 sec (10%)  0.072 sec
@@ -115 +115 @@
 (mask: 0, range: 0): 0.136 sec  0.162 sec (20%)  0.170 sec
 (mask: 1, range: 0): 0.177 sec  0.181 sec ( 2%)  0.198 sec
- 
+
 (effect of errors not tested)
- 
+
 To optmize this, use a macro and ifdef in or out the three states (errors, mask, range)
 *****************************************************************************/
@@ -320 +320 @@
 Returns
     NULL
- 
+
 using the method below, with a single loop for various options costs only a small amount and is
 much easier to debug. running 10000 tests of 1000 point vectors for the two methods gives:
- 
+
                      single
 (mask: 0, range: 0): 0.193 sec
@@ -329 +329 @@
 (mask: 0, range: 0): 0.349 sec
 (mask: 1, range: 0): 0.401 sec
- 
+
 *****************************************************************************/
 static bool vectorSampleStdev(const psVector* myVector,
@@ -408 +408 @@
 vectorClippedStats(myVector, errors, maskInput, maskValInput, stats): calculates the
 clipped stats (mean or stdev) of the input vector.
- 
+
 Inputs
     myVector
@@ -557 +557 @@
 vectorRobustStats(myVector, maskInput, maskValInput, stats): This is the new
 version of the robust stats routine.
- 
+
 XXX: MUST DO: If the errors in the input values are known, then the same
 approach is used, except that the histograms become probability density
@@ -567 +567 @@
 each input value contributes constant area. Then the robust median and
 standard deviation are estimated in the same manner as above.
- 
+
 XXX: Check for errors in psLib routines that we call.
- 
+
 XXX: Review and ensure that all memory is free'ed at premature exits.
 *****************************************************************************/
@@ -666 +666 @@
         histogram = psHistogramAlloc(min, max, numBins);
         // XXXXX we need to consider this step if errors -> variance
-        histogram = psVectorHistogram(histogram, myVector, errors, mask, maskVal);
+        if (!psVectorHistogram(histogram, myVector, errors, mask, maskVal)) {
+            psError(PS_ERR_UNKNOWN, false, "Unable to generate histogram for robust statistics.\n");
+            psFree(histogram);
+            psFree(statsMinMax);
+            return false;
+        }
         if (psTraceGetLevel("psLib.math") >= 8) {
             PS_VECTOR_PRINT_F32(histogram->bounds);
@@ -980 +985 @@
 
         psHistogram *histogram = psHistogramAlloc(min, max, numBins); // A new histogram (without outliers)
-        histogram = psVectorHistogram(histogram, myVector, errors, mask, maskVal);
+        if (!psVectorHistogram(histogram, myVector, errors, mask, maskVal)) {
+            psError(PS_ERR_UNKNOWN, false, "Unable to generate histogram for fitted statistics.\n");
+            psFree(histogram);
+            psFree(statsMinMax);
+            return false;
+        }
         if (psTraceGetLevel("psLib.math") >= 8) {
             PS_VECTOR_PRINT_F32(histogram->nums);
@@ -1174 +1184 @@
 
         psHistogram *histogram = psHistogramAlloc(min, max, numBins); // A new histogram (without outliers)
-        histogram = psVectorHistogram(histogram, myVector, errors, mask, maskVal);
+        if (!psVectorHistogram(histogram, myVector, errors, mask, maskVal)) {
+            psError(PS_ERR_UNKNOWN, false, "Unable to generate histogram for fitted statistcs v2.\n");
+            psFree(histogram);
+            psFree(statsMinMax);
+            return false;
+        }
         if (psTraceGetLevel("psLib.math") >= 8) {
             PS_VECTOR_PRINT_F32(histogram->nums);
@@ -1382 +1397 @@
 
         psHistogram *histogram = psHistogramAlloc(min, max, numBins); // A new histogram (without outliers)
-        histogram = psVectorHistogram(histogram, myVector, errors, mask, maskVal);
+        if (!psVectorHistogram(histogram, myVector, errors, mask, maskVal)) {
+            psError(PS_ERR_UNKNOWN, false, "Unable to generate histogram for fitted statistics v3.\n");
+            psFree(histogram);
+            psFree(statsMinMax);
+            return false;
+        }
         if (psTraceGetLevel("psLib.math") >= 8) {
             PS_VECTOR_PRINT_F32(histogram->nums);
@@ -1594 +1614 @@
 robustHistogram with a Gaussian of width sigma.  It returns a psVector of the
 smoothed data.
- 
+
 XXX this function is unused
 *****************************************************************************/
@@ -1753 +1773 @@
 function which calls the above private stats functions based on what bits
 were set in stats->options.
- 
+
 Inputs
     in
@@ -1761 +1781 @@
 Returns
     The stats structure.
- 
+
 XXX: Should we free stats if the asserts fail? NO; we don't own it (RHL).
 *****************************************************************************/
@@ -2107 +2127 @@
 then determines for what value x does that quadratic f(x) = yVal (the input
 parameter).
- 
+
 *****************************************************************************/
 static psF32 fitQuadraticSearchForYThenReturnX(const psVector *xVec,