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

Timestamp:

May 17, 2006, 5:10:06 PM (20 years ago)

Author:

Paul Price

Message:

Adding check for NAN data.

File:

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

trunk/psLib/src/math/psStats.c

@@ -16 +16 @@
  * use ->min and ->max (PS_STAT_USE_RANGE)
  *
- *  @version $Revision: 1.172 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2006-05-18 01:20:50 $
+ *  @version $Revision: 1.173 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-05-18 03:10:06 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -2041 +2041 @@
         // We want a copy of the mask, since we may change values
         maskU8 = psVectorCopy(NULL, mask, PS_TYPE_U8);
-    }
-
+    } else {
+        maskU8 = psVectorAlloc(in->n, PS_TYPE_U8);
+        maskU8->n = in->n;
+        memset(maskU8->data.U8, 0, in->n * sizeof(psU8));
+    }
+
+    // Mask bad data
+    for (long i = 0; i < inF32->n; i++) {
+        if (!isfinite(inF32->data.F32[i])) {
+            maskU8->data.U8[i] = 0xff;
+        }
+    }
 
     if ((stats->options & PS_STAT_USE_RANGE) && (stats->min >= stats->max)) {