Index: trunk/psModules/src/tst_pmFlatField.c
===================================================================
--- trunk/psModules/src/tst_pmFlatField.c	(revision 1802)
+++ trunk/psModules/src/tst_pmFlatField.c	(revision 1832)
@@ -5,8 +5,8 @@
  *  Detailed Description goes here.
  *
- *  @author Robert DeSonia, MHPCC
+ *  @author Ross Harman, MHPCC
  *
- *  @version $Revision: 1.2 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2004-09-14 02:29:40 $
+ *  @version $Revision: 1.3 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2004-09-20 20:36:46 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -15,10 +15,47 @@
 #include "pmFlatField.h"
 
+
+#define PRINT_MATRIX(IMAGE,TYPE)                                                                             \
+for(int i=IMAGE->numRows-1; i>-1; i--) {                                                                 \
+    for(int j=0; j<IMAGE->numCols; j++) {                                                                \
+        if(PS_IS_PSELEMTYPE_COMPLEX(IMAGE->type.type)) {                                                 \
+            printf("%f+%fi ", creal(IMAGE->data.TYPE[i][j]), cimag(IMAGE->data.TYPE[i][j]));             \
+        } else if(PS_IS_PSELEMTYPE_INT(IMAGE->type.type)) {                                              \
+            printf("%d ", (int)IMAGE->data.TYPE[i][j]);                                                  \
+        } else {                                                                                         \
+            printf("%f ", (double)IMAGE->data.TYPE[i][j]);                                               \
+        }                                                                                                \
+    }                                                                                                    \
+    printf("\n");                                                                                        \
+}                                                                                                        \
+printf("\n");
+
+
 int main(void)
 {
-    psVector* vec = psVectorAlloc(10,PS_TYPE_F32);
-    for (int i=0;i<10;i++) {
-        vec->data.F32[i] = i/9.0f;
+    psImage *inImage = psImageAlloc(3, 3, PS_TYPE_F64);
+    psImage *flatImage1 = psImageAlloc(3, 3, PS_TYPE_F64);
+    psImage *flatImage2 = psImageAlloc(3, 3, PS_TYPE_F64);
+    psReadout *inReadout  = psReadoutAlloc(0, 0, inImage);
+    psReadout *flatReadout  = psReadoutAlloc(0, 0, flatImage1);
+
+    for(int r=0;r<3;r++) {
+        for(int c=0;c<3;c++) {
+            inImage->data.F64[r][c]=8.0;
+            flatImage1->data.F64[r][c]=4.0;
+            flatImage2->data.F64[r][c]=0.0;
+        }
     }
-    //    pmFlatField(vec);
+
+    PRINT_MATRIX(inImage,F64);
+    PRINT_MATRIX(flatImage1, F64);
+    pmFlatField(inReadout, flatReadout);
+    PRINT_MATRIX(inImage,F64);
+
+    PRINT_MATRIX(inImage, F64);
+    PRINT_MATRIX(flatImage2, F64);
+    flatReadout->image = flatImage2;
+    pmFlatField(inReadout, flatReadout);
+    PRINT_MATRIX(inReadout->mask, PS_TYPE_MASK_DATA);
+    return 0;
 }