Changeset 36276 for branches/eam_branches/ipp-20130904/psModules/src/objects/pmPCMdata.c

Timestamp:

Nov 8, 2013, 3:22:59 PM (13 years ago)

Author:

eugene

Message:

working on 2d convolve issues for psphot galaxies

File:

• branches/eam_branches/ipp-20130904/psModules/src/objects/pmPCMdata.c (modified) (4 diffs)

branches/eam_branches/ipp-20130904/psModules/src/objects/pmPCMdata.c

@@ -58 +58 @@
     psFree (pcm->psfFFT);
     psFree (pcm->constraint);
-    psFree (pcm->smdata); // pre-allocated data for psImageSmooth_PreAlloc
+
+    psFree (pcm->smdata1); // pre-allocated data for psImageSmooth_PreAlloc
+    psFree (pcm->smdata2); // pre-allocated data for psImageSmooth_PreAlloc
     return;
 }
@@ -88 +90 @@
     }
 
-    pcm->smdata = NULL;
+    pcm->smdata1 = NULL;
+    pcm->smdata2 = NULL;
+
     pcm->modelConv = NULL;
     pcm->psf = NULL;
@@ -284 +288 @@
 }
 
-psImageSmoothCacheData *psImageSmoothCacheSetKernel (float *sigma, float *kappa, psImage *flux, pmModel *modelPSF) {
+// set 2 gaussians (extend to more if interesting)
+psImageSmoothCacheData *psImageSmoothCacheSetKernel (float *sigma1, float *sigma2, psImage *flux, pmModel *modelPSF) {
 
     psAssert (modelPSF, "psf model must be defined");
@@ -294 +299 @@
     
     float nsigma = 5.0;
-    *sigma = NAN;
-    *kappa = NAN;
-
+    *sigma1 = NAN;
+    *sigma2 = NAN;
+
+    // XXX need to do this more carefully
     if (modelPSF->type == modelType_GAUSS) {
         float FWHM_MAJOR = 2*modelPSF->modelRadius (modelPSF->params, 0.5*PAR[PM_PAR_I0]);
         float FWHM_MINOR = FWHM_MAJOR * (axes.minor / axes.major);
-        *sigma = 0.5 * (FWHM_MAJOR + FWHM_MINOR) / 2.35;
+        *sigma1 = 0.50 * (FWHM_MAJOR + FWHM_MINOR) / 2.35;
+        *sigma2 = 0.25 * (FWHM_MAJOR + FWHM_MINOR) / 2.35;
     }
     if (modelPSF->type == modelType_PS1_V1) {