Changeset 29506 for branches/eam_branches/ipp-20100823/psModules/src/imcombine/pmSubtractionEquation.c

Timestamp:

Oct 21, 2010, 6:14:32 AM (16 years ago)

Author:

eugene

Message:

subtract kernel(0) instead of delta function

File:

• branches/eam_branches/ipp-20100823/psModules/src/imcombine/pmSubtractionEquation.c (modified) (5 diffs)

branches/eam_branches/ipp-20100823/psModules/src/imcombine/pmSubtractionEquation.c

@@ -21 +21 @@
 
 // XXX TEST:
-# define USE_WINDOW                      // window to avoid neighbor contamination
+//# define USE_WINDOW                      // window to avoid neighbor contamination
 
 # define PENALTY false
@@ -1091 +1091 @@
         }
 
-#if 1
+#if 0
         psImage *save = psImageCopy(NULL, sumMatrix, PS_TYPE_F32);
         psFitsWriteImageSimple ("sumMatrix.fits", save, NULL);
@@ -1115 +1115 @@
         }
 
+# if (0)
         for (int i = 0; i < solution->n; i++) {
             psLogMsg("psModules.imcombine", PS_LOG_DETAIL, "Single solution %d: %lf\n", i, solution->data.F64[i]);
         }
+# endif
 
         if (!kernels->solution1) {
@@ -1169 +1171 @@
         }
 
-#if 1
+#if 0
         psImage *save = psImageCopy(NULL, sumMatrix, PS_TYPE_F32);
         psFitsWriteImageSimple ("sumMatrix.fits", save, NULL);
@@ -1187 +1189 @@
         solution = psMatrixSolveSVD(solution, sumMatrix, sumVector, NAN);
 
-#if (1)
+#if (0)
         for (int i = 0; i < solution->n; i++) {
             fprintf(stderr, "Dual solution %d: %lf\n", i, solution->data.F64[i]);