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

Timestamp:

May 23, 2007, 4:06:51 PM (19 years ago)

Author:

Paul Price

Message:

Make psStatsAlloc a macro, to deliver filename and line number information to the memory system, as most other allocators do. This is an ABI-incompatible change.

File:

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

trunk/psLib/src/math/psStats.c

@@ -13 +13 @@
  * use ->min and ->max (PS_STAT_USE_RANGE)
  *
- *  @version $Revision: 1.210 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2007-05-18 13:19:55 $
+ *  @version $Revision: 1.211 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2007-05-24 02:06:51 $
  *
  *  Copyright 2006 IfA, University of Hawaii
@@ -259 +259 @@
     for (long i = 0; i < num; i++) {
         // Check if the data is with the specified range
-        if (!isfinite(vector[i])) {
-            stats->max = stats->min = NAN;
-            psError(PS_ERR_IEEE, true, "Element %ld of vector is Inf/NaN", i);
-            return 0;
-        }
-        
+        if (!isfinite(vector[i])) {
+            stats->max = stats->min = NAN;
+            psError(PS_ERR_IEEE, true, "Element %ld of vector is Inf/NaN", i);
+            return 0;
+        }
+
         if (useRange && (vector[i] < stats->min))
             continue;
@@ -688 +688 @@
         }
         psTrace(TRACE, 6, "Data min/max is (%.2f, %.2f)\n", min, max);
-        assert (iterate < 20);
+        assert (iterate < 20);
         // fprintf (stderr, "%d @ %d : %f - %f\n", iterate, numValid, min, max);
 
@@ -1661 +1661 @@
 
 // We keep statsFree so that we can identify statistics pointers from the memblock
-static void statsFree(psStats *newStruct)
+static void statsFree(psStats *stats)
 {
     // There are non dynamic allocated items
+    return;
 }
 
@@ -1669 +1670 @@
     psStatsAlloc(): This routine must create a new psStats data structure.
 *****************************************************************************/
-psStats* psStatsAlloc(psStatsOptions options)
+psStats* p_psStatsAlloc(const char *file, unsigned int lineno, const char *func, psStatsOptions options)
 {
     psTrace(TRACE, 3,"---- %s() begin  ----\n", __func__);
-    psStats* newStruct = NULL;
-
-    newStruct = (psStats* ) psAlloc(sizeof(psStats));
-    psMemSetDeallocator(newStruct, (psFreeFunc)statsFree);
-    newStruct->sampleMean = NAN;
-    newStruct->sampleMedian = NAN;
-    newStruct->sampleStdev = NAN;
-    newStruct->sampleUQ = NAN;
-    newStruct->sampleLQ = NAN;
-    newStruct->robustMedian = NAN;
-    newStruct->robustStdev = NAN;
-    newStruct->robustUQ = NAN;
-    newStruct->robustLQ = NAN;
-    newStruct->robustN50 = -1;
-    newStruct->fittedMean = NAN;
-    newStruct->fittedStdev = NAN;
-    newStruct->fittedNfit = -1;
-    newStruct->clippedMean = NAN;
-    newStruct->clippedStdev = NAN;
-    newStruct->clippedNvalues = -1;     // XXX: This is never used
-    newStruct->clipSigma = 3.0;
-    newStruct->clipIter = 3;
-    newStruct->min = NAN;
-    newStruct->max = NAN;
-    newStruct->binsize = NAN;
-    newStruct->nSubsample = 100000;
-    newStruct->options = options;
+
+    psStats *stats = p_psAlloc(file, lineno, func, sizeof(psStats));
+    psMemSetDeallocator(stats, (psFreeFunc)statsFree);
+    stats->sampleMean = NAN;
+    stats->sampleMedian = NAN;
+    stats->sampleStdev = NAN;
+    stats->sampleUQ = NAN;
+    stats->sampleLQ = NAN;
+    stats->robustMedian = NAN;
+    stats->robustStdev = NAN;
+    stats->robustUQ = NAN;
+    stats->robustLQ = NAN;
+    stats->robustN50 = -1;
+    stats->fittedMean = NAN;
+    stats->fittedStdev = NAN;
+    stats->fittedNfit = -1;
+    stats->clippedMean = NAN;
+    stats->clippedStdev = NAN;
+    stats->clippedNvalues = -1;     // XXX: This is never used
+    stats->clipSigma = 3.0;
+    stats->clipIter = 3;
+    stats->min = NAN;
+    stats->max = NAN;
+    stats->binsize = NAN;
+    stats->nSubsample = 100000;
+    stats->options = options;
 
     psTrace(TRACE, 3, "---- %s() end  ----\n", __func__);
-    return (newStruct);
+    return stats;
 }