Changeset 13379 for trunk/psModules/src/imcombine/pmSubtraction.c

Timestamp:

May 14, 2007, 3:23:19 PM (19 years ago)

Author:

Paul Price

Message:

Adding diagnostic output: FITS files with stamp matrices or kernels

File:

• trunk/psModules/src/imcombine/pmSubtraction.c (modified) (6 diffs)

trunk/psModules/src/imcombine/pmSubtraction.c

@@ -4 +4 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.2 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2007-05-12 02:00:39 $
+ *  @version $Revision: 1.3 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2007-05-15 01:23:19 $
  *
  *  Copyright 2004-2007 Institute for Astronomy, University of Hawaii
@@ -79 +79 @@
           for (int v = -size; v <= size; v++) {
               for (int u = -size; u <= size; u++) {
-                  sum += weightFunc(kernel->kernel[v][u]) * image->data.F64[y + v][x + u];
+                  sum += weightFunc(kernel->kernel[v][u]) * image->data.F32[y + v][x + u];
               }
           }
@@ -163 +163 @@
     for (int i = 0; i < stamps->n; i++) {
         pmSubtractionStamp *stamp = stamps->data[i]; // Stamp of interest
-        psTrace("psModules.imcombine", 6, "Inspecting stamp %d\n", i);
-
         if (stamp->status == PM_SUBTRACTION_STAMP_CALCULATE) {
-            psTrace("psModules.imcombine", 6, "Calculating for stamp %d\n", i);
             psImage *stampMatrix = stamp->matrix; // Least-squares matrix for this stamp
             psVector *stampVector = stamp->vector; // Least-squares vector for this stamp
@@ -240 +237 @@
 
             stamp->status = PM_SUBTRACTION_STAMP_USED;
+
+            if (psTraceGetLevel("psModules.imcombine.equation") >= 10) {
+                psString matrixName = NULL;
+                psStringAppend(&matrixName, "matrix%d.fits", i);
+                psFits *matrixFile = psFitsOpen(matrixName, "w");
+                psFree(matrixName);
+                psFitsWriteImage(matrixFile, NULL, stampMatrix, 0, NULL);
+                psFitsClose(matrixFile);
+            }
+
         }
     }
@@ -311 +318 @@
         psError(PS_ERR_UNKNOWN, true, "Failed to solve the least-squares system.\n");
         return NULL;
+    }
+
+    if (psTraceGetLevel("psModules.imcombine") >= 7) {
+        for (int i = 0; i < solution->n; i++) {
+            printf("       Solution %d: %f\n", i, solution->data.F64[i]);
+        }
     }
 
@@ -413 +426 @@
         if (stamp->status == PM_SUBTRACTION_STAMP_USED && deviations->data.F32[i] > limit) {
             // Mask out the stamp in the image so you it's not found again
+            psTrace("psModules.imcombine", 3, "Rejecting stamp %d\n", i);
             numRejected++;
             for (int y = stamp->y - footprint; y <= stamp->y + footprint; y++) {