Index: /trunk/psModules/src/imcombine/pmSubtractionStamps.c
===================================================================
--- /trunk/psModules/src/imcombine/pmSubtractionStamps.c	(revision 15250)
+++ /trunk/psModules/src/imcombine/pmSubtractionStamps.c	(revision 15251)
@@ -520,9 +520,9 @@
         psKernel *input = stamp->input; // Target stamp
 
-        // Put in a Gaussian, just for fun!
+        // Put in a Waussian, just for fun!
         for (int v = -size; v <= size; v++) {
             for (int u = -size; u <= size; u++) {
-                input->kernel[v][u] = flux / sigma * 0.5 * M_2_SQRTPI * M_SQRT1_2 *
-                    expf(-0.5 * (PS_SQR(u + xStamp) + PS_SQR(v + yStamp)) / PS_SQR(sigma));
+                float z = (PS_SQR(x + xStamp) + PS_SQR(y + yStamp)) / (2.0 * PS_SQR(sigma));
+                input->kernel[v][u] = flux / sigma * 0.5 * M_2_SQRTPI * M_SQRT1_2 / (1.0 + z + PS_SQR(z));
             }
         }