Changeset 33468 for branches/meh_branches/ppstack_test/psModules/src/imcombine/pmSubtractionStamps.c

Timestamp:

Mar 9, 2012, 3:48:25 PM (14 years ago)

Author:

mhuber

Message:

 

File:

• branches/meh_branches/ppstack_test/psModules/src/imcombine/pmSubtractionStamps.c (modified) (3 diffs)

branches/meh_branches/ppstack_test/psModules/src/imcombine/pmSubtractionStamps.c

@@ -496 +496 @@
                                              normFrac, sysErr, skyErr);
     }
-    //MEH - stamps.dat will repeatedly be overwritten - out to trace may be btter -- see end for summary 
+    //MEH - stamps.dat will repeatedly be overwritten - out to trace may be btter but good for segfault -- see end for summary 
+    //MEH - not smart dump of sources, will segfault if empty stamp - fix?
     // XXX TEST : dump all stars in the stamps here
-    if (1) {
+    if (0) {
         FILE *f = fopen ("stamp.dat", "w");
         for (int i = 0; i < stamps->num; i++) {
@@ -1179 +1180 @@
                         float additional = image1->kernel[y][x] + image2->kernel[y][x];
                         //float additional = image2->kernel[y][x];
-//                        weight->kernel[y][x] = 1.0 / (skyErr + var->kernel[y][x] + sysErr * PS_SQR(additional));
-                        weight->kernel[y][x] = 1.0 / (skyErr + 1000.0*var->kernel[y][x] + sysErr * PS_SQR(additional));
+                        weight->kernel[y][x] = 1.0 / (skyErr + var->kernel[y][x] + sysErr * PS_SQR(additional));
+//                        weight->kernel[y][x] = 1.0 / (skyErr + 1000.0*var->kernel[y][x] + sysErr * PS_SQR(additional));
                     }
                 }
@@ -1189 +1190 @@
                 for (int y = -size; y <= size; y++) {
                     for (int x = -size; x <= size; x++) {
-//                        weight->kernel[y][x] = 1.0 / var->kernel[y][x];
-                        weight->kernel[y][x] = 1.0 / (1000.0*var->kernel[y][x]);
+                        weight->kernel[y][x] = 1.0 / var->kernel[y][x];
+//                        weight->kernel[y][x] = 1.0 / (1000.0*var->kernel[y][x]);
                     }
                 }