Changeset 15251

Timestamp:

Oct 9, 2007, 9:10:48 AM (19 years ago)

Author:

Paul Price

Message:

Using a Waussian instead of a Gaussian, so as to have some power in the wings. I believe the lack of power in the wings in the Gaussian is causing the convolved image PSF to go negative in the wings.

File:

• trunk/psModules/src/imcombine/pmSubtractionStamps.c (modified) (1 diff)

trunk/psModules/src/imcombine/pmSubtractionStamps.c

@@ -520 +520 @@
         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));
             }
         }