Changeset 7856 for trunk/psModules/src/detrend/pmFlatField.c

Timestamp:

Jul 10, 2006, 10:45:44 AM (20 years ago)

Author:

Paul Price

Message:

Fixing "special" value for integer types.

File:

• trunk/psModules/src/detrend/pmFlatField.c (modified) (4 diffs)

trunk/psModules/src/detrend/pmFlatField.c

@@ -18 +18 @@
  *  @author Ross Harman, MHPCC
  *
- *  @version $Revision: 1.10 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2006-07-10 20:42:47 $
+ *  @version $Revision: 1.11 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-07-10 20:45:44 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -90 +90 @@
 
     // Macro for all PS types
-    #define FLAT_DIVISION_CASE(TYPE) \
+    #define FLAT_DIVISION_CASE(TYPE, SPECIAL) \
 case PS_TYPE_##TYPE: \
     for (int j = 0; j < inImage->numRows; j++) { \
@@ -98 +98 @@
                     inMask->data.U8[j][i] |= PM_MASK_FLAT; \
                 } \
-                inImage->data.TYPE[j][i] = NAN; \
+                inImage->data.TYPE[j][i] = SPECIAL; \
             } else { \
                 inImage->data.TYPE[j][i] /= flatImage->data.TYPE[j][i]; \
@@ -107 +107 @@
 
     switch (inImage->type.type) {
-        FLAT_DIVISION_CASE(U8);
-        FLAT_DIVISION_CASE(U16);
-        FLAT_DIVISION_CASE(U32);
-        FLAT_DIVISION_CASE(U64);
-        FLAT_DIVISION_CASE(S8);
-        FLAT_DIVISION_CASE(S16);
-        FLAT_DIVISION_CASE(S32);
-        FLAT_DIVISION_CASE(S64);
-        FLAT_DIVISION_CASE(F32);
-        FLAT_DIVISION_CASE(F64);
+        FLAT_DIVISION_CASE(U8,  0);
+        FLAT_DIVISION_CASE(U16, 0);
+        FLAT_DIVISION_CASE(U32, 0);
+        FLAT_DIVISION_CASE(U64, 0);
+        FLAT_DIVISION_CASE(S8   0);
+        FLAT_DIVISION_CASE(S16, 0);
+        FLAT_DIVISION_CASE(S32, 0);
+        FLAT_DIVISION_CASE(S64, 0);
+        FLAT_DIVISION_CASE(F32, NAN);
+        FLAT_DIVISION_CASE(F64, NAN);
     default:
         psError(PS_ERR_BAD_PARAMETER_TYPE, true,