Index: /trunk/psModules/src/pmImageSubtract.c
===================================================================
--- /trunk/psModules/src/pmImageSubtract.c	(revision 4423)
+++ /trunk/psModules/src/pmImageSubtract.c	(revision 4424)
@@ -7,6 +7,6 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.11 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-06-28 23:23:55 $
+ *  @version $Revision: 1.12 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-29 01:17:23 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -98,5 +98,5 @@
         int spatialOrder)
 {
-    psTrace("ImageSubtract.pmSubtractionKernelsAllocPOIS", 1,
+    psTrace("ImageSubtract.pmSubtractionKernelsAllocPOIS", 3,
             "Calling pmSubtractionKernelsAllocPOIS(%d, %d)\n", size, spatialOrder);
     PS_ASSERT_INT_POSITIVE(size, NULL);
@@ -170,5 +170,5 @@
     }
 
-    psTrace("ImageSubtract.pmSubtractionKernelsAllocPOIS", 1,
+    psTrace("ImageSubtract.pmSubtractionKernelsAllocPOIS", 3,
             "Exiting pmSubtractionKernelsAllocPOIS(%d, %d)\n", size, spatialOrder);
     return(tmp);
@@ -187,5 +187,5 @@
     PS_ASSERT_VECTOR_NON_NULL(sigmas, NULL);
     PS_ASSERT_VECTOR_NON_NULL(orders, NULL);
-    psTrace("ImageSubtract.pmSubtractionKernelsAllocISIS", 1,
+    psTrace("ImageSubtract.pmSubtractionKernelsAllocISIS", 3,
             "Calling pmSubtractionKernelsAllocISIS(%d, %d, %d, %d)\n",
             sigmas->n, orders->n, size, spatialOrder);
@@ -295,5 +295,5 @@
     }
 
-    psTrace("ImageSubtract.pmSubtractionKernelsAllocISIS", 1,
+    psTrace("ImageSubtract.pmSubtractionKernelsAllocISIS", 3,
             "Exiting pmSubtractionKernelsAllocISIS(%d, %d, %d, %d)\n",
             sigmas->n, orders->n, size, spatialOrder);
@@ -322,5 +322,5 @@
                                 )
 {
-    psTrace("ImageSubtract.pmSubtractionFindStamps", 1,
+    psTrace("ImageSubtract.pmSubtractionFindStamps", 3,
             "Calling pmSubtractionFindStamps(%d, %f, %d, %d, %d)\n",
             maskVal, threshold, xNum, yNum, border);
@@ -444,5 +444,5 @@
         }
     }
-    psTrace("ImageSubtract.pmSubtractionFindStamps", 1,
+    psTrace("ImageSubtract.pmSubtractionFindStamps", 3,
             "Exiting pmSubtractionFindStamps(%d, %f, %d, %d, %d)\n",
             maskVal, threshold, xNum, yNum, border);
@@ -462,5 +462,5 @@
                                 psF32 y)
 {
-    psTrace("ImageSubtract.GenSpatialOrder", 5,
+    psTrace("ImageSubtract.GenSpatialOrder", 4,
             "Calling GenSpatialOrder(%d, %f, %f)\n", spatialOrder, x, y);
 
@@ -478,5 +478,5 @@
     }
 
-    psTrace("ImageSubtract.GenSpatialOrder", 5,
+    psTrace("ImageSubtract.GenSpatialOrder", 4,
             "Exiting GenSpatialOrder(%d, %f, %f)\n", spatialOrder, x, y);
 
@@ -497,5 +497,5 @@
 
 
-    psTrace("ImageSubtract.IsisKernelConvolve", 5,
+    psTrace("ImageSubtract.IsisKernelConvolve", 4,
             "Calling IsisKernelConvolve(%d, %d, %d)\n", kernelID, col, row);
     psS32 spatialOrder = kernels->p_spatialOrder;
@@ -529,5 +529,5 @@
     psFree(polyValues);
 
-    psTrace("ImageSubtract.IsisKernelConvolve", 5,
+    psTrace("ImageSubtract.IsisKernelConvolve", 4,
             "Exiting IsisKernelConvolve(%d, %d, %d)\n", kernelID, col, row);
     return(conv);
@@ -551,5 +551,5 @@
                                psS32 row)
 {
-    psTrace("ImageSubtract.ConvolvePixelPois", 5,
+    psTrace("ImageSubtract.ConvolvePixelPois", 4,
             "Calling ConvolvePixelPois(%d, %d)\n", col, row);
     psS32 nBF = kernels->u->n;
@@ -592,5 +592,5 @@
     }
 
-    psTrace("ImageSubtract.ConvolvePixelPois", 5,
+    psTrace("ImageSubtract.ConvolvePixelPois", 4,
             "Exiting ConvolvePixelPois(%d, %d)\n", col, row);
     return(conv);
@@ -616,5 +616,5 @@
                                psS32 row)
 {
-    psTrace("ImageSubtract.ConvolvePixelIsis", 5,
+    psTrace("ImageSubtract.ConvolvePixelIsis", 4,
             "Calling ConvolvePixelIsis(%d, %d)\n", col, row);
     psF32 background = solution->data.F64[solution->n-1];
@@ -628,5 +628,5 @@
     }
 
-    psTrace("ImageSubtract.ConvolvePixelIsis", 5,
+    psTrace("ImageSubtract.ConvolvePixelIsis", 4,
             "Exiting ConvolvePixelIsis(%d, %d)\n", col, row);
     return(conv);
@@ -643,5 +643,5 @@
                               const psSubtractionKernels *kernels)
 {
-    psTrace("ImageSubtract.ConvolveImage", 5, "Calling ConvolveImage()\n");
+    psTrace("ImageSubtract.ConvolveImage", 4, "Calling ConvolveImage()\n");
     PS_ASSERT_IMAGE_NON_NULL(input, NULL);
     PS_ASSERT_IMAGE_NON_EMPTY(input, NULL);
@@ -712,5 +712,5 @@
     }
 
-    psTrace("ImageSubtract.ConvolveImage", 5, "Exiting ConvolveImage()\n");
+    psTrace("ImageSubtract.ConvolveImage", 4, "Exiting ConvolveImage()\n");
     return convolved;
 }
@@ -731,5 +731,5 @@
                                    )
 {
-    psTrace("ImageSubtract.pmSubtractionCalculateEquation", 1,
+    psTrace("ImageSubtract.pmSubtractionCalculateEquation", 3,
             "Calling pmSubtractionCalculateEquation()\n");
     PS_ASSERT_PTR_NON_NULL(stamps, false);
@@ -793,7 +793,7 @@
     for (psS32 s = 0; s < stamps->n; s++) {
         pmStamp *stamp = (pmStamp *) stamps->data[s];
-        psTrace("pmSubtractionCalculateEquation", 5, "subCalcEqn(): stamp %d\n", s);
+        psTrace("pmSubtractionCalculateEquation", 6, "subCalcEqn(): stamp %d\n", s);
         if (stamp->status == PM_STAMP_RECALC) {
-            psTrace("pmSubtractionCalculateEquation", 5, "subCalcEqn(): stamp %d: status is PM_STAMP_RECALC.\n", s);
+            psTrace("pmSubtractionCalculateEquation", 6, "subCalcEqn(): stamp %d: status is PM_STAMP_RECALC.\n", s);
             psImage *stampMatrix = stamp->matrix;
             psVector *stampVector = stamp->vector;
@@ -816,5 +816,5 @@
             }
             PS_VECTOR_SET_F64(stampVector, 0.0);
-            psTrace("pmSubtractionCalculateEquation", 5, "subCalcEqn(): stamp %d: allocate matrix and vector.\n", s);
+            psTrace("pmSubtractionCalculateEquation", 6, "subCalcEqn(): stamp %d: allocate matrix and vector.\n", s);
 
             //
@@ -827,5 +827,5 @@
                                                   ((psF64) (stamp->y - numHalfRows)) / ((psF64) numHalfRows));
 
-            psTrace("pmSubtractionCalculateEquation", 5, "subCalcEqn(): stamp %d: generated spatial order terms.\n", s);
+            psTrace("pmSubtractionCalculateEquation", 6, "subCalcEqn(): stamp %d: generated spatial order terms.\n", s);
 
             if (kernels->type == PM_SUBTRACTION_KERNEL_POIS) {
@@ -835,5 +835,5 @@
                 for (psS32 y = stamp->y - footprint; y < stamp->y + footprint; y++) {
                     for (psS32 x = stamp->x - footprint; x < stamp->x + footprint; x++) {
-                        psTrace("pmSubtractionCalculateEquation", 5, "subCalcEqn(): pixel (%d, %d).\n", y, x);
+                        psTrace("pmSubtractionCalculateEquation", 6, "subCalcEqn(): pixel (%d, %d).\n", y, x);
 
                         // The inverse of the noise, squared.
@@ -919,5 +919,5 @@
                 for (psS32 y = stamp->y - footprint; y < stamp->y + footprint; y++) {
                     for (psS32 x = stamp->x - footprint; x < stamp->x + footprint; x++) {
-                        psTrace("pmSubtractionCalculateEquation", 5, "subCalcEqn(): pixel (%d, %d).\n", y, x);
+                        psTrace("pmSubtractionCalculateEquation", 6, "subCalcEqn(): pixel (%d, %d).\n", y, x);
                         psF32 invNoise2 = 1.0/reference->data.F32[y][x]; // The inverse of the noise, squared.
 
@@ -985,5 +985,5 @@
         }
     }
-    psTrace("ImageSubtract.pmSubtractionCalculateEquation", 1,
+    psTrace("ImageSubtract.pmSubtractionCalculateEquation", 3,
             "Exiting pmSubtractionCalculateEquation()\n");
     return(true);
@@ -999,5 +999,5 @@
                                     )
 {
-    psTrace("ImageSubtract.pmSubtractionSolveEquation", 1,
+    psTrace("ImageSubtract.pmSubtractionSolveEquation", 3,
             "Calling pmSubtractionSolveEquation()\n");
     PS_ASSERT_PTR_NON_NULL(stamps, NULL);
@@ -1067,5 +1067,5 @@
 
     // XXX: psTrace()
-    if (1) {
+    if (0) {
         PS_IMAGE_PRINT_F64(sumMatrix);
     }
@@ -1081,5 +1081,5 @@
     }
     // XXX: psTrace()
-    if (1) {
+    if (0) {
         PS_IMAGE_PRINT_F64(luMatrix);
     }
@@ -1087,5 +1087,6 @@
     solution = psMatrixLUSolve(solution, luMatrix, sumVector, permutation);
     // XXX: psTrace()
-    if (1) {
+    // XXX: should we be checking for NAN's in the solution vector?
+    if (0) {
         PS_VECTOR_PRINT_F64(solution);
     }
@@ -1104,5 +1105,5 @@
     psFree(permutation);
 
-    psTrace("ImageSubtract.pmSubtractionSolveEquation", 1,
+    psTrace("ImageSubtract.pmSubtractionSolveEquation", 3,
             "Exiting pmSubtractionSolveEquation()\n");
     return(solution);
@@ -1122,5 +1123,5 @@
                                      const psVector *solution)
 {
-    psTrace("ImageSubtract.CalculateDeviations", 5,
+    psTrace("ImageSubtract.CalculateDeviations", 4,
             "Calling CalculateDeviations()\n");
     PS_ASSERT_PTR_NON_NULL(stamps, NULL);
@@ -1203,5 +1204,5 @@
     psFree(subStamp);
 
-    psTrace("ImageSubtract.CalculateDeviations", 5,
+    psTrace("ImageSubtract.CalculateDeviations", 4,
             "Exiting CalculateDeviations()\n");
     return deviations;
@@ -1221,5 +1222,5 @@
                               )
 {
-    psTrace("ImageSubtract.pmSubtractionRejectStamps", 5,
+    psTrace("ImageSubtract.pmSubtractionRejectStamps", 3,
             "Calling pmSubtractionRejectStamps()\n");
     PS_ASSERT_PTR_NON_NULL(stamps, false);
@@ -1296,5 +1297,5 @@
 
     psFree(deviations);
-    psTrace("ImageSubtract.pmSubtractionRejectStamps", 5,
+    psTrace("ImageSubtract.pmSubtractionRejectStamps", 3,
             "Exiting pmSubtractionRejectStamps()\n");
     return(true);
@@ -1310,5 +1311,5 @@
                                  )
 {
-    psTrace("ImageSubtract.pmSubtractionKernelImage", 5,
+    psTrace("ImageSubtract.pmSubtractionKernelImage", 3,
             "Calling pmSubtractionKernelImage()\n");
     PS_ASSERT_VECTOR_NON_NULL(solution, NULL);
@@ -1384,5 +1385,5 @@
     psFree(polyValues);
 
-    psTrace("ImageSubtract.pmSubtractionKernelImage", 5,
+    psTrace("ImageSubtract.pmSubtractionKernelImage", 3,
             "Exiting pmSubtractionKernelImage()\n");
     return(out);