Changeset 13379 for trunk/psModules/src/imcombine/pmSubtractionKernels.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/pmSubtractionKernels.c (modified) (3 diffs)

trunk/psModules/src/imcombine/pmSubtractionKernels.c

@@ -29 +29 @@
     }
     float result = value;               // Result to return
-    for (int i = 2; i < exp; i++) {
+    for (int i = 2; i <= exp; i++) {
         result *= value;
     }
@@ -144 +144 @@
                         kernels->yOrder->data.S32[index] = yOrder;
                         kernels->preCalc->data[index] = psMemIncrRefCounter(preCalc);
+
+                        psTrace("psModules.imcombine", 7, "Kernel %d: %f %d %d %d %d\n", index,
+                                sigmas->data.F32[i], uOrder, vOrder, xOrder, yOrder);
                     }
                 }
+
+
+                if (psTraceGetLevel("psModules.imcombine.kernel") >= 10) {
+                    psString kernelName = NULL;
+                    psStringAppend(&kernelName, "kernel%d.fits", index);
+                    psFits *kernelFile = psFitsOpen(kernelName, "w");
+                    psFree(kernelName);
+                    psFitsWriteImage(kernelFile, NULL, preCalc->image, 0, NULL);
+                    psFitsClose(kernelFile);
+                }
+
                 psFree(preCalc);        // Drop reference
             }
@@ -156 +170 @@
            kernels->xOrder->data.S32[kernels->subIndex] == 0 &&
            kernels->yOrder->data.S32[kernels->subIndex] == 0);
+
 
     return kernels;