Index: /branches/meh_branches/ppstack_test/psModules/src/imcombine/pmSubtractionStamps.c
===================================================================
--- /branches/meh_branches/ppstack_test/psModules/src/imcombine/pmSubtractionStamps.c	(revision 33467)
+++ /branches/meh_branches/ppstack_test/psModules/src/imcombine/pmSubtractionStamps.c	(revision 33468)
@@ -496,7 +496,8 @@
                                              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,6 +1180,6 @@
                         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,6 +1190,6 @@
                 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]);
                     }
                 }
