Changeset 30344

Timestamp:

Jan 22, 2011, 12:00:07 PM (15 years ago)

Author:

eugene

Message:

add some test output blocks

File:

• branches/eam_branches/ipp-20101205/psModules/src/imcombine/pmSubtraction.c (modified) (4 diffs)

branches/eam_branches/ipp-20101205/psModules/src/imcombine/pmSubtraction.c

@@ -25 +25 @@
 #include "pmSubtractionVisual.h"
 #include "pmSubtractionThreads.h"
+
+bool psFitsWriteImageSimple (char *filename, psImage *image, psMetadata *header);
 
 #include "pmSubtraction.h"
@@ -1437 +1439 @@
     psFree(kernelErr2);
 
+    static int nOut1 = 0;
+    static int nOut2 = 0;
+
     // Calculate covariances
     // This can be fairly involved, so we only do it for a small number of instances
@@ -1450 +1455 @@
                 psKernelTruncate(kernel, covarFrac);
                 covars->data[i] = psImageCovarianceCalculate(kernel, ro1->covariance);
+                if (0) {
+                    char name[128];
+                    snprintf (name, 128, "covar.sample1.%03d.fits", nOut1);
+                    psKernel *cov = covars->data[i];
+                    psFitsWriteImageSimple (name, cov->image, NULL);
+
+                    snprintf (name, 128, "incovar.sample1.%03d.fits", nOut1);
+                    psFitsWriteImageSimple (name, ro1->covariance->image, NULL);
+
+                    snprintf (name, 128, "conv.sample1.%03d.fits", nOut1);
+                    psFitsWriteImageSimple (name, kernel->image, NULL);
+                    nOut1 ++;
+                }
                 psFree(kernel);
             }
@@ -1471 +1489 @@
                 psKernelTruncate(kernel, covarFrac);
                 covars->data[i] = psImageCovarianceCalculate(kernel, ro2->covariance);
+                if (0) {
+                    char name[128];
+                    snprintf (name, 128, "covar.sample2.%03d.fits", nOut2);
+                    psKernel *cov = covars->data[i];
+                    psFitsWriteImageSimple (name, cov->image, NULL);
+
+                    snprintf (name, 128, "incovar.sample2.%03d.fits", nOut1);
+                    psFitsWriteImageSimple (name, ro2->covariance->image, NULL);
+
+                    snprintf (name, 128, "conv.sample2.%03d.fits", nOut1);
+                    psFitsWriteImageSimple (name, kernel->image, NULL);
+                    nOut2 ++;
+                }
                 psFree(kernel);
             }