Changeset 25834 for branches/pap/psModules/src/imcombine/pmSubtractionAnalysis.c

Timestamp:

Oct 13, 2009, 11:30:52 PM (17 years ago)

Author:

Paul Price

Message:

Reworked dual convolution equation setup. Fixing up debugging (output of individual kernel basis functions). Not working yet --- the corresponding k1 and k2 terms appear anti-correllated: when k1_i is strongly positive, k2_i is strongly negative.

File:

• branches/pap/psModules/src/imcombine/pmSubtractionAnalysis.c (modified) (3 diffs)

branches/pap/psModules/src/imcombine/pmSubtractionAnalysis.c

@@ -16 +16 @@
 #define KERNEL_MOSAIC 2                 // Half-number of kernel instances in the mosaic image
 
+//#define TESTING
 
 bool pmSubtractionAnalysis(psMetadata *analysis, psMetadata *header,
@@ -117 +118 @@
 
 
-#if 0
+#ifdef TESTING
     // Generate images of the kernel components
     {
@@ -128 +129 @@
         }
         psArray *kernelImages = pmSubtractionKernelSolutions(kernels, 0.0, 0.0, false);
-        psFits *kernelFile = psFitsOpen("kernels.fits", "w");
+        psFits *kernelFile = psFitsOpen("kernels1.fits", "w");
+        (void)psFitsWriteImageCube(kernelFile, header, kernelImages, NULL);
+        psFitsClose(kernelFile);
+        psFree(kernelImages);
+        psFree(header);
+    }
+    if (kernels->solution2) {
+        psMetadata *header = psMetadataAlloc(); // Header
+        for (int i = 0; i < kernels->solution2->n; i++) {
+            psString name = NULL;       // Header keyword
+            psStringAppend(&name, "SOLN%04d", i);
+            psMetadataAddF64(header, PS_LIST_TAIL, name, 0, NULL, kernels->solution2->data.F64[i]);
+            psFree(name);
+        }
+        psArray *kernelImages = pmSubtractionKernelSolutions(kernels, 0.0, 0.0, true);
+        psFits *kernelFile = psFitsOpen("kernels2.fits", "w");
         (void)psFitsWriteImageCube(kernelFile, header, kernelImages, NULL);
         psFitsClose(kernelFile);